BUG: Managing topoBitSet as source in haloToCell. See issue #2137

src/meshTools/sets/cellSources/haloToCell/haloToCell.C

@@ -28,6 +28,7 @@ License
 #include "haloToCell.H"
 #include "polyMesh.H"
 #include "cellSet.H"
+#include "topoBitSet.H"
 #include "syncTools.H"
 #include "addToRunTimeSelectionTable.H"
 
@@ -78,15 +79,21 @@ void Foam::haloToCell::combine(topoSet& set, const bool add) const
     const labelList& faceOwn = mesh_.faceOwner();
     const labelList& faceNei = mesh_.faceNeighbour();
 
-
     // The starting set of cells
     bitSet current(cells.size());
 
-    for (const label celli : set)
+    if (isA<topoBitSet>(set))
     {
-        current.set(celli);
+        current |= refCast<const topoBitSet>(set).addressing();
     }
-
+    else 
+    {
+        for (const label celli : set)
+        {
+            current.set(celli);
+        }
+    }
+    
     // The perimeter faces of the cell set
     bitSet outsideFaces(mesh_.nFaces());