Function test_pml

tests/three_d.cpp:163–194

}

int test_pml(double eps(const vec &)) {
  double a = 10.0;

  grid_volume gv = vol3d(1.5, 1.0, 1.2, a);
  structure s(gv, eps, pml(0.401));

  master_printf("Testing pml quality...\n");
  fields f(&s);
  f.add_point_source(Ez, 0.8, 0.6, 0.0, 4.0, vec(0.751, 0.5, 0.601), 1.0);
  const double deltaT = 10.0;
  const double ttot = 3.1 * deltaT;
  double field_energy_check_time = deltaT;

  while (f.time() < f.last_source_time())
    f.step();

  double last_energy = f.field_energy();
  while (f.time() < ttot) {
    f.step();
    if (f.time() >= field_energy_check_time) {
      const double new_energy = f.field_energy();
      master_printf("Got newE/oldE of %g\n", new_energy / last_energy);
      if (new_energy > last_energy * 4e-3) {
        master_printf("Energy decaying too slowly: from %g to %g (%g)\n", last_energy, new_energy,
                      new_energy / last_energy);
        return 0;
      }
      field_energy_check_time += deltaT;
    }
  }
  return 1;
}

int test_pml_splitting(double eps(const vec &), int splitting) {
  double a = 10.0;