def METHOD_NAME(self):
    from pyamg.classical.interpolate import injection_interpolation
    A = poisson((5,), format='csr')
    splitting = np.array([1, 0, 1, 0, 1], dtype='intc')
    P = injection_interpolation(A, splitting)
    P_exact = np.array([[1., 0., 0.],
                        [0., 0., 0.],
                        [0., 1., 0.],
                        [0., 0., 0.],
                        [0., 0., 1.]])
    assert_array_almost_equal(P.toarray(), P_exact)