]> https://gitweb.dealii.org/ - dealii.git/commitdiff
VectorTools: Increase max iterations in project_matrix_free 6079/head
authorMatthias Maier <tamiko@43-1.org>
Tue, 20 Mar 2018 21:07:03 +0000 (16:07 -0500)
committerMatthias Maier <tamiko@43-1.org>
Tue, 20 Mar 2018 21:07:03 +0000 (16:07 -0500)
Here is a funny one. Running the test
  numerics/project_q_hierarchical_2
with DEAL_II_NUM_THREADS=2 OMP_NUM_THREADS=2 (as done on the tester)
worsens the reduction rate enough to run into the maximal iteration
limit. Solve this issue by slightly increasing the maximal allowed
number of iterations.

include/deal.II/numerics/vector_tools.templates.h

index e4591b66e5a666b52a89f39da593bb8af08eb4a0..ca54be3b5af6cc15f23f154fbd76240670ff2893 100644 (file)
@@ -1091,11 +1091,11 @@ namespace VectorTools
       }
 
       // now invert the matrix
-      // Allow for a maximum of 5*n steps to reduce the residual by 10^-12. n
+      // Allow for a maximum of 6*n steps to reduce the residual by 10^-12. n
       // steps may not be sufficient, since roundoff errors may accumulate for
       // badly conditioned matrices. This behavior can be observed, e.g. for
       // FE_Q_Hierarchical for degree higher than three.
-      ReductionControl     control(5.*rhs.size(), 0., 1e-12, false, false);
+      ReductionControl     control(6.*rhs.size(), 0., 1e-12, false, false);
       SolverCG<LinearAlgebra::distributed::Vector<Number> > cg(control);
       PreconditionJacobi<MatrixType> preconditioner;
       preconditioner.initialize(mass_matrix, 1.);

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.