diff --git a/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H b/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H
index 58d0d8d724799388d7f67f2d0df7ef1931adba52..e2dbf5f7675d5ce57aa29f22e789df580e067b78 100644
--- a/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H
+++ b/meshLibrary/utilities/helperClasses/geometry/quadricFitting/quadricFittingI.H
@@ -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()
diff --git a/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C b/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C
index 2ffddfde2ff6da8d2f6d0de33807f363704b5147..b9e214b9a8874f189b28eba89cb9983a3fe7aefb 100644
--- a/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C
+++ b/meshLibrary/utilities/smoothers/geometry/meshOptimizer/tetMeshOptimisation/tetMeshOptimisation.C
@@ -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;
diff --git a/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C b/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C
index 9ff1b61d619e591ed4d4c54274b952c881e16fd4..00c31ae9341aff8313cf4281a701125f9558ba1e 100644
--- a/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C
+++ b/meshLibrary/utilities/triSurfaceTools/triSurface2DCheck/triSurface2DCheck.C
@@ -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