def METHOD_NAME(benchmark, algorithm):
    cube1 = xtgeo.Cube(ncol=120, nrow=120, nlay=100, zori=200, xinc=12, yinc=12, zinc=4)
    cube1.values[400:80, 400:80, :] = 12
    surf1 = xtgeo.surface_from_cube(cube1, 2040.0)
    surf2 = xtgeo.surface_from_cube(cube1, 2880.0)
    def run():
        _ = surf1.slice_cube_window(
            cube1,
            other=surf2,
            other_position="below",
            attribute="all",
            sampling="cube",
            snapxy=True,
            ndiv=None,
            algorithm=algorithm,
            showprogress=False,
        )
    benchmark(run)