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

Clean-up of coupled updates on scalar Cholesky and ILU0 preconditioners

Browse source
This commit is contained in:
Hrvoje Jasak 2017-06-19 18:17:34 +01:00
parent dc0ad86c6f
commit bc7cc7ca9b
2 changed files with 6 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -74,7 +74,6 @@ void Foam::CholeskyPrecon::calcPreconDiag()
// Get interface coefficiens. Note: symmetric matrix // Get interface coefficiens. Note: symmetric matrix
// HJ, 19/Jun/2017 // HJ, 19/Jun/2017
const scalarField& bouCoeffs = coupleBouCoeffs_[patchI]; const scalarField& bouCoeffs = coupleBouCoeffs_[patchI];
const scalarField& intCoeffs = coupleIntCoeffs_[patchI];
forAll (fc, coeffI) forAll (fc, coeffI)
{ {
@ -184,6 +183,9 @@ void Foam::CholeskyPrecon::precondition
<< abort(FatalError); << abort(FatalError);
} }
// Note: coupled boundary updated is not needed because x is zero
// HJ and VV, 19/Jun/2017
// Diagonal block // Diagonal block
{ {
scalar* __restrict__ xPtr = x.begin(); scalar* __restrict__ xPtr = x.begin();

5
src/lduSolvers/lduPrecon/ILU0/ILU0.C
View file

@ -186,8 +186,9 @@ void Foam::ILU0::precondition
<< abort(FatalError); << abort(FatalError);
} }
// Note: coupled boundary updated is not needed because x is zero
// HJ and VV, 19/Jun/2017
// Diagonal block // Diagonal block
{ {
scalar* __restrict__ xPtr = x.begin(); scalar* __restrict__ xPtr = x.begin();