[Numpy-discussion] A bug in scipy.linalg.lu_factor?

[Nils Wagner]

On Thu, 25 Jan 2007 16:56:54 -0200
  "Paulo J. S. Silva" <pjssilva at ime.usp.br> wrote:
> Em Qui, 2007-01-25 às 19:46 +0100, Nils Wagner escreveu:
> 
>> 
>> It works if you use
>> M=num.random.rand(2,2)
>> 
>> Nils
>> 
> 
> Yes, it works for arrays but not for matrices. I thought 
>that
> scipy.linalg functions were supposed to work with 
>matrices.
> 
> Paulo
> 
> 
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at scipy.org
> http://projects.scipy.org/mailman/listinfo/numpy-discussion
  
linalg.solve works with matrix input

In [1]: import numpy as num

In [2]: import scipy.linalg as la

In [3]: M = num.mat(num.random.rand(2,2))

In [4]: b = num.random.rand(2)

In [5]: x = la.solve(M,b)

In [6]: print M*x
[[ 0.29508067  0.17152755]]

In [7]: b
Out[7]: array([ 0.29508067,  0.17152755])

Nils