zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: Do not modify the face centres when overriding cell centres.

Browse Source
This commit is contained in:
graham
2010-12-20 11:43:21 +00:00
parent 664940fa3a
commit a8727486a4
1 changed files with 17 additions and 17 deletions
Download Patch File Download Diff File Expand all files Collapse all files

34
src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C
View File

@ -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]]);
// }
}