A fast moving least squares approximation with adaptive Lagrangian mesh refinement for large scale immersed boundary simulations