Efficient Shift-and-Invert Preconditioning for Multi-GPU Accelerated Density Functional Calculations