From b84618a3c8cd5d2eaba04254b602f354effcb1ed Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Tue, 18 Oct 2016 13:05:20 +0100
Subject: [PATCH] blockMesh: Added check for curved-edge/block-edge consistency

---
 src/mesh/blockMesh/blockMesh/blockMeshCheck.C | 27 ++++++++++++++++---
 1 file changed, 23 insertions(+), 4 deletions(-)

diff --git a/src/mesh/blockMesh/blockMesh/blockMeshCheck.C b/src/mesh/blockMesh/blockMesh/blockMeshCheck.C
index 3091c16cbd..6f1fa2d906 100644
--- a/src/mesh/blockMesh/blockMesh/blockMeshCheck.C
+++ b/src/mesh/blockMesh/blockMesh/blockMeshCheck.C
@@ -29,13 +29,32 @@ License
 
 void Foam::blockMesh::check(const polyMesh& bm) const
 {
-    if (verboseOutput)
-    {
-        Info<< nl << "Check topology" << endl;
-    }
+    Info<< nl << "Check topology" << endl;
 
     bool ok = true;
 
+    // Check curved-edge/block-edge correspondence
+    const edgeList& edges = bm.edges();
+
+    forAll(edges_, cei)
+    {
+        bool found = false;
+
+        forAll(edges, ci)
+        {
+            found = edges_[cei].compare(edges[ci][0], edges[ci][1]) != 0;
+            if (found) break;
+        }
+
+        if (!found)
+        {
+            Info<< "    Curved edge " << edges_[cei]
+                << "    does not correspond to a block edge."
+                << endl;
+            ok = false;
+        }
+    }
+
     const pointField& points = bm.points();
     const faceList& faces = bm.faces();
     const cellList& cells = bm.cells();