diff --git a/src/lduSolvers/lduPrecon/CholeskyPrecon/CholeskyPrecon.C b/src/lduSolvers/lduPrecon/CholeskyPrecon/CholeskyPrecon.C
index e38815675..80e0ede30 100644
--- a/src/lduSolvers/lduPrecon/CholeskyPrecon/CholeskyPrecon.C
+++ b/src/lduSolvers/lduPrecon/CholeskyPrecon/CholeskyPrecon.C
@@ -74,7 +74,6 @@ void Foam::CholeskyPrecon::calcPreconDiag()
                 // Get interface coefficiens. Note: symmetric matrix
                 // HJ, 19/Jun/2017
                 const scalarField& bouCoeffs = coupleBouCoeffs_[patchI];
-                const scalarField& intCoeffs = coupleIntCoeffs_[patchI];
 
                 forAll (fc, coeffI)
                 {
@@ -184,6 +183,9 @@ void Foam::CholeskyPrecon::precondition
             << abort(FatalError);
     }
 
+    // Note: coupled boundary updated is not needed because x is zero
+    // HJ and VV, 19/Jun/2017
+    
     // Diagonal block
     {
         scalar* __restrict__ xPtr = x.begin();
diff --git a/src/lduSolvers/lduPrecon/ILU0/ILU0.C b/src/lduSolvers/lduPrecon/ILU0/ILU0.C
index 06ea1cb38..d4a843fe3 100644
--- a/src/lduSolvers/lduPrecon/ILU0/ILU0.C
+++ b/src/lduSolvers/lduPrecon/ILU0/ILU0.C
@@ -186,8 +186,9 @@ void Foam::ILU0::precondition
             << abort(FatalError);
     }
 
-
-
+    // Note: coupled boundary updated is not needed because x is zero
+    // HJ and VV, 19/Jun/2017
+    
     // Diagonal block
     {
         scalar* __restrict__ xPtr = x.begin();