A significant obstacle facing the approach to the continuum limit in lattice gauge theory calculations is the phenomena of critical slowing down in the fermion sector. The method of adaptive multigrid ($\alpha$-MG) methods offer a permanent solution to the superlinear growth in the cost of iterative Dirac matrix inversions. An exascale-ready implementation of $\alpha$-MG as a preconditioner is suitable for both the calculation of Dirac propagators for measurements and of force terms in HMC. Each fermion discretization offers unique challenges, and the successful formulation for Wilson-Clover cannot be applied to other discretizations in a black-box fashion. In this talk we will report on progress towards an algorithm for staggered fermions. We will also present benchmarks for an implementation in the QUDA library for GPUs.

Algorithms and Machines