On Tue, Aug 17, 2010 at 2:55 PM, Barton Willis <willisb at unk.edu> wrote:
> I'd guess this is a bug in the LU code, but the mnewton code could use
> improvement.
Well, mnewton eventually tries to compute the LU factors of
matrix([1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1])
and LU-FACTOR barfs on it.
I don't know if LU-FACTOR is generally OK with rank-deficient
matrices or what.
I agree that mnewton could use some work.
> Maybe the bug is related to the (ill-advised) way that the mnewton code
> uses an undeclared array h.
Well, that doesn't seem to have anything to do with it.
I don't see what's wrong with that anyway.
best
Robert Dodier