Inverse-distance-weighted interpolation stencil. More...

Inheritance diagram for inverseDistance:

[legend]

Collaboration diagram for inverseDistance:

[legend]

Public Member Functions
	TypeName ("inverseDistance")

	inverseDistance (const fvMesh &, const dictionary &, const bool)

virtual	~inverseDistance ()

virtual bool	update ()

virtual const labelUList &	cellTypes () const

virtual const labelUList &	interpolationCells () const

virtual const mapDistribute &	cellInterpolationMap () const

virtual const labelListList &	cellStencil () const

virtual const scalarListList &	cellInterpolationWeights () const

virtual const scalarList &	cellInterpolationWeight () const

virtual void	stencilWeights (const point &sample, const pointList &donorCcs, scalarList &weights) const

Public Member Functions inherited from cellCellStencil
	TypeName ("cellCellStencil")

	declareRunTimeSelectionTable (autoPtr, cellCellStencil, mesh,(const fvMesh &mesh, const dictionary &dict, const bool update),(mesh, dict, update))

	cellCellStencil (const fvMesh &)

virtual	~cellCellStencil ()

virtual const wordHashSet &	nonInterpolatedFields () const

virtual wordHashSet &	nonInterpolatedFields ()

bool	localStencil (const labelUList &) const

const labelIOList &	zoneID () const

template<class Type >
Foam::tmp< Foam::volScalarField >	createField (const fvMesh &mesh, const word &name, const UList< Type > &psi)

Protected Member Functions
void	markPatchesAsHoles (PstreamBuffers &pBufs, const PtrList< fvMeshSubset > &meshParts, const List< treeBoundBoxList > &patchBb, const List< labelVector > &patchDivisions, const PtrList< PackedList< 2 >> &patchParts, const label srcI, const label tgtI, labelList &allCellTypes) const

bool	betterDonor (const label destMesh, const label currentDonorMesh, const label newDonorMesh) const

void	markDonors (const globalIndex &globalCells, PstreamBuffers &pBufs, const PtrList< fvMeshSubset > &meshParts, const List< treeBoundBoxList > &meshBb, const labelList &allCellTypes, const label srcI, const label tgtI, labelListList &allStencil, labelList &allDonor) const

void	uncompactedRegionSplit (const fvMesh &mesh, const globalIndex &globalFaces, const label nZones, const labelList &zoneID, const labelList &cellTypes, const boolList &isBlockedFace, labelList &cellRegion) const

autoPtr< globalIndex >	compactedRegionSplit (const fvMesh &mesh, const globalIndex &globalFaces, labelList &cellRegion) const

void	findHoles (const globalIndex &globalCells, const fvMesh &mesh, const labelList &zoneID, const labelListList &stencil, labelList &cellTypes) const

void	seedCell (const label cellI, const scalar wantedFraction, bitSet &isFront, scalarField &fraction) const

void	walkFront (const scalar layerRelax, const labelListList &allStencil, labelList &allCellTypes, scalarField &allWeight) const

virtual void	createStencil (const globalIndex &)

Static Protected Member Functions
static label	index (const labelVector &nDivs, const labelVector &)

static labelVector	index3 (const labelVector &nDivs, const label)

static labelVector	index3 (const boundBox &bb, const labelVector &nDivs, const point &pt)

static point	position (const boundBox &bb, const labelVector &nDivs, const label boxI)

static void	fill (PackedList< 2 > &elems, const boundBox &bb, const labelVector &nDivs, const boundBox &subBb, const unsigned int val)

static bool	overlaps (const boundBox &bb, const labelVector &nDivs, const PackedList< 2 > &voxels, const treeBoundBox &subBb, const unsigned int val)

static void	markBoundaries (const fvMesh &mesh, const vector &smallVec, const boundBox &bb, const labelVector &nDivs, PackedList< 2 > &patchTypes, const labelList &cellMap, labelList &patchCellTypes)

static treeBoundBox	cellBb (const primitiveMesh &mesh, const label celli)

Static Protected Member Functions inherited from cellCellStencil
static labelList	count (const label size, const labelUList &lst)

