extrapolate intersection distance

This commit is contained in:
mattijs
2009-06-09 07:52:25 +01:00
parent 2a62ec3f39
commit 550a743b41

View File

@ -117,13 +117,13 @@ void Foam::meshRefinement::calcNeighbourData
// Extrapolate the face centre. // Extrapolate the face centre.
vector fn = faceAreas[i]; vector fn = faceAreas[i];
fn /= mag(fn)+VSMALL; fn /= mag(fn)+VSMALL;
// Normal distance from face centre to cell centre
scalar d = ((faceCentres[i] - cellCentres[i]) & fn);
label own = faceCells[i]; label own = faceCells[i];
label ownLevel = cellLevel[own]; label ownLevel = cellLevel[own];
label faceLevel = meshCutter_.getAnchorLevel(pp.start()+i); label faceLevel = meshCutter_.getAnchorLevel(pp.start()+i);
// Normal distance from face centre to cell centre
scalar d = ((faceCentres[i] - cellCentres[own]) & fn);
if (faceLevel > ownLevel) if (faceLevel > ownLevel)
{ {
// Other cell more refined. Adjust normal distance // Other cell more refined. Adjust normal distance