use correct type
This commit is contained in:
@ -254,7 +254,7 @@ void FixWallFlow::generate_velocity(int atom_i)
|
||||
double delta = gamma * flowvel;
|
||||
|
||||
const double edd = std::exp(-delta * delta) / MathConst::MY_PIS + delta * std::erf(delta);
|
||||
const double probability_threshold = 0.5f * (1.f + delta / edd);
|
||||
const double probability_threshold = 0.5 * (1.0 + delta / edd);
|
||||
|
||||
double direction = 1.0;
|
||||
|
||||
|
||||
@ -241,7 +241,7 @@ FixPressLangevin::FixPressLangevin(LAMMPS *lmp, int narg, char **arg) :
|
||||
if (iarg + 4 > narg) utils::missing_cmd_args(FLERR, "fix press/langevin temp", error);
|
||||
t_start = utils::numeric(FLERR, arg[iarg + 1], false, lmp);
|
||||
t_stop = utils::numeric(FLERR, arg[iarg + 2], false, lmp);
|
||||
seed = utils::numeric(FLERR, arg[iarg + 3], false, lmp);
|
||||
seed = utils::inumeric(FLERR, arg[iarg + 3], false, lmp);
|
||||
if (seed <= 0) error->all(FLERR, "Fix press/langevin temp seed must be > 0");
|
||||
iarg += 4;
|
||||
}
|
||||
|
||||
Reference in New Issue
Block a user