Back to index

glibc  2.9
Functions
s_csinhl.c File Reference
#include <complex.h>
#include <fenv.h>
#include <math.h>
#include "math_private.h"
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

__complex__ long double __csinhl (__complex__ long double x)

Function Documentation

__complex__ long double __csinhl ( __complex__ long double  x)

Definition at line 29 of file s_csinhl.c.

{
  __complex__ long double retval;
  int negate = signbit (__real__ x);
  int rcls = fpclassify (__real__ x);
  int icls = fpclassify (__imag__ x);

  __real__ x = fabsl (__real__ x);

  if (rcls >= FP_ZERO)
    {
      /* Real part is finite.  */
      if (icls >= FP_ZERO)
       {
         /* Imaginary part is finite.  */
         long double sinh_val = __ieee754_sinhl (__real__ x);
         long double cosh_val = __ieee754_coshl (__real__ x);
         long double sinix, cosix;

         __sincosl (__imag__ x, &sinix, &cosix);

         __real__ retval = sinh_val * cosix;
         __imag__ retval = cosh_val * sinix;

         if (negate)
           __real__ retval = -__real__ retval;
       }
      else
       {
         if (rcls == FP_ZERO)
           {
             /* Real part is 0.0.  */
             __real__ retval = __copysignl (0.0, negate ? -1.0 : 1.0);
             __imag__ retval = __nanl ("") + __nanl ("");

#ifdef FE_INVALID
             if (icls == FP_INFINITE)
              feraiseexcept (FE_INVALID);
#endif
           }
         else
           {
             __real__ retval = __nanl ("");
             __imag__ retval = __nanl ("");

#ifdef FE_INVALID
             feraiseexcept (FE_INVALID);
#endif
           }
       }
    }
  else if (rcls == FP_INFINITE)
    {
      /* Real part is infinite.  */
      if (icls == FP_ZERO)
       {
         /* Imaginary part is 0.0.  */
         __real__ retval = negate ? -HUGE_VALL : HUGE_VALL;
         __imag__ retval = __imag__ x;
       }
      else if (icls > FP_ZERO)
       {
         /* Imaginary part is finite.  */
         long double sinix, cosix;

         __sincosl (__imag__ x, &sinix, &cosix);

         __real__ retval = __copysignl (HUGE_VALL, cosix);
         __imag__ retval = __copysignl (HUGE_VALL, sinix);

         if (negate)
           __real__ retval = -__real__ retval;
       }
      else
       {
         /* The addition raises the invalid exception.  */
         __real__ retval = HUGE_VALL;
         __imag__ retval = __nanl ("") + __nanl ("");

#ifdef FE_INVALID
         if (icls == FP_INFINITE)
           feraiseexcept (FE_INVALID);
#endif
       }
    }
  else
    {
      __real__ retval = __nanl ("");
      __imag__ retval = __imag__ x == 0.0 ? __imag__ x : __nanl ("");
    }

  return retval;
}

Here is the call graph for this function: