How to find a generalized inverse of a vector in sagemath?

daviddgl
55 ●4 ●4 ●11

Given a vector $u$ of size $n\times 1 $. How to find the g-inverse of $u$?

Is there any specific command like

u.ginv()

to find the ginverse or we have define it?? If so, how to do that?

Any hint will help me lot.

Defn:

for a vector $u_{n\times 1} $ a vector $v_{1\times n}$ is said to be generalized inverse of $u$ if $$uvu=u.$$

Sort by » oldest newest most voted

This is sufficient for vectors:

V = matrix(4,1,[1,2,2,1])
V.transpose() / V.norm()^2

2429 ●1 ●22 ●55

Scipy provides methods to calculate a generalized inverse (linalg.pinv , linalg.pinv2, linalg.pinvh)

from scipy import linalg
M = matrix(4,1,(1,2,2,1))
matrix(linalg.pinv2(M))

Please start posting anonymously - your entry will be published after you log in or create a new account.