diff --git a/src/overset/oversetMesh/oversetFvPatchField/oversetFvPatchField.C b/src/overset/oversetMesh/oversetFvPatchField/oversetFvPatchField.C index cdd0a4972..baaac0d0b 100644 --- a/src/overset/oversetMesh/oversetFvPatchField/oversetFvPatchField.C +++ b/src/overset/oversetMesh/oversetFvPatchField/oversetFvPatchField.C @@ -573,7 +573,7 @@ void oversetFvPatchField::correctOffDiag // Record the coefficients into internal/boundary coeffs for AMG // coarsening - + Field& bouCoeffs = eqn.boundaryCoeffs()[this->patch().index()]; bouCoeffs.setSize(fringeFaces.size()); @@ -998,7 +998,11 @@ void oversetFvPatchField::manipulateMatrix ); eqn.setValues(holeCells, holeCellsPsi); + Field& bouCoeffs = eqn.boundaryCoeffs()[this->patch().index()]; + Field& 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 // perform overset interpolation during the solution process. Note that the // fringeUpper/Lower coefficients are collected in correctOffDiag member @@ -1017,6 +1021,9 @@ void oversetFvPatchField::manipulateMatrix const labelList& accCells = oversetPatch_.overset().acceptorCells(); eqn.setValues(accCells, accValues); + + // Kill boundary coefficients: explicit coupling + eqn.boundaryCoeffs()[this->patch().index()].setSize(0);; } }