/**
*
* Redefine secondary simulation function names to become primary names
* for systems without a Numeric Data Processor.
*
*/
#define _acos acos
#define _asin asin
#define _atan atan
#define _cos cos
#define _cosh cosh
#define _cot cot
#define _exp exp
#define _fabs fabs
#define _ldexp ldexp
#define _log log
#define _log10 log10
#define _modf modf
#define _pow pow
#define _pow2 pow2
#define _sin sin
#define _sinh sinh
#define _sqrt sqrt
#define _tan tan
#define _tanh tanh

/**
*
* Structure to hold information about math exceptions
*
*/
struct exception
   {
   int type;            /* error type */
   char *name;          /* math function name */
   double arg1, arg2;   /* function arguments */
   double retval;       /* proposed return value */
   };

/*
*
* Exception type codes, found in exception.type
*
*/
#define DOMAIN    1     /* domain error */
#define SING      2     /* singularity */
#define OVERFLOW  3     /* overflow */
#define UNDERFLOW 4     /* underflow */
#define TLOSS     5     /* total loss of significance */
#define PLOSS     6     /* partial loss of significance */

/**
*
* Error codes generated by basic arithmetic operations (+ - * /)
*
*/
#define FPEUND 1        /* underflow */
#define FPEOVF 2        /* overflow */
#define FPEZDV 3        /* zero divisor */
#define FPENAN 4        /* not a number (invalid operation) */

/**
*
* Constants
*
*/
#define PI   3.14159265358979323846
#define PID2 1.57079632679489661923   /* PI divided by 2 */
#define PID4 0.78539816339744830962   /* PI divided by 4 */
#define I_PI 0.31830988618379067154   /* Inverse of PI */
#define I_PID2 0.63661977236758134308 /* Inverse of PID2 */

#define HUGE 1.797693e308             /* huge value */
#define TINY 2.2e-308                 /* tiny value */
#define LOGHUGE 709.778               /* natural log of huge value */
#define LOGTINY -708.396              /* natural log of tiny value */

/**
*
* External declarations
*
*/
extern int _fperr;      /* floating point arithmetic error */
extern int errno;       /* UNIX error code */

extern char *ecvt();
extern short *seed48();
extern int atoi(),matherr();
extern long atol(),strtol(),lrand48(),nrand48(),mrand48(),jrand48();
extern double atof(),exp(),log(),log10(),pow(),sqrt();
extern double floor(),ceil(),fmod(),fabs(),frexp(),ldexp(),modf();
extern double sinh(),cosh(),tanh(),sin(),cos(),tan(),cot(),asin(),acos();
extern double atan(),atan2(),except();
extern double drand48(),erand48();