Integrate.C File Reference

#include "Integrate.h"
#include "GenericGraphicsInterface.h"
#include "Oges.h"
#include "CompositeGridOperators.h"
#include "SparseRep.h"
#include "UnstructuredMapping.h"
#include "Geom.h"
#include "InterpolatePoints.h"
#include "display.h"
#include "SurfaceStitcher.h"
#include "ParallelUtility.h"
#include "BoxLib.H"
#include "Box.H"
#include "BoxList.H"

Include dependency graph for Integrate.C:

Macros
#define	FOR_3D(i1, i2, i3, I1, I2, I3)
#define	FOR_3IJD(i1, i2, i3, I1, I2, I3, j1, j2, j3, J1, J2, J3)
#define	numberOfBoundaryGrids(surface, level, face)   pNumberOfBoundaryGrids[surface][(level)+(numberOfBoundaryRefinementLevels)*(face)]
#define	boundaryGrid(surface, level, face, g)   pBoundaryGrid[surface][(level)+(numberOfBoundaryRefinementLevels)*(face)][g]
#define	numberOfBoundaryGrids(surface, level, face)   pNumberOfBoundaryGrids[surface][(level)+(numberOfRefinementLevels)*(face)]
#define	boundaryGrid(surface, level, face, g)   pBoundaryGrid[surface][(level)+(numberOfRefinementLevels)*(face)][g]
#define	bWeights(surface, level, face, grid)   boundaryWeights[surface][(level)+(numberOfRefinementLevels)*(face)][grid]
#define	RXB(m, n)   rxLocal(Ib1,Ib2,Ib3,m+numberOfDimensions*(n))
#define	WEIGHTS(i0, i1, i2)   weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]
#define	WEIGHTS(i0, i1, i2)   weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]
#define	U(i0, i1, i2, i3)   up[i0+uDim0*(i1+uDim1*(i2+uDim2*(i3)))]
#define	RX(m, n)   rx(I1,I2,I3,m+numberOfDimensions*(n))
#define	RXLocal(m, n)   rxLocal(I1,I2,I3,m+numberOfDimensions*(n))
#define	GID(side, axis)   (pgid[(side)+2*(axis)])
#define	COEFF(i0, i1, i2, i3)   coeffp[i0+coeffDim0*(i1+coeffDim1*(i2+coeffDim2*(i3)))]
#define	MASK(i0, i1, i2)   maskp[(i0)+(i1)*md1+(i2)*md2]
#define	WEIGHTS(i0, i1, i2)   weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]
#define	CENTER(i0, i1, i2, i3)   centerp[i0+centerDim0*(i1+centerDim1*(i2+centerDim2*(i3)))]
#define	F(i0, i1, i2)   fp[i0+fDim0*(i1+fDim1*(i2))]
#define	JACOBIAN(i0, i1, i2)   jacobianp[i0+jacobianDim0*(i1+jacobianDim1*(i2))]
#define	CLASSIFY(i0, i1, i2)   classifyp[i0+classifyDim0*(i1+classifyDim1*(i2))]
#define	M2(m1, m2)   ((m1)+halfWidth+width*((m2)+halfWidth))
#define	M3(m1, m2, m3)   ((m1)+halfWidth+width*((m2)+halfWidth+width*((m3)+halfWidth)))
#define	WEIGHTS(i0, i1, i2)   weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]
#define	CLASSIFY(i0, i1, i2)   classifyp[i0+classifyDim0*(i1+classifyDim1*(i2))]
#define	RX(m, n)   rx(I1,I2,I3,m+numberOfDimensions*(n))
#define	KRX(m, n)   m+3*(n)
#define	FOR_3D(i1, i2, i3, I1, I2, I3)
#define	FOR_3(i1, i2, i3, I1, I2, I3)
#define	FOR_3IJD(i1, i2, i3, I1, I2, I3, j1, j2, j3, J1, J2, J3)
#define	KRX(m, n)   m+3*(n)
#define	AREA1(i1, i2, i3)
#define	AREA2(i1, i2, i3)
#define	AREA3(i1, i2, i3)
#define	KRX(m, n)   m+3*(n)

Macro Definition Documentation

#define AREA1	(		i1,
			i2,
			i3
	)

Value:

( SQRT( SQR(rx(i1,i2,i3,KRX(axis,axis1))) + \
                                SQR(rx(i1,i2,i3,KRX(axis,axis2))) + \
                                SQR(rx(i1,i2,i3,KRX(axis,axis3))) )* \
                          fabs(jacobian(i1,i2,i3))*gridSpacing(axis2)*gridSpacing(axis3) )

Referenced by Integrate::computeStitchedSurfaceWeights().

#define AREA2	(		i1,
			i2,
			i3
	)

Value:

( SQRT( SQR(rx(i1,i2,i3,KRX(axis,axis1))) + \
                                SQR(rx(i1,i2,i3,KRX(axis,axis2))) + \
                                SQR(rx(i1,i2,i3,KRX(axis,axis3))) )* \
                          fabs(jacobian(i1,i2,i3))*gridSpacing(axis1)*gridSpacing(axis3) )

