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 > Numpy-discussion@scipy.org > 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 Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion