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

Bugfix: remove bouCoeffs for explicit fringe

Browse source
This commit is contained in:
Hrvoje Jasak 2019-09-23 17:47:14 +01:00
parent 388ff4d914
commit 7d3ffd2f51
1 changed files with 8 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
src/overset/oversetMesh/oversetFvPatchField/oversetFvPatchField.C
View file

@ -573,7 +573,7 @@ void oversetFvPatchField<Type>::correctOffDiag
// Record the coefficients into internal/boundary coeffs for AMG // Record the coefficients into internal/boundary coeffs for AMG
// coarsening // coarsening
Field<Type>& bouCoeffs = eqn.boundaryCoeffs()[this->patch().index()]; Field<Type>& bouCoeffs = eqn.boundaryCoeffs()[this->patch().index()];
bouCoeffs.setSize(fringeFaces.size()); bouCoeffs.setSize(fringeFaces.size());
@ -998,7 +998,11 @@ void oversetFvPatchField<Type>::manipulateMatrix
); );
eqn.setValues(holeCells, holeCellsPsi); eqn.setValues(holeCells, holeCellsPsi);
Field<Type>& bouCoeffs = eqn.boundaryCoeffs()[this->patch().index()];
Field<Type>& intCoeffs = eqn.internalCoeffs()[this->patch().index()];
Info<< "In manipulateMatrix: " << coupledFringe_
<< " bou: " << bouCoeffs.size() << " int: " << intCoeffs.size() << endl;
// If the fringe is not coupled, set values in acceptor cells and do not // If the fringe is not coupled, set values in acceptor cells and do not
// perform overset interpolation during the solution process. Note that the // perform overset interpolation during the solution process. Note that the
// fringeUpper/Lower coefficients are collected in correctOffDiag member // fringeUpper/Lower coefficients are collected in correctOffDiag member
@ -1017,6 +1021,9 @@ void oversetFvPatchField<Type>::manipulateMatrix
const labelList& accCells = oversetPatch_.overset().acceptorCells(); const labelList& accCells = oversetPatch_.overset().acceptorCells();
eqn.setValues(accCells, accValues); eqn.setValues(accCells, accValues);
// Kill boundary coefficients: explicit coupling
eqn.boundaryCoeffs()[this->patch().index()].setSize(0);;
} }
} }