Back to index

glibc  2.9
Functions
s_ctanhl.c File Reference
#include <complex.h>
#include <fenv.h>
#include <math.h>
#include <math_ldbl_opt.h>
#include "math_private.h"

Go to the source code of this file.

Functions

__complex__ long double __ctanhl (__complex__ long double x)
 long_double_symbol (libm, __ctanhl, ctanhl)

Function Documentation

__complex__ long double __ctanhl ( __complex__ long double  x)

Definition at line 30 of file s_ctanhl.c.

{
  __complex__ long double res;

  if (!isfinite (__real__ x) || !isfinite (__imag__ x))
    {
      if (__isinfl (__real__ x))
       {
         __real__ res = __copysignl (1.0, __real__ x);
         __imag__ res = __copysignl (0.0, __imag__ x);
       }
      else if (__imag__ x == 0.0)
       {
         res = x;
       }
      else
       {
         __real__ res = __nanl ("");
         __imag__ res = __nanl ("");

#ifdef FE_INVALID
         if (__isinfl (__imag__ x))
           feraiseexcept (FE_INVALID);
#endif
       }
    }
  else
    {
      long double sin2ix, cos2ix;
      long double den;

      __sincosl (2.0 * __imag__ x, &sin2ix, &cos2ix);

      den = (__ieee754_coshl (2.0 * __real__ x) + cos2ix);

      if (den == 0.0L)
       {
         __complex__ long double ez = __cexpl (x);
         __complex__ long double emz = __cexpl (-x);

         res = (ez - emz) / (ez + emz);
       }
      else
       {
         __real__ res = __ieee754_sinhl (2.0 * __real__ x) / den;
         __imag__ res = sin2ix / den;
       }
      /* __gcc_qmul does not respect -0.0 so we need the following fixup.  */
      if ((__real__ res == 0.0) && (__real__ x == 0.0))
        __real__ res = __real__ x;

      if ((__real__ res == 0.0) && (__imag__ x == 0.0))
        __imag__ res = __imag__ x;
    }

  return res;
}

Here is the call graph for this function:

long_double_symbol ( libm  ,
__ctanhl  ,
ctanhl   
)