diff options
Diffstat (limited to 'winsup/cygwin/math/complex_internal.h')
| -rw-r--r-- | winsup/cygwin/math/complex_internal.h | 153 |
1 files changed, 153 insertions, 0 deletions
diff --git a/winsup/cygwin/math/complex_internal.h b/winsup/cygwin/math/complex_internal.h new file mode 100644 index 000000000..b230b15b8 --- /dev/null +++ b/winsup/cygwin/math/complex_internal.h @@ -0,0 +1,153 @@ +/* + This Software is provided under the Zope Public License (ZPL) Version 2.1. + + Copyright (c) 2009, 2010 by the mingw-w64 project + + See the AUTHORS file for the list of contributors to the mingw-w64 project. + + This license has been certified as open source. It has also been designated + as GPL compatible by the Free Software Foundation (FSF). + + Redistribution and use in source and binary forms, with or without + modification, are permitted provided that the following conditions are met: + + 1. Redistributions in source code must retain the accompanying copyright + notice, this list of conditions, and the following disclaimer. + 2. Redistributions in binary form must reproduce the accompanying + copyright notice, this list of conditions, and the following disclaimer + in the documentation and/or other materials provided with the + distribution. + 3. Names of the copyright holders must not be used to endorse or promote + products derived from this software without prior written permission + from the copyright holders. + 4. The right to distribute this software or to use it for any purpose does + not give you the right to use Servicemarks (sm) or Trademarks (tm) of + the copyright holders. Use of them is covered by separate agreement + with the copyright holders. + 5. If any files are modified, you must cause the modified files to carry + prominent notices stating that you changed the files and the date of + any change. + + Disclaimer + + THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS ``AS IS'' AND ANY EXPRESSED + OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES + OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO + EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY DIRECT, INDIRECT, + INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT + LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, + OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF + LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING + NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, + EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. +*/ + +#ifdef __CYGWIN__ +/* Disable __IMPORT when defining __fdlib_version. */ +#define _COMPILING_NEWLIB +#define _GNU_SOURCE +#endif +#include <math.h> +#include <complex.h> + +/* Define some PI constants for long double, as they are not defined in math.h */ +#ifndef M_PI_4l +#define M_PI_4l 0.7853981633974483096156608458198757L +#define M_PI_2l 1.5707963267948966192313216916397514L +#define M_PIl 3.1415926535897932384626433832795029L +#endif + +/* NAN builtins for gcc, as they are not part of math.h */ +#ifndef NANF +#define NANF __builtin_nanf ("") +#endif +#ifndef NANL +#define NANL __builtin_nanl ("") +#endif + +/* Some more helpers. */ +#define M_PI_3_4 (M_PI - M_PI_4) +#define M_PI_3_4l (M_PIl - M_PI_4l) + +#if defined(_NEW_COMPLEX_FLOAT) +# define __FLT_TYPE float +# define __FLT_ABI(N) N##f +# define __FLT_CST(N) N##F +# define __FLT_EPSILON __FLT_EPSILON__ +# define __FLT_NAN NANF +# define __FLT_HUGE_VAL HUGE_VALF +# define __FLT_PI M_PI +# define __FLT_PI_2 M_PI_2 +# define __FLT_PI_4 M_PI_4 +# define __FLT_PI_3_4 M_PI_3_4 +# define __FLT_MAXLOG 88.72283905206835F +# define __FLT_MINLOG -103.278929903431851103F +# define __FLT_LOGE2 0.693147180559945309F +# define __FLT_LOG10E 0.434294481903251828F +# define __FLT_REPORT(NAME) NAME "f" +#elif defined(_NEW_COMPLEX_DOUBLE) +# define __FLT_TYPE double +# define __FLT_ABI(N) N +# define __FLT_EPSILON __DBL_EPSILON__ +# define __FLT_CST(N) N +# define __FLT_NAN NAN +# define __FLT_HUGE_VAL HUGE_VAL +# define __FLT_PI M_PI +# define __FLT_PI_2 M_PI_2 +# define __FLT_PI_4 M_PI_4 +# define __FLT_PI_3_4 M_PI_3_4 +# define __FLT_MAXLOG 7.09782712893383996843E2 +# define __FLT_MINLOG -7.45133219101941108420E2 +# define __FLT_LOGE2 6.93147180559945309417E-1 +# define __FLT_LOG10E 4.34294481903251827651E-1 +# define __FLT_REPORT(NAME) NAME +#elif defined(_NEW_COMPLEX_LDOUBLE) +# define __FLT_TYPE long double +# define __FLT_ABI(N) N##l +# define __FLT_CST(N) N##L +# define __FLT_EPSILON __LDBL_EPSILON__ +# define __FLT_NAN NANL +# define __FLT_HUGE_VAL HUGE_VALL +# define __FLT_PI M_PIl +# define __FLT_PI_2 M_PI_2l +# define __FLT_PI_4 M_PI_4l +# define __FLT_PI_3_4 M_PI_3_4l +# define __FLT_MAXLOG 1.1356523406294143949492E4L +# define __FLT_MINLOG -1.1355137111933024058873E4L +# define __FLT_LOGE2 6.9314718055994530941723E-1L +# define __FLT_LOG10E 4.3429448190325182765113E-1L +# define __FLT_REPORT(NAME) NAME "l" +#else +# error "Unknown complex number type" +#endif + +#define __FLT_RPT_DOMAIN(NAME, ARG1, ARG2, RSLT) \ + errno = EDOM, \ + __mingw_raise_matherr (_DOMAIN, __FLT_REPORT(NAME), (double) (ARG1), \ + (double) (ARG2), (double) (RSLT)) +#define __FLT_RPT_ERANGE(NAME, ARG1, ARG2, RSLT, OVL) \ + errno = ERANGE, \ + __mingw_raise_matherr (((OVL) ? _OVERFLOW : _UNDERFLOW), \ + __FLT_REPORT(NAME), (double) (ARG1), \ + (double) (ARG2), (double) (RSLT)) + +#ifdef __CYGWIN__ +inline void __attribute__ ((always_inline)) +__mingw_raise_matherr (int typ, const char *name, double a1, double a2, + double rslt) +{ + if (_LIB_VERSION != _POSIX_) + { + struct exception ex; + ex.type = typ; + ex.name = (char*)name; + ex.arg1 = a1; + ex.arg2 = a2; + ex.retval = rslt; + matherr(&ex); + } +} +#define _DOMAIN DOMAIN +#define _OVERFLOW OVERFLOW +#define _UNDERFLOW UNDERFLOW +#endif |