updateRefinementNew.C File Reference

#include "Overture.h"
#include "Ogen.h"
#include "display.h"
#include "HDF_DataBase.h"
#include "ParallelUtility.h"
#include "CanInterpolate.h"

Include dependency graph for updateRefinementNew.C:

Macros
#define	FOR_3D(i1, i2, i3, I1, I2, I3)   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++)
#define	FOR_3(i1, i2, i3, I1, I2, I3)   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_3DS(i1, i2, i3, I1, I2, I3)   int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); int I1Stride=I1.getStride(), I2Stride=I2.getStride(), I3Stride=I3.getStride(); for(i3=I3Base; i3<=I3Bound; i3+=I3Stride ) for(i2=I2Base; i2<=I2Bound; i2+=I2Stride ) for(i1=I1Base; i1<=I1Bound; i1+=I1Stride )
#define	FOR_3S(i1, i2, i3, I1, I2, I3)   I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); I1Stride=I1.getStride(), I2Stride=I2.getStride(), I3Stride=I3.getStride(); for(i3=I3Base; i3<=I3Bound; i3+=I3Stride ) for(i2=I2Base; i2<=I2Bound; i2+=I2Stride ) for(i1=I1Base; i1<=I1Bound; i1+=I1Stride )
#define	FOR_3R()   I1rBase =I1r.getBase(), I2rBase =I2r.getBase(), I3rBase =I3r.getBase(); I1rBound =I1r.getBound(), I2rBound =I2r.getBound(), I3rBound =I3r.getBound(); I1rStride=I1r.getStride(),I2rStride=I2r.getStride(),I3rStride=I3r.getStride(); for(i3r=I3rBase; i3r<=I3rBound; i3r+=I3rStride) for(i2r=I2rBase; i2r<=I2rBound; i2r+=I2rStride) for(i1r=I1rBase; i1r<=I1rBound; i1r+=I1rStride)
#define	FOR_3BR()   I1bBase =I1b.getBase(), I2bBase =I2b.getBase(), I3bBase =I3b.getBase(); I1bBound=I1b.getBound(), I2bBound=I2b.getBound(), I3bBound=I3b.getBound(); I1rBase =I1r.getBase(), I2rBase =I2r.getBase(), I3rBase =I3r.getBase(); I1rBound =I1r.getBound(), I2rBound =I2r.getBound(), I3rBound =I3r.getBound(); I1rStride=I1r.getStride(),I2rStride=I2r.getStride(),I3rStride=I3r.getStride(); for(i3b=I3bBase,i3r=I3rBase; i3b<=I3bBound; i3b++,i3r+=I3rStride) for(i2b=I2bBase,i2r=I2rBase; i2b<=I2bBound; i2b++,i2r+=I2rStride) for(i1b=I1bBase,i1r=I1rBase; i1b<=I1bBound; i1b++,i1r+=I1rStride)
#define	FOR_3IJD(i1, i2, i3, I1, I2, I3, j1, j2, j3, J1, J2, J3)   int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); int J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(); 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++)
#define	MASK(i0, i1, i2)   maskp[i0+maskDim0*(i1+maskDim1*(i2))]
#define	MASKB(i0, i1, i2)   maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]
#define	MASKR(i0, i1, i2)   maskrp[i0+maskrDim0*(i1+maskrDim1*(i2))]
#define	moda(n, m)   (n)>0 ? (n)%(m) : ((m)-(n))%(m)
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	VERTEX0(i0, i1, i2)   xab[0][0]+dx[0]*(i0-iv0[0])
#define	VERTEX1(i0, i1, i2)   xab[0][1]+dx[1]*(i1-iv0[1])
#define	VERTEX2(i0, i1, i2)   xab[0][2]+dx[2]*(i2-iv0[2])
#define	MASKB(i0, i1, i2)   maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]
#define	IPBG(i0, i1)   ipbgp[i0+ipbgDim0*(i1)]
#define	interpolationPointBaseGrid(i, bg)   interpolationPointBaseGridA[g][(bg)+maximumNumberOfPossibleBaseGrids*(i)]
#define	interpoleeBaseGrid(i, bg)   interpoleeBaseGridA[g][(bg)+maximumNumberOfPossibleBaseGrids*(i)]
#define	numberOfPossibleInterpoleeBaseGrids(i)   numberOfPossibleInterpoleeBaseGridsA[g][i]
#define	interpCoords(i, bg, dir)   interpCoordsA[g][(bg)+maximumNumberOfPossibleBaseGrids*((i)+ni*(dir))]
#define	MASKB(i0, i1, i2)   maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	IPBG(i0, i1)   ipbgp[i0+ipbgDim0*(i1)]
#define	GIR(side, axis)   girp[side+2*(axis)]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	IP(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	MASKB(i0, i1, i2)   maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]
#define	MASK(i0, i1, i2)   maskp[i0+maskDim0*(i1+maskDim1*(i2))]
#define	MASKB(i0, i1, i2)   maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]

