libquadmath/math/complex.c

   1 /* GCC Quad-Precision Math Library
   2    Copyright (C) 2010, 2011 Free Software Foundation, Inc.
   3    Written by Francois-Xavier Coudert  <fxcoudert@gcc.gnu.org>
   4
   5 This file is part of the libquadmath library.
   6 Libquadmath is free software; you can redistribute it and/or
   7 modify it under the terms of the GNU Library General Public
   8 License as published by the Free Software Foundation; either
   9 version 2 of the License, or (at your option) any later version.
  10
  11 Libquadmath 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 libquadmath; see the file COPYING.LIB.  If
  18 not, write to the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
  19 Boston, MA 02110-1301, USA.  */
  20
  21 #include "quadmath-imp.h"
  22
  23 #ifdef HAVE_FENV_H
  24 # include <fenv.h>
  25 #endif
  26
  27
  28 #define REALPART(z) (__real__(z))
  29 #define IMAGPART(z) (__imag__(z))
  30 #define COMPLEX_ASSIGN(z_, r_, i_) {__real__(z_) = (r_); __imag__(z_) = (i_);}
  31
  32
  33 __float128
  34 cabsq (__complex128 z)
  35 {
  36   return hypotq (REALPART (z), IMAGPART (z));
  37 }
  38
  39
  40 __complex128
  41 cexpiq (__float128 x)
  42 {
  43   __float128 sinix, cosix;
  44   __complex128 v;
  45   sincosq (x, &sinix, &cosix);
  46   COMPLEX_ASSIGN (v, cosix, sinix);
  47   return v;
  48 }
  49
  50
  51 __float128
  52 cargq (__complex128 z)
  53 {
  54   return atan2q (IMAGPART (z), REALPART (z));
  55 }
  56
  57
  58 __complex128
  59 cpowq (__complex128 base, __complex128 power)
  60 {
  61   return cexpq (power * clogq (base));
  62 }
  63
  64
  65 __complex128
  66 ccosq (__complex128 x)
  67 {
  68   __complex128 y;
  69
  70   COMPLEX_ASSIGN (y, -IMAGPART (x), REALPART (x));
  71   return ccoshq (y);
  72 }