Activate multigrid for temperature solvers

UniMainzGeo · Dec 20, 2023 · 785ab12 · 785ab12
1 parent 5eae4a0
commit 785ab12
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 0 deletions.
diff --git a/src/JacResTemp.cpp b/src/JacResTemp.cpp
@@ -210,6 +210,12 @@ PetscErrorCode JacResCreateTempParam(JacRes *jr)
 
 	// create temperature diffusion solver
 	PetscCall(KSPCreate(PETSC_COMM_WORLD, &jr->tksp));
+
+	// enable geometric multigrid
+	PetscCall(KSPSetDM(jr->tksp, jr->DA_T));
+	PetscCall(KSPSetDMActive(jr->tksp, PETSC_FALSE));
+
+	// set options
 	PetscCall(KSPSetOptionsPrefix(jr->tksp,"ts_"));
 	PetscCall(KSPSetFromOptions(jr->tksp));
 

diff --git a/src/LaMEMLib.cpp b/src/LaMEMLib.cpp
@@ -969,6 +969,12 @@ PetscErrorCode LaMEMLibSolveTemp(LaMEMLib *lm, PetscScalar dt)
 
 	// create temperature diffusion solver
 	ierr = KSPCreate(PETSC_COMM_WORLD, &tksp); CHKERRQ(ierr);
+
+	// enable geometric multigrid
+	PetscCall(KSPSetDM(tksp, jr->DA_T));
+	PetscCall(KSPSetDMActive(tksp, PETSC_FALSE));
+
+	// set options
 	ierr = KSPSetOptionsPrefix(tksp,"its_");   CHKERRQ(ierr);
 	ierr = KSPSetFromOptions(tksp);            CHKERRQ(ierr);