Back to index

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

Go to the source code of this file.

Functions

__complex__ long double __cacoshl (__complex__ long double x)

Function Documentation

__complex__ long double __cacoshl ( __complex__ long double  x)

Definition at line 26 of file s_cacoshl.c.

{
  __complex__ long double res;
  int rcls = fpclassify (__real__ x);
  int icls = fpclassify (__imag__ x);

  if (rcls <= FP_INFINITE || icls <= FP_INFINITE)
    {
      if (icls == FP_INFINITE)
       {
         __real__ res = HUGE_VALL;

         if (rcls == FP_NAN)
           __imag__ res = __nanl ("");
         else
           __imag__ res = __copysignl ((rcls == FP_INFINITE
                                    ? (__real__ x < 0.0
                                       ? M_PIl - M_PI_4l : M_PI_4l)
                                    : M_PI_2l), __imag__ x);
       }
      else if (rcls == FP_INFINITE)
       {
         __real__ res = HUGE_VALL;

         if (icls >= FP_ZERO)
           __imag__ res = __copysignl (signbit (__real__ x) ? M_PIl : 0.0,
                                   __imag__ x);
         else
           __imag__ res = __nanl ("");
       }
      else
       {
         __real__ res = __nanl ("");
         __imag__ res = __nanl ("");
       }
    }
  else if (rcls == FP_ZERO && icls == FP_ZERO)
    {
      __real__ res = 0.0;
      __imag__ res = __copysignl (M_PI_2l, __imag__ x);
    }
  else
    {
      __complex__ long double y;

      __real__ y = (__real__ x - __imag__ x) * (__real__ x + __imag__ x) - 1.0;
      __imag__ y = 2.0 * __real__ x * __imag__ x;

      y = __csqrtl (y);

      if (__real__ x < 0.0)
       y = -y;

      __real__ y += __real__ x;
      __imag__ y += __imag__ x;

      res = __clogl (y);

      /* We have to use the positive branch.  */
      if (__real__ res < 0.0)
       res = -res;
    }

  return res;
}

Here is the call graph for this function: