[BACK]Return to s_erf.c CVS log [TXT][DIR] Up to [cvs.NetBSD.org] / src / lib / libm / src

Please note that diffs are not public domain; they are subject to the copyright notices on the relevant files.

Diff for /src/lib/libm/src/s_erf.c between version 1.5 and 1.6

version 1.5, 1994/08/10 20:32:11 version 1.6, 1994/08/18 23:06:36
Line 198  sb7  = -2.24409524465858183362e+01; /* 0
Line 198  sb7  = -2.24409524465858183362e+01; /* 0
         double x;          double x;
 #endif  #endif
 {  {
         int hx,ix,i;          int32_t hx,ix,i;
         double R,S,P,Q,s,y,z,r;          double R,S,P,Q,s,y,z,r;
         GET_HIGH_WORD(hx,x);          GET_HIGH_WORD(hx,x);
         ix = hx&0x7fffffff;          ix = hx&0x7fffffff;
         if(ix>=0x7ff00000) {            /* erf(nan)=nan */          if(ix>=0x7ff00000) {            /* erf(nan)=nan */
             i = ((unsigned)hx>>31)<<1;              i = ((u_int32_t)hx>>31)<<1;
             return (double)(1-i)+one/x; /* erf(+-inf)=+-1 */              return (double)(1-i)+one/x; /* erf(+-inf)=+-1 */
         }          }
   
Line 254  sb7  = -2.24409524465858183362e+01; /* 0
Line 254  sb7  = -2.24409524465858183362e+01; /* 0
         double x;          double x;
 #endif  #endif
 {  {
         int hx,ix;          int32_t hx,ix;
         double R,S,P,Q,s,y,z,r;          double R,S,P,Q,s,y,z,r;
         GET_HIGH_WORD(hx,x);          GET_HIGH_WORD(hx,x);
         ix = hx&0x7fffffff;          ix = hx&0x7fffffff;
         if(ix>=0x7ff00000) {                    /* erfc(nan)=nan */          if(ix>=0x7ff00000) {                    /* erfc(nan)=nan */
                                                 /* erfc(+-inf)=0,2 */                                                  /* erfc(+-inf)=0,2 */
             return (double)(((unsigned)hx>>31)<<1)+one/x;              return (double)(((u_int32_t)hx>>31)<<1)+one/x;
         }          }
   
         if(ix < 0x3feb0000) {           /* |x|<0.84375 */          if(ix < 0x3feb0000) {           /* |x|<0.84375 */

Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6

CVSweb <webmaster@jp.NetBSD.org>