diff --git a/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockILU/BlockILUSolver.C b/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockILU/BlockILUSolver.C index da84db072..823c8542e 100644 --- a/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockILU/BlockILUSolver.C +++ b/src/foam/matrices/blockLduMatrix/BlockLduSolvers/BlockILU/BlockILUSolver.C @@ -76,7 +76,7 @@ Foam::BlockILUSolver::solve Field residual(x.size()); - // Calculate residual. Note: sign of residual swapped for efficiency + // Calculate residual matrix.Amul(residual, x); forAll (b, i)