summaryrefslogtreecommitdiff
path: root/newlib/libm/mathfp/sf_tanh.c
diff options
context:
space:
mode:
authorPikalaxALT <pikalaxalt@gmail.com>2018-01-03 17:39:24 -0700
committerPikalaxALT <pikalaxalt@gmail.com>2018-01-03 17:39:24 -0700
commita6c1ed4716cf02626ea035beb6dd4a921642ba80 (patch)
treeef582c1b52819e27bdd16097ec03b69799d04ede /newlib/libm/mathfp/sf_tanh.c
parentf6c9a624fa8a6878a7fb2b02f55e4990a20feb59 (diff)
Use libc from agbcc instead of standalone newlib\nYou must have AGBCC commit 80d029caec189587f8b9294b6c8a5a489b8f5f88 in order to compile pmd_red.gbalibc
Diffstat (limited to 'newlib/libm/mathfp/sf_tanh.c')
-rw-r--r--newlib/libm/mathfp/sf_tanh.c77
1 files changed, 0 insertions, 77 deletions
diff --git a/newlib/libm/mathfp/sf_tanh.c b/newlib/libm/mathfp/sf_tanh.c
deleted file mode 100644
index 51806af..0000000
--- a/newlib/libm/mathfp/sf_tanh.c
+++ /dev/null
@@ -1,77 +0,0 @@
-
-/* @(#)z_tanhf.c 1.0 98/08/13 */
-/*****************************************************************
- * The following routines are coded directly from the algorithms
- * and coefficients given in "Software Manual for the Elementary
- * Functions" by William J. Cody, Jr. and William Waite, Prentice
- * Hall, 1980.
- *****************************************************************/
-/******************************************************************
- * Hyperbolic Tangent
- *
- * Input:
- * x - floating point value
- *
- * Output:
- * hyperbolic tangent of x
- *
- * Description:
- * This routine calculates hyperbolic tangent.
- *
- *****************************************************************/
-
-#include <float.h>
-#include "fdlibm.h"
-#include "zmath.h"
-
-static const float LN3_OVER2 = 0.5493061443;
-static const float p[] = { -0.2059432032,
- -0.0009577527 };
-static const float q[] = { 0.6178299136,
- 0.25 };
-
-float
-_DEFUN (tanhf, (float),
- float x)
-{
- float f, res, g, P, Q, R;
-
- f = fabsf (x);
-
- /* Check if the input is too big. */
- if (f > BIGX)
- res = 1.0;
-
- else if (f > LN3_OVER2)
- res = 1.0 - 2.0 / (exp (2 * f) + 1.0);
-
- /* Check if the input is too small. */
- else if (f < z_rooteps_f)
- res = f;
-
- /* Calculate the Taylor series. */
- else
- {
- g = f * f;
-
- P = p[1] * g + p[0];
- Q = (g + q[1]) * g + q[0];
- R = g * (P / Q);
-
- res = f + f * R;
- }
-
- if (x < 0.0)
- res = -res;
-
- return (res);
-}
-
-#ifdef _DOUBLE_IS_32BITS
-
-double tanh (double x)
-{
- return (double) tanhf ((float) x);
-}
-
-#endif _DOUBLE_IS_32BITS