ParallelOverlappingGridInterpolator.C File Reference

#include "ParallelOverlappingGridInterpolator.h"
#include "ParallelUtility.h"
#include "SparseArray.h"

Include dependency graph for ParallelOverlappingGridInterpolator.C:

Macros
#define	Q11(x)   (1.-(x))
#define	Q21(x)   (x)
#define	Q12(x)   .5*((x)-1.)*((x)-2.)
#define	Q22(x)   (x)*(2.-(x))
#define	Q32(x)   .5*(x)*((x)-1.)
#define	q(axis, m)   qp[axis+3*(m)]
#define	dr(axis)   drp[axis]
#define	ir(side, axis)   irp[side+2*(axis)]
#define	isCC(axis)   isCellCenteredp[axis]
#define	ilLocal(i0, i1, p)   ilap[(i1)+(numberOfDimensions+1)*(i0)]
#define	ciLocal(i0, i1, p)   ciap[(i1)+numberOfDimensions*(i0)]
#define	coeffg(i, m1, m2, m3)   coeffap[i+nil*((m1)+width[0]*((m2)+width[1]*(m3)))]
#define	q10(x)   1
#define	q20(x)   -(x)+1
#define	q21(x)   (x)
#define	q30(x)   ((x)-1)*((x)-2)/2.
#define	q31(x)   -(x)*((x)-2)
#define	q32(x)   (x)*((x)-1)/2.
#define	q40(x)   -((x)-1)*((x)-2)*((x)-3)/6.
#define	q41(x)   (x)*((x)-2)*((x)-3)/2.
#define	q42(x)   -(x)*((x)-1)*((x)-3)/2.
#define	q43(x)   (x)*((x)-1)*((x)-2)/6.
#define	q50(x)   ((x)-1)*((x)-2)*((x)-3)*((x)-4)/24.
#define	q51(x)   -(x)*((x)-2)*((x)-3)*((x)-4)/6.
#define	q52(x)   (x)*((x)-1)*((x)-3)*((x)-4)/4.
#define	q53(x)   -(x)*((x)-1)*((x)-2)*((x)-4)/6.
#define	q54(x)   (x)*((x)-1)*((x)-2)*((x)-3)/24.
#define	q60(x)   -((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)/120.
#define	q61(x)   (x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)/24.
#define	q62(x)   -(x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)/12.
#define	q63(x)   (x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)/12.
#define	q64(x)   -(x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)/24.
#define	q65(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)/120.
#define	q70(x)   ((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/720.
#define	q71(x)   -(x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/120.
#define	q72(x)   (x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/48.
#define	q73(x)   -(x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)*((x)-6)/36.
#define	q74(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)*((x)-6)/48.
#define	q75(x)   -(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-6)/120.
#define	q76(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)/720.
#define	q80(x)   -((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/5040.
#define	q81(x)   (x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/720.
#define	q82(x)   -(x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/240.
#define	q83(x)   (x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/144.
#define	q84(x)   -(x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)*((x)-6)*((x)-7)/144.
#define	q85(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-6)*((x)-7)/240.
#define	q86(x)   -(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-7)/720.
#define	q87(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/5040.
#define	q90(x)   ((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/40320.
#define	q91(x)   -(x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/5040.
#define	q92(x)   (x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/1440.
#define	q93(x)   -(x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/720.
#define	q94(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/576.
#define	q95(x)   -(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-6)*((x)-7)*((x)-8)/720.
#define	q96(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-7)*((x)-8)/1440.
#define	q97(x)   -(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-8)/5040.
#define	q98(x)   (x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/40320.
#define	INTERPOLATE_THIS_GRID(grid)   ((!onlyInterpolateSomeGrids || gridsToInterpolate(grid)) && cg.refinementLevelNumber(grid)<=maximumRefinementLevelToInterpolate )
#define	INTERPOLATE_FROM_THIS_GRID(grid)   (!onlyInterpolateFromSomeGrids || gridsToInterpolateFrom(grid))

