def METHOD_NAME():
    _z, _m = vectors("z", "m")
    actual_fn = function([_z, _m], poch(_z, _m))
    z = random_ranged(0, 5, (2,))
    m = random_ranged(0, 5, (2,))
    actual = actual_fn(z, m)
    expected = scipy_poch(z, m)
    np.testing.assert_allclose(
        actual, expected, rtol=1e-7 if config.floatX == "float64" else 1e-5
    )