From 14d08c8673a032bbb6c95ee3f463103c9d48a498 Mon Sep 17 00:00:00 2001 From: Oliver Sander <sander@igpm.rwth-aachen.de> Date: Tue, 12 May 2009 12:34:34 +0000 Subject: [PATCH] bugfix: forgot to use the special trust-region norm in all places [[Imported from SVN: r2439]] --- dune-solvers/solvers/tcgsolver.cc | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/dune-solvers/solvers/tcgsolver.cc b/dune-solvers/solvers/tcgsolver.cc index e67c3af7..a9f76d1e 100644 --- a/dune-solvers/solvers/tcgsolver.cc +++ b/dune-solvers/solvers/tcgsolver.cc @@ -82,7 +82,7 @@ void TruncatedCGSolver<MatrixType, VectorType>::solve() VectorType eta_jPlus1 = eta_j; eta_jPlus1.axpy(alpha_j, delta_j); - if (eta_jPlus1*eta_jPlus1 >= trustRegionRadius_*trustRegionRadius_) { + if (trustRegionScalarProduct(eta_jPlus1,eta_jPlus1) >= trustRegionRadius_*trustRegionRadius_) { // Compute r >= 0 such that \eta = \eta_j + \tau \delta_j satisfies || \eta || = \Delta field_type tau = positiveRoot(trustRegionScalarProduct(delta_j,delta_j), -- GitLab