fdtd-macro.h

利用VC++进行FDTD三维远近场变换

#ifndef _FDTD_MACRO_H
#define _FDTD_MACRO_H

#include "fdtd-grid.h"

/* macros that permit the "Grid" to be specified */
/* one-dimensional grid */
#define Hy1G(G,MM)     G->hy[MM]
#define Chyh1G(G,MM)   G->chyh[MM]
#define Chye1G(G,MM)   G->chye[MM]

#define Ez1G(G,MM)     G->ez[MM]
#define Ceze1G(G,MM)   G->ceze[MM]
#define Cezh1G(G,MM)   G->cezh[MM]

/* TMz grid */
#define Hx2G(G,MM,NN)   G->hx[(MM)*(SizeYG(G)-1) + NN]
#define Chxh2G(G,MM,NN) G->chxh[(MM)*(SizeYG(G)-1) + NN]
#define Chxe2G(G,MM,NN) G->chxe[(MM)*(SizeYG(G)-1) + NN]

#define Hy2G(G,MM,NN)   G->hy[(MM)*SizeYG(G) + NN]
#define Chyh2G(G,MM,NN) G->chyh[(MM)*SizeYG(G) + NN]
#define Chye2G(G,MM,NN) G->chye[(MM)*SizeYG(G) + NN]

#define Ez2G(G,MM,NN)   G->ez[(MM)*SizeYG(G) + NN]
#define Ceze2G(G,MM,NN) G->ceze[(MM)*SizeYG(G) + NN]
#define Cezh2G(G,MM,NN) G->cezh[(MM)*SizeYG(G) + NN]

/* TEz grid */
#define Ex2G(G,MM,NN)   G->ex[(MM)*SizeYG(G) + NN]
#define Cexe2G(G,MM,NN) G->cexe[(MM)*SizeYG(G) + NN]
#define Cexh2G(G,MM,NN) G->cexh[(MM)*SizeYG(G) + NN]

#define Ey2G(G,MM,NN)   G->ey[(MM)*(SizeYG(G)-1) + NN]
#define Ceye2G(G,MM,NN) G->ceye[(MM)*(SizeYG(G)-1) + NN]
#define Ceyh2G(G,MM,NN) G->ceyh[(MM)*(SizeYG(G)-1) + NN]

#define Hz2G(G,MM,NN)   G->hz[(MM)*(SizeYG(G)-1) + NN]
#define Chzh2G(G,MM,NN) G->chzh[(MM)*(SizeYG(G)-1) + NN]
#define Chze2G(G,MM,NN) G->chze[(MM)*(SizeYG(G)-1) + NN]

/* 3D grid */
#define HxG(G,MM,NN,PP)   G->hx[((MM)*(SizeYG(G)-1) + NN)*(SizeZG(G)-1) + PP]
#define ChxhG(G,MM,NN,PP) G->chxh[((MM)*(SizeYG(G)-1) + NN)*(SizeZG(G)-1) + PP]
#define ChxeG(G,MM,NN,PP) G->chxe[((MM)*(SizeYG(G)-1) + NN)*(SizeZG(G)-1) + PP]

#define HyG(G,MM,NN,PP)   G->hy[((MM)*SizeYG(G) + NN)*(SizeZG(G)-1) + PP]
#define ChyhG(G,MM,NN,PP) G->chyh[((MM)*SizeYG(G) + NN)*(SizeZG(G)-1) + PP]
#define ChyeG(G,MM,NN,PP) G->chye[((MM)*SizeYG(G) + NN)*(SizeZG(G)-1) + PP]

#define HzG(G,MM,NN,PP)   G->hz[((MM)*(SizeYG(G)-1) + NN)*SizeZG(G) + PP]
#define ChzhG(G,MM,NN,PP) G->chzh[((MM)*(SizeYG(G)-1) + NN)*SizeZG(G) + PP]
#define ChzeG(G,MM,NN,PP) G->chze[((MM)*(SizeYG(G)-1) + NN)*SizeZG(G) + PP]

#define ExG(G,MM,NN,PP)   G->ex[((MM)*SizeYG(G) + NN)*SizeZG(G) + PP]
#define CexeG(G,MM,NN,PP) G->cexe[((MM)*SizeYG(G) + NN)*SizeZG(G) + PP]
#define CexhG(G,MM,NN,PP) G->cexh[((MM)*SizeYG(G) + NN)*SizeZG(G) + PP]

#define EyG(G,MM,NN,PP)   G->ey[((MM)*(SizeYG(G)-1) + NN)*SizeZG(G) + PP]
#define CeyeG(G,MM,NN,PP) G->ceye[((MM)*(SizeYG(G)-1) + NN)*SizeZG(G) + PP]
#define CeyhG(G,MM,NN,PP) G->ceyh[((MM)*(SizeYG(G)-1) + NN)*SizeZG(G) + PP]

