message: Adjusting source functions
diff: @@ -104,7 +104,7 @@ void dEdtBCFunc(const Vector &x, double t, Vector &E);
  static double tScale_ = 1e-9;
  static double freq_ = 750.0e6;
  -static int prob_ = 0;
  +static int prob_ = -1;
  /*
  double permittivity_func(const Vector &);
  double permeability_func(const Vector &);
  @@ -440,7 +440,7 @@ int main(int argc, char *argv[])
  }
  int nsteps = SnapTimeStep(tf-t0, dtmax, dt);
  -
  + /*
  if ( nsteps > max_its )
  {
  if ( mpi.Root() )
  @@ -449,7 +449,7 @@ int main(int argc, char *argv[])
  }
  nsteps = max_its;
  }
  -
  + */
  if ( mpi.Root() )
  {
  cout << "Number of Time Steps: " << nsteps << endl;
  @@ -659,24 +659,23 @@ void dipole_pulse(const Vector &x, double t, Vector &j)
  {
  return;
  }
  + v /= h;
  double r = dp_params_[2*x.Size()+0];
  - double a = dp_params_[2*x.Size()+1];
  + double a = dp_params_[2*x.Size()+1] * tScale_;
  double b = dp_params_[2*x.Size()+2] * tScale_;
  double c = dp_params_[2*x.Size()+3] * tScale_;
  double xv = xu * v;
  - if ( h > 0.0 )
  - {
  - xu.Add(-xv/(h*h), v);
  - }
  + // Compute perpendicular vector from axis to x
  + xu.Add(-xv, v);
  double xp = xu.Norml2();
  - if ( xv >= 0.0 && xv <= h*h && xp <= r )
  + if ( xv >= 0.0 && xv <= h && xp <= r )
  {
  - j.Add(1.0/(h*M_PI*r*r), v);
  + j = v;
  }
  j *= a * (t - b) * exp(-0.5 * pow((t-b)/c, 2)) / (c * c);
  