winsup/cygwin/math/complex_internal.h

   1 /*
   2  This Software is provided under the Zope Public License (ZPL) Version 2.1.
   3
   4  Copyright (c) 2009, 2010 by the mingw-w64 project
   5
   6  See the AUTHORS file for the list of contributors to the mingw-w64 project.
   7
   8  This license has been certified as open source. It has also been designated
   9  as GPL compatible by the Free Software Foundation (FSF).
  10
  11  Redistribution and use in source and binary forms, with or without
  12  modification, are permitted provided that the following conditions are met:
  13
  14    1. Redistributions in source code must retain the accompanying copyright
  15       notice, this list of conditions, and the following disclaimer.
  16    2. Redistributions in binary form must reproduce the accompanying
  17       copyright notice, this list of conditions, and the following disclaimer
  18       in the documentation and/or other materials provided with the
  19       distribution.
  20    3. Names of the copyright holders must not be used to endorse or promote
  21       products derived from this software without prior written permission
  22       from the copyright holders.
  23    4. The right to distribute this software or to use it for any purpose does
  24       not give you the right to use Servicemarks (sm) or Trademarks (tm) of
  25       the copyright holders.  Use of them is covered by separate agreement
  26       with the copyright holders.
  27    5. If any files are modified, you must cause the modified files to carry
  28       prominent notices stating that you changed the files and the date of
  29       any change.
  30
  31  Disclaimer
  32
  33  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS ``AS IS'' AND ANY EXPRESSED
  34  OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
  35  OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO
  36  EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY DIRECT, INDIRECT,
  37  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
  38  LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA,
  39  OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
  40  LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
  41  NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE,
  42  EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  43 */
  44
  45 #ifdef __CYGWIN__
  46 /* Disable __IMPORT when defining __fdlib_version. */
  47 #define _LIBC
  48 #define _GNU_SOURCE
  49 #endif
  50 #include <math.h>
  51 #include <complex.h>
  52
  53 /* Define some PI constants for long double, as they are not defined in math.h  */
  54 #ifndef M_PI_4l
  55 #define M_PI_4l 0.7853981633974483096156608458198757L
  56 #define M_PI_2l 1.5707963267948966192313216916397514L
  57 #define M_PIl   3.1415926535897932384626433832795029L
  58 #endif
  59
  60 /* NAN builtins for gcc, as they are not part of math.h  */
  61 #ifndef NANF
  62 #define NANF __builtin_nanf ("")
  63 #endif
  64 #ifndef NANL
  65 #define NANL __builtin_nanl ("")
  66 #endif
  67
  68 /* Some more helpers.  */
  69 #define M_PI_3_4  (M_PI - M_PI_4)
  70 #define M_PI_3_4l (M_PIl - M_PI_4l)
  71
  72 #if defined(_NEW_COMPLEX_FLOAT)
  73 # define __FLT_TYPE     float
  74 # define __FLT_ABI(N)   N##f
  75 # define __FLT_CST(N)   N##F
  76 # define __FLT_EPSILON  __FLT_EPSILON__
  77 # define __FLT_NAN      NANF
  78 # define __FLT_HUGE_VAL HUGE_VALF
  79 # define __FLT_PI       M_PI
  80 # define __FLT_PI_2     M_PI_2
  81 # define __FLT_PI_4     M_PI_4
  82 # define __FLT_PI_3_4   M_PI_3_4
  83 # define __FLT_MAXLOG   88.72283905206835F
  84 # define __FLT_MINLOG   -103.278929903431851103F
  85 # define __FLT_LOGE2    0.693147180559945309F
  86 # define __FLT_LOG10E   0.434294481903251828F
  87 # define __FLT_REPORT(NAME) NAME "f"
  88 #elif defined(_NEW_COMPLEX_DOUBLE)
  89 # define __FLT_TYPE     double
  90 # define __FLT_ABI(N)   N
  91 # define __FLT_EPSILON  __DBL_EPSILON__
  92 # define __FLT_CST(N)   N
  93 # define __FLT_NAN      NAN
  94 # define __FLT_HUGE_VAL HUGE_VAL
  95 # define __FLT_PI       M_PI
  96 # define __FLT_PI_2     M_PI_2
  97 # define __FLT_PI_4     M_PI_4
  98 # define __FLT_PI_3_4   M_PI_3_4
  99 # define __FLT_MAXLOG   7.09782712893383996843E2
 100 # define __FLT_MINLOG   -7.45133219101941108420E2
 101 # define __FLT_LOGE2    6.93147180559945309417E-1
 102 # define __FLT_LOG10E   4.34294481903251827651E-1
 103 # define __FLT_REPORT(NAME)     NAME
 104 #elif defined(_NEW_COMPLEX_LDOUBLE)
 105 # define __FLT_TYPE     long double
 106 # define __FLT_ABI(N)   N##l
 107 # define __FLT_CST(N)   N##L
 108 # define __FLT_EPSILON  __LDBL_EPSILON__
 109 # define __FLT_NAN      NANL
 110 # define __FLT_HUGE_VAL HUGE_VALL
 111 # define __FLT_PI       M_PIl
 112 # define __FLT_PI_2     M_PI_2l
 113 # define __FLT_PI_4     M_PI_4l
 114 # define __FLT_PI_3_4   M_PI_3_4l
 115 # define __FLT_MAXLOG   1.1356523406294143949492E4L
 116 # define __FLT_MINLOG   -1.1355137111933024058873E4L
 117 # define __FLT_LOGE2    6.9314718055994530941723E-1L
 118 # define __FLT_LOG10E   4.3429448190325182765113E-1L
 119 # define __FLT_REPORT(NAME) NAME "l"
 120 #else
 121 # error "Unknown complex number type"
 122 #endif