Use Amesos2 Solver Package instead of just Mumps

wortiz · Apr 15, 2024 · 4ed0fc7 · 4ed0fc7
1 parent 4b42444
commit 4ed0fc7
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/sl_amesos2_interface.cpp b/src/sl_amesos2_interface.cpp
@@ -73,7 +73,7 @@ int amesos2_solve(struct GomaLinearSolverData *ams,
 
     if (solver_data->solver.is_null()) {
       Teuchos::RCP<MAT> crs_matrix = Teuchos::rcp_dynamic_cast<MAT>(tpetra_data->matrix);
-      solver_data->solver = Amesos2::create<MAT, MV>("Mumps", crs_matrix);
+      solver_data->solver = Amesos2::create<MAT, MV>(amesos2_solver, crs_matrix);
       if (amesos2_file != NULL && strlen(amesos2_file) > 0) {
         std::filesystem::path path(amesos2_file);
         Teuchos::RCP<Teuchos::ParameterList> amesos2_params;