Merge pull request #448 from sblauth/fix/petsc_destroy

Add destroy method and garbage cleanup for snes interface
sblauth · Jul 4, 2024 · 333b1a7 · 333b1a7
2 parents 9b47aa9 + 4febda3
commit 333b1a7
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/cashocs/nonlinear_solvers/snes.py b/cashocs/nonlinear_solvers/snes.py
@@ -225,6 +225,11 @@ def solve(self) -> fenics.Function:
         if converged_reason < 0:
             raise _exceptions.PETScSNESError(converged_reason)
 
+        if hasattr(PETSc, "garbage_cleanup"):
+            snes.destroy()
+            PETSc.garbage_cleanup(comm=self.comm)
+            PETSc.garbage_cleanup()
+
         return self.u