InterpolatePoints.C File Reference

#include "Overture.h"
#include "display.h"
#include "InterpolatePoints.h"
#include "MappingProjectionParameters.h"
#include "ParallelUtility.h"

Include dependency graph for InterpolatePoints.C:

Macros
#define	NRM(axis)   ( indexRange(End,axis)-indexRange(Start,axis)+1 )
#define	MODR(i, axis)
#define	IA(i0)   iap[i0]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	NUMBEROFINTERPOLATIONPOINTS(i0)   numberOfInterpolationPointsp[i0]
#define	STATUS(i0)   statusp[i0]
#define	IA0(i0)   ia0p[i0]
#define	XPROJECT(i0, i1)   xProjectp[i0+xProjectDim0*(i1)]
#define	DIMENSION(i0, i1)   dimensionp[i0+dimensionDim0*(i1)]
#define	INDEXRANGE(i0, i1)   indexRangep[i0+indexRangeDim0*(i1)]
#define	GRIDSPACING(i0)   gridSpacingp[i0]
#define	rbb(side, axis)   (prbb[(side)+2*(axis)])
#define	X(i0, i1)   xp[i0+xDim0*(i1)]
#define	RA(i0, i1)   rap[i0+raDim0*(i1)]
#define	XA(i0, i1)   xap[i0+xaDim0*(i1)]
#define	IA(i0)   iap[i0]
#define	RPROJECT(i0, i1)   rProjectp[i0+rProjectDim0*(i1)]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	IP1(i0, i1)   ip1p[i0+ip1Dim0*(i1)]
#define	DRA(i0, i1)   drap[i0+draDim0*(i1)]
#define	DR(i0, i1)   drp[i0+drDim0*(i1)]
#define	gridIndexRange(side, axis)   gridIndexRangep[(side)+2*(axis)]
#define	iRange(side, axis)   iRangep[(side)+2*(axis)]
#define	MASK(i0, i1)   maskp[i0+maskDim0*(i1)]
#define	MASK(i0, i1, i2)   maskp[i0+maskDim0*(i1+maskDim1*(i2))]
#define	IA(i0)   iap[i0]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	IP1(i0, i1)   ip1p[i0+ip1Dim0*(i1)]
#define	DRA(i0, i1)   drap[i0+draDim0*(i1)]
#define	UINTERPOLATED(i0, i1)   uInterpolatedp[i0+uInterpolatedDim0*(i1)]
#define	UG(i0, i1, i2)   ugp[i0+ugDim0*(i1+ugDim1*(i2))]
#define	UG(i0, i1, i2, i3)   ugp[i0+ugDim0*(i1+ugDim1*(i2+ugDim2*(i3)))]
#define	IA(i0)   iap[i0]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	IP1(i0, i1)   ip1p[i0+ip1Dim0*(i1)]
#define	DRA(i0, i1)   drap[i0+draDim0*(i1)]
#define	UINTERPOLATED(i0, ix, iy, iz)   uInterpolatedp[i0+ uInterpolatedDim0*( (ix)+2*((iy) +2*(iz)))]
#define	UG(i0, i1, i2)   ugp[i0+ugDim0*(i1+ugDim1*(i2))]
#define	UG(i0, i1, i2, i3)   ugp[i0+ugDim0*(i1+ugDim1*(i2+ugDim2*(i3)))]
#define	FOR_3D(i1, i2, i3, I1, I2, I3)
#define	FOR_3(i1, i2, i3, I1, I2, I3)
#define	IA(i0, i1)   iap[i0+iaDim0*(i1)]
#define	MASK(i0, i1, i2)   maskp[i0+maskDim0*(i1+maskDim1*(i2))]
#define	X(i0, i1)   xp[i0+xDim0*(i1)]
#define	UINTERPOLATED(i0, i1)   uInterpolatedp[i0+uInterpolatedDim0*(i1)]
#define	COORD0(i0, i1, i2)   (xa+dx0*(i0-i0a))
#define	COORD1(i0, i1, i2)   (ya+dy0*(i1-i1a))
#define	COORD2(i0, i1, i2)   (za+dz0*(i2-i2a))
#define	CENTER(i0, i1, i2, i3)   centerp[i0+centerDim0*(i1+centerDim1*(i2+centerDim2*(i3)))]
#define	U(i0, i1, i2, i3)   up[i0+uDim0*(i1+uDim1*(i2+uDim2*(i3)))]

