Fast Penalized Generalized Estimating Equations for Large Longitudinal Functional Datasets