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: cellMapper: map volume interpolative

Browse Source
This commit is contained in:
mattijs
2013-09-18 17:02:27 +01:00
parent a501d0ff3c
commit c5239e17eb
1 changed files with 34 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

49
src/OpenFOAM/meshes/polyMesh/mapPolyMesh/cellMapper/cellMapper.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -28,9 +28,6 @@ License
#include "polyMesh.H"
#include "mapPolyMesh.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
void Foam::cellMapper::calcAddressing() const
@ -154,6 +151,18 @@ void Foam::cellMapper::calcAddressing() const
const List<objectMap>& cfc = mpm_.cellsFromCellsMap();
const scalarField& V = mesh_.cellVolumes();
if (V.size() != sizeBeforeMapping())
{
FatalErrorIn("void cellMapper::calcAddressing() const")
<< "cellVolumes size " << V.size()
<< " is not the old number of cells " << sizeBeforeMapping()
<< ". Are your cellVolumes already mapped?"
<< " (new number of cells " << size() << ")"
<< abort(FatalError);
}
forAll(cfc, cfcI)
{
// Get addressing
@ -169,9 +178,28 @@ void Foam::cellMapper::calcAddressing() const
<< " already destination of mapping." << abort(FatalError);
}
// Map from masters, uniform weights
// Map from masters
addr[cellI] = mo;
w[cellI] = scalarList(mo.size(), 1.0/mo.size());
//- uniform weights
//w[cellI] = scalarList(mo.size(), 1.0/mo.size());
//- volume based
w[cellI].setSize(mo.size());
if (mo.size())
{
scalar sumV = 0;
forAll(mo, ci)
{
w[cellI][ci] = V[mo[ci]];
sumV += V[mo[ci]];
}
forAll(mo, ci)
{
w[cellI][ci] /= sumV;
}
}
}
@ -413,13 +441,4 @@ const Foam::labelList& Foam::cellMapper::insertedObjectLabels() const
}
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
// * * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * //
// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
// ************************************************************************* //