From b6c9f0c577a49609c1b3558d4ae6cc3a130bf549 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Wed, 15 May 2019 13:02:15 +0100
Subject: [PATCH] Tow-stage Immersed face tolerance consistency.  QUick reject
 using shortest face in edge and detailed analysis edge-by-edge. Inno Gatin

---
 .../immersedPoly/ImmersedFace.C               | 19 ++++++++++---------
 1 file changed, 10 insertions(+), 9 deletions(-)

diff --git a/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C b/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C
index 0497c6e3b..d8e922101 100644
--- a/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C
+++ b/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C
@@ -257,10 +257,11 @@ void Foam::ImmersedFace<Distance>::createSubfaces
         // than 3 points
         if (nDry < 3)
         {
-            FatalErrorInFunction
-                << "There are fewer than three points"
-                << " on the wet part of the face."
-                << abort(FatalError);
+            // The face is wet
+            isAllWet_ = true;
+            isAllDry_ = false;
+
+            drySubface_.clear();
         }
         else
         {
@@ -272,10 +273,11 @@ void Foam::ImmersedFace<Distance>::createSubfaces
 
         if (nWet < 3)
         {
-            FatalErrorInFunction
-                << "There are fewer than three points"
-                << " on the dry part of the face."
-                << abort(FatalError);
+            // The face is dry
+            isAllWet_ = false;
+            isAllDry_ = true;
+
+            wetSubface_.clear();
         }
         else
         {
@@ -342,7 +344,6 @@ Foam::ImmersedFace<Distance>::ImmersedFace
         absTol = minEdgeLength*immersedPoly::tolerance_();
     }
 
-
     // Check if all points are wet or dry, using absolute tolerance
     if (max(depth) < absTol)
     {