From 5bd79b395187b042920492a8087d4b29d9a8171c Mon Sep 17 00:00:00 2001
From: Lisa Julia Nebel <lisa_julia.nebel@tu-dresden.de>
Date: Mon, 3 Apr 2023 14:57:04 +0200
Subject: [PATCH] Print infinity norm of corrGlobal instead of the correction
on each process in Riemannian TR-solver
---
dune/gfe/riemanniantrsolver.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/riemanniantrsolver.cc b/dune/gfe/riemanniantrsolver.cc
index 87ccf8f5..fb07fcf9 100644
--- a/dune/gfe/riemanniantrsolver.cc
+++ b/dune/gfe/riemanniantrsolver.cc
@@ -546,8 +546,8 @@ void RiemannianTrustRegionSolver<Basis,TargetSpace,Assembler>::solve()
std::cout << i+1 << " trust-region steps were taken, the maximum was reached." << std::endl << "Total solver time: " << totalSolverTime << " sec., total assembly time: " << totalAssemblyTime << " sec." << std::endl;
if (solved) {
- if (this->verbosity_ == NumProc::FULL)
- std::cout << "Infinity norm of the correction: " << corr.infinity_norm() << std::endl;
+ if (this->verbosity_ == NumProc::FULL and rank==0)
+ std::cout << "Infinity norm of the correction: " << corrGlobalInfinityNorm << std::endl;
if (corrGlobalInfinityNorm < this->tolerance_ && corrGlobalInfinityNorm < trustRegion.radius()*smallestScalingParameter) {
if (this->verbosity_ == NumProc::FULL and rank==0)
--
GitLab