ENH: Adding iterative solver to S2S equations

ENH: Adding report after matrix smoothing

applications/utilities/preProcessing/viewFactorsGen/shootRays.H

@@ -4,7 +4,7 @@
 // Maximum length for dynamicList
 const label maxDynListLength
 (
-    viewFactorDict.getOrDefault<label>("maxDynListLength", 100000)
+    viewFactorDict.getOrDefault<label>("maxDynListLength", 1000000)
 );
 
 for (const int proci : Pstream::allProcs())
@@ -47,7 +47,7 @@ for (const int proci : Pstream::allProcs())
                     const vector& remA = remoteArea[j];
                     const label& remAgg = remoteAgg[j];
 
-                    const vector& d = remFc - fc;
+                    const vector d(remFc - fc);
 
                     if (((d & fA) < 0.) && ((d & remA) > 0))
                     {

applications/utilities/preProcessing/viewFactorsGen/viewFactorsGen.C

@@ -408,7 +408,7 @@ int main(int argc, char *argv[])
                 List<point> availablePoints
                 (
                     upp.faceCentres().size()
-                + upp.localPoints().size()
+                  + upp.localPoints().size()
                 );
 
                 SubList<point>
@@ -863,21 +863,6 @@ int main(int argc, char *argv[])
     IOglobalFaceFaces.write();
 
 
-    labelListIOList IOvisibleFaceFaces
-    (
-        IOobject
-        (
-            "visibleFaceFaces",
-            mesh.facesInstance(),
-            mesh,
-            IOobject::NO_READ,
-            IOobject::NO_WRITE,
-            false
-        ),
-        std::move(visibleFaceFaces)
-    );
-    IOvisibleFaceFaces.write();
-
     IOmapDistribute IOmapDist
     (
         IOobject