ENH: cvMeshSurfaceSimplify: make robust for open surfaces

applications/utilities/mesh/generation/cvMesh/cvMeshSurfaceSimplify/cvMeshSurfaceSimplify.C

@@ -173,6 +173,8 @@ int main(int argc, char *argv[])
         surfaces,
         points,
         scalarField(points.size(), sqr(GREAT)),
+        searchableSurface::OUTSIDE,     // for non-closed surfaces treat as
+                                        // outside
         nearestSurfaces,
         signedDist
     );

src/meshTools/searchableSurface/searchableSurfacesQueries.C

@@ -693,6 +693,7 @@ void Foam::searchableSurfacesQueries::signedDistance
     const labelList& surfacesToTest,
     const pointField& samples,
     const scalarField& nearestDistSqr,
+    const searchableSurface::volumeType illegalHandling,
     labelList& nearestSurfaces,
     scalarField& distance
 )
@@ -753,9 +754,32 @@ void Foam::searchableSurfacesQueries::signedDistance
             }
             else
             {
-                FatalErrorIn("signedDistance()")
+                switch (illegalHandling)
-                    << "getVolumeType failure, neither INSIDE or OUTSIDE"
+                {
-                    << exit(FatalError);
+                    case searchableSurface::OUTSIDE:
+                    {
+                        distance[pointI] = dist;
+                        break;
+                    }
+                    case searchableSurface::INSIDE:
+                    {
+                        distance[pointI] = -dist;
+                        break;
+                    }
+                    default:
+                    {
+                        FatalErrorIn("signedDistance()")
+                            << "getVolumeType failure,"
+                            << " neither INSIDE or OUTSIDE."
+                            << " point:" << surfPoints[i]
+                            << " surface:"
+                            << allSurfaces[surfacesToTest[testI]].name()
+                            << " volType:"
+                            << searchableSurface::volumeTypeNames[vT]
+                            << exit(FatalError);
+                        break;
+                    }
+                }
             }
         }
     }

src/meshTools/searchableSurface/searchableSurfacesQueries.H

@@ -184,13 +184,18 @@ public:
                 List<pointIndexHit>&
             );
 
-            //- Find signed distance to nearest surface
+            //- Find signed distance to nearest surface. Outside is positive.
+            //  illegalHandling: how to handle non-inside or outside
+            //      OUTSIDE : treat as outside
+            //      INSIDE  : treat as inside
+            //      UNKNOWN : throw fatal error
             static void signedDistance
             (
                 const PtrList<searchableSurface>& allSurfaces,
                 const labelList& surfacesToTest,
                 const pointField& samples,
                 const scalarField& nearestDistSqr,
+                const searchableSurface::volumeType illegalHandling,
                 labelList& nearestSurfaces,
                 scalarField& distance
             );