Robert Kern <robert.kern <at> gmail.com> writes: > > On 2009-06-06 00:34, Ajith Kumar wrote:
> > from numpy import * > > Please ask numpy questions on the numpy mailing list, not here: > > http://www.scipy.org/Mailing_Lists > > > a = arange(1.,10.) > > b = reshape(a, [3,3]) > > c = linalg.inv(b) > > > > And the result is > > > > [[ 3.15221191e+15 -6.30442381e+15 3.15221191e+15] > > [ -6.30442381e+15 1.26088476e+16 -6.30442381e+15] > > [ 3.15221191e+15 -6.30442381e+15 3.15221191e+15]] (gibberish) > > You have a very singular matrix (2*a[1] - a[0] == a[2]). You cannot invert it > numerically and expect sensible results. Is raising an exception (as documented (http://docs.scipy.org/doc/numpy/reference/generated/numpy.linalg.inv.html)) not a "sensible result"? -- http://mail.python.org/mailman/listinfo/python-list