Functions
int	redistribute (const intArray &u, intArray &v, const Range &P)
int	redistribute (const intArray &u, intSerialArray &v)
int	outputRefinementInfoNew (GridCollection &gc, const aString &gridFileName, const aString &fileName)

Macro Definition Documentation

#define FOR_3	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)		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_3BR

(

)

I1bBase =I1b.getBase(), I2bBase =I2b.getBase(), I3bBase =I3b.getBase(); I1bBound=I1b.getBound(), I2bBound=I2b.getBound(), I3bBound=I3b.getBound(); I1rBase =I1r.getBase(), I2rBase =I2r.getBase(), I3rBase =I3r.getBase(); I1rBound =I1r.getBound(), I2rBound =I2r.getBound(), I3rBound =I3r.getBound(); I1rStride=I1r.getStride(),I2rStride=I2r.getStride(),I3rStride=I3r.getStride(); for(i3b=I3bBase,i3r=I3rBase; i3b<=I3bBound; i3b++,i3r+=I3rStride) for(i2b=I2bBase,i2r=I2rBase; i2b<=I2bBound; i2b++,i2r+=I2rStride) for(i1b=I1bBase,i1r=I1rBase; i1b<=I1bBound; i1b++,i1r+=I1rStride)

Referenced by Ogen::markOffAxisRefinementMask(), and Ogen::updateRefinementNew().

#define FOR_3D	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)		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 Ogen::updateRefinementNew().

#define FOR_3DS	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)		int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); int I1Stride=I1.getStride(), I2Stride=I2.getStride(), I3Stride=I3.getStride(); for(i3=I3Base; i3<=I3Bound; i3+=I3Stride ) for(i2=I2Base; i2<=I2Bound; i2+=I2Stride ) for(i1=I1Base; i1<=I1Bound; i1+=I1Stride )

#define FOR_3IJD	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3,
			j1,
			j2,
			j3,
			J1,
			J2,
			J3
	)		int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); int J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(); 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 Ogen::updateRefinementNew().

#define FOR_3R

(

)

I1rBase =I1r.getBase(), I2rBase =I2r.getBase(), I3rBase =I3r.getBase(); I1rBound =I1r.getBound(), I2rBound =I2r.getBound(), I3rBound =I3r.getBound(); I1rStride=I1r.getStride(),I2rStride=I2r.getStride(),I3rStride=I3r.getStride(); for(i3r=I3rBase; i3r<=I3rBound; i3r+=I3rStride) for(i2r=I2rBase; i2r<=I2rBound; i2r+=I2rStride) for(i1r=I1rBase; i1r<=I1rBound; i1r+=I1rStride)

Referenced by Ogen::markOffAxisRefinementMask().

#define FOR_3S	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)		I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); I1Stride=I1.getStride(), I2Stride=I2.getStride(), I3Stride=I3.getStride(); for(i3=I3Base; i3<=I3Bound; i3+=I3Stride ) for(i2=I2Base; i2<=I2Bound; i2+=I2Stride ) for(i1=I1Base; i1<=I1Bound; i1+=I1Stride )

