Normal form decomposition for Gaussian-to-Gaussian superoperators