exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Removed unnecessary bPrime update on reverseSweep in BlockGaussSeidel preconditioner

Browse source
This commit is contained in:
Vuko Vukcevic 2015-09-10 10:22:13 +02:00 committed by Hrvoje Jasak
parent 22a2e59c3b
commit 2bebc90eaa
1 changed files with 2 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/foam/matrices/blockLduMatrix/BlockLduPrecons/BlockGaussSeidelPrecon/BlockGaussSeidelPrecon.C
View file

@ -204,13 +204,7 @@ void Foam::BlockGaussSeidelPrecon<Type>::BlockSweep
// Finish current x // Finish current x
curX = mult(dD[rowI], curX); curX = mult(dD[rowI], curX);
// Distribute the neighbour side using current x // No need to update bPrime on reverse sweep. VV, 10/Sep/2015.
for (curCoeff = fStart; curCoeff < fEnd; curCoeff++)
{
// lower = upper transposed
bPrime_[u[curCoeff]] -=
mult(mult.transpose(upper[curCoeff]), curX);
}
} }
} }
} }
@ -312,11 +306,7 @@ void Foam::BlockGaussSeidelPrecon<Type>::BlockSweep
// Finish current x // Finish current x
curX = mult(dD[rowI], curX); curX = mult(dD[rowI], curX);
// Distribute the neighbour side using current x // No need to update bPrime on reverse sweep. VV, 10/Sep/2015.
for (curCoeff = fStart; curCoeff < fEnd; curCoeff++)
{
bPrime_[u[curCoeff]] -= mult(lower[curCoeff], curX);
}
} }
} }
} }