Efficient evolution of unpolarized and polarized parton distributions with QCD-PEGASUS

using the symbolic moment-space solutions on a one-fits-all Mellin inversion contour. User options include the order of the evolution including the next-to-next-to-leading order in the unpolarized case, the type of the evolution including an emulation of brute-force solutions, the evolution with a fixed number n(f) of flavors or in the variable-n(f) scheme, and the evolution with a renormalization scale unequal to the factorization scale. The initial distributions are needed in a form facilitating the computation of the complex Mellin moments.