Clean-up of post-processing with topo changes

applications/utilities/postProcessing/graphics/PV3FoamReader/vtkPV3Foam/vtkPV3FoamFields.C

@@ -192,7 +192,11 @@ void Foam::vtkPV3Foam::convertPointFields
     }
 
     // Construct interpolation on the raw mesh
-    pointMesh pMesh(mesh);
+
+    // HJ, bug fix?  Point mesh handled by objectRegistry
+    // HJ, 11/Nov/2010
+    const pointMesh& pMesh = pointMesh::New(mesh);
+//     pointMesh pMesh(mesh);
 
 
     convertPointFields<scalar>

src/OpenFOAM/meshes/polyMesh/polyMeshClear.C

@@ -29,6 +29,7 @@ License
 #include "globalMeshData.H"
 #include "demandDrivenData.H"
 #include "meshObjectBase.H"
+#include "pointMesh.H"
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
@@ -93,7 +94,9 @@ void Foam::polyMesh::clearGeom()
     solutionD_ = Vector<label>::zero;
 
     // Move points all mesh objects.  HJ, 13/Oct/2010
-    meshObjectBase::allMovePoints(*this);
+    // This is a problem, because it is called in a destructor.
+    meshObjectBase::allDelete(*this);
+//     meshObjectBase::allMovePoints(*this);
 }
 
 

src/OpenFOAM/meshes/polyMesh/polyMeshIO.C

@@ -260,7 +260,6 @@ Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
             initMesh(cells);
         }
 
-
         // Even if number of patches stayed same still recalculate boundary
         // data.
 

src/OpenFOAM/meshes/polyMesh/polyMeshInitMesh.C

@@ -36,9 +36,34 @@ void Foam::polyMesh::initMesh()
             << "initialising primitiveMesh" << endl;
     }
 
+    // For backward compatibility check if the owner array is the same
+    // length as the number of allFaces and shrink to remove the -1s padding
+    if (min(owner_) < 0)
+    {
+        label nActiveFaces = 0;
+
+        forAll(owner_, faceI)
+        {
+            if (owner_[faceI] == -1)
+            {
+                break;
+            }
+            else
+            {
+                nActiveFaces++;
+            }
+        }
+
+        InfoIn("void polyMesh::initMesh()")
+            << "Truncating owner list at " << nActiveFaces
+            << " for backward compatibility" << endl;
+
+        owner_.setSize(nActiveFaces);
+    }
+
     // For backward compatibility check if the neighbour array is the same
     // length as the owner and shrink to remove the -1s padding
-    if (neighbour_.size() == owner_.size())
+    if (min(neighbour_) < 0)
     {
         label nIntFaces = 0;
 
@@ -54,6 +79,10 @@ void Foam::polyMesh::initMesh()
             }
         }
 
+        InfoIn("void polyMesh::initMesh()")
+            << "Truncating neighbour list at " << nIntFaces
+            << " for backward compatibility" << endl;
+
         neighbour_.setSize(nIntFaces);
     }
 

src/finiteVolume/interpolation/volPointInterpolation/volPointInterpolation.H

@@ -86,7 +86,8 @@ public:
 
     // Constructors
 
-        //- Constructor given fvMesh and pointMesh.
+        //- Constructor given fvMesh.  pointMesh will be created or
+        //  looked up from objectRegistry
         explicit volPointInterpolation(const fvMesh&);