On Tue, Apr 15, 2014 at 4:30 AM, Matthew Brett <[email protected]> wrote: > > It looks as though mingw-w64 is at fault, and I was confused (still > am) because of the different behavior with double and a constant: > > #include <stdio.h> > #include <math.h> > > int main() { > double z, i = -0.0; > printf("With double %f=%f, with constant=%f\n", > i, expm1(i), expm1(-0.)); > } > > gives: > > With double -0.000000=0.000000, with constant=-0.000000 > > That was ugly to track down. > > What is the right way to work round this (using the numpy version > instead of the system version I suppose)? >
The right way is to file a bug at mingw and get it fixed at the source. Additionally as this time npymath seems to be better (thats 3 bugs in npymath vs 1 in mingw on my scoreboard) one could use the mingw preprocessor define instead of HAVE_EXP1M to select this function from npymath. _______________________________________________ NumPy-Discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
