Available in R and Python, the simple analytic scheme NFWdist performs highly efficient and exact sampling of the Navarro, Frenk & White (NFW) profile as a true probability distribution function, with the only variable being the concentration.
https://ui.adsabs.harvard.edu/abs/2018RNAAS...2b..55R