It doesn't work on Windows, either. In [35]: numpy.sqrt(numpy.array([-1.0], dtype=numpy.complex192)) Out[35]: array([0.0+2.9996087e-305j], dtype=complex192)
In [36]: numpy.sqrt(numpy.array([-1.0], dtype=numpy.complex128)) Out[36]: array([ 0.+1.j]) In [37]: numpy.__version__ Out[37]: '1.0.5.dev4567' hth, L. On 1/7/08, Matthieu Brucher <[EMAIL PROTECTED]> wrote: > > i, > > I managed to reproduce your bugs on a FC6 box : > >>> import numpy as n > > >>> n.sqrt(n.array([-1.0],dtype = n.complex192)) > array([0.0+9.2747134e+492j], dtype=complex192) > > >>> n.sqrt(n.array([-1.0],dtype = n.complex128)) > array([ 0.+1.j]) > > >>> x=n.array([0.0+0.0j, 1.0+0.0j], dtype=n.complex192) > >>> x > array([0.0+0.0j, 1.0+0.0j], dtype=complex192) > > >>> -x > array([3.3621031e-4932+-1.0204727e+2057j, 1.6794099e-4932+5.5355029e-4930j], > dtype=complex192) > > >>> n.__version__ > '1.0.5.dev4494' > > Matthieu > > 2008/1/7, Matts Bjoerck < [EMAIL PROTECTED]>: > > > > Hi, > > > > I've started using complex192 for some calculations and came across two > > things > > that seems to be bugs: > > > > In [1]: sqrt(array([-1.0],dtype = complex192)) > > Out[1]: array([0.0+-6.1646549e-4703j], dtype=complex192) > > > > In [2]: sqrt(array([-1.0],dtype = complex128)) > > Out[2]: array([ 0.+1.j]) > > > > In [3]: x > > Out[3]: array([0.0+0.0j, 1.0+0.0j], dtype=complex192) > > > > In [4]: -x > > Out[4]: array([3.3621031e-4932+0.0012454e-5119j, > > 1.6794099e-4932+0.0011717e-5119j], dtype=complex192) > > > > So, sqrt and using a negation does not seem to work. Have anyone else > > came > > across this as well? > > > > using numpy 1.0.3 on ubuntu > > > > Any help appriciated > > > > /Matts > > > > _______________________________________________ > > Numpy-discussion mailing list > > [email protected] > > http://projects.scipy.org/mailman/listinfo/numpy-discussion > > > > > > -- > French PhD student > 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 > [email protected] > http://projects.scipy.org/mailman/listinfo/numpy-discussion > >
_______________________________________________ Numpy-discussion mailing list [email protected] http://projects.scipy.org/mailman/listinfo/numpy-discussion