Functions
void	allocate (int **&a, int dim0, int dim1)
void	deallocate (int **&a, int dim0, int dim1)
void	allocate (int ***&a, int dim0, int dim1, int dim2)
void	deallocate (int ***&a, int dim0, int dim1, int dim2)
void	allocate (int *****&a, int ***dim0, int dim1, int dim2, int dim3, int dim4)
void	deallocate (int *****&a, int ***dim0, int dim1, int dim2, int dim3, int dim4)
void	allocate (real *****&a, int ***dim0, int dim1, int dim2, int dim3, int dim4)
void	deallocate (real *****&a, int ***dim0, int dim1, int dim2, int dim3, int dim4)
void	allocate (int ****&a, int ***dim0, int dim1, int dim2, int dim3)
void	deallocate (int ****&a, int ***dim0, int dim1, int dim2, int dim3)

Macro Definition Documentation

#define ciLocal	(		i0,
			i1,
			p
	)		ciap[(i1)+numberOfDimensions*(i0)]

Referenced by ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), ParallelOverlappingGridInterpolator::internalInterpolate(), InterpolatePointsOnAGrid::parallelInternalInterpolate(), Ogen::updateRefinementFillInterpolationData(), and Ogen::updateRefinementFillInterpolationDataNew().

#define coeffg	(		i,
			m1,
			m2,
			m3
	)		coeffap[i+nil*((m1)+width[0]*((m2)+width[1]*(m3)))]

Referenced by Integrate::computeWeights(), ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), InterpolatePointsOnAGrid::initializeExplicitInterpolation(), and Interpolant::interpolate().

#define dr

(

axis

)

drp[axis]

Referenced by addInteriorNodes(), MappedGridOperators::applyBCaDotGradU(), MappedGridOperators::applyBCGenDiv(), MappedGridOperators::applyBCneumann(), MappedGridOperators::applyBCnormalComponent(), MappedGridOperators::applyBCnormalDerivative(), MappedGridOperators::applyBCtangentialComponent(), EllipticGridGenerator::applyBoundaryConditions(), Ogmg::assignBoundaryConditionCoefficients(), MappedGridOperators::assignCoefficientsInternal(), DataPointMapping::basicInverseS(), CompositeTopology::buildEdgeSegment(), UnstructuredMapping::buildFromAMapping(), UnstructuredMapping::buildFromARegularMapping(), InterpolatePoints::buildInterpolationInfo(), OffsetShell::buildOffsetMappings(), MappingBuilder::buildSurfacePatch(), Ogmg::checkForBetterQualityInterpolation(), Mapping::checkMapping(), DerivedFunctions::computeDerivedFunctions(), Integrate::computeStitchedSurfaceWeightsOld(), Integrate::computeSurfaceWeights(), Integrate::computeWeights(), TrimmedMappingBuilder::constructOuterBoundaryCurves(), PlotIt::contour3d(), PlotIt::contourCuts(), HyperbolicMapping::createCurveOnSurface(), Ogen::cutHoles(), Ogen::cutHolesNew(), Ogen::cutHolesNewer(), MappedGridOperators::derivativeInternal(), IntersectionMapping::determineIntersection(), HyperbolicMapping::equidistributeAndStretch(), HyperbolicMapping::findMatchingBoundaryCurve(), ApproximateGlobalInverse::findNearestGridPoint(), Mapping::findNearestGridPoint(), Ogen::getAdjustmentVectors(), ReparameterizationTransform::getBounds(), Mapping::getGrid(), Mapping::getGridSerial(), Ogen::getHoleWidth(), Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), Ogmg::getInterpolationCoordinatesNewOld(), Ogmg::getInterpolationStencil(), CrossSectionMapping::initialize(), NurbsMapping::initialize(), ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), InterpolatePointsOnAGrid::initializeExplicitInterpolation(), NurbsMapping::interpolate(), Ogen::interpolateAPoint(), interpolatePoints(), ExposedPoints::interpolatePoints(), ExactLocalInverse::inverse(), CrossSectionMapping::mapGeneral(), CrossSectionMapping::mapGeneralWithPolarSingularity(), DataPointMapping::mapS(), DataPointMapping::mapScalar(), ogen(), PlotIt::plotGrid(), PlotIt::plotStructured(), IntersectionMapping::reparameterize(), ReparameterizationTransform::scaleBounds(), ReparameterizationTransform::setBounds(), DataPointMapping::setMapping(), GridSmoother::smooth(), TrimmedMapping::triangulate(), RocketMapping::update(), GridStretcher::update(), StretchMapping::update(), SplineMapping::update(), NurbsMapping::update(), DataFormats::writePlot3d(), xInterpolate(), and xInterpolateOpt().

