def METHOD_NAME(self):
    # test abstract grid scan algorithm
    for ndim in [1, 3]:
        bounds = torch.stack([torch.zeros(ndim), torch.ones(ndim)])
        alg = GridScanAlgorithm()
        mesh = alg.create_mesh(bounds=bounds)
        # benchmark
        n_mesh = alg.n_mesh_points
        linspace_list = [
            torch.linspace(bounds.T[i][0], bounds.T[i][1], n_mesh)
            for i in range(ndim)
        ]
        xx = torch.meshgrid(*linspace_list, indexing="ij")
        benchmark_mesh = torch.stack(xx).flatten(start_dim=1).T
        assert torch.allclose(mesh, benchmark_mesh)
        assert mesh.shape == torch.Size([10**ndim, ndim])