d/freebsd-skq
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

* Update the comments that refer to erf[c](nan) and erf[c](+-inf)

Browse Source 
to use erf[c]f.
This commit is contained in:
Steve Kargl 2014-07-13 16:05:33 +00:00
parent 019ffb5df8
commit c9cb48ff4d
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/msun/src/s_erff.c
View File

@ -89,9 +89,9 @@ erff(float x)
float R,S,P,Q,s,y,z,r; float R,S,P,Q,s,y,z,r;
GET_FLOAT_WORD(hx,x); GET_FLOAT_WORD(hx,x);
ix = hx&0x7fffffff; ix = hx&0x7fffffff;
if(ix>=0x7f800000) { /* erf(nan)=nan */ if(ix>=0x7f800000) { /* erff(nan)=nan */
i = ((u_int32_t)hx>>31)<<1; i = ((u_int32_t)hx>>31)<<1;
return (float)(1-i)+one/x; /* erf(+-inf)=+-1 */ return (float)(1-i)+one/x; /* erff(+-inf)=+-1 */
} }
if(ix < 0x3f580000) { /* |x|<0.84375 */ if(ix < 0x3f580000) { /* |x|<0.84375 */
@ -136,8 +136,8 @@ erfcf(float x)
float R,S,P,Q,s,y,z,r; float R,S,P,Q,s,y,z,r;
GET_FLOAT_WORD(hx,x); GET_FLOAT_WORD(hx,x);
ix = hx&0x7fffffff; ix = hx&0x7fffffff;
if(ix>=0x7f800000) { /* erfc(nan)=nan */ if(ix>=0x7f800000) { /* erfcf(nan)=nan */
/* erfc(+-inf)=0,2 */ /* erfcf(+-inf)=0,2 */
return (float)(((u_int32_t)hx>>31)<<1)+one/x; return (float)(((u_int32_t)hx>>31)<<1)+one/x;
} }