Changeset 197d59d in mainline

Timestamp:

2021-06-18T16:06:03Z (3 years ago)

Author:

Maurizio Lombardi <mlombard@…>

Children:

ca113cf

Parents:

f250c5a

git-author:

Maurizio Lombardi <mlombard@…> (2018-11-25 16:47:56)

git-committer:

Maurizio Lombardi <mlombard@…> (2021-06-18 16:06:03)

Message:

math: sync pow() to FreeBSD 11.2

File:

• uspace/lib/math/generic/pow.c (modified) (4 diffs)

uspace/lib/math/generic/pow.c

@@ -136 +136 @@
                 if(k>20) {
                     j = ly>>(52-k);
-                    if((j<<(52-k))==(int32_t)ly) yisint = 2-(j&1);
+                    if((j<<(52-k))==ly) yisint = 2-(j&1);
                 } else if(ly==0) {
                     j = iy>>(20-k);
@@ -148 +148 @@
             if (iy==0x7ff00000) {       /* y is +-inf */
                 if(((ix-0x3ff00000)|lx)==0)
-                    return  one;        /* (-1)**+-inf is NaN */
+                    return  one;        /* (-1)**+-inf is 1 */
                 else if (ix >= 0x3ff00000)/* (|x|>1)**+-inf = inf,0 */
                     return (hy>=0)? y: zero;
@@ -157 +157 @@
                 if(hy<0) return one/x; else return x;
             }
-            if(hy==0x40000000) return x*x;   /* y is  2 */
-            if(hy==0x40080000) return x*x*x; /* y is  3 */
-            if(hy==0x40100000) {             /* y is  4 */
-                u = x*x;
-                return u*u;
-            }
-            if(hy==0x3fe00000) {             /* y is  0.5 */
+            if(hy==0x40000000) return x*x; /* y is  2 */
+            if(hy==0x3fe00000) {        /* y is  0.5 */
                 if(hx>=0)       /* x >= +0 */
-                    return sqrt(x);
+                return sqrt(x);
             }
         }
@@ -314 +309 @@
         if((j>>20)<=0) z = scalbn(z,n); /* subnormal output */
         else SET_HIGH_WORD(z,j);
-
         return s*z;
 }