Fwd: Errors in LUDecomposition class?

[boisvert@nist.gov]

Sender: Peter Coxhead <p.coxhead@cs.bham.ac.uk>
Subject: Errors in LUDecomposition class?


Dear Colleague(s),

In the LUDecomposition class, the getU() and getL() methods assume  
that m >=n, always returning L as m-by-n and U as n-by-n.  However, it  
appears from the header documentation for the class that it is  
intended to work with matrices with m < n as well.

Fix:
* the dimensions of X need to be set based on whether m >= n OR m < n
* the loops then run over X.m and X.n

Regards,
Peter


-------------------------------------
Dr Peter Coxhead
Senior Lecturer, College of Engineering & Physical Sciences
c/o School of Computer Science
The University of Birmingham
Edgbaston
Birmingham
B15 2TT

E-mail: p.coxhead@cs.bham.ac.uk
Web: http://www.cs.bham.ac.uk/~pxc/

Tel. +44 (0)121 414 4773
Fax. +44 (0)121 414 4281
-------------------------------------




----- End forwarded message -----