def test_eig_inv():
"""Test matrix pseudoinversion with setting smallest eigenvalue to zero."""
# create rank-deficient array
a = np.array([[1., 0., 1.], [0., 1., 0.], [1., 0., 1.]])
# test inversion
a_inv = np.linalg.pinv(a)
a_inv_eig = _eig_inv(a, 2)
assert_almost_equal(a_inv, a_inv_eig)
def test_eig_inv():
"""Test matrix pseudoinversion with setting smallest eigenvalue to zero."""
# create rank-deficient array
a = np.array([[1., 0., 1.], [0., 1., 0.], [1., 0., 1.]])
# test inversion
a_inv = np.linalg.pinv(a)
a_inv_eig = _eig_inv(a, 2)
assert_almost_equal(a_inv, a_inv_eig)