buildExtraLevels.C File Reference

#include "Ogmg.h"
#include "display.h"
#include "Ogen.h"
#include "ParallelUtility.h"

Include dependency graph for buildExtraLevels.C:

Macros
#define	FOR_3(i1, i2, i3, I1, I2, I3)
#define	FOR_3D(i1, i2, i3, I1, I2, I3)
#define	MASK1(i0, i1, i2)   mask1p[i0+mask1Dim0*(i1+mask1Dim1*(i2))]
#define	EIR(i0, i1)   eirp[i0+eirDim0*(i1)]
#define	X(i0, i1)   xp[i0+xDim0*(i1)]
#define	R(i0, i1)   rp[i0+rDim0*(i1)]
#define	CENTER00(i0, i1, i2)   (xab0[0][0]+dx0[0]*(i0-iv0[0]))
#define	CENTER01(i0, i1, i2)   (xab0[0][1]+dx0[1]*(i1-iv0[1]))
#define	CENTER02(i0, i1, i2)   (xab0[0][2]+dx0[2]*(i2-iv0[2]))
#define	CENTER10(i0, i1, i2)   (xab1[0][0]+dx1[0]*(i0-iv1[0]))
#define	CENTER11(i0, i1, i2)   (xab1[0][1]+dx1[1]*(i1-iv1[1]))
#define	CENTER12(i0, i1, i2)   (xab1[0][2]+dx1[2]*(i2-iv1[2]))
#define	INTERPOLEEGRID0(i0)   interpoleeGrid0p[i0]
#define	INTERPOLATIONPOINT0(i0, i1)   interpolationPoint0p[i0+interpolationPoint0Dim0*(i1)]
#define	INTERPOLATIONCOORDINATES0(i0, i1)   interpolationCoordinates0p[i0+interpolationCoordinates0Dim0*(i1)]
#define	CENTER(i0, i1, i2, i3)   centerp[i0+centerDim0*(i1+centerDim1*(i2+centerDim2*(i3)))]
#define	EIR(i0, i1)   eirp[i0+eirDim0*(i1)]
#define	INDEXRANGE0(i0, i1)   indexRange0p[i0+indexRange0Dim0*(i1)]
#define	INDEXRANGE(i0, i1)   indexRangep[i0+indexRangeDim0*(i1)]
#define	IA(i0, i1)   iap[i0+iaDim0*(i1)]
#define	INVERSEGRID(i0, i1, i2)   inverseGridp[i0+inverseGridDim0*(i1+inverseGridDim1*(i2))]
#define	INVERSECOORDINATES(i0, i1, i2, i3)   inverseCoordinatesp[i0+inverseCoordinatesDim0*(i1+inverseCoordinatesDim1*(i2+inverseCoordinatesDim2*(i3)))]
#define	IB(i0, i1)   ibp[i0+ibDim0*(i1)]
#define	INTERPOLEELOCATION(i0, i1)   interpoleeLocationp[i0+interpoleeLocationDim0*(i1)]
#define	INTERPOLEELOCATION1(i0, i1)   interpoleeLocation1p[i0+interpoleeLocation1Dim0*(i1)]
#define	INTERPOLATIONPOINT(i0, i1)   interpolationPointp[i0+interpolationPointDim0*(i1)]
#define	INTERPOLATIONPOINT1(i0, i1)   interpolationPoint1p[i0+interpolationPoint1Dim0*(i1)]
#define	INTERPOLATIONCOORDINATES(i0, i1)   interpolationCoordinatesp[i0+interpolationCoordinatesDim0*(i1)]
#define	INTERPOLATIONCOORDINATES1(i0, i1)   interpolationCoordinates1p[i0+interpolationCoordinates1Dim0*(i1)]
#define	INTERPOLEEGRID(i0)   interpoleeGridp[i0]
#define	INTERPOLEEGRID1(i0)   interpoleeGrid1p[i0]
#define	NG(i0)   ngp[i0]
#define	GRIDSTART(i0)   gridStartp[i0]
#define	VARIABLEINTERPOLATIONWIDTH(i0)   variableInterpolationWidthp[i0]
#define	VARIABLEINTERPOLATIONWIDTH1(i0)   variableInterpolationWidth1p[i0]
#define	FOR_IJ()
#define	MASK(i0, i1, i2)   maskp[i0+maskDim0*(i1+maskDim1*(i2))]
#define	CENTER0(i0, i1, i2, i3)   center0p[i0+center0Dim0*(i1+center0Dim1*(i2+center0Dim2*(i3)))]
#define	IA(i0)   iap[i0]
#define	RB(i0, i1)   rbp[i0+rbDim0*(i1)]
#define	XR(m, n)   xra[n][m]
#define	INTERPOLATIONCOORDINATES0(i0, i1)   interpolationCoordinates0p[i0+interpolationCoordinates0Dim0*(i1)]
#define	XR(m, n)   xra[n][m]
#define	indexRange(i0, i1)   indexRangep[i0+indexRangeDim0*(i1)]
#define	extendedIndexRange(i0, i1)   extendedIndexRangep[i0+extendedIndexRangeDim0*(i1)]
#define	boundaryCondition(i0, i1)   boundaryConditionp[i0+boundaryConditionDim0*(i1)]
#define	gridSpacing(i0)   gridSpacingp[i0]
#define	MASK2(i0, i1, i2)   mask2p[i0+mask2Dim0*(i1+mask2Dim1*(i2))]
#define	MASK2D(k1, k2, k3)   (MASK2(k1,k2,k3)>0 ? 1 : MASK2(k1,k2,k3)==0 ? 0 : -1)
#define	MASK2_W5A(k1, k2, k3)   MASK2(k1,k2,k3)!=0 && MASK2(k1+1,k2,k3)!=0 && MASK2(k1+2,k2,k3)!=0 && MASK2(k1+3,k2,k3)!=0 && MASK2(k1+4,k2,k3)!=0
#define	MASK2_W5B(k1, k2, k3)
#define	MASK2_W5C(k1, k2, k3)
#define	MASK2_W4A(k1, k2, k3)   MASK2(k1,k2,k3)!=0 && MASK2(k1+1,k2,k3)!=0 && MASK2(k1+2,k2,k3)!=0 && MASK2(k1+3,k2,k3)!=0
#define	MASK2_W4B(k1, k2, k3)   MASK2_W4A(k1,k2,k3) && MASK2_W4A(k1,k2+1,k3) && MASK2_W4A(k1,k2+2,k3) && MASK2_W4A(k1,k2+3,k3)
#define	MASK2_W4C(k1, k2, k3)   MASK2_W4B(k1,k2,k3) && MASK2_W4B(k1,k2,k3+1) && MASK2_W4B(k1,k2,k3+2) && MASK2_W4B(k1,k2,k3+3)

