Merge pull request #240 from Exawind/klu2_solver

CPU based runs use KLU2 rather than Basker for a Linear Solver.

src/restruct_poc/solver/solver.hpp

@@ -31,9 +31,9 @@ namespace openturbine {
 /// solver variables during the simulation
 struct Solver {
     // Define some types for the solver to make the code more readable
-    using ExecutionSpace = Kokkos::DefaultExecutionSpace;
-    using MemorySpace = ExecutionSpace::memory_space;
     using GlobalCrsMatrixType = Tpetra::CrsMatrix<>;
+    using ExecutionSpace = GlobalCrsMatrixType::execution_space;
+    using MemorySpace = GlobalCrsMatrixType::memory_space;
     using GlobalMapType = GlobalCrsMatrixType::map_type;
     using GlobalMultiVectorType = Tpetra::MultiVector<>;
     using DualViewType = GlobalMultiVectorType::dual_view_type;
@@ -306,8 +306,12 @@ struct Solver {
         b = Tpetra::createMultiVector<ScalarType>(A->getRangeMap(), 1);
         x_mv = Tpetra::createMultiVector<ScalarType>(A->getDomainMap(), 1);
 
+        const auto solver_name = (std::is_same_v<ExecutionSpace, Kokkos::DefaultHostExecutionSpace>)
+                                     ? std::string{"klu2"}
+                                     : std::string{"basker"};
+
         amesos_solver =
-            Amesos2::create<GlobalCrsMatrixType, GlobalMultiVectorType>("Basker", A, x_mv, b);
+            Amesos2::create<GlobalCrsMatrixType, GlobalMultiVectorType>(solver_name, A, x_mv, b);
         amesos_solver->symbolicFactorization();
     }
 };