Sm1 = np.linalg.pinv(S)