fdtd2.h

国外一大学教授讲章

#define _FDTD2_H    

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

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

struct Grid {
  double *ez;
  double *hy;
  int sizeX;
  int time, maxTime;
  double cdtds;
};

typedef struct Grid Grid;

/* 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);                                                     \
    }

/* macros for accessing arrays and such */
#define Hy(MM)    g->hy[MM]   
#define Ez(MM)    g->ez[MM]   
#define SizeX     g->sizeX    
#define Time      g->time     
#define MaxTime   g->maxTime  
#define Cdtds     g->cdtds    

/* function prototypes */
void gridInit2(Grid *g);      
void updateH2(Grid *g);       
void updateE2(Grid *g);