OpenCores
URL https://opencores.org/ocsvn/openrisc_me/openrisc_me/trunk

Subversion Repositories openrisc_me

[/] [openrisc/] [trunk/] [gnu-src/] [newlib-1.18.0/] [newlib/] [libm/] [mathfp/] [sf_tanh.c] - Blame information for rev 207

Details | Compare with Previous | View Log

Line No. Rev Author Line
1 207 jeremybenn
 
2
/* @(#)z_tanhf.c 1.0 98/08/13 */
3
/*****************************************************************
4
 * The following routines are coded directly from the algorithms
5
 * and coefficients given in "Software Manual for the Elementary
6
 * Functions" by William J. Cody, Jr. and William Waite, Prentice
7
 * Hall, 1980.
8
 *****************************************************************/
9
/******************************************************************
10
 * Hyperbolic Tangent
11
 *
12
 * Input:
13
 *   x - floating point value
14
 *
15
 * Output:
16
 *   hyperbolic tangent of x
17
 *
18
 * Description:
19
 *   This routine calculates hyperbolic tangent.
20
 *
21
 *****************************************************************/
22
 
23
#include <float.h>
24
#include "fdlibm.h"
25
#include "zmath.h"
26
 
27
static const float LN3_OVER2 = 0.5493061443;
28
static const float p[] = { -0.2059432032,
29
                           -0.0009577527 };
30
static const float q[] = {  0.6178299136,
31
                            0.25 };
32
 
33
float
34
_DEFUN (tanhf, (float),
35
        float x)
36
{
37
  float f, res, g, P, Q, R;
38
 
39
  f = fabsf (x);
40
 
41
  /* Check if the input is too big. */
42
  if (f > BIGX)
43
    res = 1.0;
44
 
45
  else if (f > LN3_OVER2)
46
    res = 1.0 - 2.0 / (exp (2 * f) + 1.0);
47
 
48
  /* Check if the input is too small. */
49
  else if (f < z_rooteps_f)
50
    res = f;
51
 
52
  /* Calculate the Taylor series. */
53
  else
54
    {
55
      g = f * f;
56
 
57
      P = p[1] * g + p[0];
58
      Q = (g + q[1]) * g + q[0];
59
      R = g * (P / Q);
60
 
61
      res = f + f * R;
62
    }
63
 
64
  if (x < 0.0)
65
    res = -res;
66
 
67
  return (res);
68
}
69
 
70
#ifdef _DOUBLE_IS_32BITS
71
 
72
double tanh (double x)
73
{
74
  return (double) tanhf ((float) x);
75
}
76
 
77
#endif _DOUBLE_IS_32BITS

powered by: WebSVN 2.1.0

© copyright 1999-2024 OpenCores.org, equivalent to Oliscience, all rights reserved. OpenCores®, registered trademark.