Changeset 992ffa6 in mainline for uspace/lib/math/generic/trig.c

Timestamp:

2015-09-04T06:40:20Z (9 years ago)

Author:

Jiri Svoboda <jiri@…>

Branches:

lfn, master, serial, ticket/834-toolchain-update, topic/msim-upgrade, topic/simplify-dev-export

Children:

01cdd5a

Parents:

bae1e1f

Message:

Add exp(f), log(f), pow(f). Improve precision of sin, cos.

File:

• uspace/lib/math/generic/trig.c (modified) (3 diffs)

uspace/lib/math/generic/trig.c

@@ -36 +36 @@
 #include <trig.h>
 
-#define TAYLOR_DEGREE  13
+#define TAYLOR_DEGREE_32 13
+#define TAYLOR_DEGREE_64 21
 
 /** Precomputed values for factorial (starting from 1!) */
-static float64_t factorials[TAYLOR_DEGREE] = {
+static float64_t factorials[TAYLOR_DEGREE_64] = {
         1, 2, 6, 24, 120, 720, 5040, 40320, 362880, 3628800, 39916800,
-        479001600, 6227020800
+        479001600, 6227020800.0L, 87178291200.0L, 1307674368000.0L,
+        20922789888000.0L, 355687428096000.0L, 6402373705728000.0L,
+        121645100408832000.0L, 2432902008176640000.0L, 51090942171709440000.0L
 };
 
@@ -61 +64 @@
         float64_t nom = 1;
         
-        for (unsigned int i = 0; i < TAYLOR_DEGREE; i++) {
+        for (unsigned int i = 0; i < TAYLOR_DEGREE_64; i++) {
                 nom *= arg;
                 
@@ -90 +93 @@
         float64_t nom = 1;
         
-        for (unsigned int i = 0; i < TAYLOR_DEGREE; i++) {
+        for (unsigned int i = 0; i < TAYLOR_DEGREE_64; i++) {
                 nom *= arg;