Computational Physics, Python, Linear Algebra, Math, Mathematics.
result is less accurate than the SVD method and Numpy pinv() uses the SVD. Here is an example from the same matrix A : A_plus_1 = np.linalg.inv(A.T.dot(A)).dot(A.T) A_plus_1 array([[ 0.16666667, -0.10606061, 0.03030303], [-0.16666667, 0.28787879, 0.06060606]]) In this case the result is the same as