InterfaceTransfer.C File Reference

#include "InterfaceTransfer.h"
#include "DomainSolver.h"
#include "Interface.h"
#include "ParallelUtility.h"
#include "InterpolatePointsOnAGrid.h"

Include dependency graph for InterfaceTransfer.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	gridIsOnInterface(grid)   pGridIsOnInterface[grid]
#define	q10(x)   1
	Make changes to the interface transfer parameters.
#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	IL(i0, i1)   ipp[i0+ipDim0*(i1)]
#define	VIW(i)   viw(i)
#define	IL1(i)   IL(i,0)
#define	IL2(i)   IL(i,1)
#define	IL3(i)   IL(i,2)
#define	CI(i0, i1)   cip[i0+ciDim0*(i1)]
#define	IA(i0)   iap[(i0)]
#define	cf(i, m1, m2)   coeffap[i+nil*((m1)+width[0]*((m2)))]
#define	cfs(i, m1)   coeffap[i+nil*(m1)]
#define	ID(i0, i1)   idp[i0+idDim0*(i1)]
#define	MASK(i0, i1, i2)   maskp[i0+maskDim0*(i1+maskDim1*(i2))]
#define	VS(i0, i1, i2, i3)   vsp[i0+vsDim0*(i1+vsDim1*(i2+vsDim2*(i3)))]
#define	IW5(m1)   (cf(j,0,m1)*VS(i1 ,i2+m1,i3,c1)+cf(j,1,m1)*VS(i1+1,i2+m1,i3,c1)+cf(j,2,m1)*VS(i1+2,i2+m1,i3,c1)+cf(j,3,m1)*VS(i1+3,i2+m1,i3,c1)+cf(j,4,m1)*VS(i1+4,i2+m1,i3,c1))
#define	c(i, m1, m2, m3)   coeffgp[m1+width[0]*(m2+width[1]*m3)][i]
#define	VS(i0, i1, i2, i3)   vsp[i0+vsDim0*(i1+vsDim1*(i2+vsDim2*(i3)))]
#define	c(i, m1, m2, m3)   coeffap[i+nil*((m1)+width[0]*((m2)+width[1]*(m3)))]
#define	IW5A(m1, m2)   (c(j,0,m1,m2)*VS(i1 ,i2+m1,i3+m2,c1)+c(j,1,m1,m2)*VS(i1+1,i2+m1,i3+m2,c1)+c(j,2,m1,m2)*VS(i1+2,i2+m1,i3+m2,c1)+c(j,3,m1,m2)*VS(i1+3,i2+m1,i3+m2,c1)+c(j,4,m1,m2)*VS(i1+4,i2+m1,i3+m2,c1))
#define	IW5(m2)   (IW5A(0,m2)+IW5A(1,m2)+IW5A(2,m2)+IW5A(3,m2)+IW5A(4,m2))

Macro Definition Documentation

#define c	(		i,
			m1,
			m2,
			m3
	)		coeffgp[m1+width[0]*(m2+width[1]*m3)][i]

Referenced by DomainSolver::addArtificialDissipation(), Cgcns::addConstraintEquation(), Cgad::addForcing(), DomainSolver::addGrids(), DomainSolver::advanceAdamsPredictorCorrector(), Cgad::advanceADI(), DomainSolver::advanceImplicitMultiStep(), Cgins::advanceLineSolveNew(), Cgins::advanceLineSolveOld(), DomainSolver::advanceSecondOrderSystem(), DomainSolver::advanceVariableTimeStepAdamsPredictorCorrector(), Cgasf::allSpeedImplicitTimeStep(), Cgad::applyBoundaryConditions(), Cgins::applyBoundaryConditionsForImplicitTimeStepping(), applyFourthOrderBoundaryConditionsOld(), Cgsm::assignAnnulusEigenfunction(), Cgsm::assignGaussianPulseInitialConditions(), Cgsm::assignHempInitialConditions(), DomainSolver::assignInitialConditions(), Cgins::assignPressureRHS(), Cgsm::assignTwilightZoneInitialConditions(), checkGridFunction(), DomainSolver::checkSolution(), Cgsm::computeDissipation(), computeDSIErrors(), convertToVertexCentered(), MaterialProperties::cp(), detectCollisions(), Cgins::determineErrors(), DomainSolver::determineErrors(), AdParameters::displayPdeParameters(), F90_ID(), for(), Cgasf::formAllSpeedPressureEquation(), Cgcns::formImplicitTimeSteppingMatrix(), Cgsm::getEnergy(), Cgsm::getErrors(), DomainSolver::getGridInfo(), MovingGrids::getGridVelocity(), DomainSolver::getInitialConditions(), Cgmp::getInterfaceResiduals(), Cgsm::getMaxDivAndCurl(), DomainSolver::getTimeDependentBoundaryConditions(), Cgsm::getTimeStep(), Cgad::getTimeSteppingEigenvalue(), Cgad::getUt(), Cgcns::getUt(), Cgasf::gridAccelerationBC(), Cgins::gridAccelerationBC(), MaterialProperties::h(), if(), initialConditionsFromAProfile(), DomainSolver::initializeTimeSteppingIM(), DomainSolver::initializeTimeSteppingPC(), InterfaceTransfer::internalInterpolate(), DomainSolver::jetInflow(), MaterialProperties::lnKp(), main(), mult(), DomainSolver::output(), DomainSolver::outputHeader(), DomainSolver::outputProbes(), Cgsm::outputResults(), Cgins::outputSolution(), DomainSolver::outputSolution(), DomainSolver::parabolicInflow(), ProjectVelocity::projectVelocity(), Parameters::readFromAShowFile(), DeformingBodyMotion::regenerateComponentGrids(), MovingGrids::rigidBodyMotion(), MaterialProperties::s(), Controller::saveControlSequenceData(), AdParameters::saveParametersToShowFile(), Cgad::saveShowFileComments(), Controller::saveToShowFile(), Parameters::setBoundaryConditionValues(), AdParameters::setDefaultDataForABoundaryCondition(), AdParameters::setParameters(), SmParameters::setParameters(), AdParameters::setPdeParameters(), Cgad::setPlotTitle(), Cgsm::setupGridFunctions(), DomainSolver::setupGridFunctions(), Cgasf::solveForAllSpeedPressure(), Cgasf::solveForTimeIndependentVariables(), DomainSolver::takeTimeStepAF(), DomainSolver::takeTimeStepIM(), DomainSolver::takeTimeStepPC(), trans(), MovingGrids::update(), Cgins::updatePressureEquation(), DomainSolver::updateToMatchNewGrid(), DomainSolver::userDefinedBoundaryValues(), userDefinedErrorEstimator(), DomainSolver::userDefinedGrid(), Cgad::userDefinedInitialConditions(), Cgsm::userDefinedInitialConditions(), DomainSolver::userDefinedInitialConditions(), DomainSolver::userDefinedMaterialProperties(), and Cgad::writeParameterSummary().

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

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

