Merge branch 'nextRelease' of ssh://git.code.sf.net/p/foam-extend/foam-extend-4.0 into nextRelease

src/finiteVolume/fields/fvPatchFields/derived/inletOutletTotalTemperature/inletOutletTotalTemperatureFvPatchScalarField.C

@@ -53,24 +53,6 @@ inletOutletTotalTemperatureFvPatchScalarField
 }
 
 
-Foam::inletOutletTotalTemperatureFvPatchScalarField::
-inletOutletTotalTemperatureFvPatchScalarField
-(
-    const inletOutletTotalTemperatureFvPatchScalarField& ptf,
-    const fvPatch& p,
-    const DimensionedField<scalar, volMesh>& iF,
-    const fvPatchFieldMapper& mapper
-)
-:
-    mixedFvPatchScalarField(ptf, p, iF, mapper),
-    UName_(ptf.UName_),
-    phiName_(ptf.phiName_),
-    psiName_(ptf.psiName_),
-    gamma_(ptf.gamma_),
-    T0_(ptf.T0_, mapper)
-{}
-
-
 Foam::inletOutletTotalTemperatureFvPatchScalarField::
 inletOutletTotalTemperatureFvPatchScalarField
 (
@@ -104,6 +86,24 @@ inletOutletTotalTemperatureFvPatchScalarField
 }
 
 
+Foam::inletOutletTotalTemperatureFvPatchScalarField::
+inletOutletTotalTemperatureFvPatchScalarField
+(
+    const inletOutletTotalTemperatureFvPatchScalarField& ptf,
+    const fvPatch& p,
+    const DimensionedField<scalar, volMesh>& iF,
+    const fvPatchFieldMapper& mapper
+)
+:
+    mixedFvPatchScalarField(ptf, p, iF, mapper),
+    UName_(ptf.UName_),
+    phiName_(ptf.phiName_),
+    psiName_(ptf.psiName_),
+    gamma_(ptf.gamma_),
+    T0_(ptf.T0_, mapper)
+{}
+
+
 Foam::inletOutletTotalTemperatureFvPatchScalarField::
 inletOutletTotalTemperatureFvPatchScalarField
 (

src/foam/meshes/polyMesh/polyMesh.C

@@ -227,6 +227,7 @@ Foam::polyMesh::polyMesh(const IOobject& io)
             IOobject::NO_WRITE
         )
     ),
+    syncPar_(true),  // Reading mesh from IOobject: must be valid
     clearedPrimitives_(false),
     boundary_
     (
@@ -241,7 +242,7 @@ Foam::polyMesh::polyMesh(const IOobject& io)
         ),
         *this
     ),
-    bounds_(allPoints_),
+    bounds_(allPoints_),  // Reading mesh from IOobject: syncPar
     geometricD_(Vector<label>::zero),
     solutionD_(Vector<label>::zero),
     comm_(Pstream::worldComm),
@@ -319,7 +320,6 @@ Foam::polyMesh::polyMesh(const IOobject& io)
                 "cells",
                 // Find the cells file on the basis of the faces file
                 // HJ, 8/Jul/2009
-//                 time().findInstance(meshDir(), "cells"),
                 time().findInstance(meshDir(), "faces"),
                 meshSubDir,
                 *this,
@@ -424,6 +424,7 @@ Foam::polyMesh::polyMesh
         ),
         neighbour
     ),
+    syncPar_(syncPar),
     clearedPrimitives_(false),
     boundary_
     (
@@ -585,6 +586,7 @@ Foam::polyMesh::polyMesh
         ),
         0
     ),
+    syncPar_(syncPar),
     clearedPrimitives_(false),
     boundary_
     (
@@ -1177,6 +1179,14 @@ Foam::tmp<Foam::scalarField> Foam::polyMesh::movePoints
     const pointField& newPoints
 )
 {
+    if (!syncPar_)
+    {
+        Info<< "tmp<scalarField> polyMesh::movePoints"
+            << "(const pointField&) : "
+            << "Moving the mesh for the mesh with syncPar invalidated.  "
+            << "Mesh motion should not be executed." << endl;
+    }
+
     if (debug)
     {
         Info<< "tmp<scalarField> polyMesh::movePoints(const pointField&) : "

src/foam/meshes/polyMesh/polyMesh.H

@@ -117,6 +117,9 @@ private:
             //- Face neighbour
             labelIOList neighbour_;
 
+            //- Is the mesh in a parallel sync state?
+            Switch syncPar_;
+
             //- Have the primitives been cleared
             bool clearedPrimitives_;
 
@@ -374,12 +377,27 @@ public:
                 return boundary_;
             }
 
+
+        // Mesh state and form
+
             //- Return mesh bounding box
             const boundBox& bounds() const
             {
                 return bounds_;
             }
 
+            //- Return parallel sync state
+            const Switch& syncPar() const
+            {
+                return syncPar_;
+            }
+
+            //- Return access to parallel sync state
+            Switch& syncPar()
+            {
+                return syncPar_;
+            }
+
             //- Return the vector of geometric directions in mesh.
             //  Defined according to the presence of empty and wedge patches.
             //  1 indicates unconstrained direction and -1 a constrained

src/foam/meshes/polyMesh/polyMeshFromShapeMesh.C

@@ -489,6 +489,7 @@ Foam::polyMesh::polyMesh
         ),
         0
     ),
+    syncPar_(syncPar),
     clearedPrimitives_(false),
     boundary_
     (
@@ -774,6 +775,7 @@ Foam::polyMesh::polyMesh
         ),
         0
     ),
+    syncPar_(syncPar),
     clearedPrimitives_(false),
     boundary_
     (

src/foam/meshes/polyMesh/polyMeshIO.C

@@ -97,6 +97,7 @@ Foam::polyMesh::polyMesh
         ),
         labelList(is)
     ),
+    syncPar_(syncPar),
     clearedPrimitives_(false),
     boundary_
     (