store squared cutoff as documented

src/MISC/pair_list.cpp

@@ -22,6 +22,7 @@
 #include "comm.h"
 #include "error.h"
 #include "force.h"
+#include "math_special.h"
 #include "memory.h"
 #include "text_file_reader.h"
 
@@ -31,6 +32,7 @@
 #include <map>
 
 using namespace LAMMPS_NS;
+using MathSpecial::square;
 
 enum { NONE = 0, HARM, MORSE, LJ126, QUARTIC };
 
@@ -286,7 +288,7 @@ void PairList::settings(int narg, char **arg)
           break;
         }
         if (values.has_next())
-          oneparam.cutsq = mypow(values.next_double(), 2);
+          oneparam.cutsq = square(values.next_double());
         else
           oneparam.cutsq = cut_global*cut_global;