This commit is contained in:
jtclemm
2024-11-04 09:41:02 -07:00
parent 3b3d39a8cd
commit 645d02d13b

View File

@ -231,13 +231,13 @@ void PairSPHTaitwater::coeff(int narg, char **arg)
allocate();
int ilo, ihi, jlo, jhi;
utils::bounds(FLERR,arg[0], 1, atom->ntypes, ilo, ihi, error);
utils::bounds(FLERR,arg[1], 1, atom->ntypes, jlo, jhi, error);
utils::bounds(FLERR, arg[0], 1, atom->ntypes, ilo, ihi, error);
utils::bounds(FLERR, arg[1], 1, atom->ntypes, jlo, jhi, error);
double rho0_one = utils::numeric(FLERR,arg[2],false,lmp);
double soundspeed_one = utils::numeric(FLERR,arg[3],false,lmp);
double viscosity_one = utils::numeric(FLERR,arg[4],false,lmp);
double cut_one = utils::numeric(FLERR,arg[5],false,lmp);
double rho0_one = utils::numeric(FLERR, arg[2], false, lmp);
double soundspeed_one = utils::numeric(FLERR, arg[3], false, lmp);
double viscosity_one = utils::numeric(FLERR, arg[4], false, lmp);
double cut_one = utils::numeric(FLERR, arg[5], false, lmp);
double B_one = soundspeed_one * soundspeed_one * rho0_one / 7.0;
int count = 0;
@ -254,7 +254,7 @@ void PairSPHTaitwater::coeff(int narg, char **arg)
}
if (count == 0)
error->all(FLERR,"Incorrect args for pair coefficients");
error->all(FLERR, "Incorrect args for pair coefficients");
}
/* ----------------------------------------------------------------------
@ -264,7 +264,7 @@ void PairSPHTaitwater::coeff(int narg, char **arg)
double PairSPHTaitwater::init_one(int i, int j)
{
if (setflag[i][j] == 0) {
error->all(FLERR,"All pair sph/taitwater coeffs are set");
error->all(FLERR, "All pair sph/taitwater coeffs are notset");
}
cut[j][i] = cut[i][j];