sysdeps/generic/s_ccosl.c

   1 /* Return cosine of complex long double value.
   2    Copyright (C) 1997 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
   5
   6    The GNU C Library is free software; you can redistribute it and/or
   7    modify it under the terms of the GNU Library General Public License as
   8    published by the Free Software Foundation; either version 2 of the
   9    License, or (at your option) any later version.
  10
  11    The GNU C Library is distributed in the hope that it will be useful,
  12    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14    Library General Public License for more details.
  15
  16    You should have received a copy of the GNU Library General Public
  17    License along with the GNU C Library; see the file COPYING.LIB.  If not,
  18    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  19    Boston, MA 02111-1307, USA.  */
  20
  21 #include <complex.h>
  22 #include <fenv.h>
  23 #include <math.h>
  24
  25
  26 __complex__ long double
  27 __ccosl (__complex__ long double x)
  28 {
  29   __complex__ long double res;
  30
  31   if (!isfinite (__real__ x) || __isnanl (__imag__ x))
  32     {
  33       if (__real__ x == 0.0 || __imag__ x == 0.0)
  34         {
  35           __real__ res = __nanl ("");
  36           __imag__ res = 0.0;
  37
  38 #ifdef FE_INVALID
  39           if (__isinfl (__real__ x))
  40             feraiseexcept (FE_INVALID);
  41 #endif
  42         }
  43       else if (__isinfl (__imag__ x))
  44         {
  45           __real__ res = HUGE_VALL;
  46           __imag__ res = __nanl ("");
  47
  48 #ifdef FE_INVALID
  49           if (__isinfl (__real__ x))
  50             feraiseexcept (FE_INVALID);
  51 #endif
  52         }
  53       else
  54         {
  55           __real__ res = __nanl ("");
  56           __imag__ res = __nanl ("");
  57
  58 #ifdef FE_INVALID
  59           if (isfinite (__imag__ x))
  60             feraiseexcept (FE_INVALID);
  61 #endif
  62         }
  63     }
  64   else
  65     {
  66       __complex__ long double y;
  67
  68       __real__ y = -__imag__ x;
  69       __imag__ y = __real__ x;
  70
  71       res = __ccoshl (y);
  72     }
  73
  74   return res;
  75 }
  76 weak_alias (__ccosl, ccosl)