assignBoundaryConditionsSOS.C File Reference

#include "Cgsm.h"
#include "SmParameters.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 "GridMaterialProperties.h"

Include dependency graph for assignBoundaryConditionsSOS.C:

Macros
#define	bcOptSM   EXTERN_C_NAME(bcoptsm)
#define	bcOptSmCons   EXTERN_C_NAME(bcoptsmcons)
#define	abcSolidMechanics   EXTERN_C_NAME(abcsolidmechanics)
#define	exmax   EXTERN_C_NAME(exmax)
#define	exTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-ky/(eps*cc))
#define	eyTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*( kx/(eps*cc))
#define	hzTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))
#define	exLaplacianTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(+ky*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))
#define	eyLaplacianTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-kx*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))
#define	hzLaplacianTrue(x, y, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*( -(twoPi*twoPi*(kx*kx+ky*ky) ) )
#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)))*(-ky/(eps*cc))
#define	eyTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*( kx/(eps*cc))
#define	ezTrue3d(x, y, z, t)   0
#define	hxTrue3d(x, y, z, t)   0
#define	hyTrue3d(x, y, z, t)   0
#define	hzTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))
#define	exLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(+ky*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))
#define	eyLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-kx*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))
#define	ezLaplacianTrue3d(x, y, z, t)   0
#define	hxLaplacianTrue3d(x, y, z, t)   0
#define	hyLaplacianTrue3d(x, y, z, t)   0
#define	hzLaplacianTrue3d(x, y, z, t)   sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*( -(twoPi*twoPi*(kx*kx+ky*ky) ) )
#define	pmlSolidMechanics   EXTERN_C_NAME(pmlsolidmechanics)
#define	interfaceSolidMechanics   EXTERN_C_NAME(interfacemaxwell)
#define	newInterfaceSolidMechanics   EXTERN_C_NAME(newinterfacemaxwell)
#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	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	U0(x, y, z, n, t)   (vcenter[n-uc]*(t) + rx[n-uc]*((x)-xcenter) + ry[n-uc]*((y)-ycenter))
#define	U0T(x, y, z, n, t)   (vcenter[n-uc] + rxt[n-uc]*((x)-xcenter) + ryt[n-uc]*((y)-ycenter))
#define	U0X(x, y, z, n, t)   ( rx[n-uc] )
#define	U0Y(x, y, z, n, t)   ( ry[n-uc] )
#define	dbc(s, a, side, axis)   (pdbc[(s)+2*((a)+3*((side)+2*(axis)))])
#define	addBoundaryForcing(side, axis)   (pAddBoundaryForcing[(side)+2*(axis)])
#define	bcfOffset(side, axis)   pbcfOffset[(side)+2*(axis)]
#define	VPML(m, side, axis, grid)   vpml[(m)+numberOfPMLFunctions*(side+2*(axis+3*(grid)))]
#define	WPML(m, side, axis, grid)   vpml[(m+numberOfTimeLevels)+numberOfPMLFunctions*(side+2*(axis+3*(grid)))]

Functions
void	bcOptSM (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange, real &u, const int &mask, const real &rsxy, const real &xy, const int &ndMatProp, const int &matIndex, const real &matValpc, const real &matVal, const int &boundaryCondition, const int &addBoundaryForcing, const int &interfaceType, const int &dim, const real &bcf00, const real &bcf10, const real &bcf01, const real &bcf11, const real &bcf02, const real &bcf12, const real &bcf0, const int64_t &bcfOffset, const int &ipar, const real &rpar, const int &ierr)
void	bcOptSmCons (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange, real &u, const int &mask, const real &rsxy, const real &xy, const int &ndMatProp, const int &matIndex, const real &matValpc, const real &matVal, const int &boundaryCondition, const int &addBoundaryForcing, const int &interfaceType, const int &dim, const real &bcf00, const real &bcf10, const real &bcf01, const real &bcf11, const real &bcf02, const real &bcf12, const real &bcf0, const int64_t &bcfOffset, const int &ipar, const real &rpar, const int &ierr)
void	abcSolidMechanics (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	pmlSolidMechanics (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	interfaceSolidMechanics (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange1, real &u1, const int &mask1, const real &rsxy1, const real &xy1, const int &boundaryCondition1, const int &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &gridIndexRange2, real &u2, const int &mask2, const real &rsxy2, const real &xy2, const int &boundaryCondition2, const int &ipar, const real &rpar, real &aa2, real &aa4, real &aa8, int &ipvt2, int &ipvt4, int &ipvt8, int &ierr)
void	newInterfaceSolidMechanics (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange1, real &u1, const int &mask1, const real &rsxy1, const real &xy1, const int &boundaryCondition1, const int &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &gridIndexRange2, real &u2, const int &mask2, const real &rsxy2, const real &xy2, const int &boundaryCondition2, const int &ipar, const real &rpar, int &ierr)
void	updateghostboundaries_ (realArray *&pu)
void	updateghostandperiodic_ (realMappedGridFunction *&pu)

Macro Definition Documentation

#define abcSolidMechanics   EXTERN_C_NAME(abcsolidmechanics)

#define addBoundaryForcing	(		side,
			axis
	)		(pAddBoundaryForcing[(side)+2*(axis)])

#define bcfOffset	(		side,
			axis
	)		pbcfOffset[(side)+2*(axis)]

#define bcOptSM   EXTERN_C_NAME(bcoptsm)

