Use deal.II CG solver instead of trilinos

geodynamics · Mar 22, 2024 · dcf23f4 · dcf23f4
1 parent 722dc0a
commit dcf23f4
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/source/simulator/solver.cc b/source/simulator/solver.cc
@@ -26,8 +26,7 @@
 
 #include <deal.II/base/signaling_nan.h>
 #include <deal.II/lac/solver_gmres.h>
-
-#include <deal.II/lac/trilinos_solver.h>
+#include <deal.II/lac/solver_cg.h>
 
 #include <deal.II/fe/fe_values.h>
 
@@ -276,7 +275,8 @@ namespace aspect
       {
         SolverControl solver_control(1000, src.block(1).l2_norm() * S_block_tolerance);
 
-        TrilinosWrappers::SolverCG solver(solver_control);
+        PrimitiveVectorMemory<LinearAlgebra::Vector> mem;
+        SolverCG<LinearAlgebra::Vector> solver(solver_control,mem);
 
         // Trilinos reports a breakdown
         // in case src=dst=0, even
@@ -323,7 +323,9 @@ namespace aspect
       if (do_solve_A == true)
         {
           SolverControl solver_control(10000, utmp.l2_norm() * A_block_tolerance);
-          TrilinosWrappers::SolverCG solver(solver_control);
+          PrimitiveVectorMemory<LinearAlgebra::Vector> mem;
+          SolverCG<LinearAlgebra::Vector> solver(solver_control,mem);
+
           try
             {
               dst.block(0) = 0.0;