zhyang-dev/openfoam
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

GaussSeidelSmoother: minor update

Browse Source
This commit is contained in:
Henry
2012-12-16 22:29:53 +00:00
parent f8d14b0b28
commit fae37af352
1 changed files with 13 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

34
src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C
View File

@ -107,20 +107,12 @@ void Foam::GaussSeidelSmoother::smooth
// To compensate for this, it is necessary to turn the
// sign of the contribution.
//FieldField<Field, scalar> mBouCoeffs(interfaceBouCoeffs_.size());
//
//forAll(mBouCoeffs, patchi)
//{
// if (interfaces_.set(patchi))
// {
// mBouCoeffs.set(patchi, -interfaceBouCoeffs_[patchi]);
// }
//}
FieldField<Field, scalar>& mBouCoeffs =
const_cast<FieldField<Field, scalar>&>
(
interfaceBouCoeffs_
);
forAll(mBouCoeffs, patchi)
{
if (interfaces_.set(patchi))
@ -152,35 +144,35 @@ void Foam::GaussSeidelSmoother::smooth
cmpt
);
register scalar curPsi;
register scalar psii;
register label fStart;
register label fEnd = ownStartPtr[0];
for (register label cellI=0; cellI<nCells; cellI++)
for (register label celli=0; celli<nCells; celli++)
{
// Start and end of this row
fStart = fEnd;
fEnd = ownStartPtr[cellI + 1];
fEnd = ownStartPtr[celli + 1];
// Get the accumulated neighbour side
curPsi = bPrimePtr[cellI];
psii = bPrimePtr[celli];
// Accumulate the owner product side
for (register label curFace=fStart; curFace<fEnd; curFace++)
for (register label facei=fStart; facei<fEnd; facei++)
{
curPsi -= upperPtr[curFace]*psiPtr[uPtr[curFace]];
psii -= upperPtr[facei]*psiPtr[uPtr[facei]];
}
// Finish current psi
curPsi /= diagPtr[cellI];
// Finish psi for this cell
psii /= diagPtr[celli];
// Distribute the neighbour side using current psi
for (register label curFace=fStart; curFace<fEnd; curFace++)
// Distribute the neighbour side using psi for this cell
for (register label facei=fStart; facei<fEnd; facei++)
{
bPrimePtr[uPtr[curFace]] -= lowerPtr[curFace]*curPsi;
bPrimePtr[uPtr[facei]] -= lowerPtr[facei]*psii;
}
psiPtr[cellI] = curPsi;
psiPtr[celli] = psii;
}
}