[SciPy-user] Incorrect results in linalg.solve and friends
nwagner at mecha.uni-stuttgart.de
Mon Feb 23 08:54:01 CST 2004
Pearu Peterson schrieb:
> 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
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
BTW, the cvs repository of scipy is currently not available.
> SciPy-user mailing list
> SciPy-user at scipy.net
More information about the SciPy-user