#define ilLocal	(		i0,
			i1,
			p
	)		ilap[(i1)+(numberOfDimensions+1)*(i0)]

Referenced by ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), ParallelOverlappingGridInterpolator::internalInterpolate(), and Ogen::updateRefinementFillInterpolationData().

#define INTERPOLATE_FROM_THIS_GRID

(

grid

)

(!onlyInterpolateFromSomeGrids || gridsToInterpolateFrom(grid))

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define INTERPOLATE_THIS_GRID

(

grid

)

((!onlyInterpolateSomeGrids || gridsToInterpolate(grid)) && cg.refinementLevelNumber(grid)<=maximumRefinementLevelToInterpolate )

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define ir	(		side,
			axis
	)		irp[side+2*(axis)]

Referenced by ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), InterpolatePointsOnAGrid::initializeExplicitInterpolation(), main(), GL_GraphicsInterface::saveRasterInAFile(), Ogmg::smoothInterpolationNeighbours(), and Ogen::updateRefinementFillInterpolationDataNew().

#define isCC

(

axis

)

isCellCenteredp[axis]

Referenced by ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), and InterpolatePointsOnAGrid::initializeExplicitInterpolation().

#define q	(		axis,
			m
	)		qp[axis+3*(m)]

Referenced by Ogen::adjustBoundary(), EdgeInfo::adjustOneSegmentEndPoints(), TrimmedMapping::assembleSubCurves(), CompositeTopology::buildEdgeCurves(), CanInterpolate::canInterpolate(), CanInterpolate::canInterpolateOld(), CompositeTopology::checkConsistency(), Ogen::computeInterpolationQuality(), GL_GraphicsInterface::constructor(), MappingsFromCAD::createCompositeCurve(), GL_GraphicsInterface::drawColouredSquares(), TrimmedMapping::editNurbsTrimCurve(), RotatedBox::fitBox(), GL_GraphicsInterface::generateNewDisplayList(), CompositeTopology::get(), GenericCompositeGridOperators::getInterpolationCoefficients(), GenericGraphicsInterface::getValues(), Ogen::improveQuality(), ParallelOverlappingGridInterpolator::initializeExplicitInterpolation(), InterpolatePointsOnAGrid::initializeExplicitInterpolation(), interpolateRefinements(), NurbsMapping::interpolateSurface(), NurbsMapping::interpolateVolume(), CompositeTopology::joinEdgeCurves(), main(), CompositeTopology::merge(), CompositeTopology::mergeEdgeCurves(), GL_GraphicsInterface::pickPoints(), GL_GraphicsInterface::plotLines(), GL_GraphicsInterface::plotPoints(), CompositeTopology::put(), DataFormats::readPlot3d(), GL_GraphicsInterface::setupHardCopy(), CompositeTopology::splitAndMergeEdgeCurves(), CompositeTopology::splitEdge(), Stencil::test(), CompositeTopology::update(), CompositeSurface::update(), GL_GraphicsInterface::worldToNormalizedCoordinates(), and DataFormats::writePlot3d().

#define q10

(

x

)

1

#define Q11

(

x

)

(1.-(x))

#define Q12

(

x

)

.5*((x)-1.)*((x)-2.)

#define q20

(

x

)

-(x)+1

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define Q21

(

x

)

(x)

#define q21

(

x

)

(x)

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define Q22

(

x

)

(x)*(2.-(x))

#define q30

(

x

)

((x)-1)*((x)-2)/2.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q31