Macro Definition Documentation

#define CENTER	(		i0,
			i1,
			i2,
			i3
	)		centerp[i0+centerDim0*(i1+centerDim1*(i2+centerDim2*(i3)))]

Referenced by Ogmg::buildExtraLevels(), Ogmg::buildExtraLevelsNew(), Integrate::computeWeights(), Ogmg::getInterpolationCoordinates(), InterpolatePoints::interpolateAllPoints(), InterpolatePointsOnAGrid::interpolateAllPoints(), and xInterpolateOpt().

#define COORD0	(		i0,
			i1,
			i2
	)		(xa+dx0*(i0-i0a))

#define COORD1	(		i0,
			i1,
			i2
	)		(ya+dy0*(i1-i1a))

#define COORD2	(		i0,
			i1,
			i2
	)		(za+dz0*(i2-i2a))

#define DIMENSION	(		i0,
			i1
	)		dimensionp[i0+dimensionDim0*(i1)]

Referenced by InterpolatePoints::buildInterpolationInfo().

#define DR	(		i0,
			i1
	)		drp[i0+drDim0*(i1)]

Referenced by InterpolatePoints::buildInterpolationInfo().

#define DRA	(		i0,
			i1
	)		drap[i0+draDim0*(i1)]

Referenced by InterpolatePoints::buildInterpolationInfo(), InterpolatePoints::interpolatePoints(), and InterpolatePoints::interpolationCoefficients().

#define DRA	(		i0,
			i1
	)		drap[i0+draDim0*(i1)]

#define DRA	(		i0,
			i1
	)		drap[i0+draDim0*(i1)]

#define FOR_3	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)

Value:

I1Base =I1.getBase(),   I2Base =I2.getBase(),  I3Base =I3.getBase();  \
I1Bound=I1.getBound(),  I2Bound=I2.getBound(), I3Bound=I3.getBound(); \
for(i3=I3Base; i3<=I3Bound; i3++) \
for(i2=I2Base; i2<=I2Bound; i2++) \
for(i1=I1Base; i1<=I1Bound; i1++)

#define FOR_3D	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)

Value:

int I1Base =I1.getBase(),   I2Base =I2.getBase(),  I3Base =I3.getBase();  \
int I1Bound=I1.getBound(),  I2Bound=I2.getBound(), I3Bound=I3.getBound(); \
for(i3=I3Base; i3<=I3Bound; i3++) \
for(i2=I2Base; i2<=I2Bound; i2++) \
for(i1=I1Base; i1<=I1Bound; i1++)

Referenced by InterpolatePoints::interpolateAllPoints().

#define gridIndexRange	(		side,
			axis
	)		gridIndexRangep[(side)+2*(axis)]