template<class Type >
static tmp< volScalarField >	createField (const fvMesh &mesh, const word &name, const UList< Type > &)

Protected Attributes
const dictionary	dict_

vector	smallVec_

labelList	cellTypes_

labelList	interpolationCells_

autoPtr< mapDistribute >	cellInterpolationMap_

labelListList	cellStencil_

scalarListList	cellInterpolationWeights_

volScalarField	cellInterpolationWeight_

Protected Attributes inherited from cellCellStencil
const fvMesh &	mesh_

wordHashSet	nonInterpolatedFields_

Additional Inherited Members
Public Types inherited from cellCellStencil
enum	patchCellType { OTHER = 0, PATCH = 1, OVERSET = 2 }

enum	cellType { CALCULATED = 0, INTERPOLATED = 1, HOLE = 2 }

Static Public Member Functions inherited from cellCellStencil
static autoPtr< cellCellStencil >	New (const fvMesh &, const dictionary &dict, const bool update=true)

static const labelIOList &	zoneID (const fvMesh &)

static void	globalCellCells (const globalIndex &gi, const polyMesh &mesh, const boolList &isValidDonor, const labelList &selectedCells, labelListList &cellCells, pointListList &cellCellCentres)

Static Protected Attributes inherited from cellCellStencil
static const Enum< cellType >	cellTypeNames_

Detailed Description

Inverse-distance-weighted interpolation stencil.

hole finding:

mark boundary faces on helper (voxel) mesh
mark any cell overlaying these voxels
use flood filling to find any unreachable cell Alternative is to use an octree of the boundary faces and determine directly for all cells whether we are outside. Might be slow though.

Source files

Definition at line 63 of file inverseDistanceCellCellStencil.H.

Constructor & Destructor Documentation

◆ inverseDistance()

inverseDistance	(	const fvMesh &	mesh,
		const dictionary &	dict,
		const bool	doUpdate
	)

Definition at line 1651 of file inverseDistanceCellCellStencil.C.

References Foam::expressions::patchExpr::debug, Foam::endl(), forAll, Foam::Pout, IOobject::typeHeaderOk(), update(), and zoneID.

Here is the call graph for this function:

◆ ~inverseDistance()

~inverseDistance ( )

virtual

Definition at line 1731 of file inverseDistanceCellCellStencil.C.

Member Function Documentation

◆ index()

Foam::label index	(	const labelVector &	nDivs,
		const labelVector &	ijk
	)

staticprotected

Definition at line 53 of file inverseDistanceCellCellStencil.C.

◆ index3() [1/2]

Foam::labelVector index3	(	const labelVector &	nDivs,
		const label	boxI
	)

staticprotected

Definition at line 63 of file inverseDistanceCellCellStencil.C.

References k.

◆ index3() [2/2]

Foam::labelVector index3	(	const boundBox &	bb,
		const labelVector &	nDivs,
		const point &	pt
	)

staticprotected

Definition at line 78 of file inverseDistanceCellCellStencil.C.

References boundBox::min(), and boundBox::span().

Here is the call graph for this function:

◆ position()

Foam::point position	(	const boundBox &	bb,
		const labelVector &	nDivs,
		const label	boxI
	)

staticprotected

Definition at line 97 of file inverseDistanceCellCellStencil.C.

References boundBox::min(), and boundBox::span().

Here is the call graph for this function:

◆ fill()

void fill	(	PackedList< 2 > &	elems,
		const boundBox &	bb,
		const labelVector &	nDivs,
		const boundBox &	subBb,
		const unsigned int	val
	)

staticprotected

Definition at line 114 of file inverseDistanceCellCellStencil.C.

References k, boundBox::max(), Foam::max(), boundBox::min(), and Foam::min().

Here is the call graph for this function:

◆ overlaps()

bool overlaps	(	const boundBox &	bb,
		const labelVector &	nDivs,
		const PackedList< 2 > &	voxels,
		const treeBoundBox &	subBb,
		const unsigned int	val
	)

staticprotected

