def METHOD_NAME(c3d, c2d, params):
    boundary, facet, expect = params
    utility3d.SubFunctionExtractor(c3d, c2d, boundary=boundary, elem_facet=facet).solve()
    assert numpy.allclose(c2d.dat.data_ro[:], expect)