>>>> np.log(-1 + 0j) > 3.1415926535897931j >>>> np.log(-1 - 1e-99j) > -3.1415926535897931j >>>> np.log(-1 + 1e-99j) > 3.1415926535897931j > > I'd guess this is typically harmless, but may lead to problems if your > code relies on the choice of branch at the branch cut. > > Could you check if this is really the case? If log seems to work OK, then > it's a problem with the test and not the functions.
The result is what you expected. BTW, I can confirm that the latest official MKL does not work with numpy, as it is explained on the Intel forum (http://software.intel.com/en-us/forums/intel-math-kernel-library/topic/60460). I get the i_free not defined issue. Matthieu -- French PhD student Information System Engineer Website: http://matthieu-brucher.developpez.com/ Blogs: http://matt.eifelle.com and http://blog.developpez.com/?blog=92 LinkedIn: http://www.linkedin.com/in/matthieubrucher _______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion