update2.c

国外一大学教授讲章

#include "fdtd2.h"

/* characteristic impedance of free space */
static double imp0=377.0;  

/* update magnetic field */
void updateH2(Grid *g) {   
  int mm;

  for (mm=0; mm<SizeX-1; mm++)
    Hy(mm) = Hy(mm) + (Ez(mm+1) - Ez(mm))/imp0;

  return;
}                          

/* update electric field */
void updateE2(Grid *g) {   
  int mm;

  for (mm=1; mm<SizeX-1; mm++)
    Ez(mm) = Ez(mm) + (Hy(mm) - Hy(mm-1))*imp0;

  return;
}