Functions
int	checkGrid (CompositeGrid &cg, GenericGraphicsInterface *ps=0, int debug=0)
int	displayMaskLaTeX (const intArray &mask, const aString &label=nullString, FILE *file=NULL, const DisplayParameters *displayParameters=NULL)
int Ogmg::	markGhostPoints (CompositeGrid &cg)

Macro Definition Documentation

#define boundaryCondition	(		i0,
			i1
	)		boundaryConditionp[i0+boundaryConditionDim0*(i1)]

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

#define CENTER0	(		i0,
			i1,
			i2,
			i3
	)		center0p[i0+center0Dim0*(i1+center0Dim1*(i2+center0Dim2*(i3)))]

Referenced by Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), Ogmg::getInterpolationCoordinatesNewOld(), PlotIt::grid3d(), and PlotIt::plotGrid2d().

#define CENTER00	(		i0,
			i1,
			i2
	)		(xab0[0][0]+dx0[0]*(i0-iv0[0]))

Referenced by Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), and Ogmg::getInterpolationCoordinatesNewOld().

#define CENTER01	(		i0,
			i1,
			i2
	)		(xab0[0][1]+dx0[1]*(i1-iv0[1]))

Referenced by Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), and Ogmg::getInterpolationCoordinatesNewOld().

#define CENTER02	(		i0,
			i1,
			i2
	)		(xab0[0][2]+dx0[2]*(i2-iv0[2]))

Referenced by Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), and Ogmg::getInterpolationCoordinatesNewOld().

#define CENTER10	(		i0,
			i1,
			i2
	)		(xab1[0][0]+dx1[0]*(i0-iv1[0]))

Referenced by Ogmg::buildExtraLevels(), Ogmg::buildExtraLevelsNew(), Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), and Ogmg::getInterpolationCoordinatesNewOld().

