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

Sturla Molden sturla at molden.no
Thu Sep 26 12:37:04 EDT 2013


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




More information about the SciPy-User mailing list