Re: lu decomposition without any permutations

[Kamaraju S Kusumanchi]

David Bateman wrote:

> The underlying lapack routine xGETRF doesn't allow an LU factorization
> to be performed without pivoting, and so with the existing Octave lu
> function this is not possible.. That being said, its fairly easy to
> write an LU decomposition yourself; search for "Crout" or "Gauss Jordan"
> and you'll find the algorithms and maybe even some simple code (ie about
> 50 lines) to do it.. Note that the stability of the decomposition will
> be severely impacted by the fact that you don't allow pivoting..

Thanks for all the answers. I was just being lazy to write the routine
myself. I was so confident that octave will have one of its own way of
doing lu decomposition without pivoting.

raju
-- 
Kamaraju S Kusumanchi
http://www.people.cornell.edu/pages/kk288/
http://malayamaarutham.blogspot.com/