  We present a numerical tool for solving the non-relativistic Kohn-Sham
problem for spherically-symmetric atoms. It treats the Schr\"{o}dinger equation
as an integral equation relying heavily on convolutions. The solver supports
different types of exchange-correlation functionals including screened and
long-range corrected hybrids. We implement a new method for treating range
separation based on the complementary error function kernel. The present tool
is applied in non-relativistic total energy calculations of atoms. A comparison
with ultra-precise reference data[Cinal, JOMC 58, 1571 (2020)] shows a 14-digit
agreement for Hartree-Fock results. We provide further benchmark data obtained
with 5 different exchange-correlation functionals.
