Commit 9a70121a authored by Philippose's avatar Philippose Committed by Tomislav Maric
Browse files

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

parent 2bad528b
......@@ -52,7 +52,7 @@ void quadricFitting::calculateNormalVector()
//- estimate the normal as the eigenvector associated
//- to the smallest eigenvalue
normal_ = eigenVector(mat, ev[0]);
normal_ = eigenVectors(mat, ev).x();
}
void quadricFitting::calculateCoordinateSystem()
......
......@@ -495,16 +495,16 @@ void tetMeshOptimisation::optimiseBoundaryVolumeOptimizer
if( mag(ev[2]) > (mag(ev[1]) + mag(ev[0])) )
{
//- ordinary surface vertex
vector normal = eigenVector(nt, ev[2]);
vector normal = eigenVectors(nt, ev).z();
normal /= (mag(normal)+VSMALL);
disp -= (disp & normal) * normal;
}
else if( mag(ev[1]) > 0.5 * (mag(ev[2]) + mag(ev[0])) )
{
//- this vertex is on an edge
vector normal1 = eigenVector(nt, ev[1]);
vector normal1 = eigenVectors(nt, ev).y();
normal1 /= (mag(normal1)+VSMALL);
vector normal2 = eigenVector(nt, ev[2]);
vector normal2 = eigenVectors(nt, ev).z();
normal2 /= (mag(normal2)+VSMALL);
vector eVec = normal1 ^ normal2;
......
......@@ -88,8 +88,8 @@ bool triSurface2DCheck::is2DSurface() const
//- eigenVectors spanning the plane
const vector n
(
eigenVector(covarianceMatrix_, eigenVal[1]) ^
eigenVector(covarianceMatrix_, eigenVal[2])
eigenVectors(covarianceMatrix_, eigenVal).y() ^
eigenVectors(covarianceMatrix_, eigenVal).z()
);
//- 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