Robust fast direct integral equation solver for three-dimensional quasi-periodic scattering problems with a large number of layers.