diff options
Diffstat (limited to 'Modules')
| -rw-r--r-- | Modules/cmathmodule.c | 20 | ||||
| -rw-r--r-- | Modules/mathmodule.c | 18 |
2 files changed, 35 insertions, 3 deletions
diff --git a/Modules/cmathmodule.c b/Modules/cmathmodule.c index 521d3aa..2cef27c 100644 --- a/Modules/cmathmodule.c +++ b/Modules/cmathmodule.c @@ -21,6 +21,22 @@ #define M_PI (3.141592653589793239) #endif +#ifdef SCO_ATAN2_BUG +/* + * UnixWare 7+ is known to have a bug in atan2 that will return PI instead + * of ZERO (0) if the first argument is ZERO(0). + */ +static double atan2_sco(double x, double y) +{ + if (x == 0.0) + return (double)0.0; + return atan2(x, y); +} +#define ATAN2 atan2_sco +#else +#define ATAN2 atan2 +#endif + /* First, the C functions that do the real work */ /* constants */ @@ -172,7 +188,7 @@ c_log(Py_complex x) { Py_complex r; double l = hypot(x.real,x.imag); - r.imag = atan2(x.imag, x.real); + r.imag = ATAN2(x.imag, x.real); r.real = log(l); return r; } @@ -188,7 +204,7 @@ c_log10(Py_complex x) { Py_complex r; double l = hypot(x.real,x.imag); - r.imag = atan2(x.imag, x.real)/log(10.); + r.imag = ATAN2(x.imag, x.real)/log(10.); r.real = log10(l); return r; } diff --git a/Modules/mathmodule.c b/Modules/mathmodule.c index c206ddc..7f4839a 100644 --- a/Modules/mathmodule.c +++ b/Modules/mathmodule.c @@ -31,6 +31,22 @@ extern double modf (double, double *); #define CHECK(x) /* Don't know how to check */ #endif +#ifdef SCO_ATAN2_BUG +/* + * UnixWare 7+ is known to have a bug in atan2 that will return PI instead + * of ZERO (0) if the first argument is ZERO(0). + */ +static double atan2_sco(double x, double y) +{ + if (x == 0.0) + return (double)0.0; + return atan2(x, y); +} +#define ATAN2 atan2_sco +#else +#define ATAN2 atan2 +#endif + /* Call is_error when errno != 0, and where x is the result libm * returned. is_error will usually set up an exception and return * true (1), but may return false (0) without setting up an exception. @@ -115,7 +131,7 @@ FUNC1(asin, asin, "asin(x)\n\nReturn the arc sine (measured in radians) of x.") FUNC1(atan, atan, "atan(x)\n\nReturn the arc tangent (measured in radians) of x.") -FUNC2(atan2, atan2, +FUNC2(atan2, ATAN2, "atan2(y, x)\n\nReturn the arc tangent (measured in radians) of y/x.\n" "Unlike atan(y/x), the signs of both x and y are considered.") FUNC1(ceil, ceil, |