Commit 3e2b572e authored by Richard Berger's avatar Richard Berger
Browse files

Fix special case for outside check in fix pour 

If the range between lo and hi is bigger than the extent in that dimension, in
the periodic case the value will always be inside.
parent 4a0495ae

src/GRANULAR/fix_pour.cpp

+1 −1
Original line number Diff line number Diff line
@@ -790,7 +790,7 @@ bool FixPour::outside(int dim, double value, double lo, double hi) 
  bool outside_regular_range = (value < lo || value > hi);



  if (domain->periodicity[dim]) {

    if (lo < boxlo && hi > boxhi) {

    if ((lo < boxlo && hi > boxhi) || (hi - lo) > domain->prd[dim]) {

      // value is always inside

      outside_pbc_range = false;

    } else if (lo < boxlo) {

CRA Git | Maintained and supported by SUSTech CRA and CCSE