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