Referenced by DomainSolver::advanceTrapezoidal(), InterfaceTransfer::internalInterpolate(), and Maxwell::setupDSICoefficients().

#define cfs	(		i,
			m1
	)		coeffap[i+nil*(m1)]

Referenced by InterfaceTransfer::internalInterpolate().

#define CI	(		i0,
			i1
	)		cip[i0+ciDim0*(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_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 gridIsOnInterface

(

grid

)

pGridIsOnInterface[grid]

Referenced by InterfaceTransfer::initialize().

#define IA

(

i0

)

iap[(i0)]

#define ID	(		i0,
			i1
	)		idp[i0+idDim0*(i1)]

Referenced by InterfaceTransfer::internalInterpolate().

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

#define IL1

(

i

)

IL(i,0)

Referenced by InterfaceTransfer::internalInterpolate().

#define IL2

(

i

)

IL(i,1)

Referenced by InterfaceTransfer::internalInterpolate().

#define IL3

(

i

)

IL(i,2)

Referenced by InterfaceTransfer::internalInterpolate().

#define IW5

(

m1

)

(cf(j,0,m1)*VS(i1 ,i2+m1,i3,c1)+cf(j,1,m1)*VS(i1+1,i2+m1,i3,c1)+cf(j,2,m1)*VS(i1+2,i2+m1,i3,c1)+cf(j,3,m1)*VS(i1+3,i2+m1,i3,c1)+cf(j,4,m1)*VS(i1+4,i2+m1,i3,c1))

Referenced by InterfaceTransfer::internalInterpolate().

#define IW5

(

m2

)

(IW5A(0,m2)+IW5A(1,m2)+IW5A(2,m2)+IW5A(3,m2)+IW5A(4,m2))

#define IW5A	(		m1,
			m2
	)		(c(j,0,m1,m2)*VS(i1 ,i2+m1,i3+m2,c1)+c(j,1,m1,m2)*VS(i1+1,i2+m1,i3+m2,c1)+c(j,2,m1,m2)*VS(i1+2,i2+m1,i3+m2,c1)+c(j,3,m1,m2)*VS(i1+3,i2+m1,i3+m2,c1)+c(j,4,m1,m2)*VS(i1+4,i2+m1,i3+m2,c1))

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

#define q10

(

x

)

1

Make changes to the interface transfer parameters.

#define q20

(

x

)

-(x)+1

Referenced by InterfaceTransfer::internalInterpolate().

#define q21

(

x

)

(x)

Referenced by InterfaceTransfer::internalInterpolate().

#define q30

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q31

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q32

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q40

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q41

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q42

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q43

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q50

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q51

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q52

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q53

(

x

)

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

Referenced by InterfaceTransfer::internalInterpolate().

#define q54

(

x

)

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

Referenced by InterfaceTransfer::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.

#define VIW

(

i

)

viw(i)

Referenced by InterfaceTransfer::internalInterpolate().

#define VS	(		i0,
			i1,
			i2,
			i3
	)		vsp[i0+vsDim0*(i1+vsDim1*(i2+vsDim2*(i3)))]

Referenced by InterfaceTransfer::internalInterpolate(), and DomainSolver::userDefinedInitialConditions().

#define VS	(		i0,
			i1,
			i2,
			i3
	)		vsp[i0+vsDim0*(i1+vsDim1*(i2+vsDim2*(i3)))]