git-svn-id: svn://svn.icms.temple.edu/lammps-ro/trunk@6281 f3b2605a-c512-4ea7-a41b-209d697bcdaa

This commit is contained in:
sjplimp
2011-06-01 15:03:14 +00:00
parent 2e7edca71d
commit 75a875e0fc

View File

@ -169,7 +169,7 @@ void PPPM::init()
if (typeA < 1 || typeA > atom->nangletypes ||
force->angle->setflag[typeA] == 0)
error->all("Bad TIP4P angle type for PPPM/TIP4P");
if (typeB < 1 || typeA > atom->nbondtypes ||
if (typeB < 1 || typeB > atom->nbondtypes ||
force->bond->setflag[typeB] == 0)
error->all("Bad TIP4P bond type for PPPM/TIP4P");
double theta = force->angle->equilibrium_angle(typeA);