Changeset 83932dc9 in mainline for uspace/lib/math/generic/trig.c

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.

File:

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

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);
-        }
-}
-
 /** @}
  */