  We provide a detailed exposition of our computational framework designed for
the accurate calculation of real-frequency dynamical correlation functions of
the single-impurity Anderson model (AM) in the regime of weak to intermediate
coupling. Using quantum field theory within the Keldysh formalism to directly
access the self-energy and dynamical susceptibilities in real frequencies, as
detailed in our recent publication
(https://doi.org/10.1103/PhysRevB.109.115128), the primary computational
challenge is the full three-dimensional real-frequency dependence of the
four-point vertex. Our codebase provides a fully MPI+OpenMP parallelized
implementation of the functional renormalization group (fRG) and the
self-consistent parquet equations within the parquet approximation. It
leverages vectorization to handle the additional complexity imposed by the
Keldysh formalism, using optimized data structures and highly performant
integration routines. Going beyond the results shown in the previous
publication, the code includes functionality to perform fRG calculations in the
multiloop framework, at arbitrary loop order, including self-consistent
self-energy iterations. Moreover, implementations of various regulators, such
as hybridization, interaction, frequency, and temperature are supplied.
