home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: Science / Science.zip / gmt_os2.zip / src / math / w_log10.c < prev    next >
Encoding:
C/C++ Source or Header  |  1995-03-16  |  923 b   |  43 lines
  1.  
  2. /* @(#)w_log10.c 1.3 95/01/18 */
  3. /*
  4.  * ====================================================
  5.  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
  6.  *
  7.  * Developed at SunSoft, a Sun Microsystems, Inc. business.
  8.  * Permission to use, copy, modify, and distribute this
  9.  * software is freely granted, provided that this notice 
  10.  * is preserved.
  11.  * ====================================================
  12.  */
  13.  
  14. /* 
  15.  * wrapper log10(X)
  16.  */
  17.  
  18. #include "fdlibm.h"
  19.  
  20.  
  21. #ifdef __STDC__
  22.     double log10(double x)        /* wrapper log10 */
  23. #else
  24.     double log10(x)            /* wrapper log10 */
  25.     double x;
  26. #endif
  27. {
  28. #ifdef _IEEE_LIBM
  29.     return __ieee754_log10(x);
  30. #else
  31.     double z;
  32.     z = __ieee754_log10(x);
  33.     if(_LIB_VERSION == _IEEE_ || isnan(x)) return z;
  34.     if(x<=0.0) {
  35.         if(x==0.0)
  36.             return __kernel_standard(x,x,18); /* log10(0) */
  37.         else 
  38.             return __kernel_standard(x,x,19); /* log10(x<0) */
  39.     } else
  40.         return z;
  41. #endif
  42. }
  43.