#define DBL_MEMCPY(dest,src,n) memcpy((dest),(src),(n)*sizeof(double))
#define DBL_ZERO_MEMSET(dest,n) memset((dest),0,(n)*sizeof(double))

/* Error estimation safety coefficient for methods that use step
 * doubling for error estimates. Error estimates are multiplied by
 * this constant to ensure that the error of a step is not
 * underestimated. 
 *
 * The default safety value of 8.0 ensures 90% of samples lie within
 * the error (assuming a Gaussian distribution with prior 
 * p(sigma) = 1 / sigma). Value of 1.0 conforms to equation 
 * by Ascher and Petzold (reference: Ascher, U.M., Petzold, L.R.,
 * Computer methods for ordinary differential and
 * differential-algebraic equations, SIAM, Philadelphia, 1998).
 */

#define ODEIV_ERR_SAFETY 8.0