1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
#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
|
