home *** CD-ROM | disk | FTP | other *** search
/ Scene Storm / Scene Storm - Volume 1.iso / coding / c / pdc / libsrc / math / log.c < prev    next >
Encoding:
C/C++ Source or Header  |  1990-04-07  |  1.1 KB  |  49 lines
  1. #include <stdio.h>
  2. #include <math.h>
  3. #include "pml.h"
  4. #include "cordic.h"
  5.  
  6. static char     funcname[] = "log";
  7.  
  8. extern CORDIC_Table CORDIC_Table1[1];
  9.  
  10. double
  11. log(x)
  12.     double          x;
  13. {
  14.     double          y, z;
  15.     int             k;
  16.     double          x0, y0, z0;
  17.     struct exception xcpt;
  18.  
  19.     if (x == 0.0) {
  20.         xcpt.type = SING;
  21.         xcpt.name = funcname;
  22.         xcpt.arg1 = x;
  23.         if (!matherr(&xcpt)) {
  24.             fprintf(stderr, "%s: SINGULARITY error\n", funcname);
  25.             errno = EDOM;
  26.             xcpt.retval = -MAXDOUBLE;
  27.         }
  28.     }
  29.     else if (x < 0.0) {
  30.         xcpt.type = DOMAIN;
  31.         xcpt.name = funcname;
  32.         xcpt.arg1 = x;
  33.         if (!matherr(&xcpt)) {
  34.             fprintf(stderr, "%s: DOMAIN error\n", funcname);
  35.             errno = EDOM;
  36.             xcpt.retval = -MAXDOUBLE;
  37.         }
  38.     }
  39.     else {
  40.         y = frexp(x, &k);
  41.         x0 = y + 1.0;
  42.         y0 = y - 1.0;
  43.         z0 = 0.0;
  44.         CORDIC_rotate1(0, CORDIC_Table1, &x0, &y0, &z0);
  45.         xcpt.retval = k * LN2 + 2.0 * z0;
  46.     }
  47.     return (xcpt.retval);
  48. }
  49.