| 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 | 
 #ifndef QUADMATH_H | 
 
 
 
 
 
 | 22 | 
 #define QUADMATH_H | 
 
 
 
 
 
 | 23 | 
  | 
 
 
 
 
 
 | 24 | 
 #include <stdlib.h> | 
 
 
 
 
 
 | 25 | 
  | 
 
 
 
 
 
 | 26 | 
 #ifdef __cplusplus | 
 
 
 
 
 
 | 27 | 
 extern "C" { | 
 
 
 
 
 
 | 28 | 
 #endif | 
 
 
 
 
 
 | 29 | 
  | 
 
 
 
 
 
 | 30 | 
 /* Define the complex type corresponding to __float128 | 
 
 
 
 
 
 | 31 | 
    ("_Complex __float128" is not allowed) */ | 
 
 
 
 
 
 | 32 | 
 typedef _Complex float __attribute__((mode(TC))) __complex128; | 
 
 
 
 
 
 | 33 | 
  | 
 
 
 
 
 
 | 34 | 
 #ifdef __cplusplus | 
 
 
 
 
 
 | 35 | 
 # define __quadmath_throw throw () | 
 
 
 
 
 
 | 36 | 
 # define __quadmath_nth(fct) fct throw () | 
 
 
 
 
 
 | 37 | 
 #else | 
 
 
 
 
 
 | 38 | 
 # define __quadmath_throw __attribute__((__nothrow__)) | 
 
 
 
 
 
 | 39 | 
 # define __quadmath_nth(fct) __attribute__((__nothrow__)) fct | 
 
 
 
 
 
 | 40 | 
 #endif | 
 
 
 
 
 
 | 41 | 
  | 
 
 
 
 
 
 | 42 | 
 /* Prototypes for real functions */ | 
 
 
 
 
 
 | 43 | 
 extern __float128 acosq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 44 | 
 extern __float128 acoshq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 45 | 
 extern __float128 asinq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 46 | 
 extern __float128 asinhq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 47 | 
 extern __float128 atanq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 48 | 
 extern __float128 atanhq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 49 | 
 extern __float128 atan2q (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 50 | 
 extern __float128 cbrtq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 51 | 
 extern __float128 ceilq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 52 | 
 extern __float128 copysignq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 53 | 
 extern __float128 coshq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 54 | 
 extern __float128 cosq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 55 | 
 extern __float128 erfq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 56 | 
 extern __float128 erfcq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 57 | 
 extern __float128 expq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 58 | 
 extern __float128 expm1q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 59 | 
 extern __float128 fabsq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 60 | 
 extern __float128 fdimq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 61 | 
 extern int finiteq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 62 | 
 extern __float128 floorq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 63 | 
 extern __float128 fmaq (__float128, __float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 64 | 
 extern __float128 fmaxq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 65 | 
 extern __float128 fminq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 66 | 
 extern __float128 fmodq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 67 | 
 extern __float128 frexpq (__float128, int *) __quadmath_throw; | 
 
 
 
 
 
 | 68 | 
 extern __float128 hypotq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 69 | 
 extern int isinfq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 70 | 
 extern int ilogbq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 71 | 
 extern int isnanq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 72 | 
 extern __float128 j0q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 73 | 
 extern __float128 j1q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 74 | 
 extern __float128 jnq (int, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 75 | 
 extern __float128 ldexpq (__float128, int) __quadmath_throw; | 
 
 
 
 
 
 | 76 | 
 extern __float128 lgammaq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 77 | 
 extern long long int llrintq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 78 | 
 extern long long int llroundq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 79 | 
 extern __float128 logq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 80 | 
 extern __float128 log10q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 81 | 
 extern __float128 log2q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 82 | 
 extern __float128 log1pq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 83 | 
 extern long int lrintq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 84 | 
 extern long int lroundq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 85 | 
 extern __float128 modfq (__float128, __float128 *) __quadmath_throw; | 
 
 
 
 
 
 | 86 | 
 extern __float128 nanq (const char *) __quadmath_throw; | 
 
 
 
 
 
 | 87 | 
 extern __float128 nearbyintq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 88 | 
 extern __float128 nextafterq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 89 | 
 extern __float128 powq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 90 | 
 extern __float128 remainderq (__float128, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 91 | 
 extern __float128 remquoq (__float128, __float128, int *) __quadmath_throw; | 
 
 
 
 
 
 | 92 | 
 extern __float128 rintq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 93 | 
 extern __float128 roundq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 94 | 
 extern __float128 scalblnq (__float128, long int) __quadmath_throw; | 
 
 
 
 
 
 | 95 | 
 extern __float128 scalbnq (__float128, int) __quadmath_throw; | 
 
 
 
 
 
 | 96 | 
 extern int signbitq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 97 | 
 extern void sincosq (__float128, __float128 *, __float128 *) __quadmath_throw; | 
 
 
 
 
 
 | 98 | 
 extern __float128 sinhq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 99 | 
 extern __float128 sinq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 100 | 
 extern __float128 sqrtq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 101 | 
 extern __float128 tanq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 102 | 
 extern __float128 tanhq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 103 | 
 extern __float128 tgammaq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 104 | 
 extern __float128 truncq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 105 | 
 extern __float128 y0q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 106 | 
 extern __float128 y1q (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 107 | 
 extern __float128 ynq (int, __float128) __quadmath_throw; | 
 
 
 
 
 
 | 108 | 
  | 
 
 
 
 
 
 | 109 | 
  | 
 
 
 
 
 
 | 110 | 
 /* Prototypes for complex functions */ | 
 
 
 
 
 
 | 111 | 
 extern __float128 cabsq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 112 | 
 extern __float128 cargq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 113 | 
 extern __float128 cimagq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 114 | 
 extern __float128 crealq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 115 | 
 extern __complex128 cacosq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 116 | 
 extern __complex128 cacoshq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 117 | 
 extern __complex128 casinq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 118 | 
 extern __complex128 casinhq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 119 | 
 extern __complex128 catanq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 120 | 
 extern __complex128 catanhq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 121 | 
 extern __complex128 ccosq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 122 | 
 extern __complex128 ccoshq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 123 | 
 extern __complex128 cexpq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 124 | 
 extern __complex128 cexpiq (__float128) __quadmath_throw; | 
 
 
 
 
 
 | 125 | 
 extern __complex128 clogq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 126 | 
 extern __complex128 clog10q (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 127 | 
 extern __complex128 conjq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 128 | 
 extern __complex128 cpowq (__complex128, __complex128) __quadmath_throw; | 
 
 
 
 
 
 | 129 | 
 extern __complex128 cprojq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 130 | 
 extern __complex128 csinq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 131 | 
 extern __complex128 csinhq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 132 | 
 extern __complex128 csqrtq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 133 | 
 extern __complex128 ctanq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 134 | 
 extern __complex128 ctanhq (__complex128) __quadmath_throw; | 
 
 
 
 
 
 | 135 | 
  | 
 
 
 
 
 
 | 136 | 
  | 
 
 
 
 
 
 | 137 | 
 /* Prototypes for string <-> __float128 conversion functions */ | 
 
 
 
 
 
 | 138 | 
 extern __float128 strtoflt128 (const char *, char **) __quadmath_throw; | 
 
 
 
 
 
 | 139 | 
 extern int quadmath_snprintf (char *str, size_t size, | 
 
 
 
 
 
 | 140 | 
                               const char *format, ...) __quadmath_throw; | 
 
 
 
 
 
 | 141 | 
  | 
 
 
 
 
 
 | 142 | 
  | 
 
 
 
 
 
 | 143 | 
 /* Macros */ | 
 
 
 
 
 
 | 144 | 
 #define FLT128_MAX 1.18973149535723176508575932662800702e4932Q | 
 
 
 
 
 
 | 145 | 
 #define FLT128_MIN 3.36210314311209350626267781732175260e-4932Q | 
 
 
 
 
 
 | 146 | 
 #define FLT128_EPSILON 1.92592994438723585305597794258492732e-34Q | 
 
 
 
 
 
 | 147 | 
 #define FLT128_DENORM_MIN 6.475175119438025110924438958227646552e-4966Q | 
 
 
 
 
 
 | 148 | 
 #define FLT128_MANT_DIG 113 | 
 
 
 
 
 
 | 149 | 
 #define FLT128_MIN_EXP (-16381) | 
 
 
 
 
 
 | 150 | 
 #define FLT128_MAX_EXP 16384 | 
 
 
 
 
 
 | 151 | 
 #define FLT128_DIG 33 | 
 
 
 
 
 
 | 152 | 
 #define FLT128_MIN_10_EXP (-4931) | 
 
 
 
 
 
 | 153 | 
 #define FLT128_MAX_10_EXP 4932 | 
 
 
 
 
 
 | 154 | 
  | 
 
 
 
 
 
 | 155 | 
  | 
 
 
 
 
 
 | 156 | 
 #define HUGE_VALQ __builtin_huge_valq() | 
 
 
 
 
 
 | 157 | 
 /* The following alternative is valid, but brings the warning: | 
 
 
 
 
 
 | 158 | 
    (floating constant exceeds range of ‘__float128’)  */ | 
 
 
 
 
 
 | 159 | 
 /* #define HUGE_VALQ (__extension__ 0x1.0p32767Q) */ | 
 
 
 
 
 
 | 160 | 
  | 
 
 
 
 
 
 | 161 | 
 #define M_Eq            2.7182818284590452353602874713526625Q  /* e */ | 
 
 
 
 
 
 | 162 | 
 #define M_LOG2Eq        1.4426950408889634073599246810018921Q  /* log_2 e */ | 
 
 
 
 
 
 | 163 | 
 #define M_LOG10Eq       0.4342944819032518276511289189166051Q  /* log_10 e */ | 
 
 
 
 
 
 | 164 | 
 #define M_LN2q          0.6931471805599453094172321214581766Q  /* log_e 2 */ | 
 
 
 
 
 
 | 165 | 
 #define M_LN10q         2.3025850929940456840179914546843642Q  /* log_e 10 */ | 
 
 
 
 
 
 | 166 | 
 #define M_PIq           3.1415926535897932384626433832795029Q  /* pi */ | 
 
 
 
 
 
 | 167 | 
 #define M_PI_2q         1.5707963267948966192313216916397514Q  /* pi/2 */ | 
 
 
 
 
 
 | 168 | 
 #define M_PI_4q         0.7853981633974483096156608458198757Q  /* pi/4 */ | 
 
 
 
 
 
 | 169 | 
 #define M_1_PIq         0.3183098861837906715377675267450287Q  /* 1/pi */ | 
 
 
 
 
 
 | 170 | 
 #define M_2_PIq         0.6366197723675813430755350534900574Q  /* 2/pi */ | 
 
 
 
 
 
 | 171 | 
 #define M_2_SQRTPIq     1.1283791670955125738961589031215452Q  /* 2/sqrt(pi) */ | 
 
 
 
 
 
 | 172 | 
 #define M_SQRT2q        1.4142135623730950488016887242096981Q  /* sqrt(2) */ | 
 
 
 
 
 
 | 173 | 
 #define M_SQRT1_2q      0.7071067811865475244008443621048490Q  /* 1/sqrt(2) */ | 
 
 
 
 
 
 | 174 | 
  | 
 
 
 
 
 
 | 175 | 
 #define __quadmath_extern_inline \ | 
 
 
 
 
 
 | 176 | 
   extern inline __attribute__ ((__gnu_inline__)) | 
 
 
 
 
 
 | 177 | 
  | 
 
 
 
 
 
 | 178 | 
 __quadmath_extern_inline __float128 | 
 
 
 
 
 
 | 179 | 
 __quadmath_nth (cimagq (__complex128 __z)) | 
 
 
 
 
 
 | 180 | 
 { | 
 
 
 
 
 
 | 181 | 
   return __imag__ __z; | 
 
 
 
 
 
 | 182 | 
 } | 
 
 
 
 
 
 | 183 | 
  | 
 
 
 
 
 
 | 184 | 
 __quadmath_extern_inline __float128 | 
 
 
 
 
 
 | 185 | 
 __quadmath_nth (crealq (__complex128 __z)) | 
 
 
 
 
 
 | 186 | 
 { | 
 
 
 
 
 
 | 187 | 
   return __real__ __z; | 
 
 
 
 
 
 | 188 | 
 } | 
 
 
 
 
 
 | 189 | 
  | 
 
 
 
 
 
 | 190 | 
 __quadmath_extern_inline __complex128 | 
 
 
 
 
 
 | 191 | 
 __quadmath_nth (conjq (__complex128 __z)) | 
 
 
 
 
 
 | 192 | 
 { | 
 
 
 
 
 
 | 193 | 
   return __extension__ ~__z; | 
 
 
 
 
 
 | 194 | 
 } | 
 
 
 
 
 
 | 195 | 
  | 
 
 
 
 
 
 | 196 | 
 #ifdef __cplusplus | 
 
 
 
 
 
 | 197 | 
 } | 
 
 
 
 
 
 | 198 | 
 #endif | 
 
 
 
 
 
 | 199 | 
  | 
 
 
 
 
 
 | 200 | 
 #endif |