diff --git a/applications/utilities/mesh/manipulation/objToVTK/objToVTK.C b/applications/utilities/mesh/manipulation/objToVTK/objToVTK.C index 375e38961b1625b10176931bc955a49055f055cf..ebaf8737b834153883df54981e37c98c8bfc6e2c 100644 --- a/applications/utilities/mesh/manipulation/objToVTK/objToVTK.C +++ b/applications/utilities/mesh/manipulation/objToVTK/objToVTK.C @@ -210,6 +210,14 @@ int main(int argc, char *argv[]) outFile << pt.x() << ' ' << pt.y() << ' ' << pt.z() << nl; } + outFile + << "VERTICES " << points.size() << ' ' << (2 * points.size()) << nl; + + forAll(points, i) + { + outFile << 1 << ' ' << i << nl; + } + label nItems = 0; forAll(polyLines, polyI) { diff --git a/src/regionModels/regionModel/regionModel/regionModel.C b/src/regionModels/regionModel/regionModel/regionModel.C index 5bf2951451580c3dd9a527a59088a88b5e97e3cd..ac72e9ab2477fdab5ac6e641f74c0e364d0b51de 100644 --- a/src/regionModels/regionModel/regionModel/regionModel.C +++ b/src/regionModels/regionModel/regionModel/regionModel.C @@ -103,7 +103,7 @@ void Foam::regionModels::regionModel::initialise() forAll(rbm, patchI) { const polyPatch& regionPatch = rbm[patchI]; - if (isA<mappedWallPolyPatch>(regionPatch)) + if (isA<mappedPatchBase>(regionPatch)) { if (debug) { @@ -136,7 +136,7 @@ void Foam::regionModels::regionModel::initialise() primaryPatchIDs_.transfer(primaryPatchIDs); intCoupledPatchIDs_.transfer(intCoupledPatchIDs); - if (nBoundaryFaces == 0) + if (returnReduce(nBoundaryFaces, sumOp<label>()) == 0) { WarningIn("regionModel::initialise()") << "Region model has no mapped boundary conditions - transfer "