[Numpy-discussion] Bugs using complex192

Matthieu Brucher matthieu.brucher@gmail....
Mon Jan 7 03:10:50 CST 2008


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 <Matts.Bjorck@gmail.com>:
>
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://projects.scipy.org/pipermail/numpy-discussion/attachments/20080107/e1cb67f6/attachment-0001.html 


More information about the Numpy-discussion mailing list