|
|
/* @(#)z_tanhf.c 1.0 98/08/13 */
|
/* @(#)z_tanhf.c 1.0 98/08/13 */
|
/*****************************************************************
|
/*****************************************************************
|
* The following routines are coded directly from the algorithms
|
* The following routines are coded directly from the algorithms
|
* and coefficients given in "Software Manual for the Elementary
|
* and coefficients given in "Software Manual for the Elementary
|
* Functions" by William J. Cody, Jr. and William Waite, Prentice
|
* Functions" by William J. Cody, Jr. and William Waite, Prentice
|
* Hall, 1980.
|
* Hall, 1980.
|
*****************************************************************/
|
*****************************************************************/
|
/******************************************************************
|
/******************************************************************
|
* Hyperbolic Tangent
|
* Hyperbolic Tangent
|
*
|
*
|
* Input:
|
* Input:
|
* x - floating point value
|
* x - floating point value
|
*
|
*
|
* Output:
|
* Output:
|
* hyperbolic tangent of x
|
* hyperbolic tangent of x
|
*
|
*
|
* Description:
|
* Description:
|
* This routine calculates hyperbolic tangent.
|
* This routine calculates hyperbolic tangent.
|
*
|
*
|
*****************************************************************/
|
*****************************************************************/
|
|
|
#include <float.h>
|
#include <float.h>
|
#include "fdlibm.h"
|
#include "fdlibm.h"
|
#include "zmath.h"
|
#include "zmath.h"
|
|
|
static const float LN3_OVER2 = 0.5493061443;
|
static const float LN3_OVER2 = 0.5493061443;
|
static const float p[] = { -0.2059432032,
|
static const float p[] = { -0.2059432032,
|
-0.0009577527 };
|
-0.0009577527 };
|
static const float q[] = { 0.6178299136,
|
static const float q[] = { 0.6178299136,
|
0.25 };
|
0.25 };
|
|
|
float
|
float
|
_DEFUN (tanhf, (float),
|
_DEFUN (tanhf, (float),
|
float x)
|
float x)
|
{
|
{
|
float f, res, g, P, Q, R;
|
float f, res, g, P, Q, R;
|
|
|
f = fabsf (x);
|
f = fabsf (x);
|
|
|
/* Check if the input is too big. */
|
/* Check if the input is too big. */
|
if (f > BIGX)
|
if (f > BIGX)
|
res = 1.0;
|
res = 1.0;
|
|
|
else if (f > LN3_OVER2)
|
else if (f > LN3_OVER2)
|
res = 1.0 - 2.0 / (exp (2 * f) + 1.0);
|
res = 1.0 - 2.0 / (exp (2 * f) + 1.0);
|
|
|
/* Check if the input is too small. */
|
/* Check if the input is too small. */
|
else if (f < z_rooteps_f)
|
else if (f < z_rooteps_f)
|
res = f;
|
res = f;
|
|
|
/* Calculate the Taylor series. */
|
/* Calculate the Taylor series. */
|
else
|
else
|
{
|
{
|
g = f * f;
|
g = f * f;
|
|
|
P = p[1] * g + p[0];
|
P = p[1] * g + p[0];
|
Q = (g + q[1]) * g + q[0];
|
Q = (g + q[1]) * g + q[0];
|
R = g * (P / Q);
|
R = g * (P / Q);
|
|
|
res = f + f * R;
|
res = f + f * R;
|
}
|
}
|
|
|
if (x < 0.0)
|
if (x < 0.0)
|
res = -res;
|
res = -res;
|
|
|
return (res);
|
return (res);
|
}
|
}
|
|
|
#ifdef _DOUBLE_IS_32BITS
|
#ifdef _DOUBLE_IS_32BITS
|
|
|
double tanh (double x)
|
double tanh (double x)
|
{
|
{
|
return (double) tanhf ((float) x);
|
return (double) tanhf ((float) x);
|
}
|
}
|
|
|
#endif _DOUBLE_IS_32BITS
|
#endif _DOUBLE_IS_32BITS
|
|
|
