diff --git a/src/fix_wall_lj93.cpp b/src/fix_wall_lj93.cpp index 3e2c494c0a..83c77a6caa 100644 --- a/src/fix_wall_lj93.cpp +++ b/src/fix_wall_lj93.cpp @@ -34,7 +34,8 @@ void FixWallLJ93::precompute(int m) double rinv = 1.0/cutoff[m]; double r2inv = rinv*rinv; double r4inv = r2inv*r2inv; - offset[m] = coeff3[6]*r4inv*r4inv*rinv - coeff4[6]*r2inv*rinv; + printf("AAA %d %g %g %g\n",m,coeff3[m],coeff4[m],rinv); + offset[m] = coeff3[m]*r4inv*r4inv*rinv - coeff4[m]*r2inv*rinv; } /* ---------------------------------------------------------------------- */