diff --git a/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockCG/BlockCGSolver.C b/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockCG/BlockCGSolver.C
index 375fc9a53..8f693a3f2 100644
--- a/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockCG/BlockCGSolver.C
+++ b/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockCG/BlockCGSolver.C
@@ -84,7 +84,7 @@ typename Foam::BlockSolverPerformance<Type> Foam::BlockCGSolver<Type>::solve
     Field<Type> rA(b - wA);
 
     // NOTE: Normalisation of residual per component! TU, Feb 2019
-    solverPerf.initialResidual() = cmptDivide(gSum(cmptMag(rA)),norm);
+    solverPerf.initialResidual() = cmptDivide(gSum(cmptMag(rA)), norm);
     solverPerf.finalResidual() = solverPerf.initialResidual();
 
     // Check convergence, solve if not converged