[Getfem-commits] [getfem-commits] branch master updated: Enable accessing parallel MUMPS in the scripting API linsolve function

[Konstantinos Poulios]

This is an automated email from the git hooks/post-receive script.

logari81 pushed a commit to branch master
in repository getfem.

The following commit(s) were added to refs/heads/master by this push:
     new 3f88b0cb Enable accessing parallel MUMPS in the scripting API linsolve 
function
3f88b0cb is described below

commit 3f88b0cbaf33dc600183b70dd5f69107e383627e
Author: Konstantinos Poulios <logari81@gmail.com>
AuthorDate: Sun Apr 28 10:37:44 2024 +0200

    Enable accessing parallel MUMPS in the scripting API linsolve function
---
 interface/src/gf_linsolve.cc | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/interface/src/gf_linsolve.cc b/interface/src/gf_linsolve.cc
index d3f7da65..5e07e7ca 100644
--- a/interface/src/gf_linsolve.cc
+++ b/interface/src/gf_linsolve.cc
@@ -106,7 +106,13 @@ mumps_solver(gsparse &gsp,
   garray<T> b = in.pop().to_garray(int(gsp.nrows()), T());
   garray<T> x = out.pop().create_array(b.getm(), b.getn(), T());
   gsp.to_csc();
-  gmm::MUMPS_solve(gsp.csc(T()),x,b);
+# if GETFEM_PARA_LEVEL > 1
+  double t_ref = MPI_Wtime();
+  gmm::MUMPS_distributed_matrix_solve(gsp.csc(T()), x, b);
+  if (getfem::MPI_IS_MASTER()) cout << "MUMPS solve time " << 
MPI_Wtime()-t_ref << endl;
+# else
+  gmm::MUMPS_solve(gsp.csc(T()), x, b);
+# endif
 }
 #endif