Moore–Penrose pseudoinverse?

[Jeroen Demeyer]

Does anybody know of an implementation of the Moore–Penrose
pseudoinverse (a generalised inverse for matrices) in Sage? I couldn't
find anything, but just in case I RTFM badly or it's known under some
different name, I ask here.

Jeroen.

[William Stein]

Google yields:

http://docs.scipy.org/doc/numpy-1.10.0/reference/generated/numpy.linalg.pinv.html


>
> Jeroen.
>
> --
> You received this message because you are subscribed to the Google Groups
> "sage-support" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to sage-support...@googlegroups.com.
> To post to this group, send email to sage-s...@googlegroups.com.
> Visit this group at https://groups.google.com/group/sage-support.
> For more options, visit https://groups.google.com/d/optout.



--
William (http://wstein.org)

[Jeroen Demeyer]

On 2016-03-24 15:58, William Stein wrote:
> http://docs.scipy.org/doc/numpy-1.10.0/reference/generated/numpy.linalg.pinv.html

Numpy works fine for RR or CC but not so much for other fields.

[Michael Orlitzky]

I looked once but didn't find anything. We're missing the SVD for rings
like RR, or it would be easy to fix.

[Jeroen Demeyer]

I want the pseudo-inverse over any base field. There is no reason to
restrict to RR. The SVD on the other hand is not defined over any field.

Anyway, I need the pseudo-inverse, so I'll have to implement it myself.

[Dima Pasechnik]

I recall a reference of this that even allows parameters in the matrix:

https://zbmath.org/1097.15006

[Jeroen Demeyer]

This is now implemented and ready for review at

http://trac.sagemath.org/ticket/20287