The extended Padé-Laplace method for efficient discretization of linear viscoelastic spectra