Deactivate OpenMP multithreading during MUMPS factorization (#147)

Author: julianlitz

src/DirectSolver/DirectSolver-COO-MUMPS-Give/initializeMumps.cpp

@@ -68,7 +68,7 @@ void DirectSolver_COO_MUMPS_Give::initializeMumpsSolver(DMUMPS_STRUC_C& mumps_so
     mumps_solver.ICNTL(38) = 600; // Estimates compression rate of LU factors
     mumps_solver.ICNTL(39) = 500; // Estimates compression rate of contribution blocks
     // ICNTL(40-47) Don't exist
-    mumps_solver.ICNTL(48) = 1; // Multithreading with tree parallelism
+    mumps_solver.ICNTL(48) = 0; // Multithreading with tree parallelism
     mumps_solver.ICNTL(49) = 0; // Compact workarray id%S at the end of factorization phase
     // ICNTL(50-55) Don't exist
     mumps_solver.ICNTL(56) =

src/DirectSolver/DirectSolver-COO-MUMPS-Take/initializeMumps.cpp

@@ -68,7 +68,7 @@ void DirectSolver_COO_MUMPS_Take::initializeMumpsSolver(DMUMPS_STRUC_C& mumps_so
     mumps_solver.ICNTL(38) = 600; // Estimates compression rate of LU factors
     mumps_solver.ICNTL(39) = 500; // Estimates compression rate of contribution blocks
     // ICNTL(40-47) Don't exist
-    mumps_solver.ICNTL(48) = 1; // Multithreading with tree parallelism
+    mumps_solver.ICNTL(48) = 0; // Multithreading with tree parallelism
     mumps_solver.ICNTL(49) = 0; // Compact workarray id%S at the end of factorization phase
     // ICNTL(50-55) Don't exist
     mumps_solver.ICNTL(56) =

src/ExtrapolatedSmoother/ExtrapolatedSmootherGive/initializeMumps.cpp

@@ -68,7 +68,7 @@ void ExtrapolatedSmootherGive::initializeMumpsSolver(DMUMPS_STRUC_C& mumps_solve
     mumps_solver.ICNTL(38) = 600; // Estimates compression rate of LU factors
     mumps_solver.ICNTL(39) = 500; // Estimates compression rate of contribution blocks
     // ICNTL(40-47) Don't exist
-    mumps_solver.ICNTL(48) = 1; // Multithreading with tree parallelism
+    mumps_solver.ICNTL(48) = 0; // Multithreading with tree parallelism
     mumps_solver.ICNTL(49) = 0; // Compact workarray id%S at the end of factorization phase
     // ICNTL(50-55) Don't exist
     mumps_solver.ICNTL(56) =
@@ -109,4 +109,4 @@ void ExtrapolatedSmootherGive::finalizeMumpsSolver(DMUMPS_STRUC_C& mumps_solver)
     dmumps_c(&mumps_solver);
 }
 
-#endif
+#endif

src/ExtrapolatedSmoother/ExtrapolatedSmootherTake/initializeMumps.cpp

@@ -68,7 +68,7 @@ void ExtrapolatedSmootherTake::initializeMumpsSolver(DMUMPS_STRUC_C& mumps_solve
     mumps_solver.ICNTL(38) = 600; // Estimates compression rate of LU factors
     mumps_solver.ICNTL(39) = 500; // Estimates compression rate of contribution blocks
     // ICNTL(40-47) Don't exist
-    mumps_solver.ICNTL(48) = 1; // Multithreading with tree parallelism
+    mumps_solver.ICNTL(48) = 0; // Multithreading with tree parallelism
     mumps_solver.ICNTL(49) = 0; // Compact workarray id%S at the end of factorization phase
     // ICNTL(50-55) Don't exist
     mumps_solver.ICNTL(56) =
@@ -109,4 +109,4 @@ void ExtrapolatedSmootherTake::finalizeMumpsSolver(DMUMPS_STRUC_C& mumps_solver)
     dmumps_c(&mumps_solver);
 }
 
-#endif
+#endif

src/Smoother/SmootherGive/initializeMumps.cpp

@@ -67,7 +67,7 @@ void SmootherGive::initializeMumpsSolver(DMUMPS_STRUC_C& mumps_solver, SparseMat
     mumps_solver.ICNTL(38) = 600; // Estimates compression rate of LU factors
     mumps_solver.ICNTL(39) = 500; // Estimates compression rate of contribution blocks
     // ICNTL(40-47) Don't exist
-    mumps_solver.ICNTL(48) = 1; // Multithreading with tree parallelism
+    mumps_solver.ICNTL(48) = 0; // Multithreading with tree parallelism
     mumps_solver.ICNTL(49) = 0; // Compact workarray id%S at the end of factorization phase
     // ICNTL(50-55) Don't exist
     mumps_solver.ICNTL(56) =
@@ -108,4 +108,4 @@ void SmootherGive::finalizeMumpsSolver(DMUMPS_STRUC_C& mumps_solver)
     dmumps_c(&mumps_solver);
 }
 
-#endif
+#endif

src/Smoother/SmootherTake/initializeMumps.cpp

@@ -67,7 +67,7 @@ void SmootherTake::initializeMumpsSolver(DMUMPS_STRUC_C& mumps_solver, SparseMat
     mumps_solver.ICNTL(38) = 600; // Estimates compression rate of LU factors
     mumps_solver.ICNTL(39) = 500; // Estimates compression rate of contribution blocks
     // ICNTL(40-47) Don't exist
-    mumps_solver.ICNTL(48) = 1; // Multithreading with tree parallelism
+    mumps_solver.ICNTL(48) = 0; // Multithreading with tree parallelism
     mumps_solver.ICNTL(49) = 0; // Compact workarray id%S at the end of factorization phase
     // ICNTL(50-55) Don't exist
     mumps_solver.ICNTL(56) =
@@ -108,4 +108,4 @@ void SmootherTake::finalizeMumpsSolver(DMUMPS_STRUC_C& mumps_solver)
     dmumps_c(&mumps_solver);
 }
 
-#endif
+#endif