mx/src/assignBoundaryConditions.C File Reference

#include "Maxwell.h"
#include "CompositeGridOperators.h"
#include "display.h"
#include "UnstructuredMapping.h"
#include "OGPolyFunction.h"
#include "OGTrigFunction.h"
#include "OGPulseFunction.h"
#include "RadiationBoundaryCondition.h"
#include "ParallelUtility.h"

Include dependency graph for mx/src/assignBoundaryConditions.C:

Macros
#define	bcOptMaxwell   EXTERN_C_NAME(bcoptmaxwell)
#define	abcMaxwell   EXTERN_C_NAME(abcmaxwell)
#define	bcSymmetry   EXTERN_C_NAME(bcsymmetry)
#define	exmax   EXTERN_C_NAME(exmax)
#define	adjustForIncident   EXTERN_C_NAME(adjustforincident)
#define	pmlMaxwell   EXTERN_C_NAME(pmlmaxwell)
#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	exTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[0]
#define	eyTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[1]
#define	hzTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[5]
#define	extTrue(x, y, t)   (-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[0]
#define	eytTrue(x, y, t)   (-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[1]
#define	hztTrue(x, y, t)   (-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[5]
#define	exLaplacianTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky))*pwc[0])
#define	eyLaplacianTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky))*pwc[1])
#define	hzLaplacianTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky))*pwc[5])
#define	hzGaussianPulse(xi)   exp(-betaGaussianPlaneWave*((xi)*(xi)))
#define	exGaussianPulse(xi)   hzGaussianPulse(xi)*(-ky/(eps*cc))
#define	eyGaussianPulse(xi)   hzGaussianPulse(xi)*( kx/(eps*cc))
#define	hzLaplacianGaussianPulse(xi)   ((4.*betaGaussianPlaneWave*betaGaussianPlaneWave*(kx*kx+ky*ky))*xi*xi-(2.*betaGaussianPlaneWave*(kx*kx+ky*ky)))*exp(-betaGaussianPlaneWave*((xi)*(xi)))
#define	exLaplacianGaussianPulse(xi)   hzLaplacianGaussianPulse(xi,t)*(-ky/(eps*cc))
#define	eyLaplacianGaussianPulse(xi)   hzLaplacianGaussianPulse(xi,t)*( kx/(eps*cc))
#define	exTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[0]
#define	eyTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[1]
#define	ezTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[2]
#define	extTrue3d(x, y, z, t)   (-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[0]
#define	eytTrue3d(x, y, z, t)   (-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[1]
#define	eztTrue3d(x, y, z, t)   (-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[2]
#define	hxTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[3]
#define	hyTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[4]
#define	hzTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[5]
#define	exLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[0])
#define	eyLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[1])
#define	ezLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[2])
#define	hxLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[3])
#define	hyLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[4])
#define	hzLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[5])
#define	U(i0, i1, i2, i3)   up[i0+uDim0*(i1+uDim1*(i2+uDim2*(i3)))]
#define	X(i0, i1, i2, i3)   xp[i0+xDim0*(i1+xDim1*(i2+xDim2*(i3)))]
#define	UG(i0, i1, i2, i3)   ugp[i0+ugDim0*(i1+ugDim1*(i2+ugDim2*(i3)))]
#define	PML(n, m, side, axis, grid)   vpml[(n)+numberOfTimeLevels*(m+numberOfPMLFunctions*(side+2*(axis+3*(grid))))]
#define	VPML(n, side, axis, grid)   PML(n,0,side,axis,grid)
#define	WPML(n, side, axis, grid)   PML(n,1,side,axis,grid)
#define	PML(n, m, side, axis, grid)   vpml[(n)+numberOfTimeLevels*(m+numberOfPMLFunctions*(side+2*(axis+3*(grid))))]
#define	VPML(n, side, axis, grid)   PML(n,0,side,axis,grid)
#define	WPML(n, side, axis, grid)   PML(n,1,side,axis,grid)

