zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

plane class: adding normalisation of normal vector for vector only

Browse Source 
constructor.
This commit is contained in:
graham
2010-01-28 13:56:13 +00:00
parent f026a5af0e
commit d5829988d8
1 changed files with 19 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
src/OpenFOAM/meshes/primitiveShapes/plane/plane.C
View File

@ -128,7 +128,20 @@ Foam::plane::plane(const vector& normalVector)
: :
unitVector_(normalVector), unitVector_(normalVector),
basePoint_(vector::zero) basePoint_(vector::zero)
{} {
scalar magUnitVector(mag(unitVector_));
if (magUnitVector > VSMALL)
{
unitVector_ /= magUnitVector;
}
else
{
FatalErrorIn("plane::plane(const point&, const vector&)")
<< "plane normal has zero length"
<< abort(FatalError);
}
}
// Construct from point and normal vector // Construct from point and normal vector
@ -146,7 +159,7 @@ Foam::plane::plane(const point& basePoint, const vector& normalVector)
else else
{ {
FatalErrorIn("plane::plane(const point&, const vector&)") FatalErrorIn("plane::plane(const point&, const vector&)")
<< "plane normal has got zero length" << "plane normal has zero length"
<< abort(FatalError); << abort(FatalError);
} }
} }
@ -238,7 +251,7 @@ Foam::plane::plane(Istream& is)
else else
{ {
FatalErrorIn("plane::plane(Istream& is)") FatalErrorIn("plane::plane(Istream& is)")
<< "plane normal has got zero length" << "plane normal has zero length"
<< abort(FatalError); << abort(FatalError);
} }
} }