home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 11 Util / 11-Util.zip / MAWK113.ZIP / mawk113 / matherr.c < prev    next >
C/C++ Source or Header  |  1992-03-31  |  3KB  |  189 lines
  1.  
  2. /********************************************
  3. matherr.c
  4. copyright 1991, Michael D. Brennan
  5.  
  6. This is a source file for mawk, an implementation of
  7. the AWK programming language.
  8.  
  9. Mawk is distributed without warranty under the terms of
  10. the GNU General Public License, version 2, 1991.
  11. ********************************************/
  12.  
  13. /*$Log: matherr.c,v $
  14.  * Revision 5.2  1992/03/31  16:14:44  brennan
  15.  * patch2:
  16.  * TURN_ON_FPE_TRAPS() macro
  17.  * USE_IEEEFP_H macro
  18.  *
  19.  * Revision 5.1  91/12/05  07:56:18  brennan
  20.  * 1.1 pre-release
  21.  * 
  22. */
  23.  
  24. #include  "mawk.h"
  25. #include  <math.h>
  26.  
  27. #ifdef  USE_IEEEFP_H
  28. #include <ieeefp.h>
  29. #endif
  30.  
  31. #if   FPE_TRAPS_ON
  32. #include <signal.h>
  33.  
  34. /* machine dependent changes might be needed here */
  35.  
  36. static void  fpe_catch( signal, why)
  37.   int signal, why ;
  38. {
  39.  
  40. #if   NOINFO_SIGFPE
  41.  /* some systems give no hook to find out what the exception
  42.     was -- stuff like this is why people still use fortran 
  43.  
  44.     If this fits, #define NOINFO_SIGFPE 1 in  your config.h
  45. */
  46.   rt_error("floating point exception, probably overflow") ;
  47. #else
  48.  
  49.   switch(why)
  50.   {
  51.     case FPE_ZERODIVIDE :
  52.        rt_error("division by zero") ;
  53.  
  54.     case FPE_OVERFLOW  :
  55.        rt_error("floating point overflow") ;
  56.  
  57.     default :
  58.       rt_error("floating point exception") ;
  59.   }
  60. #endif  
  61. }
  62.  
  63. void   fpe_init()
  64. { 
  65.   TURN_ON_FPE_TRAPS() ;
  66.   (void) signal(SIGFPE, fpe_catch) ; 
  67. }
  68.  
  69. #else  /* FPE_TRAPS_ON==0 */
  70.  
  71. void  fpe_init()
  72. {
  73.   TURN_OFF_FPE_TRAPS() ;
  74. }
  75. #endif
  76.  
  77. #if  HAVE_MATHERR
  78.  
  79. #if  ! FPE_TRAPS_ON
  80.  
  81. /* If we are not trapping math errors, we will shutup the library calls
  82. */
  83.  
  84. int  matherr( e )
  85.   struct exception *e ;
  86. { return 1 ; } 
  87.  
  88. #else   /* print error message and exit */
  89.  
  90. int matherr( e )
  91.   struct exception  *e ;
  92. { char *error ;
  93.  
  94.   switch( e->type )
  95.   {
  96.     case  DOMAIN :
  97.     case  SING :
  98.             error = "domain error" ;
  99.             break ;
  100.  
  101.     case  OVERFLOW :
  102.             error = "overflow" ;
  103.             break ;
  104.  
  105.     case  TLOSS :
  106.     case  PLOSS :
  107.             error = "loss of significance" ;
  108.             break ;
  109.  
  110.     case  UNDERFLOW :
  111.             e->retval = 0.0 ;
  112.             return  1 ;  /* ignore it */
  113.   }
  114.  
  115.   if ( strcmp(e->name, "atan2") == 0 )
  116.       rt_error("atan2(%g,%g) : %s" ,
  117.          e->arg1, e->arg2, error ) ;
  118.   else
  119.       rt_error("%s(%g) : %s" , e->name, e->arg1, error) ;
  120.  
  121.   /* won't get here */
  122.   return 0 ;
  123. }
  124. #endif   /* FPE_TRAPS */
  125.  
  126. #endif   /*  HAVE_MATHERR */
  127.  
  128.  
  129. /* this is how one gets the libm calls to do the right
  130. thing on bsd43_vax
  131. */
  132.  
  133. #ifdef   BSD43_VAX
  134.  
  135. #include <errno.h>
  136.  
  137. double infnan( arg )
  138.   int arg ;
  139. {
  140.   switch(arg)
  141.   {
  142.     case  ERANGE : errno = ERANGE ; return HUGE ;
  143.     case -ERANGE : errno = EDOM ; return -HUGE ;
  144.     default :  errno = EDOM ; 
  145.   }
  146.   return 0.0 ;
  147. }
  148.  
  149. #endif  /* BSD43_VAX */
  150.  
  151. /* This routine is for XENIX-68K 2.3A.
  152.     Error check routine to be called after fp arithmetic.
  153. */
  154.  
  155. #if SW_FP_CHECK
  156. /* Definitions of bit values in iserr() return value */
  157.  
  158. #define OVFLOW        2
  159. #define UFLOW        4
  160. #define ZERODIV        8
  161. #define OVFLFIX        32
  162. #define INFNAN        64
  163.  
  164. void
  165. fpcheck()
  166. {
  167.     register int fperrval ;
  168.     char *errdesc ;
  169.  
  170.     if ((fperrval = iserr()) == 0)
  171.         return ;    /* no error */
  172.  
  173.     errdesc = (char *) 0 ;
  174.  
  175.     if (fperrval & INFNAN)
  176.         errdesc = "arg is infinity or NAN" ;
  177.     else if (fperrval & ZERODIV)
  178.         errdesc = "division by zero" ;
  179.     else if (fperrval & OVFLOW)
  180.         errdesc = "overflow" ;
  181.     else if (fperrval & UFLOW)
  182.         ;        /* ignored */
  183.  
  184.     if (errdesc)
  185.         rt_error("%s", errdesc) ;
  186. }
  187.  
  188. #endif
  189.