def METHOD_NAME():
    """
    This test is testing the validity of trapezoidal rule for block structure implementation.
    Here we are using the _trapezoidal for validity check, because it is already tested in the
    above test function.
    Returns
    -------
    """
    blocksize = 30
    ghost_idx = np.array([14, 15])
    input_vector = np.random.randn(3, blocksize)
    correct_vector = np.hstack(
        (
            _trapezoidal(input_vector[..., : ghost_idx[0]]),
            np.zeros((3, 1)),
            _trapezoidal(input_vector[..., ghost_idx[1] + 1 :]),
        )
    )
    test_vector = _trapezoidal_for_block_structure(input_vector, ghost_idx)
    assert_allclose(test_vector, correct_vector, atol=Tolerance.atol())