Back to index

glibc  2.9
Functions
s_csinhf.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__ float __csinhf (__complex__ float x)

Function Documentation

__complex__ float __csinhf ( __complex__ float  x)

Definition at line 29 of file s_csinhf.c.

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

  __real__ x = fabsf (__real__ x);

  if (rcls >= FP_ZERO)
    {
      /* Real part is finite.  */
      if (icls >= FP_ZERO)
       {
         /* Imaginary part is finite.  */
         float sinh_val = __ieee754_sinhf (__real__ x);
         float cosh_val = __ieee754_coshf (__real__ x);
         float sinix, cosix;

         __sincosf (__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 = __copysignf (0.0, negate ? -1.0 : 1.0);
             __imag__ retval = __nanf ("") + __nanf ("");

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

#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_VALF : HUGE_VALF;
         __imag__ retval = __imag__ x;
       }
      else if (icls > FP_ZERO)
       {
         /* Imaginary part is finite.  */
         float sinix, cosix;

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

         __real__ retval = __copysignf (HUGE_VALF, cosix);
         __imag__ retval = __copysignf (HUGE_VALF, sinix);

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

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

  return retval;
}

Here is the call graph for this function: