ENH: Changing tolerance overlap. Taking out filter by small overlap

This commit is contained in:
sergio
2019-08-19 13:38:36 -07:00
committed by Andrew Heather
parent 31bfd40001
commit 8fe96e5dee

View File

@ -49,7 +49,7 @@ namespace cellCellStencils
}
Foam::scalar
Foam::cellCellStencils::cellVolumeWeight::defaultOverlapTolerance_ = 1e-9;
Foam::cellCellStencils::cellVolumeWeight::defaultOverlapTolerance_ = 1e-6;
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
@ -953,7 +953,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update()
// Interpolate. Check if enough overlap
scalar v = mesh_.V()[cellI];
scalar overlapVol = sum(allWeights[cellI]);
if (overlapVol > (1.0-overlapTolerance_)*v)
if (overlapVol > overlapTolerance_*v)
{
allCellTypes[cellI] = INTERPOLATED;
}
@ -971,7 +971,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update()
}
}
/*
// Knock out cell with insufficient interpolation weights
forAll(allCellTypes, cellI)
{
@ -989,7 +989,7 @@ bool Foam::cellCellStencils::cellVolumeWeight::update()
}
}
}
*/
if (debug)
{
tmp<volScalarField> tfld