diff --git a/docs/julia/heater_3d.jl b/docs/julia/heater_3d.jl
index 90da2a00..2f64f4e5 100644
--- a/docs/julia/heater_3d.jl
+++ b/docs/julia/heater_3d.jl
@@ -133,11 +133,9 @@ GridapPETSc.with(args = split(options)) do
         temperature(T0) * 0,
         2e-6,
         2e-4,
-        solver = PETScLinearSolver(),
+        #solver = PETScLinearSolver(),
     )
 
-    stophere
-
     #createpvd("poisson_transient_solution") do pvd
     #    for (uₕ, t) in uₕₜ
     #        pvd[t] = createvtk(