Functions
void	bcOptMaxwell (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &ndf1a, const int &ndf1b, const int &ndf2a, const int &ndf2b, const int &ndf3a, const int &ndf3b, const int &gid, const int &dimension, const real &u, const real &f, const int &mask, const real &rsxy, const real &xy, const int &bc, const int &boundaryCondition, const int &ipar, const real &rpar, int &ierr)
void	abcMaxwell (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &ndf1a, const int &ndf1b, const int &ndf2a, const int &ndf2b, const int &ndf3a, const int &ndf3b, const int &gid, const real &u, const real &un, const real &f, const int &mask, const real &rsxy, const real &xy, const int &bc, const int &boundaryCondition, const int &ipar, const real &rpar, int &ierr)
void	bcSymmetry (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gid, const real &u, const int &mask, const real &rsxy, const real &xy, const int &bc, const int &boundaryCondition, const int &ipar, const real &rpar, int &ierr)
void	adjustForIncident (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gid, const real &um, const real &u, const real &un, const int &mask, const real &rsxy, const real &xy, const real &initialConditionBoundingBox, const int &bc, const int &boundaryCondition, const int &ipar, const real &rpar, int &ierr)
void	exmax (double &Ez, double &Bx, double &By, const int &nsources, const double &xs, const double &ys, const double &tau, const double &var, const double &amp, const double &a, const double &x, const double &y, const double &time)
void	pmlMaxwell (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &ndf1a, const int &ndf1b, const int &ndf2a, const int &ndf2b, const int &ndf3a, const int &ndf3b, const int &gid, const real &um, const real &u, const real &un, const int &ndra1a, const int &ndra1b, const int &ndra2a, const int &ndra2b, const int &ndra3a, const int &ndra3b, const real &vram, const real &vra, const real &vran, const real &wram, const real &wra, const real &wran, const int &ndrb1a, const int &ndrb1b, const int &ndrb2a, const int &ndrb2b, const int &ndrb3a, const int &ndrb3b, const real &vrbm, const real &vrb, const real &vrbn, const real &wrbm, const real &wrb, const real &wrbn, const int &ndsa1a, const int &ndsa1b, const int &ndsa2a, const int &ndsa2b, const int &ndsa3a, const int &ndsa3b, const real &vsam, const real &vsa, const real &vsan, const real &wsam, const real &wsa, const real &wsan, const int &ndsb1a, const int &ndsb1b, const int &ndsb2a, const int &ndsb2b, const int &ndsb3a, const int &ndsb3b, const real &vsbm, const real &vsb, const real &vsbn, const real &wsbm, const real &wsb, const real &wsbn, const int &ndta1a, const int &ndta1b, const int &ndta2a, const int &ndta2b, const int &ndta3a, const int &ndta3b, const real &vtam, const real &vta, const real &vtan, const real &wtam, const real &wta, const real &wtan, const int &ndtb1a, const int &ndtb1b, const int &ndtb2a, const int &ndtb2b, const int &ndtb3a, const int &ndtb3b, const real &vtbm, const real &vtb, const real &vtbn, const real &wtbm, const real &wtb, const real &wtbn, const real &f, const int &mask, const real &rsxy, const real &xy, const int &bc, const int &boundaryCondition, const int &ipar, const real &rpar, int &ierr)
void	updateghostboundaries_ (realArray *&pu)
void	updateghostandperiodic_ (realMappedGridFunction *&pu)

Macro Definition Documentation

#define abcMaxwell   EXTERN_C_NAME(abcmaxwell)

Referenced by Maxwell::assignBoundaryConditions().

#define adjustForIncident   EXTERN_C_NAME(adjustforincident)

Referenced by Maxwell::assignBoundaryConditions().

#define bcOptMaxwell   EXTERN_C_NAME(bcoptmaxwell)

Referenced by Maxwell::assignBoundaryConditions().

#define bcSymmetry   EXTERN_C_NAME(bcsymmetry)

Referenced by Maxwell::assignBoundaryConditions().

#define exGaussianPulse

(

xi

)

hzGaussianPulse(xi)*(-ky/(eps*cc))

#define exLaplacianGaussianPulse

(

xi

)

hzLaplacianGaussianPulse(xi,t)*(-ky/(eps*cc))

