From 89c0655809998f498a172c4018f46337cb7d3087 Mon Sep 17 00:00:00 2001
From: ckadding <cody.addington@gmail.com>
Date: Tue, 2 Oct 2018 09:46:53 -0400
Subject: [PATCH] Add input validation

---
 src/USER-MISC/compute_pressure_cylinder.cpp | 12 +++++++++++-
 1 file changed, 11 insertions(+), 1 deletion(-)

diff --git a/src/USER-MISC/compute_pressure_cylinder.cpp b/src/USER-MISC/compute_pressure_cylinder.cpp
index 4c3b227842..a40e105b5a 100644
--- a/src/USER-MISC/compute_pressure_cylinder.cpp
+++ b/src/USER-MISC/compute_pressure_cylinder.cpp
@@ -26,6 +26,7 @@
 #include "memory.h"
 #include "error.h"
 #include "citeme.h"
+#include "domain.h"
 
 using namespace LAMMPS_NS;
 
@@ -56,10 +57,19 @@ ComputePressureCyl::ComputePressureCyl(LAMMPS *lmp, int narg, char **arg) :
   zhi=force->numeric(FLERR,arg[4]);
   Rmax=force->numeric(FLERR,arg[5]);
   bin_width=force->numeric(FLERR,arg[6]);
+     
+  if (bin_width<0.0 || bin_width<Rmax) 
+    error->all(FLERR,"Illegal compute pressure/cylinder command");
+  if (zhi<zlo || (zhi-zlo)<bin_width)
+    error->all(FLERR,"Illegal compute pressure/cylinder command");
+  if (zhi>domain->boxhi[2] || zlo<domain->boxlo[2])
+    error->all(FLERR,"Illegal compute pressure/cylinder command");
 
   nbins=(int)(Rmax/bin_width);
-
   nzbins=(int)((zhi-zlo)/bin_width);
+     
+  if (nbins<1 or nzbins<1) 
+    error->all(FLERR,"Illegal compute pressure/cylinder command");
 
   array_flag=1;
   vector_flag=0;
-- 
GitLab