Re: [Matlisp-users] pseudoinverse

[Michael A. K. <ma...@ll...>]

Tunc,

> When M<N, should the matrix have rank=M or does it
> not matter?

1.  Good question.  The answer is, PINV as written will fail or
    produce a numerically unstable result due to taking the reciprocal
    of a very small singular value.

2.  The solution is straight forward, let all singular values < MACH-TOL
    equal zero.  ==> rank(A) = r <= min(M,N).  And use just those 
    r elements to form PINV.

3.  I have an all day meeting today, but will make these changes in the
    next couple of days.  ....BTW what is the name of the "*MACH" in 
    the matlisp core....

more later


mike