#define exLaplacianTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky))*pwc[0])

#define exLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[0])

#define exmax   EXTERN_C_NAME(exmax)

#define exTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[0]

Referenced by Maxwell::assignBoundaryConditions().

#define exTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[0]

Referenced by Maxwell::assignBoundaryConditions().

#define extTrue	(		x,
			y,
			t
	)		(-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[0]

Referenced by Maxwell::assignBoundaryConditions().

#define extTrue3d	(		x,
			y,
			z,
			t
	)		(-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[0]

Referenced by Maxwell::assignBoundaryConditions().

#define eyGaussianPulse

(

xi

)

hzGaussianPulse(xi)*( kx/(eps*cc))

#define eyLaplacianGaussianPulse

(

xi

)

hzLaplacianGaussianPulse(xi,t)*( kx/(eps*cc))

#define eyLaplacianTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky))*pwc[1])

#define eyLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[1])

#define eyTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[1]

Referenced by Maxwell::assignBoundaryConditions().

#define eyTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[1]

Referenced by Maxwell::assignBoundaryConditions().

#define eytTrue	(		x,
			y,
			t
	)		(-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[1]

Referenced by Maxwell::assignBoundaryConditions().

#define eytTrue3d	(		x,
			y,
			z,
			t
	)		(-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[1]

Referenced by Maxwell::assignBoundaryConditions().

#define ezLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[2])

#define ezTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[2]

Referenced by Maxwell::assignBoundaryConditions().

#define eztTrue3d	(		x,
			y,
			z,
			t
	)		(-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[2]

Referenced by Maxwell::assignBoundaryConditions().

#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 hxLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[3])

#define hxTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[3]

Referenced by Maxwell::assignBoundaryConditions().

#define hyLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[4])

#define hyTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[4]

Referenced by Maxwell::assignBoundaryConditions().

#define hzGaussianPulse

(

xi

)

exp(-betaGaussianPlaneWave*((xi)*(xi)))

Referenced by Maxwell::assignBoundaryConditions().

#define hzLaplacianGaussianPulse

(

xi

)

((4.*betaGaussianPlaneWave*betaGaussianPlaneWave*(kx*kx+ky*ky))*xi*xi-(2.*betaGaussianPlaneWave*(kx*kx+ky*ky)))*exp(-betaGaussianPlaneWave*((xi)*(xi)))

#define hzLaplacianTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky))*pwc[5])

#define hzLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-(twoPi*twoPi*(kx*kx+ky*ky+kz*kz))*pwc[5])

#define hzTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[5]

Referenced by Maxwell::assignBoundaryConditions().

#define hzTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*pwc[5]

Referenced by Maxwell::assignBoundaryConditions().

#define hztTrue	(		x,
			y,
			t
	)		(-twoPi*cc)*cos(twoPi*(kx*(x)+ky*(y)-cc*(t)))*pwc[5]

Referenced by Maxwell::assignBoundaryConditions().

#define PML	(		n,
			m,
			side,
			axis,
			grid
	)		vpml[(n)+numberOfTimeLevels*(m+numberOfPMLFunctions*(side+2*(axis+3*(grid))))]

Referenced by Maxwell::assignBoundaryConditions().

#define PML	(		n,
			m,
			side,
			axis,
			grid
	)		vpml[(n)+numberOfTimeLevels*(m+numberOfPMLFunctions*(side+2*(axis+3*(grid))))]

#define pmlMaxwell   EXTERN_C_NAME(pmlmaxwell)

Referenced by Maxwell::assignBoundaryConditions().

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

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

Referenced by Maxwell::assignBoundaryConditions(), Maxwell::assignInitialConditions(), and Maxwell::getErrors().

#define VPML	(		n,
			side,
			axis,
			grid
	)		PML(n,0,side,axis,grid)

#define VPML	(		n,
			side,
			axis,
			grid
	)		PML(n,0,side,axis,grid)

#define WPML	(		n,
			side,
			axis,
			grid
	)		PML(n,1,side,axis,grid)

