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 } 73