Practical solutions to the evolution equations for double-parton distributions (DPDs) that are accurate but feasible in terms of computing resources are a challenging task in the study of double-parton scattering phenomenology in perturbative QCD. Due to the higher dimensionality of DPDs with respect to single PDFs, a simple extension of the available methods for PDF evolution is computationally prohibitive. We present an efficient method to implement the fast solution of DGLAP equations and

more »

... erpolation for PDFs and DPDs, based on the use of Chebyshev polynomials. For PDF evolution, our method allows for a higher numerical accuracy using a considerably smaller number of grid points with respect to existing methods. The DPD evolution is realized using an affordable number of grid points, and importantly allows for two independent renormalization scales for the two partons.