From d25102187a50d6abe220536d21f8f520325992ae Mon Sep 17 00:00:00 2001
From: Dominik Christ <d.christ@wikki.co.uk>
Date: Mon, 7 Jan 2013 17:55:06 +0000
Subject: [PATCH] Preconditioning bugfix by Alexander Monakov

---
 src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C b/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
index eb8066528..cae48e437 100644
--- a/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
+++ b/src/lduSolvers/lduSolver/bicgStabSolver/bicgStabSolver.C
@@ -172,8 +172,9 @@ Foam::lduSolverPerformance Foam::bicgStabSolver::solve
                 s[i] = r[i] - alpha*v[i];
             }
 
-            // Execute preconditioning transpose
-            preconPtr_->preconditionT(sh, s, cmpt);
+            // Execute preconditioning
+            // Bug fix, Alexander Monakov, 11/Jul/2012
+            preconPtr_->precondition(sh, s, cmpt);
             matrix_.Amul(t, sh, coupleBouCoeffs_, interfaces_, cmpt);
             omega = gSumProd(t, s)/gSumProd(t, t);