Towards exascale simulations of the magnetic universe