(

x

)

-(x)*((x)-2)

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define Q32

(

x

)

.5*(x)*((x)-1.)

#define q32

(

x

)

(x)*((x)-1)/2.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q40

(

x

)

-((x)-1)*((x)-2)*((x)-3)/6.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q41

(

x

)

(x)*((x)-2)*((x)-3)/2.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q42

(

x

)

-(x)*((x)-1)*((x)-3)/2.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q43

(

x

)

(x)*((x)-1)*((x)-2)/6.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q50

(

x

)

((x)-1)*((x)-2)*((x)-3)*((x)-4)/24.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q51

(

x

)

-(x)*((x)-2)*((x)-3)*((x)-4)/6.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q52

(

x

)

(x)*((x)-1)*((x)-3)*((x)-4)/4.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q53

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-4)/6.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q54

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)/24.

Referenced by ParallelOverlappingGridInterpolator::internalInterpolate().

#define q60

(

x

)

-((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)/120.

#define q61

(

x

)

(x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)/24.

#define q62

(

x

)

-(x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)/12.

#define q63

(

x

)

(x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)/12.

#define q64

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)/24.

#define q65

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)/120.

#define q70

(

x

)

((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/720.

#define q71

(

x

)

-(x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/120.

#define q72

(

x

)

(x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/48.

#define q73

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)*((x)-6)/36.

#define q74

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)*((x)-6)/48.

#define q75

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-6)/120.

#define q76

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)/720.

#define q80

(

x

)

-((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/5040.

#define q81

(

x

)

(x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/720.

#define q82

(

x

)

-(x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/240.

#define q83

(

x

)

(x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/144.

#define q84

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)*((x)-6)*((x)-7)/144.

#define q85

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-6)*((x)-7)/240.

#define q86

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-7)/720.

#define q87

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)/5040.

#define q90

(

x

)

((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/40320.

#define q91

(

x

)

-(x)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/5040.

#define q92

(

x

)

(x)*((x)-1)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/1440.

#define q93

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-4)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/720.

#define q94

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-5)*((x)-6)*((x)-7)*((x)-8)/576.

#define q95

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-6)*((x)-7)*((x)-8)/720.

#define q96

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-7)*((x)-8)/1440.

#define q97

(

x

)

-(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-8)/5040.

#define q98

(

x

)

(x)*((x)-1)*((x)-2)*((x)-3)*((x)-4)*((x)-5)*((x)-6)*((x)-7)/40320.

Function Documentation

void allocate	(	int **&	a,
		int	dim0,
		int	dim1
	)

References max().

Referenced by HarwellEquationSolver::solve(), YaleEquationSolver::solve(), and SlapEquationSolver::solve().

void allocate	(	int ***&	a,
		int	dim0,
		int	dim1,
		int	dim2
	)

References i2, and max().

void allocate	(	int *****&	a,
		int ***	dim0,
		int	dim1,
		int	dim2,
		int	dim3,
		int	dim4
	)

References i2, and max().

void allocate	(	real *****&	a,
		int ***	dim0,
		int	dim1,
		int	dim2,
		int	dim3,
		int	dim4
	)

References i2, and max().

void allocate	(	int ****&	a,
		int ***	dim0,
		int	dim1,
		int	dim2,
		int	dim3
	)

References i2, and max().

void deallocate	(	int **&	a,
		int	dim0,
		int	dim1
	)

References NULL.

void deallocate	(	int ***&	a,
		int	dim0,
		int	dim1,
		int	dim2
	)

References i2, and NULL.

void deallocate	(	int *****&	a,
		int ***	dim0,
		int	dim1,
		int	dim2,
		int	dim3,
		int	dim4
	)

References i2, and NULL.

void deallocate	(	real *****&	a,
		int ***	dim0,
		int	dim1,
		int	dim2,
		int	dim3,
		int	dim4
	)

References i2, and NULL.

void deallocate	(	int ****&	a,
		int ***	dim0,
		int	dim1,
		int	dim2,
		int	dim3
	)

References i2, and NULL.