[Numpy-discussion] how to use lu decomposition in lapack library?

[Piotr Luszczek]

You would have to be more specific. If you just want to
solve a system of linear equations with matrix 'a' and
right-hand-side 'b':

import numarray.linear_algebra as LA

x = LA.solve_linear_equations(a, b)

Unlike LAPACK, the above will leave your 'a' untouched.
So, if you have another right hand side 'b1' and the
same matrix 'a' you'll have to pay the cost of factorization
all over again.

There is a way around it but I don't know what you really need.

Piotr

meng at are.berkeley.edu wrote:
> Hi there-
> 
> Can someone help me on this? Thanks!
> 
> Best,
> Xiangyi
> 
> Xiangyi Meng
> Department of Agricultural and Resource Economics
> Room 303, Giannini Hall #3310
> University of California, Berkeley
> Berkeley, CA 94720-3310
> Tel:  (510) 643-4124
> Fax: (510) 643-8911
> Email: meng at are.berkeley.edu
>          
> 
> 
> -------------------------------------------------------
> This SF.Net email is sponsored by:
> Power Architecture Resource Center: Free content, downloads, discussions,
> and more. http://solutions.newsforge.com/ibmarch.tmpl
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion