Changeset 777 for soft/giet_vm/giet_libs/math/e_pow.c

Timestamp:

Feb 4, 2016, 6:25:22 PM (9 years ago)

Author:

meunier

Message:

• Ajout de quelques fonction dans la lib math
• Déplacement de certaines fonctions de stdlib vers string

File:

• soft/giet_vm/giet_libs/math/e_pow.c (modified) (7 diffs)

soft/giet_vm/giet_libs/math/e_pow.c

@@ -103 +103 @@
         int32_t i,j,k,yisint,n;
         int32_t hx,hy,ix,iy;
-        u_int32_t lx,ly;
+        uint32_t lx,ly;
 
         EXTRACT_WORDS(hx,lx,x);
@@ -180 +180 @@
 
     /* (x<0)**(non-int) is NaN */
-        if(((((u_int32_t)hx>>31)-1)|yisint)==0) return (x-x)/(x-x);
+        if(((((uint32_t)hx>>31)-1)|yisint)==0) return (x-x)/(x-x);
 
     /* |y| is huge */
@@ -251 +251 @@
 
         s = one; /* s (sign of result -ve**odd) = -1 else = 1 */
-        if(((((u_int32_t)hx>>31)-1)|(yisint-1))==0)
+        if(((((uint32_t)hx>>31)-1)|(yisint-1))==0)
             s = -one;/* (-ve)**(odd int) */
 
@@ -312 +312 @@
 {
         double z = __ieee754_pow(x, y);
-        if (isnan(y))
+        if (isnan(y)) {
                 return z;
+    }
         if (isnan(x)) {
                 if (y == 0.0) {
@@ -322 +323 @@
         if (x == 0.0) {
                 if (y == 0.0) {
-                        return 0; /* pow(0.0,0.0) */
+            return 0; /* pow(0.0,0.0) */
                 }
                 if (isfinite(y) && y < 0.0) {
@@ -332 +333 @@
                 if (isfinite(x) && isfinite(y)) {
                         if (isnan(z)) {
-                                return 0.0/0.0; /* pow neg**non-int */
+                                return 0.0 / 0.0; /* pow neg**non-int */
                         }
                         {
                                 double ret = huge;
                                 y *= 0.5;
-                                if(x<zero&&rint(y)!=y) ret = -huge;
+                                if (x < zero && rint(y) !=y ) {
+                    ret = -huge;
+                }
                                 return ret;
                         }
@@ -347 +350 @@
         return z;
 }
+
+/*
+ * Pseudo exp function (QM)
+ * There probably is a better implementation
+ */
+double exp(double x)
+{
+    const double e = 2.71828182846;
+    return pow(e, x);
+}
+