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

Revert to original CholeskyPrecon on coupled interfaces

Browse source
This commit is contained in:
Hrvoje Jasak 2017-07-07 12:52:44 +01:00
parent a18c6f4a23
commit d52d673767
1 changed files with 1 additions and 33 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

34
src/lduSolvers/lduPrecon/CholeskyPrecon/CholeskyPrecon.C
View file

@ -52,40 +52,8 @@ namespace Foam
void Foam::CholeskyPrecon::calcPreconDiag() void Foam::CholeskyPrecon::calcPreconDiag()
{ {
// Precondition the diagonal // Precondition the diagonal
if (matrix_.symmetric())
// Do coupled interfaces
// Note: ordering of ILU coefficients on the coupled boundary is
// out of sequence. Formally, this can be fixed by visiting the coupled
// boundary from the cell loop, but HJ thinks this does not
// make a difference
// HJ and VV, 19/Jun/2017
if (!matrix_.diagonal())
{ {
// Do coupled interfaces
forAll (interfaces_, patchI)
{
if (interfaces_.set(patchI))
{
// Get face-cells addressing
const unallocLabelList& fc =
interfaces_[patchI].coupledInterface().faceCells();
// Get interface coefficiens. Note: symmetric matrix
// HJ, 19/Jun/2017
const scalarField& bouCoeffs = coupleBouCoeffs_[patchI];
forAll (fc, coeffI)
{
// Note: sign fixed by HJ, 19/Jun/2017
preconDiag_[fc[coeffI]] +=
sqr(bouCoeffs[coeffI])/preconDiag_[fc[coeffI]];
}
}
}
// Do core matrix
const unallocLabelList& upperAddr = matrix_.lduAddr().upperAddr(); const unallocLabelList& upperAddr = matrix_.lduAddr().upperAddr();
const unallocLabelList& lowerAddr = matrix_.lduAddr().lowerAddr(); const unallocLabelList& lowerAddr = matrix_.lduAddr().lowerAddr();