From a8727486a4ff628e0b31d0b1741cb6601a18544d Mon Sep 17 00:00:00 2001 From: graham Date: Mon, 20 Dec 2010 11:43:21 +0000 Subject: [PATCH] ENH: Do not modify the face centres when overriding cell centres. --- .../primitiveMeshCellCentresAndVols.C | 34 +++++++++---------- 1 file changed, 17 insertions(+), 17 deletions(-) diff --git a/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C b/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C index 016ba98a13..499deee588 100644 --- a/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C +++ b/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C @@ -206,28 +206,28 @@ void Foam::primitiveMesh::overrideCellCentres // Set internal face centres to the midpoint of the cell-centre delta vector - if (debug) - { - Pout<< "void Foam::primitiveMesh::overrideCellCentres" - << "(const vectorField& newCellCtrs) const : " - << "overriding internal face centres." << endl; - } + // if (debug) + // { + // Pout<< "void Foam::primitiveMesh::overrideCellCentres" + // << "(const vectorField& newCellCtrs) const : " + // << "overriding internal face centres." << endl; + // } - deleteDemandDrivenData(faceCentresPtr_); - deleteDemandDrivenData(faceAreasPtr_); + // deleteDemandDrivenData(faceCentresPtr_); + // deleteDemandDrivenData(faceAreasPtr_); - calcFaceCentresAndAreas(); + // calcFaceCentresAndAreas(); - vectorField& fCtrs = *faceCentresPtr_; + // vectorField& fCtrs = *faceCentresPtr_; - const vectorField& C = cellCentres(); - const labelUList& owner = faceOwner(); - const labelUList& neighbour = faceNeighbour(); + // const vectorField& C = cellCentres(); + // const labelUList& owner = faceOwner(); + // const labelUList& neighbour = faceNeighbour(); - forAll(neighbour, faceI) - { - fCtrs[faceI] = 0.5*(C[neighbour[faceI]] + C[owner[faceI]]); - } + // forAll(neighbour, faceI) + // { + // fCtrs[faceI] = 0.5*(C[neighbour[faceI]] + C[owner[faceI]]); + // } }