fdtd3.h

国外一大学教授讲章

#define _FDTD3_H 

#ifndef _STDIO_H
#include <stdio.h>
#endif

#ifndef _STDLIB_H
#include <stdlib.h>
#endif

struct Grid {
  double *ez, *ceze, *cezh; 
  double *hy, *chyh, *chye; 
  int sizeX;
  int time, maxTime;
  double cdtds;
};

typedef struct Grid Grid;

/* macros for accessing arrays and such */

#define Hy(MM)    g->hy[MM]
#define Chyh(MM)  g->chyh[MM]
#define Chye(MM)  g->chye[MM]  

#define Ez(MM)    g->ez[MM]
#define Ceze(MM)  g->ceze[MM]
#define Cezh(MM)  g->cezh[MM]  

#define SizeX     g->sizeX
#define Time      g->time
#define MaxTime   g->maxTime
#define Cdtds     g->cdtds

/* memory allocation macro */
#define ALLOC_1D(name,nx,type)                                      \
    name = (type *)calloc(nx,sizeof(type));                         \
    if (!name) {                                                    \
      perror("ALLOC_1D");                                           \
      fprintf(stderr,                                               \
              "Allocation failed for " #name ". Terminating...\n"); \
      exit(-1);                                                     \
    }

/* Function prototypes */ 
void abc(Grid *g);

void updateE3(Grid *g);
void updateH3(Grid *g);

void gridInit3(Grid *g);

void snapshotInit(Grid *g);
void snapshot(Grid *g);

void tfsfInit(Grid *g);
void tfsfUpdate(Grid *g);