Dynamical potentials appear in many advanced electronic-structure methods, including self-energies from many-body perturbation theory, dynamical mean-field theory, electronic-transport formulations, and many embedding approaches. Here, we propose a novel treatment for the frequency dependence, introducing an algorithmic inversion method that can be applied to dynamical potentials expanded as sum-over-poles. This approach allows for an exact solution of Dyson-like equations at all frequencies via a mapping to a matrix diagonalization, and provides simultaneously frequency-dependent (spectral) and frequency-integrated (thermodynamic) properties of the Dyson-inverted propagators. The transformation to a sum-over-poles is performed introducing n-th order generalized Lorentzians as an improved basis set to represent the spectral function of a propagator. Numerical results for the homogeneous electron gas at the G0W0 level are provided to argue for the accuracy and efficiency of such unified approach. In this record, we provide the input and outputs of the AGWX suite (see article) used in the work. Also, we add the plots of the obtained Green's function for all densities, k-points, and frequencies.