Method do

numpy/linalg/tests/test_linalg.py:906–912

(self, a, b, tags)

                LinalgGeneralizedNonsquareTestCase):

    def do(self, a, b, tags):
        a_ginv = linalg.pinv(a)
        # `a @ a_ginv == I` does not hold if a is singular
        dot = matmul
        result = dot(dot(a, a_ginv), a)
        assert_almost_equal(result, a, single_decimal=5, double_decimal=11)
        assert_(consistent_subclass(a_ginv, a))


class TestPinv(PinvCases):