Scalable smoothing strategies for a geometric multigrid method for the immersed boundary equations