The naive discretization of nonlocal operators leads to matrices with significant density, as compared to classical PDEs. This makes the efficient solution of nonlocal models a challenging task. In this presentation, we will discuss ongoing research into efficient hierarchical matrix assembly and geometric and algebraic multigrid preconditioners that are suitable for nonlocal models.