LeastSquaresGrad: Corrected handling of boundary values in flattening

src/finiteVolume/finiteVolume/gradSchemes/LeastSquaresGrad/LeastSquaresGrad.C

@@ -101,7 +101,11 @@ Foam::fv::LeastSquaresGrad<Type, Stencil>::calcGrad
     forAll(vtf.boundaryField(), patchi)
     {
         const fvPatchField<Type>& ptf = vtf.boundaryField()[patchi];
-        label nCompact = ptf.patch().start();
+
+        label nCompact =
+            ptf.patch().start()
+          - mesh.nInternalFaces()
+          + mesh.nCells();
 
         forAll(ptf, i)
         {

src/finiteVolume/finiteVolume/gradSchemes/LeastSquaresGrad/LeastSquaresGrads.C

@@ -24,6 +24,7 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "LeastSquaresGrad.H"
+//#include "centredCFCCellToCellStencilObject.H"
 #include "centredCPCCellToCellStencilObject.H"
 #include "centredCECCellToCellStencilObject.H"
 
@@ -33,6 +34,12 @@ namespace Foam
 {
 namespace fv
 {
+    // makeLeastSquaresGradScheme
+    // (
+    //     faceCellsLeastSquares,
+    //     centredCFCCellToCellStencilObject
+    // )
+
     makeLeastSquaresGradScheme
     (
         pointCellsLeastSquares,