Definition at line 264 of file inverseDistanceCellCellStencil.C.

References k, boundBox::max(), Foam::max(), boundBox::min(), and Foam::min().

Here is the call graph for this function:

◆ markBoundaries()

void markBoundaries	(	const fvMesh &	mesh,
		const vector &	smallVec,
		const boundBox &	bb,
		const labelVector &	nDivs,
		PackedList< 2 > &	patchTypes,
		const labelList &	cellMap,
		labelList &	patchCellTypes
	)

staticprotected

Definition at line 152 of file inverseDistanceCellCellStencil.C.

References fvPatch::faceCells(), forAll, boundBox::max(), mesh, boundBox::min(), boundBox::overlaps(), fvPatch::patch(), patchTypes(), and PrimitivePatch< FaceList, PointField >::points().

Here is the call graph for this function:

◆ cellBb()

Foam::treeBoundBox cellBb	(	const primitiveMesh &	mesh,
		const label	celli
	)

staticprotected

Definition at line 230 of file inverseDistanceCellCellStencil.C.

References cells, f(), forAll, boundBox::max(), Foam::max(), mesh, boundBox::min(), Foam::min(), p, and points.

Here is the call graph for this function:

◆ markPatchesAsHoles()

void markPatchesAsHoles	(	PstreamBuffers &	pBufs,
		const PtrList< fvMeshSubset > &	meshParts,
		const List< treeBoundBoxList > &	patchBb,
		const List< labelVector > &	patchDivisions,
		const PtrList< PackedList< 2 >> &	patchParts,
		const label	srcI,
		const label	tgtI,
		labelList &	allCellTypes
	)		const

protected

Definition at line 308 of file inverseDistanceCellCellStencil.C.

References PstreamBuffers::clear(), Foam::exit(), Foam::FatalError, FatalErrorInFunction, PstreamBuffers::finishedSends(), forAll, boundBox::max(), meshParts(), boundBox::min(), os(), and treeBoundBox::overlaps().

Here is the call graph for this function:

◆ betterDonor()

bool betterDonor	(	const label	destMesh,
		const label	currentDonorMesh,
		const label	newDonorMesh
	)		const

protected

Definition at line 434 of file inverseDistanceCellCellStencil.C.

References Foam::mag().

Here is the call graph for this function:

◆ markDonors()

void markDonors	(	const globalIndex &	globalCells,
		PstreamBuffers &	pBufs,
		const PtrList< fvMeshSubset > &	meshParts,
		const List< treeBoundBoxList > &	meshBb,
		const labelList &	allCellTypes,
		const label	srcI,
		const label	tgtI,
		labelListList &	allStencil,
		labelList &	allDonor
	)		const

protected

Definition at line 476 of file inverseDistanceCellCellStencil.C.

References Foam::abort(), List::append(), DynamicList::append(), primitiveMesh::cellCentres(), PstreamBuffers::clear(), Foam::FatalError, FatalErrorInFunction, polyMesh::findCell(), PstreamBuffers::finishedSends(), forAll, boundBox::max(), meshBb(), meshParts(), boundBox::min(), primitiveMesh::nCells(), os(), treeBoundBox::overlaps(), samples(), List::setSize(), polyMesh::tetBasePtIs(), and globalIndex::toGlobal().

Here is the call graph for this function:

◆ uncompactedRegionSplit()

void uncompactedRegionSplit	(	const fvMesh &	mesh,
		const globalIndex &	globalFaces,
		const label	nZones,
		const labelList &	zoneID,
		const labelList &	cellTypes,
		const boolList &	isBlockedFace,
		labelList &	cellRegion
	)		const

protected

◆ compactedRegionSplit()

autoPtr<globalIndex> compactedRegionSplit	(	const fvMesh &	mesh,
		const globalIndex &	globalFaces,
		labelList &	cellRegion
	)		const

protected

◆ findHoles()

void findHoles	(	const globalIndex &	globalCells,
		const fvMesh &	mesh,
		const labelList &	zoneID,
		const labelListList &	stencil,
		labelList &	cellTypes
	)		const