Changeset 83932dc9 in mainline

Timestamp:

2018-09-12T16:51:07Z (6 years ago)

Author:

Jiří Zárevúcky <jiri.zarevucky@…>

Branches:

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

Children:

b2aaaa0

Parents:

ec9aadd

git-author:

Jiří Zárevúcky <jiri.zarevucky@…> (2018-09-12 16:50:02)

git-committer:

Jiří Zárevúcky <jiri.zarevucky@…> (2018-09-12 16:51:07)

Message:

Split up libmath functions into more files

so that they don't conflict with fdlibm during static linking.

Location:

uspace/lib/math

Files:

• Makefile (modified) (1 diff)
• generic/cos.c (added)
• generic/fmod.c (modified) (1 diff)
• generic/fmodf.c (added)
• generic/internal.h (added)
• generic/sin.c (added)
• generic/sincos.c (added)
• generic/trig.c (modified) (6 diffs)

uspace/lib/math/Makefile

@@ -39 +39 @@
         generic/fabs.c \
         generic/fmod.c \
+        generic/fmodf.c \
         generic/nearbyint.c \
         generic/round.c \
         generic/trig.c \
+        generic/sin.c \
+        generic/cos.c \
+        generic/sincos.c \
         generic/trunc.c
 

uspace/lib/math/generic/fmod.c

@@ -35 +35 @@
 #include <math.h>
 
-/** Remainder function (32-bit floating point)
- *
- * Calculate the modulo of dividend by divisor.
- *
- * This is a very basic implementation that uses
- * division and multiplication (instead of exact
- * arithmetics). Thus the result might be very
- * imprecise (depending on the magnitude of the
- * arguments).
- *
- * @param dividend Dividend.
- * @param divisor  Divisor.
- *
- * @return Modulo.
- *
- */
-float fmodf(float dividend, float divisor)
-{
-        // FIXME: replace with exact arithmetics
-
-        float quotient = truncf(dividend / divisor);
-
-        return (dividend - quotient * divisor);
-}
-
 /** Remainder function (64-bit floating point)
  *

uspace/lib/math/generic/trig.c

@@ -35 +35 @@
 
 #include <math.h>
+#include "internal.h"
 
 #define TAYLOR_DEGREE_32 13
@@ -175 +176 @@
  *
  */
-static float base_sin_32(float arg)
+float __math_base_sin_32(float arg)
 {
         unsigned int period = arg / (M_PI / 4);
@@ -208 +209 @@
  *
  */
-static double base_sin_64(double arg)
+double __math_base_sin_64(double arg)
 {
         unsigned int period = arg / (M_PI / 4);
@@ -241 +242 @@
  *
  */
-static float base_cos_32(float arg)
+float __math_base_cos_32(float arg)
 {
         unsigned int period = arg / (M_PI / 4);
@@ -274 +275 @@
  *
  */
-static double base_cos_64(double arg)
+double __math_base_cos_64(double arg)
 {
         unsigned int period = arg / (M_PI / 4);
@@ -295 +296 @@
 }
 
-/** Sine (32-bit floating point)
- *
- * Compute sine value.
- *
- * @param arg Sine argument.
- *
- * @return Sine value.
- *
- */
-float sinf(float arg)
-{
-        float base_arg = fmodf(arg, 2 * M_PI);
-
-        if (base_arg < 0)
-                return -base_sin_32(-base_arg);
-
-        return base_sin_32(base_arg);
-}
-
-/** Sine (64-bit floating point)
- *
- * Compute sine value.
- *
- * @param arg Sine argument.
- *
- * @return Sine value.
- *
- */
-double sin(double arg)
-{
-        double base_arg = fmod(arg, 2 * M_PI);
-
-        if (base_arg < 0)
-                return -base_sin_64(-base_arg);
-
-        return base_sin_64(base_arg);
-}
-
-/** Cosine (32-bit floating point)
- *
- * Compute cosine value.
- *
- * @param arg Cosine argument.
- *
- * @return Cosine value.
- *
- */
-float cosf(float arg)
-{
-        float base_arg = fmodf(arg, 2 * M_PI);
-
-        if (base_arg < 0)
-                return base_cos_32(-base_arg);
-
-        return base_cos_32(base_arg);
-}
-
-/** Cosine (64-bit floating point)
- *
- * Compute cosine value.
- *
- * @param arg Cosine argument.
- *
- * @return Cosine value.
- *
- */
-double cos(double arg)
-{
-        double base_arg = fmod(arg, 2 * M_PI);
-
-        if (base_arg < 0)
-                return base_cos_64(-base_arg);
-
-        return base_cos_64(base_arg);
-}
-
-/**
- * Computes sine and cosine at the same time, which might be more efficient than
- * computing each separately.
- *
- * @param x  Input value.
- * @param s  Output sine value, *s = sinf(x).
- * @param c  Output cosine value, *c = cosf(x).
- */
-void sincosf(float x, float *s, float *c)
-{
-        float base_arg = fmodf(x, 2 * M_PI);
-
-        if (base_arg < 0) {
-                *s = -base_sin_32(-base_arg);
-                *c = base_cos_32(-base_arg);
-        } else {
-                *s = base_sin_32(base_arg);
-                *c = base_cos_32(base_arg);
-        }
-}
-
-/**
- * Computes sine and cosine at the same time, which might be more efficient than
- * computing each separately.
- *
- * @param x  Input value.
- * @param s  Output sine value, *s = sin(x).
- * @param c  Output cosine value, *c = cos(x).
- */
-void sincos(double x, double *s, double *c)
-{
-        double base_arg = fmod(x, 2 * M_PI);
-
-        if (base_arg < 0) {
-                *s = -base_sin_64(-base_arg);
-                *c = base_cos_64(-base_arg);
-        } else {
-                *s = base_sin_64(base_arg);
-                *c = base_cos_64(base_arg);
-        }
-}
-
 /** @}
  */