Commit 559637f4 authored by Tim Mattox's avatar Tim Mattox
Browse files

USER-DPD Bugfix: reset_dt() is not called when I thought it should be called.

Move the calculation of dtsqrt inside FixShardlow::initial_integrate()
parent e9fed809
Loading
Loading
Loading
Loading
+2 −7
Original line number Diff line number Diff line
@@ -169,13 +169,6 @@ void FixShardlow::min_setup_pre_exchange()

/* ---------------------------------------------------------------------- */

void FixShardlow::reset_dt()
{
  dtsqrt = sqrt(update->dt);
}

/* ---------------------------------------------------------------------- */

void FixShardlow::setup(int vflag)
{
  bool fixShardlow = false;
@@ -452,6 +445,8 @@ void FixShardlow::initial_integrate(int vflag)
  inum = list->inum;
  ilist = list->ilist;

  dtsqrt = sqrt(update->dt);

  //Loop over all 14 directions (8 stages)
  for (airnum = 1; airnum <=8; airnum++){

+0 −2
Original line number Diff line number Diff line
@@ -40,8 +40,6 @@ class FixShardlow : public Fix {
  void copy_arrays(int, int, int);
  void set_arrays(int);

  void reset_dt();

  int pack_border(int, int *, double *);
  int unpack_border(int, int, double *);
  int unpack_exchange(int, double *);