Referenced by Ogen::adjustBoundary(), Ogen::adjustBoundarySerial(), EllipticGridGenerator::applyBoundaryConditions(), Integrate::buildAdaptiveMeshRefinementSurfaceWeights(), Ogmg::buildExtraLevels(), Ogmg::buildExtraLevelsNew(), InterpolatePoints::buildInterpolationInfo(), MappingBuilder::buildSurfacePatch(), Ogen::changeParameters(), DerivedFunctions::computeDerivedFunctions(), computeMaximumSpeed(), AdvancingFront::computeTransformationAtPoint(), Integrate::computeWeights(), PlotIt::contour1d(), PlotIt::contour2d(), PlotIt::contour3d(), PlotIt::contourOpt2d(), EllipticGridGenerator::defineBoundaryControlFunction(), EllipticGridGenerator::defineSurfaceControlFunction(), EllipticGridGenerator::determineBoundarySpacing(), Mapping::display(), dudt(), EllipticGridGenerator::estimateUnderRelaxationCoefficients(), fileOutput(), EllipticGridGenerator::generateGrid(), Ugen::generateSpacingControlMesh(), Ugen::generateSpacingControlMeshForSurface(), Mapping::getBoundingBox(), EllipticGridGenerator::getControlFunctions(), PlotIt::getGridBounds(), Mapping::getGridDimensions(), PlotIt::grid3d(), interpolateFromControlFunction(), Ogen::interpolateMixedBoundary(), EllipticGridGenerator::jacobi(), EllipticGridGenerator::lineSmoother(), main(), EllipticGridGenerator::multigridVcycle(), CompositeGrid::numberOfPossibleMultigridLevels(), ogen(), doubleMappedGridFunction::periodicUpdate(), floatMappedGridFunction::periodicUpdate(), intMappedGridFunction::periodicUpdate(), PlotIt::plotGrid(), PlotIt::plotGrid2d(), PlotIt::plotGridBoundaries(), PlotIt::plotStructured(), MappedGridData::put(), Mapping::put(), EllipticGridGenerator::redBlack(), EllipticGridGenerator::restrictMovement(), FileOutput::saveData(), Ogen::setGridParameters(), EllipticGridGenerator::smoothJacobi(), EllipticGridGenerator::startingGrid(), PlotIt::streamLines2d(), PlotIt::streamLines3d(), EllipticGridGenerator::stretchTheGrid(), PlotIt::surfaceGrid3d(), Integrate::surfaceIntegral(), Interpolant::testInterpolation(), GridStretcher::update(), EllipticGridGenerator::update(), EllipticGridGenerator::updateRightHandSideWithFASCorrection(), PlotIt::userDefinedOutput(), EllipticGridGenerator::weightFunction(), xInterpolate(), and xInterpolateOpt().

#define GRIDSPACING

(

i0

)

gridSpacingp[i0]

Referenced by InterpolatePoints::buildInterpolationInfo().

#define IA

(

i0

)

iap[i0]

#define IA

(

i0

)

iap[i0]

#define IA

(

i0

)

iap[i0]

#define IA

(

i0

)

iap[i0]

#define IA	(		i0,
			i1
	)		iap[i0+iaDim0*(i1)]

#define IA0

(

i0

)

ia0p[i0]

Referenced by InterpolatePoints::buildInterpolationInfo(), InterpolatePointsOnAGrid::buildInterpolationInfo(), and Ogmg::getInterpolationCoordinatesNew().

#define INDEXRANGE	(		i0,
			i1
	)		indexRangep[i0+indexRangeDim0*(i1)]

Referenced by Ogmg::buildExtraLevels(), Ogmg::buildExtraLevelsNew(), and InterpolatePoints::buildInterpolationInfo().

#define IP	(		i0,
			i1
	)		ipp[i0+ipDim0*(i1)]

#define IP	(		i0,
			i1
	)		ipp[i0+ipDim0*(i1)]

#define IP	(		i0,
			i1
	)		ipp[i0+ipDim0*(i1)]

#define IP	(		i0,
			i1
	)		ipp[i0+ipDim0*(i1)]

#define IP1	(		i0,
			i1
	)		ip1p[i0+ip1Dim0*(i1)]

Referenced by InterpolatePoints::buildInterpolationInfo(), and InterpolatePoints::interpolatePoints().

#define IP1	(		i0,
			i1
	)		ip1p[i0+ip1Dim0*(i1)]

#define IP1	(		i0,
			i1
	)		ip1p[i0+ip1Dim0*(i1)]

#define iRange	(		side,
			axis
	)		iRangep[(side)+2*(axis)]

Referenced by InterpolatePoints::buildInterpolationInfo().

#define MASK	(		i0,
			i1
	)		maskp[i0+maskDim0*(i1)]

#define MASK	(		i0,
			i1,
			i2
	)		maskp[i0+maskDim0*(i1+maskDim1*(i2))]

#define MASK	(		i0,
			i1,
			i2
	)		maskp[i0+maskDim0*(i1+maskDim1*(i2))]

#define MODR	(		i,
			axis
	)

Value:

( \
  ( (i-indexRange(Start,axis)+NRM(axis)) % NRM(axis)) \
      +indexRange(Start,axis) \
                           )

Referenced by InterpolatePoints::buildInterpolationInfo(), and xInterpolateOpt().

#define NRM

(

axis

)

( indexRange(End,axis)-indexRange(Start,axis)+1 )

#define NUMBEROFINTERPOLATIONPOINTS

(

i0

)

numberOfInterpolationPointsp[i0]

Referenced by InterpolatePoints::buildInterpolationInfo(), and InterpolatePointsOnAGrid::buildInterpolationInfo().

#define RA	(		i0,
			i1
	)		rap[i0+raDim0*(i1)]

#define rbb	(		side,
			axis
	)		(prbb[(side)+2*(axis)])

Referenced by InterpolatePoints::buildInterpolationInfo(), and InterpolatePointsOnAGrid::buildInterpolationInfo().

#define RPROJECT	(		i0,
			i1
	)		rProjectp[i0+rProjectDim0*(i1)]

Referenced by InterpolatePoints::buildInterpolationInfo(), and InterpolatePointsOnAGrid::buildInterpolationInfo().

#define STATUS

(

i0

)

statusp[i0]

Referenced by InterpolatePoints::buildInterpolationInfo(), InterpolatePointsOnAGrid::buildInterpolationInfo(), and SurfaceStitcher::enlargeGapWidth().

#define U	(		i0,
			i1,
			i2,
			i3
	)		up[i0+uDim0*(i1+uDim1*(i2+uDim2*(i3)))]

#define UG	(		i0,
			i1,
			i2
	)		ugp[i0+ugDim0*(i1+ugDim1*(i2))]

#define UG	(		i0,
			i1,
			i2,
			i3
	)		ugp[i0+ugDim0*(i1+ugDim1*(i2+ugDim2*(i3)))]

#define UG	(		i0,
			i1,
			i2
	)		ugp[i0+ugDim0*(i1+ugDim1*(i2))]

#define UG	(		i0,
			i1,
			i2,
			i3
	)		ugp[i0+ugDim0*(i1+ugDim1*(i2+ugDim2*(i3)))]

#define UINTERPOLATED	(		i0,
			i1
	)		uInterpolatedp[i0+uInterpolatedDim0*(i1)]

Referenced by InterpolatePoints::interpolateAllPoints(), InterpolatePointsOnAGrid::interpolateAllPoints(), InterpolatePoints::interpolatePoints(), and InterpolatePoints::interpolationCoefficients().

#define UINTERPOLATED	(		i0,
			ix,
			iy,
			iz
	)		uInterpolatedp[i0+ uInterpolatedDim0*( (ix)+2*((iy) +2*(iz)))]

#define UINTERPOLATED	(		i0,
			i1
	)		uInterpolatedp[i0+uInterpolatedDim0*(i1)]

#define X	(		i0,
			i1
	)		xp[i0+xDim0*(i1)]

#define X	(		i0,
			i1
	)		xp[i0+xDim0*(i1)]

#define XA	(		i0,
			i1
	)		xap[i0+xaDim0*(i1)]

#define XPROJECT	(		i0,
			i1
	)		xProjectp[i0+xProjectDim0*(i1)]

Referenced by InterpolatePoints::buildInterpolationInfo().