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  }