From 940ff62548798bca06bc953bef5e21d0f86c204b Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Mon, 31 Jan 2011 08:24:25 +0000
Subject: [PATCH] make the cgsolver not print any output -- at a paranoia check
to make sure the residual really is zero afterwards
[[Imported from SVN: r6900]]
---
.../gfe/coupling/rodcontinuumsteklovpoincarestep.hh | 13 +++++++++++--
1 file changed, 11 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
index 72902e61..89eb4f7a 100644
--- a/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
+++ b/dune/gfe/coupling/rodcontinuumsteklovpoincarestep.hh
@@ -687,14 +687,23 @@ linearizedRodNeumannToDirichletMap(const std::string& rodName,
Dune::SeqILU0<MatrixType,RodCorrectionType,RodCorrectionType> ilu0(stiffnessMatrix,1.0);
// A preconditioned conjugate-gradient solver
- Dune::CGSolver<RodCorrectionType> cg(op,ilu0,1E-4,100,2);
+ Dune::CGSolver<RodCorrectionType> cg(op,ilu0,1E-4,100,0);
// Object storing some statistics about the solving process
Dune::InverseOperatorResult statistics;
// Solve!
+#ifndef NDEBUG
+ RodCorrectionType residual = rhs;
+#endif
cg.apply(x, rhs, statistics);
-
+
+#ifndef NDEBUG
+ // paranoia
+ stiffnessMatrix.mmv(x,residual);
+ assert(residual.infinity_norm() < 1e-4);
+#endif
+
///////////////////////////////////////////////////////////////////////////////////////
// Extract the solution at the interface boundaries
///////////////////////////////////////////////////////////////////////////////////////
--
GitLab