/************************************************************************
 *                                                                      *
 *                              N O T I C E                             *
 *                                                                      *
 *                      Copyright Abandoned, 1987, Fred Fish            *
 *                                                                      *
 *      This previously copyrighted work has been placed into the       *
 *      public domain by the author (Fred Fish) and may be freely used  *
 *      for any purpose, private or commercial.  I would appreciate     *
 *      it, as a courtesy, if this notice is left in all copies and     *
 *      derivative works.  Thank you, and enjoy...                      *
 *                                                                      *
 *      The author makes no warranty of any kind with respect to this   *
 *      product and explicitly disclaims any implied warranties of      *
 *      merchantability or fitness for any particular purpose.          *
 *                                                                      *
 ************************************************************************
 */

/*
 *      This file gets included with all of the floating point math
 *      library routines when they are compiled.  Note that
 *      this is the proper place to put machine dependencies
 *      whenever possible.
 *
 *      It should be pointed out that for simplicity's sake, the
 *      environment parameters are defined as floating point constants,
 *      rather than octal or hexadecimal initializations of allocated
 *      storage areas.  This means that the range of allowed numbers
 *      may not exactly match the hardware's capabilities.  For example,
 *      if the maximum positive double precision floating point number
 *      is EXACTLY 1.11...E100 and the constant "MAXDOUBLE is
 *      defined to be 1.11E100 then the numbers between 1.11E100 and
 *      1.11...E100 are considered to be undefined.  For most
 *      applications, this will cause no problems.
 *
 *      An alternate method is to allocate a global static "double" variable,
 *      say "maxdouble", and use a union declaration and initialization
 *      to initialize it with the proper bits for the EXACT maximum value.
 *      This was not done because the only compilers available to the
 *      author did not fully support union initialization features.
 *
 */

#ifndef NO_DBUG
#    include <dbug.h>
#else
#    define DBUG_ENTER(a1)
#    define DBUG_RETURN(a1) return(a1)
#    define DBUG_VOID_RETURN return
#    define DBUG_EXECUTE(keyword,a1)
#    define DBUG_2(keyword,format)
#    define DBUG_3(keyword,format,a1)
#    define DBUG_4(keyword,format,a1,a2)
#    define DBUG_5(keyword,format,a1,a2,a3)
#    define DBUG_PUSH(a1)
#    define DBUG_POP()
#    define DBUG_PROCESS(a1)
#    define DBUG_FILE (stderr)
#endif

#include <errno.h>

/*
 *      MAXDOUBLE       =>      Maximum double precision number
 *      MINDOUBLE       =>      Minimum double precision number
 *      DMAXEXP         =>      Maximum exponent of a double
 *      DMINEXP         =>      Minimum exponent of a double
 */
 
#ifndef MINDOUBLE
#define MINDOUBLE       (1.0/MAXDOUBLE) 
#endif 

#define LOG2_MAXDOUBLE  (DMAXEXP + 1)
#define LOG2_MINDOUBLE  (DMINEXP - 1)
#define LOGE_MAXDOUBLE  (LOG2_MAXDOUBLE / LOG2E)
#define LOGE_MINDOUBLE  (LOG2_MINDOUBLE / LOG2E)

/*
 *      The following are hacks which should be fixed when I understand all
 *      the issues a little better.   |tanh(TANH_MAXARG)| = 1.0
 */
#define TANH_MAXARG 16
#define SQRT_MAXDOUBLE 1.304380e19

#define PI              3.14159265358979323846
#define TWOPI           (2.0 * PI)
#define HALFPI          (PI / 2.0)
#define FOURTHPI        (PI / 4.0)
#define SIXTHPI         (PI / 6.0)
#define LOG2E           1.4426950408889634074   /* Log to base 2 of e */
#define LOG10E          0.4342944819032518276
#define SQRT2           1.4142135623730950488
#define SQRT3           1.7320508075688772935
#define LN2             0.6931471805599453094
#define LNSQRT2         0.3465735902799726547


/*
 *      MC68000 HARDWARE DEPENDENCIES
 *
 *              cc -DIEEE       =>      uses IEEE floating point format
 *
 */

#define TRUE 1                  /* This really should be in stdio.h */
#define FALSE 0                 /* This too */
#define VOID void

struct exception {
        int type;
        char *name;
        double arg1;
        double arg2;
        double retval;
};

 
#define DOMAIN          1
#define SING            2
#define OVERFLOW        3
#define UNDERFLOW       4
#define TLOSS           5
#define PLOSS           6

#define MAX_ERRORS              8

#define LOG                             1
#define COUNT                   2
#define CONTINUE                4