#define WPML	(		n,
			side,
			axis,
			grid
	)		PML(n,1,side,axis,grid)

#define X	(		i0,
			i1,
			i2,
			i3
	)		xp[i0+xDim0*(i1+xDim1*(i2+xDim2*(i3)))]

Function Documentation

void abcMaxwell	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	ndf1a,
		const int &	ndf1b,
		const int &	ndf2a,
		const int &	ndf2b,
		const int &	ndf3a,
		const int &	ndf3b,
		const int &	gid,
		const real &	u,
		const real &	un,
		const real &	f,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const int &	bc,
		const int &	boundaryCondition,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void adjustForIncident	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gid,
		const real &	um,
		const real &	u,
		const real &	un,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const real &	initialConditionBoundingBox,
		const int &	bc,
		const int &	boundaryCondition,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void bcOptMaxwell	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	ndf1a,
		const int &	ndf1b,
		const int &	ndf2a,
		const int &	ndf2b,
		const int &	ndf3a,
		const int &	ndf3b,
		const int &	gid,
		const int &	dimension,
		const real &	u,
		const real &	f,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const int &	bc,
		const int &	boundaryCondition,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void bcSymmetry	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gid,
		const real &	u,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const int &	bc,
		const int &	boundaryCondition,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void exmax	(	double &	Ez,
		double &	Bx,
		double &	By,
		const int &	nsources,
		const double &	xs,
		const double &	ys,
		const double &	tau,
		const double &	var,
		const double &	amp,
		const double &	a,
		const double &	x,
		const double &	y,
		const double &	time
	)

void pmlMaxwell	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	ndf1a,
		const int &	ndf1b,
		const int &	ndf2a,
		const int &	ndf2b,
		const int &	ndf3a,
		const int &	ndf3b,
		const int &	gid,
		const real &	um,
		const real &	u,
		const real &	un,
		const int &	ndra1a,
		const int &	ndra1b,
		const int &	ndra2a,
		const int &	ndra2b,
		const int &	ndra3a,
		const int &	ndra3b,
		const real &	vram,
		const real &	vra,
		const real &	vran,
		const real &	wram,
		const real &	wra,
		const real &	wran,
		const int &	ndrb1a,
		const int &	ndrb1b,
		const int &	ndrb2a,
		const int &	ndrb2b,
		const int &	ndrb3a,
		const int &	ndrb3b,
		const real &	vrbm,
		const real &	vrb,
		const real &	vrbn,
		const real &	wrbm,
		const real &	wrb,
		const real &	wrbn,
		const int &	ndsa1a,
		const int &	ndsa1b,
		const int &	ndsa2a,
		const int &	ndsa2b,
		const int &	ndsa3a,
		const int &	ndsa3b,
		const real &	vsam,
		const real &	vsa,
		const real &	vsan,
		const real &	wsam,
		const real &	wsa,
		const real &	wsan,
		const int &	ndsb1a,
		const int &	ndsb1b,
		const int &	ndsb2a,
		const int &	ndsb2b,
		const int &	ndsb3a,
		const int &	ndsb3b,
		const real &	vsbm,
		const real &	vsb,
		const real &	vsbn,
		const real &	wsbm,
		const real &	wsb,
		const real &	wsbn,
		const int &	ndta1a,
		const int &	ndta1b,
		const int &	ndta2a,
		const int &	ndta2b,
		const int &	ndta3a,
		const int &	ndta3b,
		const real &	vtam,
		const real &	vta,
		const real &	vtan,
		const real &	wtam,
		const real &	wta,
		const real &	wtan,
		const int &	ndtb1a,
		const int &	ndtb1b,
		const int &	ndtb2a,
		const int &	ndtb2b,
		const int &	ndtb3a,
		const int &	ndtb3b,
		const real &	vtbm,
		const real &	vtb,
		const real &	vtbn,
		const real &	wtbm,
		const real &	wtb,
		const real &	wtbn,
		const real &	f,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const int &	bc,
		const int &	boundaryCondition,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void updateghostandperiodic_

(

realMappedGridFunction *&

pu

)

void updateghostboundaries_

(

realArray *&

pu

)