diff --git a/docs/literate/tutorials/dam_break.jl b/docs/literate/tutorials/dam_break.jl
index 630c8bc1..b311c4e5 100644
--- a/docs/literate/tutorials/dam_break.jl
+++ b/docs/literate/tutorials/dam_break.jl
@@ -261,7 +261,7 @@ function variational_multiscale_method(state)
             ## Solve nonlinear system
             U = zeros(ndofs(dofmap)) # Initialize nodal dispacement and pressure with zero
             solved = Tesserae.newton!(U, U->residual(U,state), U->J;
-                                      linsolve=(x,A,b)->ldiv!(x,A,b), backtracking=true, verbose=true)
+                                      linsolve=(x,A,b)->ldiv!(x,A,b), backtracking=true, verbose=true, rtol=1e-12)
             solved && break
         end