[SciPy-User] [cython-users] bug report + feature request (was: scipy cblas return value)

[Oleksandr Huziy]

2013/9/26 Sturla Molden <sturla at molden.no>

>
> On Sep 26, 2013, at 3:09 PM, Radim Řehůřek <me at radimrehurek.com> wrote:
>
> > I guess I was too verbose, let me rephrase:
> >
> > 1. BUG REPORT
> > `scipy.linalg.blas.sdot` gives wrong results on mac:
> >
> > >>> scipy.linalg.blas.sdot(np.array([ 10.], dtype=np.float32),
> np.array([ 0.01], dtype=np.float32))
> > -0.0
> >
>
> Let me submit in evidence to the contrary:
>
> In [20]:
> scipy.linalg.blas.sdot(np.array([10.],dtype=np.float32),np.array([0.01],dtype=np.float32))
> Out[20]: 0.09999999403953552
>
> In [21]: sys.platform
> Out[21]: 'darwin'
>
> It is your BLAS library that makes the error, not SciPy.
>
>
> Sturla
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>




Confirming Sturla's evidence

In [9]:
linalg.blas.sdot(np.array([10.],dtype=np.float32),np.array([0.01],dtype=np.float32))
Out[9]: 0.09999999403953552

In [10]: sys.platform
Out[10]: 'darwin'

In [11]: scipy.__version__
Out[11]: '0.12.0'

-- 
Sasha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20130926/0b83205a/attachment.html>