From c7121e7df251b3b01a312c3ee054f5cd0f8c71bf Mon Sep 17 00:00:00 2001
From: "Bernhard F.W. Gschaider" <bgschaid@ice-sf.at>
Date: Thu, 18 Jul 2013 23:59:38 +0200
Subject: [PATCH] BUGFIX: For the diagonal the size of the upper Matrix
 coefficients was used which led to an 'out of range'-error

--HG--
branch : bgschaid/minorAdditionsBranch
---
 .../BlockLduMatrix/BlockLduMatrixATmul.C      | 20 +++++++++----------
 1 file changed, 10 insertions(+), 10 deletions(-)

diff --git a/src/OpenFOAM/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrixATmul.C b/src/OpenFOAM/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrixATmul.C
index 1e4612ace..298c6415d 100644
--- a/src/OpenFOAM/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrixATmul.C
+++ b/src/OpenFOAM/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrixATmul.C
@@ -42,7 +42,7 @@ void Foam::BlockLduMatrix<Type>::Amul
 
     // Initialise the update of coupled interfaces
     initInterfaces(coupleUpper_, Ax, x);
-    
+
     AmulCore(Ax, x);
 
     // Update coupled interfaces
@@ -64,13 +64,13 @@ void Foam::BlockLduMatrix<Type>::AmulCore
     const unallocLabelList& u = lduAddr().upperAddr();
     const unallocLabelList& l = lduAddr().lowerAddr();
 
-    
+
     const TypeCoeffField& Diag = this->diag();
     const TypeCoeffField& Upper = this->upper();
 
     // Create multiplication function object
     typename BlockCoeff<Type>::multiply mult;
-    
+
     // AmulCore must be additive to account for initialisation step
     // in ldu interfaces.  HJ, 6/Nov/2007
     // Fixed by IC 19/Oct/2011
@@ -78,7 +78,7 @@ void Foam::BlockLduMatrix<Type>::AmulCore
     {
         const scalarTypeField& activeDiag = Diag.asScalar();
 
-        for (register label cellI = 0; cellI < u.size(); cellI++)
+        for (register label cellI = 0; cellI < activeDiag.size(); cellI++)
         {
             Ax[cellI] += mult(activeDiag[cellI], x[cellI]);
         }
@@ -87,7 +87,7 @@ void Foam::BlockLduMatrix<Type>::AmulCore
     {
         const linearTypeField& activeDiag = Diag.asLinear();
 
-        for (register label cellI = 0; cellI < u.size(); cellI++)
+        for (register label cellI = 0; cellI < activeDiag.size(); cellI++)
         {
             Ax[cellI] += mult(activeDiag[cellI], x[cellI]);
         }
@@ -96,7 +96,7 @@ void Foam::BlockLduMatrix<Type>::AmulCore
     {
         const squareTypeField& activeDiag = Diag.asSquare();
 
-        for (register label cellI = 0; cellI < u.size(); cellI++)
+        for (register label cellI = 0; cellI < activeDiag.size(); cellI++)
         {
             Ax[cellI] += mult(activeDiag[cellI], x[cellI]);
         }
@@ -248,7 +248,7 @@ void Foam::BlockLduMatrix<Type>::TmulCore
     {
         const scalarTypeField& activeDiag = Diag.asScalar();
 
-        for (register label cellI = 0; cellI < u.size(); cellI++)
+        for (register label cellI = 0; cellI < activeDiag.size(); cellI++)
         {
             Tx[cellI] += mult(activeDiag[cellI], x[cellI]);
         }
@@ -257,7 +257,7 @@ void Foam::BlockLduMatrix<Type>::TmulCore
     {
         const linearTypeField& activeDiag = Diag.asLinear();
 
-        for (register label cellI = 0; cellI < u.size(); cellI++)
+        for (register label cellI = 0; cellI < activeDiag.size(); cellI++)
         {
             Tx[cellI] += mult(activeDiag[cellI], x[cellI]);
         }
@@ -266,12 +266,12 @@ void Foam::BlockLduMatrix<Type>::TmulCore
     {
         const squareTypeField& activeDiag = Diag.asSquare();
 
-        for (register label cellI = 0; cellI < u.size(); cellI++)
+        for (register label cellI = 0; cellI < activeDiag.size(); cellI++)
         {
             Tx[cellI] += mult(activeDiag[cellI].T(), x[cellI]);
         }
     }
-    
+
     // Upper multiplication
 
     if (Upper.activeType() == blockCoeffBase::SCALAR)