From 3e480df8a7e22a4ee8adff4b1d86a1de8d32b9b7 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 9 Jun 2014 11:09:51 +0100
Subject: [PATCH 1/5] BUG: Cloud: force parallel construction of tetBasePtIs

---
 src/lagrangian/basic/Cloud/Cloud.C | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/src/lagrangian/basic/Cloud/Cloud.C b/src/lagrangian/basic/Cloud/Cloud.C
index 3773de66..8c365799 100644
--- a/src/lagrangian/basic/Cloud/Cloud.C
+++ b/src/lagrangian/basic/Cloud/Cloud.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -413,6 +413,12 @@ void Foam::Cloud<ParticleType>::autoMap
 //    polyMesh_.clearCellTree();
     cellWallFacesPtr_.clear();
 
+    // Ask for the tetBasePtIs to trigger all processors to build
+    // them, otherwise, if some processors have no particles then
+    // there is a comms mismatch.
+    polyMesh_.tetBasePtIs();
+
+
     forAllIter(typename Cloud<ParticleType>, *this, pIter)
     {
         ParticleType& p = pIter();

From f3559174840fecc731f28373e72618a318ebd370 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 9 Jun 2014 11:12:14 +0100
Subject: [PATCH 2/5] STYLE: surfaceCheck: use boundBox non-reduce constructor

---
 applications/utilities/surface/surfaceCheck/surfaceCheck.C | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

diff --git a/applications/utilities/surface/surfaceCheck/surfaceCheck.C b/applications/utilities/surface/surfaceCheck/surfaceCheck.C
index db6b1b02..fa3fc28c 100644
--- a/applications/utilities/surface/surfaceCheck/surfaceCheck.C
+++ b/applications/utilities/surface/surfaceCheck/surfaceCheck.C
@@ -31,12 +31,10 @@ Description
 
 #include "triangle.H"
 #include "triSurface.H"
-#include "triSurfaceTools.H"
 #include "triSurfaceSearch.H"
 #include "argList.H"
 #include "OFstream.H"
 #include "OBJstream.H"
-#include "surfaceIntersection.H"
 #include "SortableList.H"
 #include "PatchTools.H"
 #include "vtkSurfaceWriter.H"
@@ -223,7 +221,7 @@ int main(int argc, char *argv[])
     // write bounding box corners
     if (args.optionFound("blockMesh"))
     {
-        pointField cornerPts(boundBox(surf.points()).points());
+        pointField cornerPts(boundBox(surf.points(), false).points());
 
         Info<<"// blockMeshDict info" << nl
             <<"vertices\n(" << nl;

From 420997e125612db1a5a1e4d49133c9858878c5a6 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 9 Jun 2014 12:51:01 +0100
Subject: [PATCH 3/5] BUG: splitMeshRegions: #1282

---
 .../splitMeshRegions/splitMeshRegions.C          | 16 +++++++++++-----
 1 file changed, 11 insertions(+), 5 deletions(-)

diff --git a/applications/utilities/mesh/manipulation/splitMeshRegions/splitMeshRegions.C b/applications/utilities/mesh/manipulation/splitMeshRegions/splitMeshRegions.C
index 6d7b9fe6..aec7160e 100644
--- a/applications/utilities/mesh/manipulation/splitMeshRegions/splitMeshRegions.C
+++ b/applications/utilities/mesh/manipulation/splitMeshRegions/splitMeshRegions.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -496,10 +496,16 @@ void getInterfaceSizes
             }
             interfaceSizes[nInterfaces] = infoIter();
 
-            Map<label> zoneAndInterface;
-            zoneAndInterface.insert(zoneID, nInterfaces);
-            regionsToInterface.insert(e, zoneAndInterface);
-
+            if (regionsToInterface.found(e))
+            {
+                regionsToInterface[e].insert(zoneID, nInterfaces);
+            }
+            else
+            {
+                Map<label> zoneAndInterface;
+                zoneAndInterface.insert(zoneID, nInterfaces);
+                regionsToInterface.insert(e, zoneAndInterface);
+            }
             nInterfaces++;
         }
     }

From c70020edaebc0d66ebea82b55ca4786d4c05f49c Mon Sep 17 00:00:00 2001
From: william <william>
Date: Wed, 11 Jun 2014 11:22:33 +0100
Subject: [PATCH 4/5] BUG: Mantis 1325: Johnson Jackson BCs evaluating before
 dependent sub model data is created

---
 ...JohnsonJacksonParticleSlipFvPatchVectorField.C | 15 ++++-----------
 ...ohnsonJacksonParticleThetaFvPatchScalarField.C | 15 ++++-----------
 2 files changed, 8 insertions(+), 22 deletions(-)

diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C
index 104201d4..78b6ab7e 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C
@@ -102,17 +102,10 @@ JohnsonJacksonParticleSlipFvPatchVectorField
             << abort(FatalError);
     }
 
-    if (dict.found("value"))
-    {
-        fvPatchVectorField::operator=
-        (
-            vectorField("value", dict, p.size())
-        );
-    }
-    else
-    {
-        partialSlipFvPatchVectorField::evaluate();
-    }
+    fvPatchVectorField::operator=
+    (
+        vectorField("value", dict, p.size())
+    );
 }
 
 
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C
index 4369c4ad..1bea5199 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C
@@ -129,17 +129,10 @@ JohnsonJacksonParticleThetaFvPatchScalarField
             << abort(FatalError);
     }
 
-    if (dict.found("value"))
-    {
-        fvPatchScalarField::operator=
-        (
-            scalarField("value", dict, p.size())
-        );
-    }
-    else
-    {
-        evaluate();
-    }
+    fvPatchScalarField::operator=
+    (
+        scalarField("value", dict, p.size())
+    );
 }
 
 

From ddafa9836ad9fdbb12c9ee18916f55d419596a24 Mon Sep 17 00:00:00 2001
From: william <william>
Date: Wed, 11 Jun 2014 11:31:12 +0100
Subject: [PATCH 5/5] BUG: Mantis 1325: laplacian name in kinetic theory model
 incorrectly formatted

---
 .../kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
index f158b243..4fdc9a5f 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
@@ -466,7 +466,7 @@ void Foam::RASModels::kineticTheoryModel::correct()
               + fvm::div(alphaRhoPhi, Theta_)
               - fvc::Sp(fvc::ddt(alpha, rho) + fvc::div(alphaRhoPhi), Theta_)
             )
-          - fvm::laplacian(kappa_, Theta_, "laplacian(kappa, Theta)")
+          - fvm::laplacian(kappa_, Theta_, "laplacian(kappa,Theta)")
          ==
             fvm::SuSp(-((PsCoeff*I) && gradU), Theta_)
           + (tau && gradU)