From e9648a1635f3ee50f9e0e732b8e7aa05cf66726c Mon Sep 17 00:00:00 2001 From: Peter Boyle Date: Tue, 14 Jun 2022 23:40:04 -0400 Subject: [PATCH] Useful periodic print. CG convergence bound is remarkably accurate on low eigenvalue in numerical tests --- Grid/algorithms/iterative/ConjugateGradient.h | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/Grid/algorithms/iterative/ConjugateGradient.h b/Grid/algorithms/iterative/ConjugateGradient.h index 14f3d306..c20e267d 100644 --- a/Grid/algorithms/iterative/ConjugateGradient.h +++ b/Grid/algorithms/iterative/ConjugateGradient.h @@ -152,8 +152,13 @@ public: LinearCombTimer.Stop(); LinalgTimer.Stop(); - std::cout << GridLogIterative << "ConjugateGradient: Iteration " << k + if ( (k % 500) == 0 ) { + std::cout << GridLogMessage << "ConjugateGradient: Iteration " << k << " residual " << sqrt(cp/ssq) << " target " << Tolerance << std::endl; + } else { + std::cout << GridLogIterative << "ConjugateGradient: Iteration " << k + << " residual " << sqrt(cp/ssq) << " target " << Tolerance << std::endl; + } // Stopping condition if (cp <= rsq) {