LAPACK Archives

[Lapack] psuedo inverse algorithm


I have recently published an algorithm that computes the Moore-Penrose 
pseudoinverse, and does so in exact arithmetic. It accomplishes this with 2mn^2 
multiplications (dimensions with m>n), thus is comparable to a single matrix 
multiply. I believe your function makes use of singular value decomposition. 
Mine does not achieve the singular values, which is why it works quicker for 
just the pseudoinverse.

Here is the (thus far only) link to my algorithm:

I believe that this could enhance your wonderful software, how may I help this 
happen? Do I need to submit publication into a journal, or maybe write some 
code? (I have my code in python currently, but the algorithm is quite simple 
really and could easily be ported-- at least it would be easy for me to do.)

-------------- next part --------------
An HTML attachment was scrubbed...

<Prev in Thread] Current Thread [Next in Thread>

For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or