[SciPy-user] Incorrect results in linalg.solve and friends

Nils Wagner nwagner at mecha.uni-stuttgart.de
Mon Feb 23 08:54:01 CST 2004


Pearu Peterson schrieb:
> 
> Hi,
> 
> This is heads-up for linalg crunchers. Incorrect results appear when
> the following condition holds:
> 
>   an input array is a *slice* of a complex array.
> 
> > for i in arange(0,n):
> 
> For example, changing
> 
> >    r = R[:,i]
> 
> to `r = R[:,i].copy()` will give correct results.
> 
> >    X[:,i] = linalg.solve(A,r)
> 
> This gives a workaround to avoid this problem until it gets resolved.
> Also, this is not due to some bug in Atlas but rather of copy_ND_array in
> fortranobject.c or of PyArray_CopyFromObject, but most probably of
> copy_ND_array.
> 
> Pearu

Pearu,

Afaik fortranobject.c belongs to your f2py project. Therefore, it would
be very kind
of you, if you put a new F2PY release on your website.

Thanks in advance

                     Nils
 
BTW, the cvs repository of scipy is currently not available.
> 
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-user



More information about the SciPy-user mailing list