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

Removed bPrime update on reverseSweep in symGaussSeidel preconditioner

Browse source
This commit is contained in:
Vuko Vukcevic 2015-05-20 12:26:21 +02:00 committed by Hrvoje Jasak
parent 5c260e9115
commit 239cb5d376
1 changed files with 2 additions and 11 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
src/lduSolvers/lduPrecon/symGaussSeidelPrecon/symGaussSeidelPrecon.C
View file

@ -203,11 +203,7 @@ void Foam::symGaussSeidelPrecon::precondition
// Finish current x // Finish current x
curX /= diagPtr[rowI]; curX /= diagPtr[rowI];
// Distribute the neighbour side using current x // No need to update bPrime on reverse sweep. VV, 20/May/2015.
for (register label curCoeff = fStart; curCoeff < fEnd; curCoeff++)
{
bPrimePtr[uPtr[curCoeff]] -= lowerPtr[curCoeff]*curX;
}
} }
} }
} }
@ -325,12 +321,7 @@ void Foam::symGaussSeidelPrecon::preconditionT
// Finish current x // Finish current x
curX /= diagPtr[rowI]; curX /= diagPtr[rowI];
// Distribute the neighbour side using current x // No need to update bPrime on reverse sweep. VV, 20/May/2015.
for (register label curCoeff = fStart; curCoeff < fEnd; curCoeff++)
{
// Transpose multiplication. HJ, 10/Jul/2007
bPrimePtr[uPtr[curCoeff]] -= upperPtr[curCoeff]*curX;
}
} }
} }
} }