diff --git a/src/ASPHERE/pair_gayberne.cpp b/src/ASPHERE/pair_gayberne.cpp
index 4aae7177dae1c01ccd883413b39f0d8122a0c339..bdff7a5cd684e7481832aaae8ed7cdea3db4408e 100644
--- a/src/ASPHERE/pair_gayberne.cpp
+++ b/src/ASPHERE/pair_gayberne.cpp
@@ -297,8 +297,8 @@ void PairGayBerne::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/ASPHERE/pair_line_lj.cpp b/src/ASPHERE/pair_line_lj.cpp
index aa3493ef47c1511424d284c236d628c654631943..4e3df473a385bc92d635ae66406613bd061382dc 100644
--- a/src/ASPHERE/pair_line_lj.cpp
+++ b/src/ASPHERE/pair_line_lj.cpp
@@ -371,8 +371,8 @@ void PairLineLJ::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double size_itype = force->numeric(FLERR,arg[2]);
double size_jtype = force->numeric(FLERR,arg[3]);
diff --git a/src/ASPHERE/pair_resquared.cpp b/src/ASPHERE/pair_resquared.cpp
index 870ec95740eb0d0f314653953ec6cff584b0eebb..172516aa49021691de908ca9b24cfcff3db156eb 100644
--- a/src/ASPHERE/pair_resquared.cpp
+++ b/src/ASPHERE/pair_resquared.cpp
@@ -269,8 +269,8 @@ void PairRESquared::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/ASPHERE/pair_tri_lj.cpp b/src/ASPHERE/pair_tri_lj.cpp
index ef5123dc11b3f3da24bb0ab0f58908eda9a3cbc2..773ad2d6a3737707d53112462818d4ebec33aea9 100644
--- a/src/ASPHERE/pair_tri_lj.cpp
+++ b/src/ASPHERE/pair_tri_lj.cpp
@@ -442,8 +442,8 @@ void PairTriLJ::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/BODY/pair_body.cpp b/src/BODY/pair_body.cpp
index 86fea7e859c43bdad3f51192db3261398d34e29e..2a9edb37cc5fc365e0d9391440950535286a57ce 100644
--- a/src/BODY/pair_body.cpp
+++ b/src/BODY/pair_body.cpp
@@ -388,8 +388,8 @@ void PairBody::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/CLASS2/angle_class2.cpp b/src/CLASS2/angle_class2.cpp
index ceeba556787ad08471d158c8139410912c919642..0315e6601f65033e9185f19eba39ec2eab9c70f6 100644
--- a/src/CLASS2/angle_class2.cpp
+++ b/src/CLASS2/angle_class2.cpp
@@ -271,7 +271,7 @@ void AngleClass2::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
int count = 0;
diff --git a/src/CLASS2/bond_class2.cpp b/src/CLASS2/bond_class2.cpp
index b87deed5a65b41217c973c8c19a571060fac2e51..f358720e7f77b074820e80c3e4cf4bea72b43427 100644
--- a/src/CLASS2/bond_class2.cpp
+++ b/src/CLASS2/bond_class2.cpp
@@ -132,7 +132,7 @@ void BondClass2::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double r0_one = force->numeric(FLERR,arg[1]);
double k2_one = force->numeric(FLERR,arg[2]);
diff --git a/src/CLASS2/dihedral_class2.cpp b/src/CLASS2/dihedral_class2.cpp
index d18d75b155fdf038aec42508b5fe3dfefa9fa31e..78be10fbcac7554c7686702d1de1ea6a636096f5 100644
--- a/src/CLASS2/dihedral_class2.cpp
+++ b/src/CLASS2/dihedral_class2.cpp
@@ -640,7 +640,7 @@ void DihedralClass2::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
int count = 0;
diff --git a/src/CLASS2/improper_class2.cpp b/src/CLASS2/improper_class2.cpp
index 6c599276ca1b27c1fa2c168c53ed4aef45ae7cc2..618e20f6c38a320103642a6d5e9426eb5c17fc45 100644
--- a/src/CLASS2/improper_class2.cpp
+++ b/src/CLASS2/improper_class2.cpp
@@ -529,7 +529,7 @@ void ImproperClass2::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
int count = 0;
diff --git a/src/CLASS2/pair_lj_class2.cpp b/src/CLASS2/pair_lj_class2.cpp
index 0974769b9c7c3464b3306b488f9e329e8c480a2b..ee61aaae1fdcc905bf4d8ad6ed811547f29afc25 100644
--- a/src/CLASS2/pair_lj_class2.cpp
+++ b/src/CLASS2/pair_lj_class2.cpp
@@ -189,8 +189,8 @@ void PairLJClass2::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/CLASS2/pair_lj_class2_coul_cut.cpp b/src/CLASS2/pair_lj_class2_coul_cut.cpp
index c7c473868947b11ade979bb966625ef773bfd08d..45f0ccfe2788069aa47d5473195545127aef4d63 100644
--- a/src/CLASS2/pair_lj_class2_coul_cut.cpp
+++ b/src/CLASS2/pair_lj_class2_coul_cut.cpp
@@ -220,8 +220,8 @@ void PairLJClass2CoulCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/CLASS2/pair_lj_class2_coul_long.cpp b/src/CLASS2/pair_lj_class2_coul_long.cpp
index 62ffdab6e1226a7edfbdfcabf173c2751045dd59..b58094713f79ff181c72504b5497ce91831d44cd 100644
--- a/src/CLASS2/pair_lj_class2_coul_long.cpp
+++ b/src/CLASS2/pair_lj_class2_coul_long.cpp
@@ -256,8 +256,8 @@ void PairLJClass2CoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/COLLOID/pair_brownian.cpp b/src/COLLOID/pair_brownian.cpp
index c1cc523d3262c1885a1dcb485a907b8f36094f5c..84fda485a036b0f93a1997335be4f73f2c78f5f8 100644
--- a/src/COLLOID/pair_brownian.cpp
+++ b/src/COLLOID/pair_brownian.cpp
@@ -423,8 +423,8 @@ void PairBrownian::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_inner_one = cut_inner_global;
double cut_one = cut_global;
diff --git a/src/COLLOID/pair_colloid.cpp b/src/COLLOID/pair_colloid.cpp
index b5d83233ab9ca1713609d6a2096428fdfae08d53..440d6f9d4f8f86da23ab16e97ffb3f693b20d563 100644
--- a/src/COLLOID/pair_colloid.cpp
+++ b/src/COLLOID/pair_colloid.cpp
@@ -272,8 +272,8 @@ void PairColloid::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a12_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/COLLOID/pair_lubricate.cpp b/src/COLLOID/pair_lubricate.cpp
index ea398c340a9467711ad95fb19cab7573b5eae20e..71e08f3f1993501268402885d61530e422a60c9e 100644
--- a/src/COLLOID/pair_lubricate.cpp
+++ b/src/COLLOID/pair_lubricate.cpp
@@ -509,8 +509,8 @@ void PairLubricate::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_inner_one = cut_inner_global;
double cut_one = cut_global;
diff --git a/src/COLLOID/pair_lubricateU.cpp b/src/COLLOID/pair_lubricateU.cpp
index 214bc3c2f3ed68e71fb87d2d234087582d6165f0..eafa57973ca135404c90fdbe086d28acc93f38f1 100644
--- a/src/COLLOID/pair_lubricateU.cpp
+++ b/src/COLLOID/pair_lubricateU.cpp
@@ -1739,8 +1739,8 @@ void PairLubricateU::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_inner_one = cut_inner_global;
double cut_one = cut_global;
diff --git a/src/DIPOLE/pair_lj_cut_dipole_cut.cpp b/src/DIPOLE/pair_lj_cut_dipole_cut.cpp
index 53c977f876ede3268d7d75946a166042843e0148..c57eb09e52e2f8fca83b4dac4aefe09fb8af7d1a 100644
--- a/src/DIPOLE/pair_lj_cut_dipole_cut.cpp
+++ b/src/DIPOLE/pair_lj_cut_dipole_cut.cpp
@@ -326,8 +326,8 @@ void PairLJCutDipoleCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/DIPOLE/pair_lj_cut_dipole_long.cpp b/src/DIPOLE/pair_lj_cut_dipole_long.cpp
index 3f44579b8f265615ec61e1dd9dcb1dcfcaced593..ae85b55ff7063089c10067e04ce4eea1e9b834c4 100644
--- a/src/DIPOLE/pair_lj_cut_dipole_long.cpp
+++ b/src/DIPOLE/pair_lj_cut_dipole_long.cpp
@@ -376,8 +376,8 @@ void PairLJCutDipoleLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/DIPOLE/pair_lj_long_dipole_long.cpp b/src/DIPOLE/pair_lj_long_dipole_long.cpp
index 248865ef728ca04fb658a0929397d11094d6c31c..ef865b66cd50efd89d75a8c5b5d5e24770d99d07 100644
--- a/src/DIPOLE/pair_lj_long_dipole_long.cpp
+++ b/src/DIPOLE/pair_lj_long_dipole_long.cpp
@@ -194,8 +194,8 @@ void PairLJLongDipoleLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/GRANULAR/pair_gran_hooke_history.cpp b/src/GRANULAR/pair_gran_hooke_history.cpp
index 2652381c966c8c3a7c93ae164ce35681a9b8a114..edd12b4b2720ef4df54c46598b6e6efa2524e226 100644
--- a/src/GRANULAR/pair_gran_hooke_history.cpp
+++ b/src/GRANULAR/pair_gran_hooke_history.cpp
@@ -371,8 +371,8 @@ void PairGranHookeHistory::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/KIM/pair_kim.cpp b/src/KIM/pair_kim.cpp
index 9140b73f3b6ef324e875d55c9c5fcb5c30058043..b35e90511f64bfa954f16ddd1de84a02f6df1c16 100644
--- a/src/KIM/pair_kim.cpp
+++ b/src/KIM/pair_kim.cpp
@@ -343,8 +343,8 @@ void PairKIM::coeff(int narg, char **arg)
error->all(FLERR,"Incorrect args for pair coefficients");
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
// read args that map atom species to KIM elements
// lmps_map_species_to_unique[i] =
diff --git a/src/KOKKOS/pair_table_kokkos.cpp b/src/KOKKOS/pair_table_kokkos.cpp
index fec6512a331d474935529a0542cdaf0ee33a68e2..278c5b0a2fd40d1f6739306738d3c4fccdc3f77d 100644
--- a/src/KOKKOS/pair_table_kokkos.cpp
+++ b/src/KOKKOS/pair_table_kokkos.cpp
@@ -488,8 +488,8 @@ void PairTableKokkos<DeviceType>::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/KSPACE/pair_born_coul_long.cpp b/src/KSPACE/pair_born_coul_long.cpp
index 15bcd1081b892a610383b8d33c5ef9e1953c698c..14d43f4c6383fcb1dcc82e9fbbb89f6823f62ba7 100644
--- a/src/KSPACE/pair_born_coul_long.cpp
+++ b/src/KSPACE/pair_born_coul_long.cpp
@@ -265,8 +265,8 @@ void PairBornCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/KSPACE/pair_buck_coul_long.cpp b/src/KSPACE/pair_buck_coul_long.cpp
index 683e5bf0045e17f979e9d44e489576ec1fe762aa..9cd8485e5ccbcc849e26543f5fbcbea516317ac4 100644
--- a/src/KSPACE/pair_buck_coul_long.cpp
+++ b/src/KSPACE/pair_buck_coul_long.cpp
@@ -256,8 +256,8 @@ void PairBuckCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/KSPACE/pair_buck_long_coul_long.cpp b/src/KSPACE/pair_buck_long_coul_long.cpp
index 52d250fbe180c8d9284d8f8b6eaa804bab777e1f..7b42bb47d8539ee937163ae9301aea62f00a3e25 100644
--- a/src/KSPACE/pair_buck_long_coul_long.cpp
+++ b/src/KSPACE/pair_buck_long_coul_long.cpp
@@ -197,8 +197,8 @@ void PairBuckLongCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(*(arg++),atom->ntypes,ilo,ihi);
- force->bounds(*(arg++),atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,*(arg++),atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,*(arg++),atom->ntypes,jlo,jhi);
double buck_a_one = force->numeric(FLERR,*(arg++));
double buck_rho_one = force->numeric(FLERR,*(arg++));
diff --git a/src/KSPACE/pair_coul_long.cpp b/src/KSPACE/pair_coul_long.cpp
index c1cd14523e59b50b204f6b41821a2a460c635a53..988694153afa952e71a403f859a0106860496a5a 100644
--- a/src/KSPACE/pair_coul_long.cpp
+++ b/src/KSPACE/pair_coul_long.cpp
@@ -215,8 +215,8 @@ void PairCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/KSPACE/pair_lj_charmm_coul_long.cpp b/src/KSPACE/pair_lj_charmm_coul_long.cpp
index 7fad9c5eef45dc79b4ee47a49671da08a173c2e9..bd020a439a67b95c99820bc95a229999eed5a55d 100644
--- a/src/KSPACE/pair_lj_charmm_coul_long.cpp
+++ b/src/KSPACE/pair_lj_charmm_coul_long.cpp
@@ -650,8 +650,8 @@ void PairLJCharmmCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/KSPACE/pair_lj_cut_coul_long.cpp b/src/KSPACE/pair_lj_cut_coul_long.cpp
index 866db6ba659120ab6b17da7a864bea10088c13b9..0d0649fe6a2d4b757b864335f2ec896b58684023 100644
--- a/src/KSPACE/pair_lj_cut_coul_long.cpp
+++ b/src/KSPACE/pair_lj_cut_coul_long.cpp
@@ -624,8 +624,8 @@ void PairLJCutCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/KSPACE/pair_lj_long_coul_long.cpp b/src/KSPACE/pair_lj_long_coul_long.cpp
index 3675ea76dc36668ff2865f61a59d5c26483c03aa..3b4ed79786fcaff010c820f756eb1df2a05af1d2 100644
--- a/src/KSPACE/pair_lj_long_coul_long.cpp
+++ b/src/KSPACE/pair_lj_long_coul_long.cpp
@@ -192,8 +192,8 @@ void PairLJLongCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/MANYBODY/pair_eam.cpp b/src/MANYBODY/pair_eam.cpp
index 53fbef60e85af97803706f52679675cf8c9d872c..880b8f62fce5bb8d33c0869be1810e5689213f86 100644
--- a/src/MANYBODY/pair_eam.cpp
+++ b/src/MANYBODY/pair_eam.cpp
@@ -366,8 +366,8 @@ void PairEAM::coeff(int narg, char **arg)
// parse pair of atom types
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
// read funcfl file if hasn't already been read
// store filename in Funcfl data struct
diff --git a/src/MC/pair_dsmc.cpp b/src/MC/pair_dsmc.cpp
index 7aa0237a91ffea889079984f185de897567c2a81..344faf87f62bff5da684c9b77f694877782513e4 100644
--- a/src/MC/pair_dsmc.cpp
+++ b/src/MC/pair_dsmc.cpp
@@ -245,8 +245,8 @@ void PairDSMC::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double sigma_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MISC/compute_ti.cpp b/src/MISC/compute_ti.cpp
index 02f16085057bd82c3dc29f0789907e3a586ab0ce..eca136c9359c388566ac8342552df6f61a0cb98d 100644
--- a/src/MISC/compute_ti.cpp
+++ b/src/MISC/compute_ti.cpp
@@ -79,7 +79,7 @@ ComputeTI::ComputeTI(LAMMPS *lmp, int narg, char **arg) :
int n = strlen(arg[iarg]) + 1;
pstyle[nterms] = new char[n];
strcpy(pstyle[nterms],arg[iarg]);
- force->bounds(arg[iarg+1],atom->ntypes,ilo[nterms],ihi[nterms]);
+ force->bounds(FLERR,arg[iarg+1],atom->ntypes,ilo[nterms],ihi[nterms]);
iarg += 1;
if (strstr(arg[iarg+1],"v_") == arg[iarg+1]) {
diff --git a/src/MISC/pair_nm_cut.cpp b/src/MISC/pair_nm_cut.cpp
index f0c9806fb8183722a0188e2bbb6028117ca9da9c..467be1b7bedb369c86d93bced6f0f28e365f8bcb 100644
--- a/src/MISC/pair_nm_cut.cpp
+++ b/src/MISC/pair_nm_cut.cpp
@@ -203,8 +203,8 @@ void PairNMCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double e0_one = force->numeric(FLERR,arg[2]);
double r0_one = force->numeric(FLERR,arg[3]);
diff --git a/src/MISC/pair_nm_cut_coul_cut.cpp b/src/MISC/pair_nm_cut_coul_cut.cpp
index 61326ee1eb1d3f947414edba1929eb889774ebfc..86fa09f176700239016bfca52e83f4e7859c0349 100644
--- a/src/MISC/pair_nm_cut_coul_cut.cpp
+++ b/src/MISC/pair_nm_cut_coul_cut.cpp
@@ -232,8 +232,8 @@ void PairNMCutCoulCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double e0_one = force->numeric(FLERR,arg[2]);
double r0_one = force->numeric(FLERR,arg[3]);
diff --git a/src/MISC/pair_nm_cut_coul_long.cpp b/src/MISC/pair_nm_cut_coul_long.cpp
index 0cd12ac0358d08db8b55ad25a9e62c500b3307ec..c186d19539da00ec8448b1d7e8e5e1724ccb7b93 100644
--- a/src/MISC/pair_nm_cut_coul_long.cpp
+++ b/src/MISC/pair_nm_cut_coul_long.cpp
@@ -271,8 +271,8 @@ void PairNMCutCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double e0_one = force->numeric(FLERR,arg[2]);
double r0_one = force->numeric(FLERR,arg[3]);
diff --git a/src/MOLECULE/angle_charmm.cpp b/src/MOLECULE/angle_charmm.cpp
index c2d44b26dde801fbf23be044c62dd2beeb5b1f7f..c7284add8ccfa765b96a70e27d861b454de67cc2 100644
--- a/src/MOLECULE/angle_charmm.cpp
+++ b/src/MOLECULE/angle_charmm.cpp
@@ -196,7 +196,7 @@ void AngleCharmm::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double theta0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/angle_cosine.cpp b/src/MOLECULE/angle_cosine.cpp
index 135d4bbf8f103edbb936ef4c2eb099f90c52daa0..f064a34f20d17aa02f07c90ef152ab87246d8357 100644
--- a/src/MOLECULE/angle_cosine.cpp
+++ b/src/MOLECULE/angle_cosine.cpp
@@ -156,7 +156,7 @@ void AngleCosine::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
diff --git a/src/MOLECULE/angle_cosine_periodic.cpp b/src/MOLECULE/angle_cosine_periodic.cpp
index dbd93b3b7b2da57486766f1fb4353f793f108481..77fb335b785c769c6494012e2e5798b29ca44555 100644
--- a/src/MOLECULE/angle_cosine_periodic.cpp
+++ b/src/MOLECULE/angle_cosine_periodic.cpp
@@ -201,7 +201,7 @@ void AngleCosinePeriodic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double c_one = force->numeric(FLERR,arg[1]);
int b_one = force->inumeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/angle_cosine_squared.cpp b/src/MOLECULE/angle_cosine_squared.cpp
index 6d073f156eb1035d7f3bf1bd3456160e257e7c21..bc38d76dcdb4fded224f9994d08f756c8532255b 100644
--- a/src/MOLECULE/angle_cosine_squared.cpp
+++ b/src/MOLECULE/angle_cosine_squared.cpp
@@ -168,7 +168,7 @@ void AngleCosineSquared::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double theta0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/angle_harmonic.cpp b/src/MOLECULE/angle_harmonic.cpp
index cffc3182ddf8e20374ea3649c4d634a8b7fa0fbf..7e8b57b04fad7748eeaf0ce0f1ee75a743992425 100644
--- a/src/MOLECULE/angle_harmonic.cpp
+++ b/src/MOLECULE/angle_harmonic.cpp
@@ -168,7 +168,7 @@ void AngleHarmonic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double theta0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/angle_table.cpp b/src/MOLECULE/angle_table.cpp
index 9dd9925725f0540503554d8cef23051b9c6df3b4..4d9007adb7263014f9ed1c795e7bf40afa274683 100644
--- a/src/MOLECULE/angle_table.cpp
+++ b/src/MOLECULE/angle_table.cpp
@@ -215,7 +215,7 @@ void AngleTable::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/MOLECULE/bond_fene.cpp b/src/MOLECULE/bond_fene.cpp
index 66c5c08174f42c69c7ad7cc2ec18c0194046da5b..5d58f2683a0a4969663776323761d101d50bb31b 100644
--- a/src/MOLECULE/bond_fene.cpp
+++ b/src/MOLECULE/bond_fene.cpp
@@ -154,7 +154,7 @@ void BondFENE::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double r0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/bond_fene_expand.cpp b/src/MOLECULE/bond_fene_expand.cpp
index f65632cc3f3400df7555ab1c9aa8838a8ddc432c..033f8d6502ef91d2fcb53b9b25a972bfa01eb391 100644
--- a/src/MOLECULE/bond_fene_expand.cpp
+++ b/src/MOLECULE/bond_fene_expand.cpp
@@ -160,7 +160,7 @@ void BondFENEExpand::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double r0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/bond_harmonic.cpp b/src/MOLECULE/bond_harmonic.cpp
index 8c5663bf79c1a97f70c857568979b6dbd1ff8572..f164a51de4309d57cff2d06a90ffa6b67c203639 100644
--- a/src/MOLECULE/bond_harmonic.cpp
+++ b/src/MOLECULE/bond_harmonic.cpp
@@ -121,7 +121,7 @@ void BondHarmonic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double r0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/bond_morse.cpp b/src/MOLECULE/bond_morse.cpp
index da2717899a670b60f206d8809585c40c84da011b..3204a0ca6183b003e0ec012083a8cd817f6e6f2c 100644
--- a/src/MOLECULE/bond_morse.cpp
+++ b/src/MOLECULE/bond_morse.cpp
@@ -126,7 +126,7 @@ void BondMorse::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double d0_one = force->numeric(FLERR,arg[1]);
double alpha_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/bond_nonlinear.cpp b/src/MOLECULE/bond_nonlinear.cpp
index 9da2993c54b3b048cd2de44595e05429ad03dd0d..edd9dcd4e5c9e8e8be20a2ea85c1ee5b3fb13a3a 100644
--- a/src/MOLECULE/bond_nonlinear.cpp
+++ b/src/MOLECULE/bond_nonlinear.cpp
@@ -123,7 +123,7 @@ void BondNonlinear::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double epsilon_one = force->numeric(FLERR,arg[1]);
double r0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/bond_quartic.cpp b/src/MOLECULE/bond_quartic.cpp
index c92048b2ca354349adb4690516514fb17a5fbcc1..2a3e395a77f81f9c92b59f6137d0cfbdc6c13a26 100644
--- a/src/MOLECULE/bond_quartic.cpp
+++ b/src/MOLECULE/bond_quartic.cpp
@@ -203,7 +203,7 @@ void BondQuartic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double b1_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/bond_table.cpp b/src/MOLECULE/bond_table.cpp
index 858bc83a3665fb0b324221b80b99034e91430efa..38cbe7e406823af76b118c6eb0529515e8d7971d 100644
--- a/src/MOLECULE/bond_table.cpp
+++ b/src/MOLECULE/bond_table.cpp
@@ -168,7 +168,7 @@ void BondTable::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/MOLECULE/dihedral_charmm.cpp b/src/MOLECULE/dihedral_charmm.cpp
index 8b6909f1ad3df52c3f6ce18a61f9280ced338771..b9d1c440d44cbb7ef5acb06a219ef1007191be9b 100644
--- a/src/MOLECULE/dihedral_charmm.cpp
+++ b/src/MOLECULE/dihedral_charmm.cpp
@@ -330,7 +330,7 @@ void DihedralCharmm::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
// require integer values of shift for backwards compatibility
// arbitrary phase angle shift could be allowed, but would break
diff --git a/src/MOLECULE/dihedral_harmonic.cpp b/src/MOLECULE/dihedral_harmonic.cpp
index f9ea31ac7cfcb521064151eec30444abfa7f3ca2..82c5fe31507d667359a208a00bb43056af1eef6a 100644
--- a/src/MOLECULE/dihedral_harmonic.cpp
+++ b/src/MOLECULE/dihedral_harmonic.cpp
@@ -276,7 +276,7 @@ void DihedralHarmonic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
int sign_one = force->inumeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/dihedral_helix.cpp b/src/MOLECULE/dihedral_helix.cpp
index b0565e29fb090fc57f54fba7cc9af2e107a9a278..43794bd6aabfacca837f8b9a57e10452b46f2a70 100644
--- a/src/MOLECULE/dihedral_helix.cpp
+++ b/src/MOLECULE/dihedral_helix.cpp
@@ -284,7 +284,7 @@ void DihedralHelix::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
double aphi_one = force->numeric(FLERR,arg[1]);
double bphi_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/dihedral_multi_harmonic.cpp b/src/MOLECULE/dihedral_multi_harmonic.cpp
index 0842595fba01c500aeb8ea9e308b22ed1e640c24..3a4be46fc7876360636365724968ca6fc007a710 100644
--- a/src/MOLECULE/dihedral_multi_harmonic.cpp
+++ b/src/MOLECULE/dihedral_multi_harmonic.cpp
@@ -273,7 +273,7 @@ void DihedralMultiHarmonic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
double a1_one = force->numeric(FLERR,arg[1]);
double a2_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/dihedral_opls.cpp b/src/MOLECULE/dihedral_opls.cpp
index 609fac74bdd21c05d2e2b179eb9fecb6f0327b55..c3a7ee6aa7f9a4388f058b39587b4e695702020b 100644
--- a/src/MOLECULE/dihedral_opls.cpp
+++ b/src/MOLECULE/dihedral_opls.cpp
@@ -289,7 +289,7 @@ void DihedralOPLS::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
double k1_one = force->numeric(FLERR,arg[1]);
double k2_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/improper_cvff.cpp b/src/MOLECULE/improper_cvff.cpp
index 28b3612c7838ca497356b5cab562f2d4dc24f445..45740d7880f0b57409d8de9c73872d7176861cf7 100644
--- a/src/MOLECULE/improper_cvff.cpp
+++ b/src/MOLECULE/improper_cvff.cpp
@@ -298,7 +298,7 @@ void ImproperCvff::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
int sign_one = force->inumeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/improper_harmonic.cpp b/src/MOLECULE/improper_harmonic.cpp
index 56ddce6f5bc772b1fc99a28f547441f22a8ea69d..c482055c76f8c57b4acb7e1dee7af71066663d11 100644
--- a/src/MOLECULE/improper_harmonic.cpp
+++ b/src/MOLECULE/improper_harmonic.cpp
@@ -240,7 +240,7 @@ void ImproperHarmonic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double chi_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/improper_umbrella.cpp b/src/MOLECULE/improper_umbrella.cpp
index a2f176e732d1e05482ba7e0d1754eb1174e0c139..8a9b0b44345711174b26c84799ac13b7793c5081 100644
--- a/src/MOLECULE/improper_umbrella.cpp
+++ b/src/MOLECULE/improper_umbrella.cpp
@@ -278,7 +278,7 @@ void ImproperUmbrella::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double w_one = force->numeric(FLERR,arg[2]);
diff --git a/src/MOLECULE/pair_hbond_dreiding_lj.cpp b/src/MOLECULE/pair_hbond_dreiding_lj.cpp
index af54c9bdfa736b7d51cb4bd6476ff3c28b20292d..ecb4883cb364b07afbb56208e62c73afe91c50b2 100644
--- a/src/MOLECULE/pair_hbond_dreiding_lj.cpp
+++ b/src/MOLECULE/pair_hbond_dreiding_lj.cpp
@@ -321,9 +321,9 @@ void PairHbondDreidingLJ::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi,klo,khi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
- force->bounds(arg[2],atom->ntypes,klo,khi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[2],atom->ntypes,klo,khi);
int donor_flag;
if (strcmp(arg[3],"i") == 0) donor_flag = 0;
diff --git a/src/MOLECULE/pair_hbond_dreiding_morse.cpp b/src/MOLECULE/pair_hbond_dreiding_morse.cpp
index ab793405dde71defb98c792fdfe9ef26e1252a32..5c61ea4a45316c50933e9dc37ccfbf4f7e71aa0c 100644
--- a/src/MOLECULE/pair_hbond_dreiding_morse.cpp
+++ b/src/MOLECULE/pair_hbond_dreiding_morse.cpp
@@ -245,9 +245,9 @@ void PairHbondDreidingMorse::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi,klo,khi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
- force->bounds(arg[2],atom->ntypes,klo,khi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[2],atom->ntypes,klo,khi);
int donor_flag;
if (strcmp(arg[3],"i") == 0) donor_flag = 0;
diff --git a/src/MOLECULE/pair_lj_charmm_coul_charmm.cpp b/src/MOLECULE/pair_lj_charmm_coul_charmm.cpp
index f12bc8f3b833acbde215c404dc12a486888e7ae1..3700de00b0fe380a96abaf358ab0e06f0416ab5a 100644
--- a/src/MOLECULE/pair_lj_charmm_coul_charmm.cpp
+++ b/src/MOLECULE/pair_lj_charmm_coul_charmm.cpp
@@ -247,8 +247,8 @@ void PairLJCharmmCoulCharmm::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/MOLECULE/pair_lj_cut_tip4p_cut.cpp b/src/MOLECULE/pair_lj_cut_tip4p_cut.cpp
index 0f5a1f68d11565f7e86463c4959eb886397c04cc..a9e00e80b40cbf73faf80628194f2bba121e4504 100644
--- a/src/MOLECULE/pair_lj_cut_tip4p_cut.cpp
+++ b/src/MOLECULE/pair_lj_cut_tip4p_cut.cpp
@@ -447,8 +447,8 @@ void PairLJCutTIP4PCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/MOLECULE/pair_tip4p_cut.cpp b/src/MOLECULE/pair_tip4p_cut.cpp
index 47aa7730474d9f81004b63ee525b2e78656ac00f..dd48637f23037f5c47174684b79aeb5bafcf6d9c 100644
--- a/src/MOLECULE/pair_tip4p_cut.cpp
+++ b/src/MOLECULE/pair_tip4p_cut.cpp
@@ -387,8 +387,8 @@ void PairTIP4PCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/PERI/pair_peri_eps.cpp b/src/PERI/pair_peri_eps.cpp
index ecf6f0a3fa8d9c40bf6fdd0326fd488fab8c2f63..b5807c0e3c19c2c148e6b8f43033bb36a145722d 100644
--- a/src/PERI/pair_peri_eps.cpp
+++ b/src/PERI/pair_peri_eps.cpp
@@ -450,8 +450,8 @@ void PairPeriEPS::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double bulkmodulus_one = atof(arg[2]);
double shearmodulus_one = atof(arg[3]);
diff --git a/src/PERI/pair_peri_lps.cpp b/src/PERI/pair_peri_lps.cpp
index 7b2ccec92cce19f3e29c84669eb2a7dab94a96e9..cd88b41825d3b15125116f9ed0d85640c5c76e86 100644
--- a/src/PERI/pair_peri_lps.cpp
+++ b/src/PERI/pair_peri_lps.cpp
@@ -375,8 +375,8 @@ void PairPeriLPS::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double bulkmodulus_one = force->numeric(FLERR,arg[2]);
double shearmodulus_one = force->numeric(FLERR,arg[3]);
diff --git a/src/PERI/pair_peri_pmb.cpp b/src/PERI/pair_peri_pmb.cpp
index bc971cb2442ff175a9e01433701fea438dd2a974..9e598289bd01842bb255a5501f5163fad9488d0e 100644
--- a/src/PERI/pair_peri_pmb.cpp
+++ b/src/PERI/pair_peri_pmb.cpp
@@ -312,8 +312,8 @@ void PairPeriPMB::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double kspring_one = force->numeric(FLERR,arg[2]);
double cut_one = force->numeric(FLERR,arg[3]);
diff --git a/src/PERI/pair_peri_ves.cpp b/src/PERI/pair_peri_ves.cpp
index c34fb87e20f2f0f6df8009d7c1199694396f240e..87213b5b59c10de830ae0e131335f5f99779b0f9 100644
--- a/src/PERI/pair_peri_ves.cpp
+++ b/src/PERI/pair_peri_ves.cpp
@@ -426,8 +426,8 @@ void PairPeriVES::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double bulkmodulus_one = atof(arg[2]);
double shearmodulus_one = atof(arg[3]);
diff --git a/src/RIGID/fix_rigid.cpp b/src/RIGID/fix_rigid.cpp
index 2f2f84570e0b438566618e12dd4397350d011e61..8c1ecdacc9a2b1047fc21647ea3b598af57ce959 100644
--- a/src/RIGID/fix_rigid.cpp
+++ b/src/RIGID/fix_rigid.cpp
@@ -282,7 +282,7 @@ FixRigid::FixRigid(LAMMPS *lmp, int narg, char **arg) :
if (iarg+5 > narg) error->all(FLERR,"Illegal fix rigid command");
int mlo,mhi;
- force->bounds(arg[iarg+1],nbody,mlo,mhi);
+ force->bounds(FLERR,arg[iarg+1],nbody,mlo,mhi);
double xflag,yflag,zflag;
if (strcmp(arg[iarg+2],"off") == 0) xflag = 0.0;
@@ -313,7 +313,7 @@ FixRigid::FixRigid(LAMMPS *lmp, int narg, char **arg) :
if (iarg+5 > narg) error->all(FLERR,"Illegal fix rigid command");
int mlo,mhi;
- force->bounds(arg[iarg+1],nbody,mlo,mhi);
+ force->bounds(FLERR,arg[iarg+1],nbody,mlo,mhi);
double xflag,yflag,zflag;
if (strcmp(arg[iarg+2],"off") == 0) xflag = 0.0;
diff --git a/src/USER-AWPMD/pair_awpmd_cut.cpp b/src/USER-AWPMD/pair_awpmd_cut.cpp
index 29feb37ca15eeea7f20701cad24a45fa21642aa2..cd89c3984d536b9373e51c0ee868dac30b68e3d5 100644
--- a/src/USER-AWPMD/pair_awpmd_cut.cpp
+++ b/src/USER-AWPMD/pair_awpmd_cut.cpp
@@ -494,8 +494,8 @@ void PairAWPMDCut::coeff(int narg, char **arg)
}
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_one = cut_global;
if (narg == 3) cut_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-CG-CMM/angle_sdk.cpp b/src/USER-CG-CMM/angle_sdk.cpp
index 31136e5cd93e9392b8bf59fbd93790be5232779f..cc54985991008fccc7ec04081f45931a7df745b4 100644
--- a/src/USER-CG-CMM/angle_sdk.cpp
+++ b/src/USER-CG-CMM/angle_sdk.cpp
@@ -240,7 +240,7 @@ void AngleSDK::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double theta0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-CG-CMM/pair_lj_sdk.cpp b/src/USER-CG-CMM/pair_lj_sdk.cpp
index 1ab93447f4bd6d17f8b56d157d4e3c555df83262..665f188ce9d740d50a86fdc879f91420c15875f6 100644
--- a/src/USER-CG-CMM/pair_lj_sdk.cpp
+++ b/src/USER-CG-CMM/pair_lj_sdk.cpp
@@ -263,8 +263,8 @@ void PairLJSDK::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int lj_type_one = find_lj_type(arg[2],lj_type_list);
if (lj_type_one == LJ_NOT_SET)
diff --git a/src/USER-CG-CMM/pair_lj_sdk_coul_long.cpp b/src/USER-CG-CMM/pair_lj_sdk_coul_long.cpp
index 1549acfc283ae4a5004a5d506417cb791b90273b..5e4a0db31c721b19f262c99396f6915cd7c8d789 100644
--- a/src/USER-CG-CMM/pair_lj_sdk_coul_long.cpp
+++ b/src/USER-CG-CMM/pair_lj_sdk_coul_long.cpp
@@ -324,8 +324,8 @@ void PairLJSDKCoulLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int lj_type_one = find_lj_type(arg[2],lj_type_list);
if (lj_type_one == LJ_NOT_SET)
diff --git a/src/USER-DPD/pair_dpd_fdt.cpp b/src/USER-DPD/pair_dpd_fdt.cpp
index d151b77c8565a82c28fc4d81bad688097b6fe00b..3b5804ff8ec0fb44d70cdfcf9446c778310a6806 100644
--- a/src/USER-DPD/pair_dpd_fdt.cpp
+++ b/src/USER-DPD/pair_dpd_fdt.cpp
@@ -280,8 +280,8 @@ void PairDPDfdt::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a0_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-DPD/pair_dpd_fdt_energy.cpp b/src/USER-DPD/pair_dpd_fdt_energy.cpp
index 902ca005256829733278855d204b53ab4baeeb85..84fa352a45548c9b655da83e3e13f26b1ca93702 100644
--- a/src/USER-DPD/pair_dpd_fdt_energy.cpp
+++ b/src/USER-DPD/pair_dpd_fdt_energy.cpp
@@ -364,8 +364,8 @@ void PairDPDfdtEnergy::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a0_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-DPD/pair_exp6_rx.cpp b/src/USER-DPD/pair_exp6_rx.cpp
index ad2dbf14657287953de18426381b5184996765e7..9af28026ae166e2c984e66cdb647c0d6d3aca814 100644
--- a/src/USER-DPD/pair_exp6_rx.cpp
+++ b/src/USER-DPD/pair_exp6_rx.cpp
@@ -558,8 +558,8 @@ void PairExp6rx::coeff(int narg, char **arg)
int ilo,ihi,jlo,jhi;
int n;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
nspecies = atom->nspecies_dpd;
if(nspecies==0) error->all(FLERR,"There are no rx species specified.");
diff --git a/src/USER-DPD/pair_multi_lucy.cpp b/src/USER-DPD/pair_multi_lucy.cpp
index a063eff1958435b90cc45160a1b597f6840a37a2..a4ea8cad229946605eb8b904842cefa47b822bd4 100644
--- a/src/USER-DPD/pair_multi_lucy.cpp
+++ b/src/USER-DPD/pair_multi_lucy.cpp
@@ -267,8 +267,8 @@ void PairMultiLucy::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/USER-DPD/pair_multi_lucy_rx.cpp b/src/USER-DPD/pair_multi_lucy_rx.cpp
index cf10a70756f663295ebbace3121cb39ab7c5abe2..586e58908da4292a3b0d0e041e2bc80163922611 100644
--- a/src/USER-DPD/pair_multi_lucy_rx.cpp
+++ b/src/USER-DPD/pair_multi_lucy_rx.cpp
@@ -353,8 +353,8 @@ void PairMultiLucyRX::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/USER-DPD/pair_table_rx.cpp b/src/USER-DPD/pair_table_rx.cpp
index 74b9222a0b7f290fd86a8d46dc781eb6dd9df82b..902d0e5bb4397144599637f8c03c977b7784ba6c 100644
--- a/src/USER-DPD/pair_table_rx.cpp
+++ b/src/USER-DPD/pair_table_rx.cpp
@@ -326,8 +326,8 @@ void PairTableRX::coeff(int narg, char **arg)
if (!rx_flag) error->all(FLERR,"PairTableRX requires a fix rx command.");
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/USER-DRUDE/pair_lj_cut_thole_long.cpp b/src/USER-DRUDE/pair_lj_cut_thole_long.cpp
index efa11c6c5e5d00322f0b95a4b2a54b1ce31c2058..726baae2400c077075159c2b95cc666cf70124c3 100644
--- a/src/USER-DRUDE/pair_lj_cut_thole_long.cpp
+++ b/src/USER-DRUDE/pair_lj_cut_thole_long.cpp
@@ -314,8 +314,8 @@ void PairLJCutTholeLong::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-DRUDE/pair_thole.cpp b/src/USER-DRUDE/pair_thole.cpp
index 13b1265a7edfc3109533aa8cb1d835693f9eabf6..0ed94ebbc919b3b5be5d691caa15d86edede5989 100644
--- a/src/USER-DRUDE/pair_thole.cpp
+++ b/src/USER-DRUDE/pair_thole.cpp
@@ -218,8 +218,8 @@ void PairThole::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double polar_one = force->numeric(FLERR,arg[2]);
double thole_one = thole_global;
diff --git a/src/USER-EFF/pair_eff_cut.cpp b/src/USER-EFF/pair_eff_cut.cpp
index b12e1f43b3e28741290881c1035a09cf4db0d950..66f59c86c31aee1fe3b1d6ffc0f4bc352b16000e 100644
--- a/src/USER-EFF/pair_eff_cut.cpp
+++ b/src/USER-EFF/pair_eff_cut.cpp
@@ -891,8 +891,8 @@ void PairEffCut::coeff(int narg, char **arg)
if ((strcmp(arg[0],"*") == 0) || (strcmp(arg[1],"*") == 0)) {
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_one = cut_global;
if (narg == 3) cut_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-FEP/compute_fep.cpp b/src/USER-FEP/compute_fep.cpp
index d3c3bf9409ec0946761ac5112633f174431df4ad..bfa3bf7e2c91648cb08b74cca10e394da9a4a9e9 100644
--- a/src/USER-FEP/compute_fep.cpp
+++ b/src/USER-FEP/compute_fep.cpp
@@ -95,9 +95,9 @@ ComputeFEP::ComputeFEP(LAMMPS *lmp, int narg, char **arg) :
n = strlen(arg[iarg+2]) + 1;
perturb[npert].pparam = new char[n];
strcpy(perturb[npert].pparam,arg[iarg+2]);
- force->bounds(arg[iarg+3],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+3],atom->ntypes,
perturb[npert].ilo,perturb[npert].ihi);
- force->bounds(arg[iarg+4],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+4],atom->ntypes,
perturb[npert].jlo,perturb[npert].jhi);
if (strstr(arg[iarg+5],"v_") == arg[iarg+5]) {
n = strlen(&arg[iarg+5][2]) + 1;
@@ -112,7 +112,7 @@ ComputeFEP::ComputeFEP(LAMMPS *lmp, int narg, char **arg) :
perturb[npert].aparam = CHARGE;
chgflag = 1;
} else error->all(FLERR,"Illegal atom argument in compute fep");
- force->bounds(arg[iarg+2],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+2],atom->ntypes,
perturb[npert].ilo,perturb[npert].ihi);
if (strstr(arg[iarg+3],"v_") == arg[iarg+3]) {
int n = strlen(&arg[iarg+3][2]) + 1;
diff --git a/src/USER-FEP/fix_adapt_fep.cpp b/src/USER-FEP/fix_adapt_fep.cpp
index 3e0e950cb777ebcde5869a65aec3b8e9e3926ac5..1b7fc559420f2ee6bd932d872c50524cebe208d9 100644
--- a/src/USER-FEP/fix_adapt_fep.cpp
+++ b/src/USER-FEP/fix_adapt_fep.cpp
@@ -95,9 +95,9 @@ FixAdaptFEP::FixAdaptFEP(LAMMPS *lmp, int narg, char **arg) :
n = strlen(arg[iarg+2]) + 1;
adapt[nadapt].pparam = new char[n];
strcpy(adapt[nadapt].pparam,arg[iarg+2]);
- force->bounds(arg[iarg+3],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+3],atom->ntypes,
adapt[nadapt].ilo,adapt[nadapt].ihi);
- force->bounds(arg[iarg+4],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+4],atom->ntypes,
adapt[nadapt].jlo,adapt[nadapt].jhi);
if (strstr(arg[iarg+5],"v_") == arg[iarg+5]) {
n = strlen(&arg[iarg+5][2]) + 1;
@@ -126,7 +126,7 @@ FixAdaptFEP::FixAdaptFEP(LAMMPS *lmp, int narg, char **arg) :
adapt[nadapt].aparam = CHARGE;
chgflag = 1;
} else error->all(FLERR,"Illegal fix adapt/fep command");
- force->bounds(arg[iarg+2],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+2],atom->ntypes,
adapt[nadapt].ilo,adapt[nadapt].ihi);
if (strstr(arg[iarg+3],"v_") == arg[iarg+3]) {
int n = strlen(&arg[iarg+3][2]) + 1;
diff --git a/src/USER-FEP/pair_coul_cut_soft.cpp b/src/USER-FEP/pair_coul_cut_soft.cpp
index b82abc1ed165ee2d6bc5eae94a64c0bef66fd2a6..2c675c607f4fee76bc93a8830a640e4c6f4a0941 100644
--- a/src/USER-FEP/pair_coul_cut_soft.cpp
+++ b/src/USER-FEP/pair_coul_cut_soft.cpp
@@ -184,8 +184,8 @@ void PairCoulCutSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double lambda_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-FEP/pair_coul_long_soft.cpp b/src/USER-FEP/pair_coul_long_soft.cpp
index d20962810497339d133a56a595390d5985783e5d..3d24a997ddf0162af51b1a89a93dc0e525d4cf30 100644
--- a/src/USER-FEP/pair_coul_long_soft.cpp
+++ b/src/USER-FEP/pair_coul_long_soft.cpp
@@ -206,8 +206,8 @@ void PairCoulLongSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double lambda_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-FEP/pair_lj_charmm_coul_long_soft.cpp b/src/USER-FEP/pair_lj_charmm_coul_long_soft.cpp
index d2a9a04ba61eddab1029a2e8f8d0a15d80719e65..81b82e97743c2c8c088c6b036cea040d15be6df4 100644
--- a/src/USER-FEP/pair_lj_charmm_coul_long_soft.cpp
+++ b/src/USER-FEP/pair_lj_charmm_coul_long_soft.cpp
@@ -651,8 +651,8 @@ void PairLJCharmmCoulLongSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-FEP/pair_lj_cut_coul_cut_soft.cpp b/src/USER-FEP/pair_lj_cut_coul_cut_soft.cpp
index 631d608afb03a4945f65d0c3cc0ffd9c8f17c89e..16da07a657214b2b3be32010ea058dfe4a029d3a 100644
--- a/src/USER-FEP/pair_lj_cut_coul_cut_soft.cpp
+++ b/src/USER-FEP/pair_lj_cut_coul_cut_soft.cpp
@@ -230,8 +230,8 @@ void PairLJCutCoulCutSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-FEP/pair_lj_cut_coul_long_soft.cpp b/src/USER-FEP/pair_lj_cut_coul_long_soft.cpp
index 1b3c2b5eceba08f9d2af7bdd2907e9c037264ee3..6636e727159be3c3a85c69a5e57a48d08305b72c 100644
--- a/src/USER-FEP/pair_lj_cut_coul_long_soft.cpp
+++ b/src/USER-FEP/pair_lj_cut_coul_long_soft.cpp
@@ -598,8 +598,8 @@ void PairLJCutCoulLongSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-FEP/pair_lj_cut_soft.cpp b/src/USER-FEP/pair_lj_cut_soft.cpp
index 7f04158c1775adc23c2eafd46b9aba65085066ca..3798b279369e56a98cf5584b5eb6c3483e3af7e3 100644
--- a/src/USER-FEP/pair_lj_cut_soft.cpp
+++ b/src/USER-FEP/pair_lj_cut_soft.cpp
@@ -478,8 +478,8 @@ void PairLJCutSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-FEP/pair_morse_soft.cpp b/src/USER-FEP/pair_morse_soft.cpp
index 877fa74d69d43c419035869d0c066c4919c35df7..6c86d8916fb1611c001f7f13c8fff46fe5390b37 100644
--- a/src/USER-FEP/pair_morse_soft.cpp
+++ b/src/USER-FEP/pair_morse_soft.cpp
@@ -178,8 +178,8 @@ void PairMorseSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double d0_one = force->numeric(FLERR,arg[2]);
double alpha_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-LB/fix_lb_rigid_pc_sphere.cpp b/src/USER-LB/fix_lb_rigid_pc_sphere.cpp
index 54637ce349f1fa2ae834a96a7fdc8558913fecc8..b5383fe1a6743dcb987e8912367ab9afb7c1675f 100644
--- a/src/USER-LB/fix_lb_rigid_pc_sphere.cpp
+++ b/src/USER-LB/fix_lb_rigid_pc_sphere.cpp
@@ -235,7 +235,7 @@ FixLbRigidPCSphere::FixLbRigidPCSphere(LAMMPS *lmp, int narg, char **arg) :
if (iarg+5 > narg) error->all(FLERR,"Illegal fix lb/rigid/pc/sphere command");
int mlo,mhi;
- force->bounds(arg[iarg+1],nbody,mlo,mhi);
+ force->bounds(FLERR,arg[iarg+1],nbody,mlo,mhi);
double xflag,yflag,zflag;
if (strcmp(arg[iarg+2],"off") == 0) xflag = 0.0;
@@ -262,7 +262,7 @@ FixLbRigidPCSphere::FixLbRigidPCSphere(LAMMPS *lmp, int narg, char **arg) :
if (iarg+5 > narg) error->all(FLERR,"Illegal fix lb/rigid/pc/sphere command");
int mlo,mhi;
- force->bounds(arg[iarg+1],nbody,mlo,mhi);
+ force->bounds(FLERR,arg[iarg+1],nbody,mlo,mhi);
double xflag,yflag,zflag;
if (strcmp(arg[iarg+2],"off") == 0) xflag = 0.0;
diff --git a/src/USER-MISC/angle_cosine_shift.cpp b/src/USER-MISC/angle_cosine_shift.cpp
index 8e7c67fd70b1dd9e20ef74873c062c7581bd4ce5..66f5c82c8445185b07e618d77fcace93f03cba32 100644
--- a/src/USER-MISC/angle_cosine_shift.cpp
+++ b/src/USER-MISC/angle_cosine_shift.cpp
@@ -172,7 +172,7 @@ void AngleCosineShift::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double umin = force->numeric(FLERR,arg[1]);
double theta0 = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/angle_cosine_shift_exp.cpp b/src/USER-MISC/angle_cosine_shift_exp.cpp
index e19b9e20efaf5ed67c8d68dcb90a3221122b77c6..6d11c6894e0b7d1734b6393b582b8373c82f2c2d 100644
--- a/src/USER-MISC/angle_cosine_shift_exp.cpp
+++ b/src/USER-MISC/angle_cosine_shift_exp.cpp
@@ -194,7 +194,7 @@ void AngleCosineShiftExp::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double umin_ = force->numeric(FLERR,arg[1]);
double theta0_ = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/angle_dipole.cpp b/src/USER-MISC/angle_dipole.cpp
index 52e94ad3f9cc976fbb30bbf827072083a32026ef..990096ae7fce5cf63ccc9ed1797804f4f018b924 100644
--- a/src/USER-MISC/angle_dipole.cpp
+++ b/src/USER-MISC/angle_dipole.cpp
@@ -156,7 +156,7 @@ void AngleDipole::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double gamma0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/angle_fourier.cpp b/src/USER-MISC/angle_fourier.cpp
index dee6ed1301be486e968a85ef3851ded0410a0ebf..c4dbf3f08b4343b0ea3de63c1ae461706a3550d4 100644
--- a/src/USER-MISC/angle_fourier.cpp
+++ b/src/USER-MISC/angle_fourier.cpp
@@ -173,7 +173,7 @@ void AngleFourier::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double C0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/angle_fourier_simple.cpp b/src/USER-MISC/angle_fourier_simple.cpp
index 63991f320e1e6cb675dd61d4df93b2a4420218e4..bd4322fe78ca77e385d4c3668b1594408f9a55b9 100644
--- a/src/USER-MISC/angle_fourier_simple.cpp
+++ b/src/USER-MISC/angle_fourier_simple.cpp
@@ -188,7 +188,7 @@ void AngleFourierSimple::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double C_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/angle_quartic.cpp b/src/USER-MISC/angle_quartic.cpp
index a1a3e4cf3f4d0871a3b5afd6501373e9c48fe905..f9801edb0e5cb84e811a088689ecb72e4454f6e9 100644
--- a/src/USER-MISC/angle_quartic.cpp
+++ b/src/USER-MISC/angle_quartic.cpp
@@ -183,7 +183,7 @@ void AngleQuartic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double theta0_one = force->numeric(FLERR,arg[1]);
double k2_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/bond_harmonic_shift.cpp b/src/USER-MISC/bond_harmonic_shift.cpp
index 6f1b67e2459ae226ef252b20514f6e007506feab..c4ee00cb27164138fd77025423fc23fbf3cca21b 100644
--- a/src/USER-MISC/bond_harmonic_shift.cpp
+++ b/src/USER-MISC/bond_harmonic_shift.cpp
@@ -129,7 +129,7 @@ void BondHarmonicShift::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double Umin = force->numeric(FLERR,arg[1]); // energy at minimum
double r0_one = force->numeric(FLERR,arg[2]); // position of minimum
diff --git a/src/USER-MISC/bond_harmonic_shift_cut.cpp b/src/USER-MISC/bond_harmonic_shift_cut.cpp
index 1a1be169990a42f697d27722d296532d5a276a55..26fd32c1b596162444a759b7adceeb0ab749f66b 100644
--- a/src/USER-MISC/bond_harmonic_shift_cut.cpp
+++ b/src/USER-MISC/bond_harmonic_shift_cut.cpp
@@ -131,7 +131,7 @@ void BondHarmonicShiftCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double Umin = force->numeric(FLERR,arg[1]); // energy at minimum
double r0_one = force->numeric(FLERR,arg[2]); // position of minimum
diff --git a/src/USER-MISC/dihedral_cosine_shift_exp.cpp b/src/USER-MISC/dihedral_cosine_shift_exp.cpp
index 8706b5fb1a18670a68e915240789329b43c524ca..85ff4b27191e3b6d19bce335781f7bd715fbfff1 100644
--- a/src/USER-MISC/dihedral_cosine_shift_exp.cpp
+++ b/src/USER-MISC/dihedral_cosine_shift_exp.cpp
@@ -274,7 +274,7 @@ void DihedralCosineShiftExp::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
double umin_ = force->numeric(FLERR,arg[1]);
double theta0_ = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/dihedral_fourier.cpp b/src/USER-MISC/dihedral_fourier.cpp
index 93e818183da81490c1829c2497d99d03f6a81c10..86ab8894d8133fb3272cee4e7e4366eaa152aa9a 100644
--- a/src/USER-MISC/dihedral_fourier.cpp
+++ b/src/USER-MISC/dihedral_fourier.cpp
@@ -297,7 +297,7 @@ void DihedralFourier::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
// require integer values of shift for backwards compatibility
// arbitrary phase angle shift could be allowed, but would break
diff --git a/src/USER-MISC/dihedral_nharmonic.cpp b/src/USER-MISC/dihedral_nharmonic.cpp
index 67b363f5526528940074d3b868230a19d16b34c7..ba8af7e65b3b08e94507527c1390c0a39d14adfd 100644
--- a/src/USER-MISC/dihedral_nharmonic.cpp
+++ b/src/USER-MISC/dihedral_nharmonic.cpp
@@ -279,7 +279,7 @@ void DihedralNHarmonic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/USER-MISC/dihedral_quadratic.cpp b/src/USER-MISC/dihedral_quadratic.cpp
index 3a5eac6fa0b549f66138e5c68fbbb62d068be494..ac261153b4b8d4194144b7dfc52952682945fbfa 100644
--- a/src/USER-MISC/dihedral_quadratic.cpp
+++ b/src/USER-MISC/dihedral_quadratic.cpp
@@ -285,7 +285,7 @@ void DihedralQuadratic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double phi0_one= force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/dihedral_spherical.cpp b/src/USER-MISC/dihedral_spherical.cpp
index e6b431b63a338b6553f7e2637fd49653630f70f9..c3345453c9a9e76668bc593ab28debf31faf608f 100644
--- a/src/USER-MISC/dihedral_spherical.cpp
+++ b/src/USER-MISC/dihedral_spherical.cpp
@@ -681,7 +681,7 @@ void DihedralSpherical::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
int nterms_one = force->inumeric(FLERR,arg[1]);
diff --git a/src/USER-MISC/dihedral_table.cpp b/src/USER-MISC/dihedral_table.cpp
index 4e3b05bc4b4a751f10978a1b6c2051ea158761e3..f84814cc5e232408fc7cc3eb78d8b18610433e67 100644
--- a/src/USER-MISC/dihedral_table.cpp
+++ b/src/USER-MISC/dihedral_table.cpp
@@ -813,7 +813,7 @@ void DihedralTable::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/USER-MISC/improper_cossq.cpp b/src/USER-MISC/improper_cossq.cpp
index 77033381e98a3ab39bfbfb7fb5ea0a52b8728507..223df4a76d064abcfead348518d17f1034bbd1a9 100644
--- a/src/USER-MISC/improper_cossq.cpp
+++ b/src/USER-MISC/improper_cossq.cpp
@@ -271,7 +271,7 @@ void ImproperCossq::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double chi_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/improper_distance.cpp b/src/USER-MISC/improper_distance.cpp
index a98005a2474d3c474df01179c9856b4480df3266..ce6441d6746e3d4946bb8e21eaff27c39844cbe5 100644
--- a/src/USER-MISC/improper_distance.cpp
+++ b/src/USER-MISC/improper_distance.cpp
@@ -213,7 +213,7 @@ void ImproperDistance::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double chi_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/improper_fourier.cpp b/src/USER-MISC/improper_fourier.cpp
index ea5f78662e56d5cd5725899470545c64fb87222e..b139fdad569dc14d67560301c1550fb867204d16 100644
--- a/src/USER-MISC/improper_fourier.cpp
+++ b/src/USER-MISC/improper_fourier.cpp
@@ -282,7 +282,7 @@ void ImproperFourier::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double C0_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/improper_ring.cpp b/src/USER-MISC/improper_ring.cpp
index 4134b4b172449da5e49edee4c9d2e930daf83fa0..5b5341f5061b3cf113aed66fb668121aaa037b83 100644
--- a/src/USER-MISC/improper_ring.cpp
+++ b/src/USER-MISC/improper_ring.cpp
@@ -291,7 +291,7 @@ void ImproperRing ::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
double k_one = force->numeric(FLERR,arg[1]);
double chi_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-MISC/pair_buck_mdf.cpp b/src/USER-MISC/pair_buck_mdf.cpp
index 96f49e42dd42dd9b040d24ec291d53db593b5b0f..3a433b16bd1daaca6179f825a23247d84939a93c 100644
--- a/src/USER-MISC/pair_buck_mdf.cpp
+++ b/src/USER-MISC/pair_buck_mdf.cpp
@@ -212,8 +212,8 @@ void PairBuckMDF::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_coul_diel.cpp b/src/USER-MISC/pair_coul_diel.cpp
index 1913e48b67d705c2eba026725e7637aa41b6039e..a732ace1a0ab04d836e8c48cfc7687056283a33b 100644
--- a/src/USER-MISC/pair_coul_diel.cpp
+++ b/src/USER-MISC/pair_coul_diel.cpp
@@ -183,8 +183,8 @@ void PairCoulDiel::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
eps_s = force->numeric(FLERR,arg[2]);
double rme_one =force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_gauss_cut.cpp b/src/USER-MISC/pair_gauss_cut.cpp
index f9cdcb96e68782a8f4c593689df63eee9c7d4953..f44b1bbd2d0d16a2fcc249d33ea3d1ed675de768 100644
--- a/src/USER-MISC/pair_gauss_cut.cpp
+++ b/src/USER-MISC/pair_gauss_cut.cpp
@@ -190,8 +190,8 @@ void PairGaussCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double hgauss_one = force->numeric(FLERR,arg[2]);
double rmh_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_lennard_mdf.cpp b/src/USER-MISC/pair_lennard_mdf.cpp
index 81c4435963f4bc2e5e7b8886d45e6e091e5d03c9..3a81955199d1e19c8b7506b8e3f0405250b3b6a9 100644
--- a/src/USER-MISC/pair_lennard_mdf.cpp
+++ b/src/USER-MISC/pair_lennard_mdf.cpp
@@ -216,8 +216,8 @@ void PairLJ_AB_MDF::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double aparm_one = force->numeric(FLERR,arg[2]);
double bparm_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_list.cpp b/src/USER-MISC/pair_list.cpp
index 5f053926318d771715730ca0a6a5962ed6aa17f8..3adbe8b69d948097c275cd11504079c5a77bfa94 100644
--- a/src/USER-MISC/pair_list.cpp
+++ b/src/USER-MISC/pair_list.cpp
@@ -349,8 +349,8 @@ void PairList::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/USER-MISC/pair_lj_mdf.cpp b/src/USER-MISC/pair_lj_mdf.cpp
index 7c3b8d586e0155cdb8122c0b1c9262d6a9c75a01..3b52cf0b86f84227581daecab051d73899d9d5e1 100644
--- a/src/USER-MISC/pair_lj_mdf.cpp
+++ b/src/USER-MISC/pair_lj_mdf.cpp
@@ -216,8 +216,8 @@ void PairLJMDF::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_lj_sf.cpp b/src/USER-MISC/pair_lj_sf.cpp
index 69a7efcfde4d49bee3241e2b0dda5c91b4ea2537..32f45ff48d181b80e856a7b5d515ce85110fa740 100644
--- a/src/USER-MISC/pair_lj_sf.cpp
+++ b/src/USER-MISC/pair_lj_sf.cpp
@@ -197,8 +197,8 @@ void PairLJShiftedForce::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_lj_sf_dipole_sf.cpp b/src/USER-MISC/pair_lj_sf_dipole_sf.cpp
index b88f72bac0b7670e7acf315d952a3c07e7681c71..33f10f2f12ad51889ded6c00c1cf247d88e89d33 100644
--- a/src/USER-MISC/pair_lj_sf_dipole_sf.cpp
+++ b/src/USER-MISC/pair_lj_sf_dipole_sf.cpp
@@ -361,8 +361,8 @@ void PairLJSFDipoleSF::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_morse_smooth_linear.cpp b/src/USER-MISC/pair_morse_smooth_linear.cpp
index 844d0113bea3d8bd14502e0a8625eb8b1aee61df..ea33510b585334508b738834429fb0b674bf9909 100644
--- a/src/USER-MISC/pair_morse_smooth_linear.cpp
+++ b/src/USER-MISC/pair_morse_smooth_linear.cpp
@@ -186,8 +186,8 @@ void PairMorseSmoothLinear::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double d0_one = force->numeric(FLERR,arg[2]);
double alpha_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-MISC/pair_srp.cpp b/src/USER-MISC/pair_srp.cpp
index 5e254b74b7433592ea379b1d38e581dc64fc0765..18ea4dc332ed971477b83cb98b0d8ee6674eadad 100644
--- a/src/USER-MISC/pair_srp.cpp
+++ b/src/USER-MISC/pair_srp.cpp
@@ -425,8 +425,8 @@ void PairSRP::coeff(int narg, char **arg)
// set ij bond-bond cutoffs
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], bptype, ilo, ihi);
- force->bounds(arg[1], bptype, jlo, jhi);
+ force->bounds(FLERR,arg[0], bptype, ilo, ihi);
+ force->bounds(FLERR,arg[1], bptype, jlo, jhi);
double a0_one = force->numeric(FLERR,arg[2]);
double cut_one = cut_global;
diff --git a/src/USER-SMD/pair_smd_hertz.cpp b/src/USER-SMD/pair_smd_hertz.cpp
index 0bcf2c8a3e5612dcd1fdcf76474ca816340dafd2..3c7d7cffb9012139f7f526b7a421da35863a9fb5 100644
--- a/src/USER-SMD/pair_smd_hertz.cpp
+++ b/src/USER-SMD/pair_smd_hertz.cpp
@@ -266,8 +266,8 @@ void PairHertz::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double bulkmodulus_one = atof(arg[2]);
diff --git a/src/USER-SMD/pair_smd_triangulated_surface.cpp b/src/USER-SMD/pair_smd_triangulated_surface.cpp
index f5db6acd8e97bb5bed5c333ab3304a12d53502d8..8410f2ec05d62e8a6949f851fa27385854df7e2f 100644
--- a/src/USER-SMD/pair_smd_triangulated_surface.cpp
+++ b/src/USER-SMD/pair_smd_triangulated_surface.cpp
@@ -350,8 +350,8 @@ void PairTriSurf::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double bulkmodulus_one = atof(arg[2]);
diff --git a/src/USER-SPH/pair_sph_heatconduction.cpp b/src/USER-SPH/pair_sph_heatconduction.cpp
index 0b3716658cdb46662fe0e42aa71bccd86928aac8..02e667e437a2d2ff96c1f85600273954dff5af62 100644
--- a/src/USER-SPH/pair_sph_heatconduction.cpp
+++ b/src/USER-SPH/pair_sph_heatconduction.cpp
@@ -172,8 +172,8 @@ void PairSPHHeatConduction::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double alpha_one = force->numeric(FLERR,arg[2]);
double cut_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-SPH/pair_sph_idealgas.cpp b/src/USER-SPH/pair_sph_idealgas.cpp
index caea8eb3e4bfe489dc1893070af15524ee985728..2ba0d9bfbe6a970bcffd36e17aa412e7d9df0e4f 100644
--- a/src/USER-SPH/pair_sph_idealgas.cpp
+++ b/src/USER-SPH/pair_sph_idealgas.cpp
@@ -214,8 +214,8 @@ void PairSPHIdealGas::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double viscosity_one = force->numeric(FLERR,arg[2]);
double cut_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-SPH/pair_sph_lj.cpp b/src/USER-SPH/pair_sph_lj.cpp
index f0f3e3df339c55f3fbd2eedc4566d93c19e59725..98752ac126aa6385c3bb1bc4f2df5b91d34d73bd 100644
--- a/src/USER-SPH/pair_sph_lj.cpp
+++ b/src/USER-SPH/pair_sph_lj.cpp
@@ -222,8 +222,8 @@ void PairSPHLJ::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double viscosity_one = force->numeric(FLERR,arg[2]);
double cut_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-SPH/pair_sph_rhosum.cpp b/src/USER-SPH/pair_sph_rhosum.cpp
index 4b7b028ad272ab812a945f52529c938ecdc260ad..4bff0ec27c6e75e33e904627d06c6336581a3627 100644
--- a/src/USER-SPH/pair_sph_rhosum.cpp
+++ b/src/USER-SPH/pair_sph_rhosum.cpp
@@ -243,8 +243,8 @@ void PairSPHRhoSum::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double cut_one = force->numeric(FLERR,arg[2]);
diff --git a/src/USER-SPH/pair_sph_taitwater.cpp b/src/USER-SPH/pair_sph_taitwater.cpp
index 846ec118936fee79dfb85cc3df9582fa2583b205..b1887c6164ac1e0d9528c8fec71d864a447d5e2c 100644
--- a/src/USER-SPH/pair_sph_taitwater.cpp
+++ b/src/USER-SPH/pair_sph_taitwater.cpp
@@ -243,8 +243,8 @@ void PairSPHTaitwater::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double rho0_one = force->numeric(FLERR,arg[2]);
double soundspeed_one = force->numeric(FLERR,arg[3]);
diff --git a/src/USER-SPH/pair_sph_taitwater_morris.cpp b/src/USER-SPH/pair_sph_taitwater_morris.cpp
index d6a94f51cd2c4da83d6141f46c2f90bd6b1c62c7..76b538e46bd38fb0037a412ae2203adaa15f0d98 100644
--- a/src/USER-SPH/pair_sph_taitwater_morris.cpp
+++ b/src/USER-SPH/pair_sph_taitwater_morris.cpp
@@ -243,8 +243,8 @@ void PairSPHTaitwaterMorris::coeff(int narg, char **arg) {
allocate();
int ilo, ihi, jlo, jhi;
- force->bounds(arg[0], atom->ntypes, ilo, ihi);
- force->bounds(arg[1], atom->ntypes, jlo, jhi);
+ force->bounds(FLERR,arg[0], atom->ntypes, ilo, ihi);
+ force->bounds(FLERR,arg[1], atom->ntypes, jlo, jhi);
double rho0_one = force->numeric(FLERR,arg[2]);
double soundspeed_one = force->numeric(FLERR,arg[3]);
diff --git a/src/angle_hybrid.cpp b/src/angle_hybrid.cpp
index 42b6860b6cc9d8f9fa5c58fb58391beb424b7ba1..a477d7f8f6dcdba430122fa47bcb43d2eee45c26 100644
--- a/src/angle_hybrid.cpp
+++ b/src/angle_hybrid.cpp
@@ -245,7 +245,7 @@ void AngleHybrid::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
// 2nd arg = angle sub-style name
// allow for "none" or "skip" as valid sub-style name
diff --git a/src/angle_zero.cpp b/src/angle_zero.cpp
index 375d5ad81375fb3b36e03e9a1899097d55a36544..ba84295ca76a997a956991c35e6b0b194f925cdc 100644
--- a/src/angle_zero.cpp
+++ b/src/angle_zero.cpp
@@ -88,7 +88,7 @@ void AngleZero::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nangletypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nangletypes,ilo,ihi);
double theta0_one = 0.0;
if (coeffflag && (narg == 2))
diff --git a/src/bond_hybrid.cpp b/src/bond_hybrid.cpp
index 519a1a300825ef5f9d122c4892576cb4ce4648e5..1244c2822bb27ae64d756b3c3f66e1e510c18a6e 100644
--- a/src/bond_hybrid.cpp
+++ b/src/bond_hybrid.cpp
@@ -243,7 +243,7 @@ void BondHybrid::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
// 2nd arg = bond sub-style name
// allow for "none" as valid sub-style name
diff --git a/src/bond_zero.cpp b/src/bond_zero.cpp
index 0a354a758eb18e9bf2fe2dc97fe7152c18d31466..f38b7754f826e297e1f820e8dcb22f6f6315bc16 100644
--- a/src/bond_zero.cpp
+++ b/src/bond_zero.cpp
@@ -86,7 +86,7 @@ void BondZero::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nbondtypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nbondtypes,ilo,ihi);
double r0_one = 0.0;
if (coeffflag && (narg == 2))
diff --git a/src/comm.cpp b/src/comm.cpp
index a47807b410638076f14816359c3110ccf9ae013f..b558b3fd8c0e0d4b37e4fb3b0b87a5e332089105 100644
--- a/src/comm.cpp
+++ b/src/comm.cpp
@@ -287,7 +287,7 @@ void Comm::modify_params(int narg, char **arg)
for (i=0; i < ntypes+1; ++i)
cutusermulti[i] = -1.0;
}
- force->bounds(arg[iarg+1],ntypes,nlo,nhi,1);
+ force->bounds(FLERR,arg[iarg+1],ntypes,nlo,nhi,1);
cut = force->numeric(FLERR,arg[iarg+2]);
cutghostuser = MAX(cutghostuser,cut);
if (cut < 0.0)
diff --git a/src/compute_coord_atom.cpp b/src/compute_coord_atom.cpp
index 5747c869dc6fc565bd14ebcfff28ab652b8cdac1..21744dcc93839b4e1baf264c554b77de8dfe13b2 100644
--- a/src/compute_coord_atom.cpp
+++ b/src/compute_coord_atom.cpp
@@ -53,7 +53,7 @@ ComputeCoordAtom::ComputeCoordAtom(LAMMPS *lmp, int narg, char **arg) :
ncol = 0;
int iarg = 4;
while (iarg < narg) {
- force->bounds(arg[iarg],ntypes,typelo[ncol],typehi[ncol]);
+ force->bounds(FLERR,arg[iarg],ntypes,typelo[ncol],typehi[ncol]);
if (typelo[ncol] > typehi[ncol])
error->all(FLERR,"Illegal compute coord/atom command");
ncol++;
diff --git a/src/compute_rdf.cpp b/src/compute_rdf.cpp
index 86245080178d4f19853df130e934ec6edd0cdd43..08863940689bd44ced0cd1427f52936f3f62fae6 100644
--- a/src/compute_rdf.cpp
+++ b/src/compute_rdf.cpp
@@ -72,8 +72,8 @@ ComputeRDF::ComputeRDF(LAMMPS *lmp, int narg, char **arg) :
npairs = 0;
int iarg = 4;
while (iarg < narg) {
- force->bounds(arg[iarg],atom->ntypes,ilo[npairs],ihi[npairs]);
- force->bounds(arg[iarg+1],atom->ntypes,jlo[npairs],jhi[npairs]);
+ force->bounds(FLERR,arg[iarg],atom->ntypes,ilo[npairs],ihi[npairs]);
+ force->bounds(FLERR,arg[iarg+1],atom->ntypes,jlo[npairs],jhi[npairs]);
if (ilo[npairs] > ihi[npairs] || jlo[npairs] > jhi[npairs])
error->all(FLERR,"Illegal compute rdf command");
npairs++;
diff --git a/src/delete_bonds.cpp b/src/delete_bonds.cpp
index 3871641e0ec3ebc452f62e246007a7704d73ecb6..97eb1c09a857246f3236be089f73c1cb46008325 100644
--- a/src/delete_bonds.cpp
+++ b/src/delete_bonds.cpp
@@ -74,7 +74,7 @@ void DeleteBonds::command(int narg, char **arg)
else error->all(FLERR,"Illegal delete_bonds command");
// setup list of types (atom,bond,etc) to consider
- // use force->bounds() to allow setting of range of types
+ // use force->bounds(FLERR,) to allow setting of range of types
// range can be 0 to ntypes inclusive
int *tlist = NULL;
@@ -93,7 +93,7 @@ void DeleteBonds::command(int narg, char **arg)
tlist = new int[n+1];
for (int i = 0; i <= n; i++) tlist[i] = 0;
int nlo,nhi;
- force->bounds(arg[2],n,nlo,nhi,0);
+ force->bounds(FLERR,arg[2],n,nlo,nhi,0);
for (int i = nlo; i <= nhi; i++) tlist[i] = 1;
iarg++;
diff --git a/src/dihedral_hybrid.cpp b/src/dihedral_hybrid.cpp
index ed73351b3aef365bacb5de8246f778394dedd553..0ae396b887c9f7150582bbbcae047fb7dd235dc7 100644
--- a/src/dihedral_hybrid.cpp
+++ b/src/dihedral_hybrid.cpp
@@ -247,7 +247,7 @@ void DihedralHybrid::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
// 2nd arg = dihedral sub-style name
// allow for "none" or "skip" as valid sub-style name
diff --git a/src/dihedral_zero.cpp b/src/dihedral_zero.cpp
index 1e4a4ab3ae0b3ca53c0fada325411034e7f231d2..9b0c569fcac8743bf3124bcfe1359aeb35f07e40 100644
--- a/src/dihedral_zero.cpp
+++ b/src/dihedral_zero.cpp
@@ -84,7 +84,7 @@ void DihedralZero::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ndihedraltypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ndihedraltypes,ilo,ihi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/dump_image.cpp b/src/dump_image.cpp
index fd3dbed9e916a8093453ecce3224da24006f15f1..6d971d016aacbb732032f311c1748bfd4e5af708 100644
--- a/src/dump_image.cpp
+++ b/src/dump_image.cpp
@@ -1241,7 +1241,7 @@ int DumpImage::modify_param(int narg, char **arg)
if (strcmp(arg[0],"acolor") == 0) {
if (narg < 3) error->all(FLERR,"Illegal dump_modify command");
int nlo,nhi;
- force->bounds(arg[1],atom->ntypes,nlo,nhi);
+ force->bounds(FLERR,arg[1],atom->ntypes,nlo,nhi);
// ptrs = list of ncount colornames separated by '/'
@@ -1275,7 +1275,7 @@ int DumpImage::modify_param(int narg, char **arg)
if (strcmp(arg[0],"adiam") == 0) {
if (narg < 3) error->all(FLERR,"Illegal dump_modify command");
int nlo,nhi;
- force->bounds(arg[1],atom->ntypes,nlo,nhi);
+ force->bounds(FLERR,arg[1],atom->ntypes,nlo,nhi);
double diam = force->numeric(FLERR,arg[2]);
if (diam <= 0.0) error->all(FLERR,"Illegal dump_modify command");
for (int i = nlo; i <= nhi; i++) diamtype[i] = diam;
@@ -1304,7 +1304,7 @@ int DumpImage::modify_param(int narg, char **arg)
if (atom->nbondtypes == 0)
error->all(FLERR,"Dump modify bcolor not allowed with no bond types");
int nlo,nhi;
- force->bounds(arg[1],atom->nbondtypes,nlo,nhi);
+ force->bounds(FLERR,arg[1],atom->nbondtypes,nlo,nhi);
// ptrs = list of ncount colornames separated by '/'
@@ -1340,7 +1340,7 @@ int DumpImage::modify_param(int narg, char **arg)
if (atom->nbondtypes == 0)
error->all(FLERR,"Dump modify bdiam not allowed with no bond types");
int nlo,nhi;
- force->bounds(arg[1],atom->ntypes,nlo,nhi);
+ force->bounds(FLERR,arg[1],atom->ntypes,nlo,nhi);
double diam = force->numeric(FLERR,arg[2]);
if (diam <= 0.0) error->all(FLERR,"Illegal dump_modify command");
for (int i = nlo; i <= nhi; i++) bdiamtype[i] = diam;
diff --git a/src/fix_adapt.cpp b/src/fix_adapt.cpp
index 2a3f7db1f1affb779a6702af60c8bcabdd5b42b0..4c7eb5f218e4a7d0a1e83996ca10cbe6d74c8b07 100644
--- a/src/fix_adapt.cpp
+++ b/src/fix_adapt.cpp
@@ -92,9 +92,9 @@ nadapt(0), id_fix_diam(NULL), id_fix_chg(NULL), adapt(NULL)
adapt[nadapt].pparam = new char[n];
adapt[nadapt].pair = NULL;
strcpy(adapt[nadapt].pparam,arg[iarg+2]);
- force->bounds(arg[iarg+3],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+3],atom->ntypes,
adapt[nadapt].ilo,adapt[nadapt].ihi);
- force->bounds(arg[iarg+4],atom->ntypes,
+ force->bounds(FLERR,arg[iarg+4],atom->ntypes,
adapt[nadapt].jlo,adapt[nadapt].jhi);
if (strstr(arg[iarg+5],"v_") == arg[iarg+5]) {
n = strlen(&arg[iarg+5][2]) + 1;
diff --git a/src/force.cpp b/src/force.cpp
index dbc71b49a9b5b35959623bce9cfb6ccab3f84225..3154139d54a5f5e6116ee9f87a762e7965f04ae6 100644
--- a/src/force.cpp
+++ b/src/force.cpp
@@ -835,7 +835,8 @@ void Force::set_special(int narg, char **arg)
return nlo,nhi
------------------------------------------------------------------------- */
-void Force::bounds(char *str, int nmax, int &nlo, int &nhi, int nmin)
+void Force::bounds(const char *file, int line, char *str,
+ int nmax, int &nlo, int &nhi, int nmin)
{
char *ptr = strchr(str,'*');
@@ -856,7 +857,7 @@ void Force::bounds(char *str, int nmax, int &nlo, int &nhi, int nmin)
}
if (nlo < nmin || nhi > nmax || nlo > nhi)
- error->all(FLERR,"Numeric index is out of bounds");
+ error->all(file,line,"Numeric index is out of bounds");
}
/* ----------------------------------------------------------------------
@@ -868,8 +869,8 @@ void Force::bounds(char *str, int nmax, int &nlo, int &nhi, int nmin)
return nlo,nhi
------------------------------------------------------------------------- */
-void Force::boundsbig(char *str, bigint nmax, bigint &nlo, bigint &nhi,
- bigint nmin)
+void Force::boundsbig(const char *file, int line, char *str,
+ bigint nmax, bigint &nlo, bigint &nhi, bigint nmin)
{
char *ptr = strchr(str,'*');
@@ -890,7 +891,7 @@ void Force::boundsbig(char *str, bigint nmax, bigint &nlo, bigint &nhi,
}
if (nlo < nmin || nhi > nmax || nlo > nhi)
- error->all(FLERR,"Numeric index is out of bounds");
+ error->all(file,line,"Numeric index is out of bounds");
}
/* ----------------------------------------------------------------------
diff --git a/src/force.h b/src/force.h
index 78788826a4f851bcff0ad19ad6f55f50d87e933a..f2d9abc7dd235db067bab708669e0ae5cd633f08 100644
--- a/src/force.h
+++ b/src/force.h
@@ -125,8 +125,8 @@ class Force : protected Pointers {
void store_style(char *&, const char *, int);
void set_special(int, char **);
- void bounds(char *, int, int &, int &, int nmin=1);
- void boundsbig(char *, bigint, bigint &, bigint &, bigint nmin=1);
+ void bounds(const char *, int, char *, int, int &, int &, int nmin=1);
+ void boundsbig(const char *, int, char *, bigint, bigint &, bigint &, bigint nmin=1);
double numeric(const char *, int, char *);
int inumeric(const char *, int, char *);
bigint bnumeric(const char *, int, char *);
diff --git a/src/improper_hybrid.cpp b/src/improper_hybrid.cpp
index d819365b74103a43b8ff3445d350a9f386eb6191..78ee69569e0e5a55efc044238d0034ed6b5c6c54 100644
--- a/src/improper_hybrid.cpp
+++ b/src/improper_hybrid.cpp
@@ -256,7 +256,7 @@ void ImproperHybrid::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
// 2nd arg = improper sub-style name
// allow for "none" as valid sub-style name
diff --git a/src/improper_zero.cpp b/src/improper_zero.cpp
index 958f6621624b198e0e2150d130988bb57aef3f54..556e427ff8b9203112c8cb53f475ceab4bb6379a 100644
--- a/src/improper_zero.cpp
+++ b/src/improper_zero.cpp
@@ -84,7 +84,7 @@ void ImproperZero::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->nimpropertypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->nimpropertypes,ilo,ihi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/input.cpp b/src/input.cpp
index 905e7f0590438ad068f347b4b30cd87d5a2102b0..d1343ee38dbffaceeac83a9c39de3dba2586e53a 100644
--- a/src/input.cpp
+++ b/src/input.cpp
@@ -683,7 +683,7 @@ int Input::expand_args(int narg, char **arg, int mode, char **&earg)
if (expandflag) {
*ptr2 = '\0';
- force->bounds(ptr1+1,nmax,nlo,nhi);
+ force->bounds(FLERR,ptr1+1,nmax,nlo,nhi);
*ptr2 = ']';
if (newarg+nhi-nlo+1 > maxarg) {
maxarg += nhi-nlo+1;
@@ -1130,7 +1130,7 @@ void Input::partition()
else error->all(FLERR,"Illegal partition command");
int ilo,ihi;
- force->bounds(arg[1],universe->nworlds,ilo,ihi);
+ force->bounds(FLERR,arg[1],universe->nworlds,ilo,ihi);
// copy original line to copy, since will use strtok() on it
// ptr = start of 4th word
diff --git a/src/pair_beck.cpp b/src/pair_beck.cpp
index 0e05afb43ba758a6e35926a06e091b4726e9d0d6..e3e8b0c5cc3d679e4ca7fbdb2e96c4d2b8ca6299 100644
--- a/src/pair_beck.cpp
+++ b/src/pair_beck.cpp
@@ -199,8 +199,8 @@ void PairBeck::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double AA_one = force->numeric(FLERR,arg[2]);
double BB_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_born.cpp b/src/pair_born.cpp
index d925af1030beeb1567331de7c568a6260153e854..5fc26e2529f11552cf2085f627f4a4627dd74f44 100644
--- a/src/pair_born.cpp
+++ b/src/pair_born.cpp
@@ -200,8 +200,8 @@ void PairBorn::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_born_coul_wolf.cpp b/src/pair_born_coul_wolf.cpp
index 48066fc43e11ff2d4ad95f242d470df3a55c7d11..31c0cc715c661ef625ac7300c0436456045071e0 100644
--- a/src/pair_born_coul_wolf.cpp
+++ b/src/pair_born_coul_wolf.cpp
@@ -245,8 +245,8 @@ void PairBornCoulWolf::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_buck.cpp b/src/pair_buck.cpp
index baaea7ccdd2beb95ec8915f908ebccba3f7c7ffe..78a5321cf04877099f9c0deeffb93f67bece213e 100644
--- a/src/pair_buck.cpp
+++ b/src/pair_buck.cpp
@@ -191,8 +191,8 @@ void PairBuck::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_buck_coul_cut.cpp b/src/pair_buck_coul_cut.cpp
index 06aeea3083e04724a333cf8e2908245880f1c28d..7c948f58a8e6c323e448d94c3a0df9eec557eaa8 100644
--- a/src/pair_buck_coul_cut.cpp
+++ b/src/pair_buck_coul_cut.cpp
@@ -224,8 +224,8 @@ void PairBuckCoulCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double rho_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_coul_cut.cpp b/src/pair_coul_cut.cpp
index 4f21ba5ac1a8bbdeb0f0b0a66e458e55b9177dbe..fec592bb191e079c2e5bca23179622a00328ca82 100644
--- a/src/pair_coul_cut.cpp
+++ b/src/pair_coul_cut.cpp
@@ -171,8 +171,8 @@ void PairCoulCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_one = cut_global;
if (narg == 3) cut_one = force->numeric(FLERR,arg[2]);
diff --git a/src/pair_coul_dsf.cpp b/src/pair_coul_dsf.cpp
index a72b8d72d34075e2bc6c2ef11c9344d93e8dc432..3595f9cbee8d8b64369b5eff1c3a77ee28df8c4f 100644
--- a/src/pair_coul_dsf.cpp
+++ b/src/pair_coul_dsf.cpp
@@ -185,8 +185,8 @@ void PairCoulDSF::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/pair_coul_wolf.cpp b/src/pair_coul_wolf.cpp
index 320eed3e9fac5b95e307fb1877820332b97c25d4..1a7f4772608c98aeb8da19bb7f273e8fe64a18de 100644
--- a/src/pair_coul_wolf.cpp
+++ b/src/pair_coul_wolf.cpp
@@ -187,8 +187,8 @@ void PairCoulWolf::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int count = 0;
for (int i = ilo; i <= ihi; i++) {
diff --git a/src/pair_dpd.cpp b/src/pair_dpd.cpp
index 99166918c77e6b111bb45cb0d57672d91674358a..b5b959f85b62a3d57ef5c8de1dacc94cdaf84a66 100644
--- a/src/pair_dpd.cpp
+++ b/src/pair_dpd.cpp
@@ -223,8 +223,8 @@ void PairDPD::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a0_one = force->numeric(FLERR,arg[2]);
double gamma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_dpd_tstat.cpp b/src/pair_dpd_tstat.cpp
index 2c51bb8d850b5a1edc5addd1ad7b01f75e2011c7..6d8f75d95d61a5e7b852d4198b61a4d38404cf88 100644
--- a/src/pair_dpd_tstat.cpp
+++ b/src/pair_dpd_tstat.cpp
@@ -175,8 +175,8 @@ void PairDPDTstat::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a0_one = 0.0;
double gamma_one = force->numeric(FLERR,arg[2]);
diff --git a/src/pair_gauss.cpp b/src/pair_gauss.cpp
index 1b922a215fba430a308459632e627c80b2797537..c7b77c32709af9e5a854fc671bfa6143dc5ac5d3 100644
--- a/src/pair_gauss.cpp
+++ b/src/pair_gauss.cpp
@@ -189,8 +189,8 @@ void PairGauss::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
double b_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_hybrid.cpp b/src/pair_hybrid.cpp
index 90b847339e55063dda99d70c26b576164cef84fd..3f4392d02ab5768f04c54d1f3dcd5bdb746d589d 100644
--- a/src/pair_hybrid.cpp
+++ b/src/pair_hybrid.cpp
@@ -365,8 +365,8 @@ void PairHybrid::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
// 3rd arg = pair sub-style name
// 4th arg = pair sub-style index if name used multiple times
diff --git a/src/pair_hybrid_overlay.cpp b/src/pair_hybrid_overlay.cpp
index 30c1d7d54ee5c08d6b66c6731c5d03bb33e05c3e..6ffd6a19d6fbcd93495f6a538f1fb1ec4be4f602 100644
--- a/src/pair_hybrid_overlay.cpp
+++ b/src/pair_hybrid_overlay.cpp
@@ -37,8 +37,8 @@ void PairHybridOverlay::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
// 3rd arg = pair sub-style name
// 4th arg = pair sub-style index if name used multiple times
diff --git a/src/pair_lj96_cut.cpp b/src/pair_lj96_cut.cpp
index cba6fdbdb45780db25a9bcfc81bb3fe251de5e66..1579a9fbf4524de1ace40a2d8c2f5afce9f21fb1 100644
--- a/src/pair_lj96_cut.cpp
+++ b/src/pair_lj96_cut.cpp
@@ -458,8 +458,8 @@ void PairLJ96Cut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_cubic.cpp b/src/pair_lj_cubic.cpp
index 68b664fe210dc242bb50fd7427b9fbff12d79970..633c12019ea312d294d1e10084427800da6772a9 100644
--- a/src/pair_lj_cubic.cpp
+++ b/src/pair_lj_cubic.cpp
@@ -200,8 +200,8 @@ void PairLJCubic::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_cut.cpp b/src/pair_lj_cut.cpp
index ab1710b2fc190bfd9e8efefac22d6b637361dbd8..b8ea4a72b671d4d92b0b0dbf9b020dc2d6f81afe 100644
--- a/src/pair_lj_cut.cpp
+++ b/src/pair_lj_cut.cpp
@@ -452,8 +452,8 @@ void PairLJCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_cut_coul_cut.cpp b/src/pair_lj_cut_coul_cut.cpp
index 79d3fa12a843c6cccfc77aae12371c24d0a5c86e..0d2bff3c9f8c96941215bb16d3b3a83dfb60a118 100644
--- a/src/pair_lj_cut_coul_cut.cpp
+++ b/src/pair_lj_cut_coul_cut.cpp
@@ -217,8 +217,8 @@ void PairLJCutCoulCut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_cut_coul_dsf.cpp b/src/pair_lj_cut_coul_dsf.cpp
index d80c7e74c05a569c6d56ec1f9088ef528b6c1a6a..538336d8e648f4964dd8db46f1e20f200dd135c1 100644
--- a/src/pair_lj_cut_coul_dsf.cpp
+++ b/src/pair_lj_cut_coul_dsf.cpp
@@ -241,8 +241,8 @@ void PairLJCutCoulDSF::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_expand.cpp b/src/pair_lj_expand.cpp
index 237ee91cce83a5e284f8aaec8f67c242f151b902..90f1ae0df29abbcb81e50cc7c4f622bdc3ce20f1 100644
--- a/src/pair_lj_expand.cpp
+++ b/src/pair_lj_expand.cpp
@@ -194,8 +194,8 @@ void PairLJExpand::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_gromacs.cpp b/src/pair_lj_gromacs.cpp
index d5b24b733ec83d6d139b6986c4b7e973144d05ef..bb0a6e647e8f76cbe892c6f61af042c97596825b 100644
--- a/src/pair_lj_gromacs.cpp
+++ b/src/pair_lj_gromacs.cpp
@@ -223,8 +223,8 @@ void PairLJGromacs::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_gromacs_coul_gromacs.cpp b/src/pair_lj_gromacs_coul_gromacs.cpp
index b68a557dcb2750883e297a87e5f0eb615cbfaea7..26e433b38113a6ef282c7034f1ddcdb9f1b5539e 100644
--- a/src/pair_lj_gromacs_coul_gromacs.cpp
+++ b/src/pair_lj_gromacs_coul_gromacs.cpp
@@ -247,8 +247,8 @@ void PairLJGromacsCoulGromacs::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_smooth.cpp b/src/pair_lj_smooth.cpp
index b73545f57b509df0d172d20b1c1b86062e59def5..c59b35aebf87287ce0c6b6afdb05cafac9f4f4a0 100644
--- a/src/pair_lj_smooth.cpp
+++ b/src/pair_lj_smooth.cpp
@@ -225,8 +225,8 @@ void PairLJSmooth::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_lj_smooth_linear.cpp b/src/pair_lj_smooth_linear.cpp
index fa41bcae4b5b7a8c5693b0de7a9a02daeb1ea861..189475aa71bd62e2598890cbac33db48c39a3496 100644
--- a/src/pair_lj_smooth_linear.cpp
+++ b/src/pair_lj_smooth_linear.cpp
@@ -192,8 +192,8 @@ void PairLJSmoothLinear::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_mie_cut.cpp b/src/pair_mie_cut.cpp
index 7cb67b9437b3d11618d826945517936119320b21..b79c6c9bf92c98b662658a6dbf64f0df7c1ef4d1 100644
--- a/src/pair_mie_cut.cpp
+++ b/src/pair_mie_cut.cpp
@@ -463,8 +463,8 @@ void PairMIECut::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double epsilon_one = force->numeric(FLERR,arg[2]);
double sigma_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_morse.cpp b/src/pair_morse.cpp
index 4e2d47be5c4ab92fc68ce9a729ac4753af4290d2..503a882595a22d9bd3843f9aa2af5849730c4029 100644
--- a/src/pair_morse.cpp
+++ b/src/pair_morse.cpp
@@ -180,8 +180,8 @@ void PairMorse::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double d0_one = force->numeric(FLERR,arg[2]);
double alpha_one = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_soft.cpp b/src/pair_soft.cpp
index d9b9bb819ba55b922322cd03b3d032a89d806839..8ffd1393077f82e0e9db88fe5f7a65519b536961 100644
--- a/src/pair_soft.cpp
+++ b/src/pair_soft.cpp
@@ -174,8 +174,8 @@ void PairSoft::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double prefactor_one = force->numeric(FLERR,arg[2]);
diff --git a/src/pair_table.cpp b/src/pair_table.cpp
index f83e499a23cbe8ce8abe487040b9083963f9ef6a..c4bc3e7dd2096a379874618317187399e66cb1e7 100644
--- a/src/pair_table.cpp
+++ b/src/pair_table.cpp
@@ -263,8 +263,8 @@ void PairTable::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
int me;
MPI_Comm_rank(world,&me);
diff --git a/src/pair_yukawa.cpp b/src/pair_yukawa.cpp
index 2cc309d949f2bb4fbe0731d509f17f0336f34a5c..a38e2aa8808aa1a9f1ada33654d2328783c2d7f7 100644
--- a/src/pair_yukawa.cpp
+++ b/src/pair_yukawa.cpp
@@ -178,8 +178,8 @@ void PairYukawa::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double a_one = force->numeric(FLERR,arg[2]);
diff --git a/src/pair_zbl.cpp b/src/pair_zbl.cpp
index 69622009826f1cfabcdd7cf3bc355644500b7eab..86c6c64d77feef8fc0b013925c6c54e82e1e6556 100644
--- a/src/pair_zbl.cpp
+++ b/src/pair_zbl.cpp
@@ -211,10 +211,10 @@ void PairZBL::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
int jlo,jhi;
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
z_one = force->numeric(FLERR,arg[2]);
z_two = force->numeric(FLERR,arg[3]);
diff --git a/src/pair_zero.cpp b/src/pair_zero.cpp
index 787ed29935b7a1e162ac02da903fe5e52a3581ad..ab6d713cc4aa433524eff1484a38b7dd4eba5d31 100644
--- a/src/pair_zero.cpp
+++ b/src/pair_zero.cpp
@@ -111,8 +111,8 @@ void PairZero::coeff(int narg, char **arg)
if (!allocated) allocate();
int ilo,ihi,jlo,jhi;
- force->bounds(arg[0],atom->ntypes,ilo,ihi);
- force->bounds(arg[1],atom->ntypes,jlo,jhi);
+ force->bounds(FLERR,arg[0],atom->ntypes,ilo,ihi);
+ force->bounds(FLERR,arg[1],atom->ntypes,jlo,jhi);
double cut_one = cut_global;
if (coeffflag && (narg == 3)) cut_one = force->numeric(FLERR,arg[2]);
diff --git a/src/set.cpp b/src/set.cpp
index 56fbf42e543707528c8925b46f9ec33e16e8063b..82d0bce129b229e73c189d6183f3b44599eb4de6 100644
--- a/src/set.cpp
+++ b/src/set.cpp
@@ -498,7 +498,7 @@ void Set::selection(int n)
if (atom->tag_enable == 0)
error->all(FLERR,"Cannot use set atom with no atom IDs defined");
bigint nlobig,nhibig;
- force->boundsbig(id,MAXTAGINT,nlobig,nhibig);
+ force->boundsbig(FLERR,id,MAXTAGINT,nlobig,nhibig);
tagint *tag = atom->tag;
for (int i = 0; i < n; i++)
@@ -509,7 +509,7 @@ void Set::selection(int n)
if (atom->molecule_flag == 0)
error->all(FLERR,"Cannot use set mol with no molecule IDs defined");
bigint nlobig,nhibig;
- force->boundsbig(id,MAXTAGINT,nlobig,nhibig);
+ force->boundsbig(FLERR,id,MAXTAGINT,nlobig,nhibig);
tagint *molecule = atom->molecule;
for (int i = 0; i < n; i++)
@@ -517,7 +517,7 @@ void Set::selection(int n)
else select[i] = 0;
} else if (style == TYPE_SELECT) {
- force->bounds(id,atom->ntypes,nlo,nhi);
+ force->bounds(FLERR,id,atom->ntypes,nlo,nhi);
int *type = atom->type;
for (int i = 0; i < n; i++)