#define EzG(G,MM,NN,PP)   G->ez[((MM)*SizeYG(G) + NN)*(SizeZG(G)-1) + PP]
#define CezeG(G,MM,NN,PP) G->ceze[((MM)*SizeYG(G) + NN)*(SizeZG(G)-1) + PP]
#define CezhG(G,MM,NN,PP) G->cezh[((MM)*SizeYG(G) + NN)*(SizeZG(G)-1) + PP]

#define SizeXG(G)      G->sizeX
#define SizeYG(G)      G->sizeY
#define SizeZG(G)      G->sizeZ
#define TimeG(G)       G->time
#define MaxTimeG(G)    G->maxTime
#define CdtdsG(G)      G->cdtds
#define TypeG(G)       G->type

/* macros that assume the "Grid" is "g" */
/* one-dimensional grid */
#define Hy1(MM)     Hy1G(g,MM)
#define Chyh1(MM)   Chyh1G(g,MM)   
#define Chye1(MM)   Chye1G(g,MM)   

#define Ez1(MM)     Ez1G(g,MM)     
#define Ceze1(MM)   Ceze1G(g,MM)   
#define Cezh1(MM)   Cezh1G(g,MM)   

/* TMz grid */
#define Hx2(MM,NN)   Hx2G(g,MM,NN)   
#define Chxh2(MM,NN) Chxh2G(g,MM,NN) 
#define Chxe2(MM,NN) Chxe2G(g,MM,NN) 

#define Hy2(MM,NN)   Hy2G(g,MM,NN)   
#define Chyh2(MM,NN) Chyh2G(g,MM,NN) 
#define Chye2(MM,NN) Chye2G(g,MM,NN) 

#define Ez2(MM,NN)   Ez2G(g,MM,NN)   
#define Ceze2(MM,NN) Ceze2G(g,MM,NN) 
#define Cezh2(MM,NN) Cezh2G(g,MM,NN) 

/* TEz grid */
#define Hz2(MM,NN)   Hz2G(g,MM,NN)   
#define Chzh2(MM,NN) Chzh2G(g,MM,NN) 
#define Chze2(MM,NN) Chze2G(g,MM,NN) 

#define Ex2(MM,NN)   Ex2G(g,MM,NN)   
#define Cexe2(MM,NN) Cexe2G(g,MM,NN) 
#define Cexh2(MM,NN) Cexh2G(g,MM,NN) 

#define Ey2(MM,NN)   Ey2G(g,MM,NN)   
#define Ceye2(MM,NN) Ceye2G(g,MM,NN) 
#define Ceyh2(MM,NN) Ceyh2G(g,MM,NN) 

/* 3D grid */
#define Hx(MM,NN,PP)   HxG(g,MM,NN,PP)   
#define Chxh(MM,NN,PP) ChxhG(g,MM,NN,PP) 
#define Chxe(MM,NN,PP) ChxeG(g,MM,NN,PP) 

#define Hy(MM,NN,PP)   HyG(g,MM,NN,PP)   
#define Chyh(MM,NN,PP) ChyhG(g,MM,NN,PP) 
#define Chye(MM,NN,PP) ChyeG(g,MM,NN,PP) 

#define Hz(MM,NN,PP)   HzG(g,MM,NN,PP)   
#define Chzh(MM,NN,PP) ChzhG(g,MM,NN,PP) 
#define Chze(MM,NN,PP) ChzeG(g,MM,NN,PP) 

#define Ex(MM,NN,PP)   ExG(g,MM,NN,PP)   
#define Cexe(MM,NN,PP) CexeG(g,MM,NN,PP) 
#define Cexh(MM,NN,PP) CexhG(g,MM,NN,PP) 

#define Ey(MM,NN,PP)   EyG(g,MM,NN,PP)   
#define Ceye(MM,NN,PP) CeyeG(g,MM,NN,PP) 
#define Ceyh(MM,NN,PP) CeyhG(g,MM,NN,PP) 

#define Ez(MM,NN,PP)   EzG(g,MM,NN,PP)
#define Ceze(MM,NN,PP) CezeG(g,MM,NN,PP)
#define Cezh(MM,NN,PP) CezhG(g,MM,NN,PP) 

#define SizeX       SizeXG(g)
#define SizeY       SizeYG(g)
#define SizeZ       SizeZG(g)
#define Time        TimeG(g)
#define MaxTime     MaxTimeG(g)   
#define Cdtds       CdtdsG(g)
#define Type        TypeG(g)

#endif