#define CENTER11	(		i0,
			i1,
			i2
	)		(xab1[0][1]+dx1[1]*(i1-iv1[1]))

Referenced by Ogmg::buildExtraLevels(), Ogmg::buildExtraLevelsNew(), Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), and Ogmg::getInterpolationCoordinatesNewOld().

#define CENTER12	(		i0,
			i1,
			i2
	)		(xab1[0][2]+dx1[2]*(i2-iv1[2]))

Referenced by Ogmg::buildExtraLevels(), Ogmg::buildExtraLevelsNew(), Ogmg::getInterpolationCoordinates(), Ogmg::getInterpolationCoordinatesNew(), and Ogmg::getInterpolationCoordinatesNewOld().

#define EIR	(		i0,
			i1
	)		eirp[i0+eirDim0*(i1)]

Referenced by Ogmg::buildExtraLevels().

#define EIR	(		i0,
			i1
	)		eirp[i0+eirDim0*(i1)]

#define extendedIndexRange	(		i0,
			i1
	)		extendedIndexRangep[i0+extendedIndexRangeDim0*(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++ )

Referenced by Ogmg::buildExtraLevels().

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

Value:

int I1Base,I2Base,I3Base;\
  int I1Bound,I2Bound,I3Bound;\
  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++ )

Referenced by Ogmg::buildExtraLevels().

#define FOR_IJ

(

)

Value:

for( i3=I3Base, j3=J3Base; i3<=I3Bound; i3++,j3++ ) \
  for( i2=I2Base, j2=J2Base; i2<=I2Bound; i2++,j2++ ) \
  for( i1=I1Base, j1=J1Base; i1<=I1Bound; i1++,j1++ )

Referenced by markGhostPoints().

#define gridSpacing

(

i0

)

gridSpacingp[i0]

#define GRIDSTART

(

i0

)

gridStartp[i0]

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

#define IA

(

i0

)

iap[i0]

#define IB	(		i0,
			i1
	)		ibp[i0+ibDim0*(i1)]

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

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

#define INDEXRANGE0	(		i0,
			i1
	)		indexRange0p[i0+indexRange0Dim0*(i1)]

Referenced by Ogmg::buildExtraLevels().

#define INTERPOLATIONCOORDINATES	(		i0,
			i1
	)		interpolationCoordinatesp[i0+interpolationCoordinatesDim0*(i1)]

#define INTERPOLATIONCOORDINATES0	(		i0,
			i1
	)		interpolationCoordinates0p[i0+interpolationCoordinates0Dim0*(i1)]

Referenced by Ogmg::buildExtraLevels(), and Ogmg::getInterpolationCoordinates().

#define INTERPOLATIONCOORDINATES0	(		i0,
			i1
	)		interpolationCoordinates0p[i0+interpolationCoordinates0Dim0*(i1)]

#define INTERPOLATIONCOORDINATES1	(		i0,
			i1
	)		interpolationCoordinates1p[i0+interpolationCoordinates1Dim0*(i1)]

#define INTERPOLATIONPOINT	(		i0,
			i1
	)		interpolationPointp[i0+interpolationPointDim0*(i1)]

#define INTERPOLATIONPOINT0	(		i0,
			i1
	)		interpolationPoint0p[i0+interpolationPoint0Dim0*(i1)]

Referenced by Ogmg::buildExtraLevels(), and Ogmg::getInterpolationCoordinates().

#define INTERPOLATIONPOINT1	(		i0,
			i1
	)		interpolationPoint1p[i0+interpolationPoint1Dim0*(i1)]

#define INTERPOLEEGRID

(

i0

)

interpoleeGridp[i0]

#define INTERPOLEEGRID0

(

i0

)

interpoleeGrid0p[i0]

Referenced by Ogmg::buildExtraLevels(), and Ogmg::getInterpolationCoordinates().

#define INTERPOLEEGRID1

(

i0

)

interpoleeGrid1p[i0]

#define INTERPOLEELOCATION	(		i0,
			i1
	)		interpoleeLocationp[i0+interpoleeLocationDim0*(i1)]

#define INTERPOLEELOCATION1	(		i0,
			i1
	)		interpoleeLocation1p[i0+interpoleeLocation1Dim0*(i1)]

#define INVERSECOORDINATES	(		i0,
			i1,
			i2,
			i3
	)		inverseCoordinatesp[i0+inverseCoordinatesDim0*(i1+inverseCoordinatesDim1*(i2+inverseCoordinatesDim2*(i3)))]

