Merge branch 'olesenm'

applications/utilities/mesh/conversion/polyDualMesh/Make/files

@@ -1,4 +1,3 @@
-polyDualMesh.C
 polyDualMeshApp.C
 
 EXE = $(FOAM_APPBIN)/polyDualMesh

applications/utilities/mesh/conversion/polyDualMesh/Make/options

@@ -1,5 +1,6 @@
 EXE_INC = \
-    -I$(LIB_SRC)/meshTools/lnInclude
+    -I$(LIB_SRC)/meshTools/lnInclude \
+    -I$(LIB_SRC)/conversion/lnInclude
 
 EXE_LIBS = \
-    -lmeshTools
+    -lmeshTools -lconversion

applications/utilities/parallelProcessing/decomposePar/decomposePar.C

@@ -50,6 +50,14 @@ Usage
     Remove any existing @a processor subdirectories before decomposing the
     geometry.
 
+    @param -lazy \n
+    Only decompose the geometry if the number of domains has changed from a
+    previous decomposition. No @a processor subdirectories will be removed
+    unless the @a -force option is also specified. This option can be used
+    to avoid redundant geometry decomposition (eg, in scripts), but should
+    be used with caution when the underlying (serial) geometry or the
+    decomposition method etc. have been changed between decompositions.
+
 \*---------------------------------------------------------------------------*/
 
 #include "OSspecific.H"
@@ -80,6 +88,7 @@ int main(int argc, char *argv[])
     argList::validOptions.insert("fields", "");
     argList::validOptions.insert("filterPatches", "");
     argList::validOptions.insert("force", "");
+    argList::validOptions.insert("lazy", "");
 
 #   include "setRootCase.H"
 
@@ -88,6 +97,7 @@ int main(int argc, char *argv[])
     bool decomposeFieldsOnly(args.options().found("fields"));
     bool filterPatches(args.options().found("filterPatches"));
     bool forceOverwrite(args.options().found("force"));
+    bool lazyDecomposition(args.options().found("lazy"));
 
 #   include "createTime.H"
 
@@ -100,27 +110,61 @@ int main(int argc, char *argv[])
         ++nProcs;
     }
 
-    // Check for previously decomposed case first
+    // get requested numberOfSubdomains
+    label nDomains = 0;
+    {
+        IOdictionary decompDict
+        (
+            IOobject
+            (
+                "decomposeParDict",
+                runTime.time().system(),
+                runTime,
+                IOobject::MUST_READ,
+                IOobject::NO_WRITE,
+                false
+            )
+        );
+
+        decompDict.lookup("numberOfSubdomains") >> nDomains;
+    }
+
     if (decomposeFieldsOnly)
     {
-        if (!nProcs)
+        // Sanity check on previously decomposed case
+        if (nProcs != nDomains)
         {
             FatalErrorIn(args.executable())
-                << "Specifying -fields requires a decomposed geometry!"
+                << "Specified -fields, but the case was decomposed with "
+                << nProcs << " domains"
+                << nl
+                << "instead of " << nDomains
+                << " domains as specified in decomposeParDict"
                 << nl
                 << exit(FatalError);
         }
     }
-    else
+    else if (nProcs)
     {
-        if (nProcs)
+        bool procDirsProblem = true;
+
+        if (lazyDecomposition && nProcs == nDomains)
         {
+            // we can reuse the decomposition
+            decomposeFieldsOnly = true;
+            procDirsProblem = false;
+            forceOverwrite = false;
+
+            Info<< "Using existing processor directories" << nl;
+        }
+
         if (forceOverwrite)
         {
             Info<< "Removing " << nProcs
                 << " existing processor directories" << endl;
 
             // remove existing processor dirs
+            // reverse order to avoid gaps if someone interrupts the process
             for (label procI = nProcs-1; procI >= 0; --procI)
             {
                 fileName procDir
@@ -130,19 +174,22 @@ int main(int argc, char *argv[])
 
                 rmDir(procDir);
             }
+
+            procDirsProblem = false;
         }
-            else
+
+        if (procDirsProblem)
         {
             FatalErrorIn(args.executable())
-                    << "Case is already decomposed, "
+                << "Case is already decomposed with " << nProcs
-                        "use the -force option or manually remove" << nl
+                << " domains, use the -force option or manually" << nl
-                    << "processor directories before decomposing. e.g.," << nl
+                << "remove processor directories before decomposing. e.g.,"
+                << nl
                 << "    rm -rf " << runTime.path().c_str() << "/processor*"
                 << nl
                 << exit(FatalError);
         }
     }
-    }
 
     Info<< "Create mesh" << endl;
     domainDecomposition mesh

applications/utilities/parallelProcessing/decomposePar/decomposeParDict

@@ -1,33 +1,31 @@
-// Mesh decomposition control dictionary
+/*-------------------------------*- C++ -*---------------------------------*\
-
+|    =========                                                              |
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+|    \\      /     OpenFOAM                                                 |
-
+|     \\    /                                                               |
+|      \\  /       The Open Source CFD Toolbox                              |
+|       \\/                                        http://www.OpenFOAM.org  |
+\*-------------------------------------------------------------------------*/
 FoamFile
 {
-    version 0.5;
+    version     2.0;
     format      ascii;
-
-    root     "ROOT";
-    case     "CASE";
-    instance "system";
-    local    "";
-
     class       dictionary;
-
+    note        "mesh decomposition control dictionary";
-    object decompositionDict;
+    location    "system";
+    object      decomposeParDict;
 }
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 numberOfSubdomains  4;
 
-//preservePatches (inlet);
+// preservePatches (inlet);
-//preserveFaceZones (heater solid1 solid3);
+// preserveFaceZones (heater solid1 solid3);
 
 method          simple;
-//method               hierarchical;
+// method          hierarchical;
-//method               metis;
+// method          metis;
-//method               manual;
+// method          manual;
 
 simpleCoeffs
 {
@@ -44,13 +42,15 @@ hierarchicalCoeffs
 
 metisCoeffs
 {
-    //processorWeights 
+ /*
-    //(
+    processorWeights
-    //    1
+    (
-    //    1
+        1
-    //    1
+        1
-    //    1
+        1
-    //);
+        1
+    );
+  */
 }
 
 manualCoeffs

src/Allwmake

@@ -31,6 +31,7 @@ wmake libso randomProcesses
 ( cd turbulenceModels && ./Allwmake )
 ( cd lagrangian && ./Allwmake )
 ( cd postProcessing && ./Allwmake )
+( cd conversion && ./Allwmake )
 
 wmake libso autoMesh
 wmake libso errorEstimation

src/conversion/Allwmake

@@ -0,0 +1,4 @@
+#!/bin/sh
+set -x
+
+wmake libso

src/conversion/Make/files

@@ -0,0 +1,3 @@
+polyDualMesh/polyDualMesh.C
+
+LIB = $(FOAM_LIBBIN)/libconversion

src/conversion/Make/options

@@ -0,0 +1,5 @@
+EXE_INC = \
+    -I$(LIB_SRC)/meshTools/lnInclude
+
+LIB_LIBS = \
+    -lmeshTools