Referenced by Cgsm::assignBoundaryConditionsSOS().

#define bcOptSmCons   EXTERN_C_NAME(bcoptsmcons)

Referenced by Cgsm::assignBoundaryConditionsSOS().

#define dbc	(		s,
			a,
			side,
			axis
	)		(pdbc[(s)+2*((a)+3*((side)+2*(axis)))])

#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)))*(+ky*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))

#define exLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(+ky*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))

#define exmax   EXTERN_C_NAME(exmax)

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

#define exTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-ky/(eps*cc))

#define exTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*(-ky/(eps*cc))

#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)))*(-kx*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))

#define eyLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*(-kx*(twoPi*twoPi*(kx*kx+ky*ky))/(eps*cc))

#define eyTrue	(		x,
			y,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*( kx/(eps*cc))

#define eyTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)+kz*(z)-cc*(t)))*( kx/(eps*cc))

#define ezLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		0

#define ezTrue3d	(		x,
			y,
			z,
			t
	)		0

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

#define hxTrue3d	(		x,
			y,
			z,
			t
	)		0

#define hyLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		0

#define hyTrue3d	(		x,
			y,
			z,
			t
	)		0

#define hzGaussianPulse

(

xi

)

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

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

#define hzLaplacianTrue3d	(		x,
			y,
			z,
			t
	)		sin(twoPi*(kx*(x)+ky*(y)-cc*(t)))*( -(twoPi*twoPi*(kx*kx+ky*ky) ) )

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

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

#define interfaceSolidMechanics   EXTERN_C_NAME(interfacemaxwell)

#define newInterfaceSolidMechanics   EXTERN_C_NAME(newinterfacemaxwell)

#define pmlSolidMechanics   EXTERN_C_NAME(pmlsolidmechanics)

Referenced by Cgsm::assignBoundaryConditionsSOS().

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

#define U0	(		x,
			y,
			z,
			n,
			t
	)		(vcenter[n-uc]*(t) + rx[n-uc]*((x)-xcenter) + ry[n-uc]*((y)-ycenter))

#define U0T	(		x,
			y,
			z,
			n,
			t
	)		(vcenter[n-uc] + rxt[n-uc]*((x)-xcenter) + ryt[n-uc]*((y)-ycenter))

#define U0X	(		x,
			y,
			z,
			n,
			t
	)		( rx[n-uc] )

#define U0Y	(		x,
			y,
			z,
			n,
			t
	)		( ry[n-uc] )

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

Referenced by Maxwell::assignBoundaryConditions(), and Cgsm::assignBoundaryConditionsSOS().

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

Referenced by Maxwell::assignBoundaryConditions(), and Cgsm::assignBoundaryConditionsSOS().

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

Function Documentation

void abcSolidMechanics	(	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 bcOptSM	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gridIndexRange,
		real &	u,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const int &	ndMatProp,
		const int &	matIndex,
		const real &	matValpc,
		const real &	matVal,
		const int &	boundaryCondition,
		const int &	addBoundaryForcing,
		const int &	interfaceType,
		const int &	dim,
		const real &	bcf00,
		const real &	bcf10,
		const real &	bcf01,
		const real &	bcf11,
		const real &	bcf02,
		const real &	bcf12,
		const real &	bcf0,
		const int64_t &	bcfOffset,
		const int &	ipar,
		const real &	rpar,
		const int &	ierr
	)

void bcOptSmCons	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gridIndexRange,
		real &	u,
		const int &	mask,
		const real &	rsxy,
		const real &	xy,
		const int &	ndMatProp,
		const int &	matIndex,
		const real &	matValpc,
		const real &	matVal,
		const int &	boundaryCondition,
		const int &	addBoundaryForcing,
		const int &	interfaceType,
		const int &	dim,
		const real &	bcf00,
		const real &	bcf10,
		const real &	bcf01,
		const real &	bcf11,
		const real &	bcf02,
		const real &	bcf12,
		const real &	bcf0,
		const int64_t &	bcfOffset,
		const int &	ipar,
		const real &	rpar,
		const int &	ierr
	)

void interfaceSolidMechanics	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gridIndexRange1,
		real &	u1,
		const int &	mask1,
		const real &	rsxy1,
		const real &	xy1,
		const int &	boundaryCondition1,
		const int &	md1a,
		const int &	md1b,
		const int &	md2a,
		const int &	md2b,
		const int &	md3a,
		const int &	md3b,
		const int &	gridIndexRange2,
		real &	u2,
		const int &	mask2,
		const real &	rsxy2,
		const real &	xy2,
		const int &	boundaryCondition2,
		const int &	ipar,
		const real &	rpar,
		real &	aa2,
		real &	aa4,
		real &	aa8,
		int &	ipvt2,
		int &	ipvt4,
		int &	ipvt8,
		int &	ierr
	)

void newInterfaceSolidMechanics	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gridIndexRange1,
		real &	u1,
		const int &	mask1,
		const real &	rsxy1,
		const real &	xy1,
		const int &	boundaryCondition1,
		const int &	md1a,
		const int &	md1b,
		const int &	md2a,
		const int &	md2b,
		const int &	md3a,
		const int &	md3b,
		const int &	gridIndexRange2,
		real &	u2,
		const int &	mask2,
		const real &	rsxy2,
		const real &	xy2,
		const int &	boundaryCondition2,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void pmlSolidMechanics	(	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

)