Just apply Jacobi when computing facet areas

thetisproject · Oct 28, 2021 · 327cac1 · 327cac1
1 parent 4a3d90d
commit 327cac1
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/thetis/utility.py b/thetis/utility.py
@@ -563,7 +563,13 @@ def get_facet_areas(mesh):
     facet_areas = Function(HDivTrace, name="Facet areas")
     mass_term = v('+')*u('+')*dS + v*u*ds
     rhs = v('+')*FacetArea(mesh)*dS + v*FacetArea(mesh)*ds
-    solve(mass_term == rhs, facet_areas)
+    sp = {
+        "mat_type": "matfree",
+        "snes_type": "ksponly",
+        "ksp_type": "preonly",
+        "pc_type": "jacobi",
+    }
+    solve(mass_term == rhs, facet_areas, solver_parameters=sp)
     return facet_areas