diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinement.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinement.C
index 2bf1272938db8b1fd5977ca5a1b32fa5250de4bc..cb31a5c87882973e70a285ee8aaa7e96022d9eff 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinement.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinement.C
@@ -2754,9 +2754,9 @@ const
         label nMasterFaces = 0;
         forAll(isMasterFace, i)
         {
-            if (isMasterFace[i])
+            if (isMasterFace.test(i))
             {
-                nMasterFaces++;
+                ++nMasterFaces;
             }
         }
 
@@ -2764,9 +2764,9 @@ const
         label nMasterPoints = 0;
         forAll(isMeshMasterPoint, i)
         {
-            if (isMeshMasterPoint[i])
+            if (isMeshMasterPoint.test(i))
             {
-                nMasterPoints++;
+                ++nMasterPoints;
             }
         }
 
diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
index 34ae55b76e2cef086647935208568fac9e540be3..03279701c4c4f81174b40280a2c81398257a6cb6 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
@@ -2333,20 +2333,20 @@ void Foam::meshRefinement::getIntersections
             )
             {
                 namedSurfaceIndex[faceI] = surface2[i];
-                posOrientation[faceI] = ((area&normal2[i]) > 0);
+                posOrientation.set(faceI, ((area&normal2[i]) > 0));
                 nSurfFaces[surface2[i]]++;
             }
             else
             {
                 namedSurfaceIndex[faceI] = surface1[i];
-                posOrientation[faceI] = ((area&normal1[i]) > 0);
+                posOrientation.set(faceI, ((area&normal1[i]) > 0));
                 nSurfFaces[surface1[i]]++;
             }
         }
         else if (surface2[i] != -1)
         {
             namedSurfaceIndex[faceI] = surface2[i];
-            posOrientation[faceI] = ((area&normal2[i]) > 0);
+            posOrientation.set(faceI, ((area&normal2[i]) > 0));
             nSurfFaces[surface2[i]]++;
         }
     }
diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
index 9a6fe089dfa611a29d4a7a14179a3d1aaa2ab141..714aad7bc93f768f8eb51bab853c37c278bc05ab 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
@@ -1025,9 +1025,9 @@ Foam::labelList Foam::meshRefinement::markFacesOnProblemCells
                         else
                         {
                             facePatch[facei] = nearestAdaptPatch[facei];
-                            if (isMasterFace[facei])
+                            if (isMasterFace.test(facei))
                             {
-                                nBaffleFaces++;
+                                ++nBaffleFaces;
                             }
 
                             // Do NOT update boundary data since this would grow
diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementRefine.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementRefine.C
index 9954987315aa98aa3bab3b1563f3dc2b825e0c78..7643ab1822c93e07b363927b6698c29fb400339f 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementRefine.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementRefine.C
@@ -261,7 +261,7 @@ Foam::labelList Foam::meshRefinement::getChangedFaces
 
         forAll(changedFace, faceI)
         {
-            if (changedFace[faceI] && isMasterFace[faceI])
+            if (changedFace[faceI] && isMasterFace.test(faceI))
             {
                 nMasterChanged++;
             }
diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
index 989426026588dd1aba2854a332971481d86ff091..e109e96d304e5e92b9e4302c1f5984e4272b9b45 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappyLayerDriver.C
@@ -634,7 +634,7 @@ void Foam::snappyLayerDriver::handleNonManifolds
             (
                 pp.edgeFaces()[edgei].size() == 1
              && edgeGlobalFaces[edgei].size() == 1
-             && isCoupledEdge[meshEdgei]
+             && isCoupledEdge.test(meshEdgei)
             )
             {
                 // Edge of patch but no continuation across processor.
@@ -4336,7 +4336,7 @@ void Foam::snappyLayerDriver::addLayers
             )
             {
                 label oldFacei = map.faceMap()[facei];
-                if (oldFacei != -1 && oldBaffleFace[oldFacei])
+                if (oldFacei != -1 && oldBaffleFace.test(oldFacei))
                 {
                     const face& f = mesh.faces()[facei];
                     forAll(f, fp)
@@ -4680,7 +4680,7 @@ void Foam::snappyLayerDriver::doLayers
                 facei++
             )
             {
-                if (intOrCoupled[facei] && isExtrudedZoneFace[facei])
+                if (intOrCoupled[facei] && isExtrudedZoneFace.test(facei))
                 {
                     faceZoneOnCoupledFace = true;
                     break;
diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C
index c7a165672a948ce8fab64d8ec16cf144569b3973..e25b4185fa67c570635c3bfb157e2ef8f46353af 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriver.C
@@ -183,7 +183,7 @@ Foam::tmp<Foam::pointField> Foam::snappySnapDriver::smoothInternalDisplacement
         label ownLevel = cellLevel[mesh.faceOwner()[facei]];
         label neiLevel = cellLevel[mesh.faceNeighbour()[facei]];
 
-        if (!isFront[facei] && ownLevel != neiLevel)
+        if (!isFront.test(facei) && ownLevel != neiLevel)
         {
             const face& f = mesh.faces()[facei];
             isMovingPoint.set(f);
@@ -200,7 +200,7 @@ Foam::tmp<Foam::pointField> Foam::snappySnapDriver::smoothInternalDisplacement
         label ownLevel = cellLevel[mesh.faceOwner()[facei]];
         label neiLevel = neiCellLevel[facei-mesh.nInternalFaces()];
 
-        if (!isFront[facei] && ownLevel != neiLevel)
+        if (!isFront.test(facei) && ownLevel != neiLevel)
         {
             const face& f = mesh.faces()[facei];
             isMovingPoint.set(f);
@@ -239,7 +239,7 @@ Foam::tmp<Foam::pointField> Foam::snappySnapDriver::smoothInternalDisplacement
 
     forAll(isMovingPoint, pointi)
     {
-        if (isMovingPoint[pointi])
+        if (isMovingPoint.test(pointi))
         {
             const labelList& pCells = mesh.pointCells(pointi);
 
diff --git a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C
index b4c2dc996a521e66e490865723ab12a3afb8b33a..e7f391fdc4ede11800565a824c1f0ef3f0b96959 100644
--- a/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C
+++ b/src/mesh/snappyHexMesh/snappyHexMeshDriver/snappySnapDriverFeature.C
@@ -271,7 +271,7 @@ void Foam::snappySnapDriver::calcNearestFace
                 << exit(FatalError);
         }
         const faceZone& fZone = mesh.faceZones()[zonei];
-        bitSet isZonedFace(mesh.nFaces(), fZone);
+        const bitSet isZonedFace(mesh.nFaces(), fZone);
 
         DynamicList<label> ppFaces(fZone.size());
         DynamicList<label> meshFaces(fZone.size());
@@ -451,7 +451,7 @@ void Foam::snappySnapDriver::calcNearestFacePointProperties
     List<List<point>>& pointFaceSurfNormals,
     List<List<point>>& pointFaceDisp,
     List<List<point>>& pointFaceCentres,
-    List<labelList>&    pointFacePatchID
+    List<labelList>&  pointFacePatchID
 ) const
 {
     const fvMesh& mesh = meshRefiner_.mesh();
@@ -603,7 +603,7 @@ void Foam::snappySnapDriver::calcNearestFacePointProperties
                 {
                     label pointi = meshToPatchPoint[f[fp]];
 
-                    if (pointi != -1 && isBoundaryPoint[pointi])
+                    if (pointi != -1 && isBoundaryPoint.test(pointi))
                     {
                         List<point>& pNormals = pointFaceSurfNormals[pointi];
                         List<point>& pDisp = pointFaceDisp[pointi];
@@ -1440,13 +1440,13 @@ void Foam::snappySnapDriver::releasePointsNextToMultiPatch
             pointFacePatchID[pointi],
             pointFaceCentres[pointi]
         );
-        isMultiPatchPoint[pointi] = multiPatchPt.hit();
+        isMultiPatchPoint.set(pointi, multiPatchPt.hit());
     }
 
     // 2. Make sure multi-patch points are also attracted
     forAll(isMultiPatchPoint, pointi)
     {
-        if (isMultiPatchPoint[pointi])
+        if (isMultiPatchPoint.test(pointi))
         {
             if
             (
@@ -1484,11 +1484,11 @@ void Foam::snappySnapDriver::releasePointsNextToMultiPatch
             label pointi = f[fp];
             if
             (
-                isMultiPatchPoint[pointi]
+                isMultiPatchPoint.test(pointi)
              && patchConstraints[pointi].first() > 1
             )
             {
-                nMultiPatchPoints++;
+                ++nMultiPatchPoints;
             }
         }
 
@@ -1499,7 +1499,7 @@ void Foam::snappySnapDriver::releasePointsNextToMultiPatch
                 label pointi = f[fp];
                 if
                 (
-                   !isMultiPatchPoint[pointi]
+                   !isMultiPatchPoint.test(pointi)
                  && patchConstraints[pointi].first() > 1
                 )
                 {
@@ -2870,7 +2870,7 @@ void Foam::snappySnapDriver::determineBaffleFeatures
         if (efn.size() == 2 && (efn[0]&efn[1]) < baffleFeatureCos)
         {
             isBaffleEdge.set(edgei);
-            nBaffleEdges++;
+            ++nBaffleEdges;
             const edge& e = pp.edges()[edgei];
             pointStatus[e[0]] = 0;
             pointStatus[e[1]] = 0;
diff --git a/src/mesh/snappyHexMesh/trackedParticle/trackedParticle.H b/src/mesh/snappyHexMesh/trackedParticle/trackedParticle.H
index 2a30d95915695df7e0a3ea9137d9c3ca6223c93a..59db354d5f65936f5d53458fce4666988e0f3668 100644
--- a/src/mesh/snappyHexMesh/trackedParticle/trackedParticle.H
+++ b/src/mesh/snappyHexMesh/trackedParticle/trackedParticle.H
@@ -36,6 +36,7 @@ SourceFiles
 #ifndef trackedParticle_H
 #define trackedParticle_H
 
+#include "bitSet.H"
 #include "particle.H"
 #include "autoPtr.H"
 
@@ -44,11 +45,9 @@ SourceFiles
 namespace Foam
 {
 
-class trackedParticleCloud;
-
-
-// Forward declaration of friend functions and operators
+// Forward declarations
 
+class trackedParticleCloud;
 class trackedParticle;
 
 Ostream& operator<<(Ostream&, const trackedParticle&);