/* Some constant values --  */

#include "mconf.h"
/* Many of these values should ideally come from <float.h> or
   <values.h>, which should be included in mconf.h if required */

#ifdef DBL_EPSILON
double MACHEP =  DBL_EPSILON;
#else
double MACHEP =  1.11022302462515654042E-16;   /* 2**-53 */
#endif

#if defined DBL_MIN
double UFLOWTHRESH = DBL_MIN;
#elif defined MINDOUBLE
double UFLOWTHRESH = MINDOUBLE;
#else
double UFLOWTHRESH =  2.22507385850720138309E-308; /* 2**-1022 */
#endif

#ifdef DBL_MAX_10_EXP
double MAXLOG = DBL_MAX_10_EXP;
#else
double MAXLOG =  7.08396418532264106224E2;     /* log 2**1022 */
#endif

#ifdef DBL_MIN_10_EXP
double MINLOG = DBL_MIN_10_EXP;
#else
double MINLOG = -7.08396418532264106224E2;     /* log 2**-1022 */
#endif

#if defined MAXDOUBLE
double MAXNUM =  MAXDOUBLE;
#elif defined DBL_MAX
double MAXNUM =  DBL_MAX;
#else
double MAXNUM =  1.79769313486231570815E308;    /* 2**1024*(1-MACHEP) */
#endif

#ifdef M_PI
#ifndef PI
double PI     =  M_PI;
#endif
double PIO2   =  M_PI/2;
double PIO4   =  M_PI/4;
double THPIO4 =  0.75*M_PI;
double TWOOPI =  2/M_PI;
#else
#ifndef PI
double PI     =  3.14159265358979323846;       /* pi */
#endif
double PIO2   =  1.57079632679489661923;       /* pi/2 */
double PIO4   =  7.85398163397448309616E-1;    /* pi/4 */
double THPIO4 =  2.35619449019234492885;       /* 3*pi/4 */
double TWOOPI =  6.36619772367581343075535E-1; /* 2/pi */
#endif

#ifdef M_SQRT2
double SQRT2  =  M_SQRT2;       /* sqrt(2) */
double SQRTH  =  M_SQRT2/2;    /* sqrt(2)/2 */
#else
double SQRT2  =  1.41421356237309504880;       /* sqrt(2) */
double SQRTH  =  7.07106781186547524401E-1;    /* sqrt(2)/2 */
#endif

double SQ2OPI =  7.9788456080286535587989E-1;  /* sqrt( 2/pi ) */

#ifdef M_LN2
double LOGE2  =  M_LN2;
double LOGSQ2 =  M_LN2/2;    /* log(2)/2 */
double LOG2E  =  1/M_LN2;     /* 1/log(2) */
#else
double LOGE2  =  6.93147180559945309417E-1;    /* log(2) */
double LOGSQ2 =  3.46573590279972654709E-1;    /* log(2)/2 */
double LOG2E  =  1.4426950408889634073599;     /* 1/log(2) */
#endif