Commit 3426919c authored by Philippose's avatar Philippose

Porting to OpenFOAM-v1606+ ... change in the eigenvector API

parent fbad999d
...@@ -52,7 +52,7 @@ void quadricFitting::calculateNormalVector() ...@@ -52,7 +52,7 @@ void quadricFitting::calculateNormalVector()
//- estimate the normal as the eigenvector associated //- estimate the normal as the eigenvector associated
//- to the smallest eigenvalue //- to the smallest eigenvalue
normal_ = eigenVector(mat, ev[0]); normal_ = eigenVectors(mat, ev).x();
} }
void quadricFitting::calculateCoordinateSystem() void quadricFitting::calculateCoordinateSystem()
......
...@@ -495,16 +495,16 @@ void tetMeshOptimisation::optimiseBoundaryVolumeOptimizer ...@@ -495,16 +495,16 @@ void tetMeshOptimisation::optimiseBoundaryVolumeOptimizer
if( mag(ev[2]) > (mag(ev[1]) + mag(ev[0])) ) if( mag(ev[2]) > (mag(ev[1]) + mag(ev[0])) )
{ {
//- ordinary surface vertex //- ordinary surface vertex
vector normal = eigenVector(nt, ev[2]); vector normal = eigenVectors(nt, ev).z();
normal /= (mag(normal)+VSMALL); normal /= (mag(normal)+VSMALL);
disp -= (disp & normal) * normal; disp -= (disp & normal) * normal;
} }
else if( mag(ev[1]) > 0.5 * (mag(ev[2]) + mag(ev[0])) ) else if( mag(ev[1]) > 0.5 * (mag(ev[2]) + mag(ev[0])) )
{ {
//- this vertex is on an edge //- this vertex is on an edge
vector normal1 = eigenVector(nt, ev[1]); vector normal1 = eigenVectors(nt, ev).y();
normal1 /= (mag(normal1)+VSMALL); normal1 /= (mag(normal1)+VSMALL);
vector normal2 = eigenVector(nt, ev[2]); vector normal2 = eigenVectors(nt, ev).z();
normal2 /= (mag(normal2)+VSMALL); normal2 /= (mag(normal2)+VSMALL);
vector eVec = normal1 ^ normal2; vector eVec = normal1 ^ normal2;
......
...@@ -88,8 +88,8 @@ bool triSurface2DCheck::is2DSurface() const ...@@ -88,8 +88,8 @@ bool triSurface2DCheck::is2DSurface() const
//- eigenVectors spanning the plane //- eigenVectors spanning the plane
const vector n const vector n
( (
eigenVector(covarianceMatrix_, eigenVal[1]) ^ eigenVectors(covarianceMatrix_, eigenVal).y() ^
eigenVector(covarianceMatrix_, eigenVal[2]) eigenVectors(covarianceMatrix_, eigenVal).z()
); );
//- check if the plane is in the x-y plane of the coordinate system //- check if the plane is in the x-y plane of the coordinate system
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment