From a3c239af848bb09b27b78a8df6e39fd0b9dbfaee Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 28 Feb 2019 15:09:32 +0000
Subject: [PATCH] BUG: snappyHexMesh: deal with non-compact regions. Fixes
 #1221.

---
 .../triSurfaceSearch/triSurfaceRegionSearch.C | 22 ++++++++++---------
 1 file changed, 12 insertions(+), 10 deletions(-)

diff --git a/src/meshTools/triSurface/triSurfaceSearch/triSurfaceRegionSearch.C b/src/meshTools/triSurface/triSurfaceSearch/triSurfaceRegionSearch.C
index 933c3960fa..faeb953d65 100644
--- a/src/meshTools/triSurface/triSurfaceSearch/triSurfaceRegionSearch.C
+++ b/src/meshTools/triSurface/triSurfaceSearch/triSurfaceRegionSearch.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2015 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2015,2019 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
                             | Copyright (C) 2011-2016 OpenFOAM Foundation
@@ -74,15 +74,20 @@ Foam::triSurfaceRegionSearch::treeByRegion() const
 {
     if (treeByRegion_.empty())
     {
-        Map<label> regionSizes;
+        label maxRegion = -1;
         forAll(surface(), fI)
         {
             const label regionI = surface()[fI].region();
-
-            regionSizes(regionI)++;
+            maxRegion = max(maxRegion, regionI);
         }
+        const label nRegions = maxRegion+1;
 
-        label nRegions = regionSizes.size();
+        labelList nFacesInRegions(nRegions, 0);
+        forAll(surface(), fI)
+        {
+            const label regionI = surface()[fI].region();
+            nFacesInRegions[regionI]++;
+        }
 
         indirectRegionPatches_.setSize(nRegions);
         treeByRegion_.setSize(nRegions);
@@ -91,15 +96,12 @@ Foam::triSurfaceRegionSearch::treeByRegion() const
 
         forAll(regionsAddressing, regionI)
         {
-            regionsAddressing[regionI] = labelList(regionSizes[regionI], -1);
+            regionsAddressing[regionI].setSize(nFacesInRegions[regionI]);
         }
-
-        labelList nFacesInRegions(nRegions, Zero);
-
+        nFacesInRegions = Zero;
         forAll(surface(), fI)
         {
             const label regionI = surface()[fI].region();
-
             regionsAddressing[regionI][nFacesInRegions[regionI]++] = fI;
         }