Referenced by Integrate::computeStitchedSurfaceWeights().

#define AREA3	(		i1,
			i2,
			i3
	)

Value:

( SQRT( SQR(rx(i1,i2,i3,KRX(axis,axis1))) + \
                                SQR(rx(i1,i2,i3,KRX(axis,axis2))) + \
                                SQR(rx(i1,i2,i3,KRX(axis,axis3))) )* \
                          fabs(jacobian(i1,i2,i3))*gridSpacing(axis1)*gridSpacing(axis2) )

Referenced by Integrate::computeStitchedSurfaceWeights().

#define boundaryGrid	(		surface,
			level,
			face,
			g
	)		pBoundaryGrid[surface][(level)+(numberOfBoundaryRefinementLevels)*(face)][g]

Referenced by Integrate::buildAdaptiveMeshRefinementSurfaceWeights(), Integrate::getFace(), and Integrate::surfaceIntegral().

#define boundaryGrid	(		surface,
			level,
			face,
			g
	)		pBoundaryGrid[surface][(level)+(numberOfRefinementLevels)*(face)][g]

#define bWeights	(		surface,
			level,
			face,
			grid
	)		boundaryWeights[surface][(level)+(numberOfRefinementLevels)*(face)][grid]

Referenced by Integrate::buildAdaptiveMeshRefinementSurfaceWeights(), and Integrate::surfaceIntegral().

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

#define CLASSIFY	(		i0,
			i1,
			i2
	)		classifyp[i0+classifyDim0*(i1+classifyDim1*(i2))]

#define CLASSIFY	(		i0,
			i1,
			i2
	)		classifyp[i0+classifyDim0*(i1+classifyDim1*(i2))]

#define COEFF	(		i0,
			i1,
			i2,
			i3
	)		coeffp[i0+coeffDim0*(i1+coeffDim1*(i2+coeffDim2*(i3)))]

#define F	(		i0,
			i1,
			i2
	)		fp[i0+fDim0*(i1+fDim1*(i2))]

#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++)

#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++)

#define FOR_3IJD	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3,
			j1,
			j2,
			j3,
			J1,
			J2,
			J3
	)

Value:

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 FOR_3IJD	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3,
			j1,
			j2,
			j3,
			J1,
			J2,
			J3
	)

Value:

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 GID	(		side,
			axis
	)		(pgid[(side)+2*(axis)])

Referenced by Integrate::computeWeights().

#define JACOBIAN	(		i0,
			i1,
			i2
	)		jacobianp[i0+jacobianDim0*(i1+jacobianDim1*(i2))]

Referenced by Integrate::computeWeights().

#define KRX	(		m,
			n
	)		m+3*(n)

Referenced by Integrate::computeStitchedSurfaceWeightsOld().

#define KRX	(		m,
			n
	)		m+3*(n)

#define KRX	(		m,
			n
	)		m+3*(n)

#define M2	(		m1,
			m2
	)		((m1)+halfWidth+width*((m2)+halfWidth))

#define M3	(		m1,
			m2,
			m3
	)		((m1)+halfWidth+width*((m2)+halfWidth+width*((m3)+halfWidth)))

#define MASK	(		i0,
			i1,
			i2
	)		maskp[(i0)+(i1)*md1+(i2)*md2]

#define numberOfBoundaryGrids	(		surface,
			level,
			face
	)		pNumberOfBoundaryGrids[surface][(level)+(numberOfBoundaryRefinementLevels)*(face)]

Referenced by Integrate::buildAdaptiveMeshRefinementSurfaceWeights(), Integrate::getFace(), and Integrate::surfaceIntegral().

#define numberOfBoundaryGrids	(		surface,
			level,
			face
	)		pNumberOfBoundaryGrids[surface][(level)+(numberOfRefinementLevels)*(face)]

#define RX	(		m,
			n
	)		rx(I1,I2,I3,m+numberOfDimensions*(n))

#define RX	(		m,
			n
	)		rx(I1,I2,I3,m+numberOfDimensions*(n))

#define RXB	(		m,
			n
	)		rxLocal(Ib1,Ib2,Ib3,m+numberOfDimensions*(n))

Referenced by Integrate::buildAdaptiveMeshRefinementSurfaceWeights().

#define RXLocal	(		m,
			n
	)		rxLocal(I1,I2,I3,m+numberOfDimensions*(n))

Referenced by Integrate::computeSurfaceWeights(), and Integrate::computeWeights().

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

#define WEIGHTS	(		i0,
			i1,
			i2
	)		weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]

Referenced by Integrate::computeWeights(), Integrate::volume(), and Integrate::volumeIntegral().

#define WEIGHTS	(		i0,
			i1,
			i2
	)		weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]

#define WEIGHTS	(		i0,
			i1,
			i2
	)		weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]

#define WEIGHTS	(		i0,
			i1,
			i2
	)		weightsp[i0+weightsDim0*(i1+weightsDim1*(i2))]