From 1bad64ac6a3710a22abb7432266f48240c57ff40 Mon Sep 17 00:00:00 2001 From: Daniel Richtmann Date: Fri, 27 Oct 2017 14:35:04 +0200 Subject: [PATCH] Some formatting --- lib/algorithms/iterative/MinimalResidual.h | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/lib/algorithms/iterative/MinimalResidual.h b/lib/algorithms/iterative/MinimalResidual.h index 87dc5941..beda875b 100644 --- a/lib/algorithms/iterative/MinimalResidual.h +++ b/lib/algorithms/iterative/MinimalResidual.h @@ -72,7 +72,7 @@ template class MinimalResidual : public OperatorFunction { RealD cp = norm2(r); // Cp = |r[0]|^2 // 2 Nc Ns flops // flopcount.addSiteFlops(4*Nc*Ns, s); - if(cp <= rsd_sq) { + if (cp <= rsd_sq) { return; } @@ -85,7 +85,7 @@ template class MinimalResidual : public OperatorFunction { SolverTimer.Start(); int k; - for(k = 1; k <= MaxIterations; k++) { // a[k-1] := < M.r[k-1], r[k-1] >/ < M.r[k-1], M.r[k-1] > + for (k = 1; k <= MaxIterations; k++) { // a[k-1] := < M.r[k-1], r[k-1] >/ < M.r[k-1], M.r[k-1] > MatrixTimer.Start(); Linop.Op(r, Mr); // Mr = M * r // flopcount.addFlops(M.nFlops()); @@ -114,7 +114,7 @@ template class MinimalResidual : public OperatorFunction { std::cout << GridLogDebug << "a = " << a << " c = " << c << " d = " << d << std::endl; // Stopping condition - if(cp <= rsd_sq) { + if (cp <= rsd_sq) { SolverTimer.Stop(); Linop.Op(psi, Mr); @@ -146,7 +146,7 @@ template class MinimalResidual : public OperatorFunction { std::cout << GridLogMessage << "MinimalResidual did NOT converge" << std::endl; - if(ErrorOnNoConverge) + if (ErrorOnNoConverge) assert(0); IterationsToComplete = k;