advanceStructured.C File Reference

#include "Maxwell.h"
#include "display.h"
#include "CompositeGridOperators.h"
#include "ParallelUtility.h"
#include "ParallelGridUtility.h"

Include dependency graph for advanceStructured.C:

Macros
#define	advMaxwell   EXTERN_C_NAME(advmaxwell)
#define	mxFilter   EXTERN_C_NAME(mxfilter)
#define	FD4_2D(u, i1, i2, i3, c)
#define	FD4_3D(u, i1, i2, i3, c)
#define	FOR_3(i1, i2, i3, I1, I2, I3)   for( i3=I3.getBase(); i3<=I3.getBound(); i3++ ) for( i2=I2.getBase(); i2<=I2.getBound(); i2++ ) for( i1=I1.getBase(); i1<=I1.getBound(); i1++ )
#define	FN(m)   fn[m+numberOfFunctions*(grid)]
#define	FD4A_2D(u, i1, i2, i3, c)
#define	FD4A_3D(u, i1, i2, i3, c)
#define	FD4V_2D(u, i1, i2, i3, c)
#define	FD4V_3D(u, i1, i2, i3, c)

Functions
void	advMaxwell (const int &nd, const int &n1a, const int &n1b, const int &n2a, const int &n2b, const int &n3a, const int &n3b, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &nd4a, const int &nd4b, const int &mask, const real &rx, const real &um, const real &u, real &un, const real &f, const real &ut1, const real &ut2, const real &ut3, const real &ut4, const real &ut5, const real &ut6, const real &ut7, const int &bc, const real &dis, const real &varDis, const int &ipar, const real &rpar, int &ierr)
void	mxFilter (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange, const real &u, const real &d, const int &mask, const int &boundaryCondition, const int &ipar, const real &rpar, int &ierr)

Macro Definition Documentation

#define advMaxwell   EXTERN_C_NAME(advmaxwell)

Referenced by Maxwell::advanceNFDTD().

#define FD4_2D	(		u,
			i1,
			i2,
			i3,
			c
	)

Value:

(    -( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c) )   \
        +4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c) ) \
       -12.*u(i1,i2,i3,c) )

Referenced by Maxwell::computeDissipation().

#define FD4_3D	(		u,
			i1,
			i2,
			i3,
			c
	)

Value:

(    -( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c)+u(i1,i2,i3-2,c)+u(i1,i2,i3+2,c) )   \
        +4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c)+u(i1,i2,i3-1,c)+u(i1,i2,i3+1,c) ) \
       -18.*u(i1,i2,i3,c) )

Referenced by Maxwell::computeDissipation().

#define FD4A_2D	(		u,
			i1,
			i2,
			i3,
			c
	)

Value:

(     ( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c) )   \
        -4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c) ) \
       +12.*u(i1,i2,i3,c) )

Referenced by Maxwell::addFilter().

#define FD4A_3D	(		u,
			i1,
			i2,
			i3,
			c
	)

Value:

(     ( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c)+u(i1,i2,i3-2,c)+u(i1,i2,i3+2,c) )   \
        -4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c)+u(i1,i2,i3-1,c)+u(i1,i2,i3+1,c) ) \
       +18.*u(i1,i2,i3,c) )

Referenced by Maxwell::addFilter().

#define FD4V_2D	(		u,
			i1,
			i2,
			i3,
			c
	)

Value:

(     ax*( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c) )  \
           +ay*( u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c) )   \
          -ax4*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c) ) \
          -ay4*( u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c) ) \
       + axy12*u(i1,i2,i3,c) )

#define FD4V_3D	(		u,
			i1,
			i2,
			i3,
			c
	)

Value:

(     ax*( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c) )\
           +ay*( u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c) )\
           +az*( u(i1,i2,i3-2,c)+u(i1,i2,i3+2,c) ) \
          -ax4*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c) ) \
          -ay4*( u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c) ) \
          -az4*( u(i1,i2,i3-1,c)+u(i1,i2,i3+1,c) ) \
          +axyz18*u(i1,i2,i3,c) )

#define FN

(

m

)

fn[m+numberOfFunctions*(grid)]

Referenced by Maxwell::advanceNFDTD().

#define FOR_3	(		i1,
			i2,
			i3,
			I1,
			I2,
			I3
	)		for( i3=I3.getBase(); i3<=I3.getBound(); i3++ ) for( i2=I2.getBase(); i2<=I2.getBound(); i2++ ) for( i1=I1.getBase(); i1<=I1.getBound(); i1++ )

#define mxFilter   EXTERN_C_NAME(mxfilter)

Referenced by Maxwell::addFilter().

Function Documentation

void advMaxwell	(	const int &	nd,
		const int &	n1a,
		const int &	n1b,
		const int &	n2a,
		const int &	n2b,
		const int &	n3a,
		const int &	n3b,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	nd4a,
		const int &	nd4b,
		const int &	mask,
		const real &	rx,
		const real &	um,
		const real &	u,
		real &	un,
		const real &	f,
		const real &	ut1,
		const real &	ut2,
		const real &	ut3,
		const real &	ut4,
		const real &	ut5,
		const real &	ut6,
		const real &	ut7,
		const int &	bc,
		const real &	dis,
		const real &	varDis,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)

void mxFilter	(	const int &	nd,
		const int &	nd1a,
		const int &	nd1b,
		const int &	nd2a,
		const int &	nd2b,
		const int &	nd3a,
		const int &	nd3b,
		const int &	gridIndexRange,
		const real &	u,
		const real &	d,
		const int &	mask,
		const int &	boundaryCondition,
		const int &	ipar,
		const real &	rpar,
		int &	ierr
	)