From 629999a45eb95c010da36c49f864d6f5aa0b6fd5 Mon Sep 17 00:00:00 2001
From: Helge Gehring <42973196+HelgeGehring@users.noreply.github.com>
Date: Wed, 18 Oct 2023 19:37:55 -0700
Subject: [PATCH] use LU solver

---
 docs/julia/heater_3d.jl | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

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(