fix

library/NekMeshUtils/Octree/Octree.cpp

@@ -185,9 +185,21 @@ NekDouble Octree::Query(Array<OneD, NekDouble> loc)
             found = true;
         }
     }
+
     return min(n->GetDelta(),tmp);
 }
 
+NekDouble Octree::GetMinDelta()
+{
+    NekDouble tmp = numeric_limits<double>::max();
+
+    for(int i = 0; i < m_lsources.size(); i++)
+    {
+        tmp = min(m_lsources[i].delta,tmp);
+    }
+    return min(m_minDelta,tmp);
+}
+
 void Octree::WriteOctree(string nm)
 {
     MeshSharedPtr oct = boost::shared_ptr<Mesh>(new Mesh());

library/NekMeshUtils/Octree/Octree.h

@@ -132,10 +132,7 @@ public:
      *
      * @return miminum delta in octree
      */
-    NekDouble GetMinDelta()
-    {
-        return m_minDelta;
-    }
+    NekDouble GetMinDelta();
 
     /**
      * @brief sets the parameters used for curvature sampling

library/NekMeshUtils/SurfaceMeshing/FaceMesh.cpp

@@ -918,17 +918,17 @@ bool FaceMesh::Validate()
 
         int numValid = 0;
 
-        if (r[0] < triDelta[0] && r[2] < triDelta[0])
+        if (r[0] < (triDelta[0] + triDelta[1]) / 2.0 *1.41)
         {
             numValid++;
         }
 
-        if (r[1] < triDelta[1] && r[0] < triDelta[1])
+        if (r[1] < (triDelta[1] + triDelta[2]) / 2.0 *1.41)
         {
             numValid++;
         }
 
-        if (r[2] < triDelta[2] && r[1] < triDelta[2])
+        if (r[2] < (triDelta[2] + triDelta[0]) / 2.0 *1.41)
         {
             numValid++;
         }

utilities/NekMesh/Tests/MeshGen/cylinder.mcf

@@ -7,9 +7,9 @@
         </INFORMATION>
 
         <PARAMETERS>
-            <P PARAM="MinDelta" VALUE="0.005"    />
+            <P PARAM="MinDelta" VALUE="0.003"    />
             <P PARAM="MaxDelta" VALUE="0.2"     />
-            <P PARAM="EPS"      VALUE="0.09"     />
+            <P PARAM="EPS"      VALUE="0.05"     />
 
             <P PARAM="Order"    VALUE="4"       />
         </PARAMETERS>