diff --git a/src/USER-OMP/pair_reaxc_omp.cpp b/src/USER-OMP/pair_reaxc_omp.cpp index f2413de430e57d0a99f2ae3d73c8cbd7692c63a0..d4448ad868c3ec84763d56d8a753fb5adf2353a6 100644 --- a/src/USER-OMP/pair_reaxc_omp.cpp +++ b/src/USER-OMP/pair_reaxc_omp.cpp @@ -352,6 +352,9 @@ void PairReaxCOMP::init_style( ) neighbor->requests[irequest]->ghost = 1; cutmax = MAX3(control->nonb_cut, control->hbond_cut, control->bond_cut); + if ((cutmax < 2.0*control->bond_cut) && (comm->me == 0)) + error->warning(FLERR,"Total cutoff < 2*bond cutoff. May need to use an " + "increased neighbor list skin."); for( int i = 0; i < LIST_N; ++i ) lists[i].allocated = 0;