This version only works on non-singular matrices.
----- Original Message ----- From: Ric Sherlock <[email protected]> Date: Saturday, March 12, 2011 18:26 Subject: [Jprogramming] LU decomposition essay To: Programming forum <[email protected]> > I was investigating the LU Decomposition essay on the J wiki ( > http://www.jsoftware.com/jwiki/Essays/LU%20Decomposition ) > > The verb presented there seems to fail (index error) for > singular matricies. > LU 1 1 1 , 2 2 2 ,: 3 3 3 > |index error: LU > | p;(=1);p {"1 A[p=.C.(n- > 1);~.0,(0~:,A)i.1 > My understanding is that LU decomposition should work with any square > matrix (singular or not). > Is this a known downside of this particular algorithm? ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