Referenced by Ogmg::buildExtraLevels().

#define INVERSEGRID	(		i0,
			i1,
			i2
	)		inverseGridp[i0+inverseGridDim0*(i1+inverseGridDim1*(i2))]

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

#define MASK1	(		i0,
			i1,
			i2
	)		mask1p[i0+mask1Dim0*(i1+mask1Dim1*(i2))]

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

#define MASK2	(		i0,
			i1,
			i2
	)		mask2p[i0+mask2Dim0*(i1+mask2Dim1*(i2))]

#define MASK2_W4A	(		k1,
			k2,
			k3
	)		MASK2(k1,k2,k3)!=0 && MASK2(k1+1,k2,k3)!=0 && MASK2(k1+2,k2,k3)!=0 && MASK2(k1+3,k2,k3)!=0

#define MASK2_W4B	(		k1,
			k2,
			k3
	)		MASK2_W4A(k1,k2,k3) && MASK2_W4A(k1,k2+1,k3) && MASK2_W4A(k1,k2+2,k3) && MASK2_W4A(k1,k2+3,k3)

#define MASK2_W4C	(		k1,
			k2,
			k3
	)		MASK2_W4B(k1,k2,k3) && MASK2_W4B(k1,k2,k3+1) && MASK2_W4B(k1,k2,k3+2) && MASK2_W4B(k1,k2,k3+3)

#define MASK2_W5A	(		k1,
			k2,
			k3
	)		MASK2(k1,k2,k3)!=0 && MASK2(k1+1,k2,k3)!=0 && MASK2(k1+2,k2,k3)!=0 && MASK2(k1+3,k2,k3)!=0 && MASK2(k1+4,k2,k3)!=0

#define MASK2_W5B	(		k1,
			k2,
			k3
	)

Value:

MASK2_W5A(k1,k2,k3) && MASK2_W5A(k1,k2+1,k3) && MASK2_W5A(k1,k2+2,k3) && \
  MASK2_W5A(k1,k2+3,k3) && MASK2_W5A(k1,k2+4,k3)

#define MASK2_W5C	(		k1,
			k2,
			k3
	)

Value:

MASK2_W5B(k1,k2,k3) && MASK2_W5B(k1,k2,k3+1) && MASK2_W5B(k1,k2,k3+2) && \
  MASK2_W5B(k1,k2,k3+3) && MASK2_W5B(k1,k2,k3+4)

#define MASK2D	(		k1,
			k2,
			k3
	)		(MASK2(k1,k2,k3)>0 ? 1 : MASK2(k1,k2,k3)==0 ? 0 : -1)

#define NG

(

i0

)

ngp[i0]

#define R	(		i0,
			i1
	)		rp[i0+rDim0*(i1)]

#define RB	(		i0,
			i1
	)		rbp[i0+rbDim0*(i1)]

#define VARIABLEINTERPOLATIONWIDTH

(

i0

)

variableInterpolationWidthp[i0]

#define VARIABLEINTERPOLATIONWIDTH1

(

i0

)

variableInterpolationWidth1p[i0]

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

#define XR	(		m,
			n
	)		xra[n][m]

#define XR	(		m,
			n
	)		xra[n][m]

Function Documentation

int checkGrid	(	CompositeGrid &	cg,
		GenericGraphicsInterface *	ps = 0,
		int	debug = 0
	)

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

int displayMaskLaTeX	(	const intArray &	mask,
		const aString &	label = nullString,
		FILE *	file = NULL,
		const DisplayParameters *	displayParameters = NULL
	)

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

int Ogmg:: markGhostPoints

(

CompositeGrid &

cg

)

References MappedGrid::boundaryCondition(), MappedGrid::boundaryFlag(), MappedGrid::dimension(), End, MappedGrid::extendedIndexRange(), FOR_IJ, getBoundaryIndex(), getGhostIndex(), ParallelUtility::getLocalArrayBounds(), getLocalArrayWithGhostBoundaries(), MappedGrid::gridIndexRange(), i2, MappedGrid::ISghostPoint, mask, MASK, MappedGrid::mask(), MappedGrid::mixedPhysicalInterpolationBoundary, GenericGridCollection::numberOfComponentGrids(), intMappedGridFunction::periodicUpdate(), MappedGrid::physicalBoundary, and Start.