Fixed bug in a call to memset

splinterpolator.h

@@ -644,7 +644,7 @@ double Splinterpolator<T>::value_and_derivatives_at(const double       *coord,
 						    double             *dval) 
 const
 {
-  if (should_be_zero(coord)) { memset(dval,n_nonzero(deriv)*sizeof(double),0); return(0.0); }
+  if (should_be_zero(coord)) { memset(dval,0,n_nonzero(deriv)*sizeof(double)); return(0.0); }
 
   double       iwgt[8], jwgt[8], kwgt[8], lwgt[8], mwgt[8];
   double       *wgts[] = {iwgt, jwgt, kwgt, lwgt, mwgt};