Referenced by Ogen::updateRefinementNew().

#define GIR	(		side,
			axis
	)		girp[side+2*(axis)]

#define interpCoords	(		i,
			bg,
			dir
	)		interpCoordsA[g][(bg)+maximumNumberOfPossibleBaseGrids*((i)+ni*(dir))]

#define interpolationPointBaseGrid	(		i,
			bg
	)		interpolationPointBaseGridA[g][(bg)+maximumNumberOfPossibleBaseGrids*(i)]

#define interpoleeBaseGrid	(		i,
			bg
	)		interpoleeBaseGridA[g][(bg)+maximumNumberOfPossibleBaseGrids*(i)]

#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 IPBG	(		i0,
			i1
	)		ipbgp[i0+ipbgDim0*(i1)]

#define IPBG	(		i0,
			i1
	)		ipbgp[i0+ipbgDim0*(i1)]

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

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

#define MASKB	(		i0,
			i1,
			i2
	)		maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]

#define MASKB	(		i0,
			i1,
			i2
	)		maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]

#define MASKB	(		i0,
			i1,
			i2
	)		maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]

#define MASKB	(		i0,
			i1,
			i2
	)		maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]

#define MASKB	(		i0,
			i1,
			i2
	)		maskbp[i0+maskbDim0*(i1+maskbDim1*(i2))]

#define MASKR	(		i0,
			i1,
			i2
	)		maskrp[i0+maskrDim0*(i1+maskrDim1*(i2))]

#define moda	(		n,
			m
	)		(n)>0 ? (n)%(m) : ((m)-(n))%(m)

#define numberOfPossibleInterpoleeBaseGrids

(

i

)

numberOfPossibleInterpoleeBaseGridsA[g][i]

#define VERTEX0	(		i0,
			i1,
			i2
	)		xab[0][0]+dx[0]*(i0-iv0[0])

#define VERTEX1	(		i0,
			i1,
			i2
	)		xab[0][1]+dx[1]*(i1-iv0[1])

#define VERTEX2	(		i0,
			i1,
			i2
	)		xab[0][2]+dx[2]*(i2-iv0[2])

Function Documentation

int outputRefinementInfoNew	(	GridCollection &	gc,
		const aString &	gridFileName,
		const aString &	fileName
	)

References assert, aString, GenericGridCollection::baseGridNumber, c, GenericGridCollectionData::computedGeometry, CompositeGrid::deleteRefinementLevels(), CompositeGrid::destroy(), CompositeGrid::EVERYTHING, MappedGrid::EVERYTHING, CompositeGrid::get(), GridCollection::getClassName(), MappedGrid::gridIndexRange(), HDF_DataBase::mount(), GenericDataBase::normalMode, GenericDataBase::noStreamMode, GenericGridCollection::numberOfComponentGrids(), GridCollection::numberOfDimensions(), GenericGridCollection::numberOfGrids(), CompositeGrid::numberOfInterpolationPoints, HDF_DataBase::put(), CompositeGrid::put(), GridCollection::refinementFactor, GenericGridCollection::refinementLevelNumber, HDF_DataBase::setMode(), CompositeGrid::THEinterpolationCoordinates, CompositeGrid::THEinterpolationPoint, CompositeGrid::THEinterpoleeGrid, CompositeGrid::THEinterpoleeLocation, GridCollection::THEmask, MappedGrid::THEmask, HDF_DataBase::unmount(), and CompositeGrid::updateReferences().

Referenced by Ogen::updateRefinementFillInterpolationData(), Ogen::updateRefinementFillInterpolationDataNew(), Ogen::updateRefinementNew(), and Ogen::updateRefinementNewer().

int redistribute	(	const intArray &	u,
		intArray &	v,
		const Range &	P
	)

References NULL.

Referenced by Ogen::updateRefinementNew().

int redistribute	(	const intArray &	u,
		intSerialArray &	v
	)

References CopyArray::copyArray(), max(), and IndexBox::setBounds().