Skip to content
Snippets Groups Projects
Commit 14d08c86 authored by Oliver Sander's avatar Oliver Sander Committed by sander@PCPOOL.MI.FU-BERLIN.DE
Browse files

bugfix: forgot to use the special trust-region norm in all places 

[[Imported from SVN: r2439]]
parent 384ea5f7
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune-solvers/solvers/tcgsolver.cc
+ 1
− 1
View file @ 14d08c86
......@@ -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),
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment