]> https://gitweb.dealii.org/ - dealii.git/commitdiff
SolverGMRES: Fix a test case 16757/head
authorMartin Kronbichler <martin.kronbichler@rub.de>
Sat, 16 Mar 2024 09:07:48 +0000 (10:07 +0100)
committerMartin Kronbichler <martin.kronbichler@rub.de>
Sat, 16 Mar 2024 09:07:48 +0000 (10:07 +0100)
tests/petsc/deal_solver_03.cc
tests/trilinos/deal_solver_03.cc
tests/trilinos/deal_solver_03.output

index f01274a84264f1f60741311d540bafe098ef86a5..057b9e5675c11b12a660e17d1292951e22b2d032 100644 (file)
@@ -64,9 +64,10 @@ main(int argc, char **argv)
     f = 1.;
     A.compress(VectorOperation::insert);
 
-    GrowingVectorMemory<PETScWrappers::MPI::Vector> mem;
-    SolverGMRES<PETScWrappers::MPI::Vector>         solver(control, mem);
-    PreconditionIdentity                            preconditioner;
+    GrowingVectorMemory<PETScWrappers::MPI::Vector>         mem;
+    SolverGMRES<PETScWrappers::MPI::Vector>::AdditionalData data(28);
+    SolverGMRES<PETScWrappers::MPI::Vector> solver(control, mem, data);
+    PreconditionIdentity                    preconditioner;
     deallog << "Solver type: " << typeid(solver).name() << std::endl;
     check_solver_within_range(solver.solve(A, u, f, preconditioner),
                               control.last_step(),
index ffa1065eb8793450c581922cca5511c93e63bf30..cc22fcc5e7d58e9defa5755c0432aa3bc163c164 100644 (file)
@@ -74,15 +74,16 @@ main(int argc, char **argv)
     f.compress(VectorOperation::insert);
     u.compress(VectorOperation::insert);
 
-    GrowingVectorMemory<TrilinosWrappers::MPI::Vector> mem;
-    SolverGMRES<TrilinosWrappers::MPI::Vector>         solver(control, mem);
-    PreconditionIdentity                               preconditioner;
+    GrowingVectorMemory<TrilinosWrappers::MPI::Vector>         mem;
+    SolverGMRES<TrilinosWrappers::MPI::Vector>::AdditionalData data(28);
+    SolverGMRES<TrilinosWrappers::MPI::Vector> solver(control, mem, data);
+    PreconditionIdentity                       preconditioner;
 
     deallog << "Solver type: " << typeid(solver).name() << std::endl;
 
     check_solver_within_range(solver.solve(A, u, f, preconditioner),
                               control.last_step(),
-                              71,
+                              74,
                               76);
   }
 }
index 3e4091158b3e03f47589f7267e0909f4c3f406fb..5b6f34323b5067354d98a2ccde4355d5f9ce520d 100644 (file)
@@ -1,4 +1,4 @@
 
 DEAL::Size 32 Unknowns 961
 DEAL::Solver type: N6dealii11SolverGMRESINS_16TrilinosWrappers3MPI6VectorEEE
-DEAL::Solver stopped within 71 - 76 iterations
+DEAL::Solver stopped within 74 - 76 iterations

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.