diff --git a/src/finiteVolume/interpolation/interpolation/interpolationPoint/interpolationPoint.H b/src/finiteVolume/interpolation/interpolation/interpolationPoint/interpolationPoint.H index 49ee13e308e38c03c06ad55d7525f739182c0751..a3dc106eb0295b42d5761d4f60ec8871c23c69aa 100644 --- a/src/finiteVolume/interpolation/interpolation/interpolationPoint/interpolationPoint.H +++ b/src/finiteVolume/interpolation/interpolation/interpolationPoint/interpolationPoint.H @@ -60,7 +60,7 @@ protected: public: //- Runtime type information - TypeName("cellPoint"); + TypeName("point"); // Constructors diff --git a/src/finiteVolume/interpolation/interpolation/interpolationPoint/pointMVCWeight.C b/src/finiteVolume/interpolation/interpolation/interpolationPoint/pointMVCWeight.C index 2945198e46d488a515e886916db16b995989dd0b..42e201c6ab6d9b52e9680da4008993c7b8b1b93e 100644 --- a/src/finiteVolume/interpolation/interpolation/interpolationPoint/pointMVCWeight.C +++ b/src/finiteVolume/interpolation/interpolation/interpolationPoint/pointMVCWeight.C @@ -119,7 +119,15 @@ void Foam::pointMVCWeight::calcWeights //Pout<< " uj:" << u[j] << " ujPlus1:" << u[jPlus1] << endl; vector temp = u[j] ^ u[jPlus1]; - temp /= mag(temp); + + scalar magTemp = mag(temp); + + if (magTemp < VSMALL) + { + continue; + } + + temp /= magTemp; //Pout<< " uj:" << u[j] << " ujPlus1:" << u[jPlus1] // << " temp:" << temp << endl;