We present a thermal velocity sampling method for calculating Doppler-broadened atomic spectra, which more efficiently reaches convergence than regular weighted sampling. The uses equal-population of the 1D distribution, “inverse transform” cumulative distribution function, and is broadly applicable to normal distributions. also discuss efficiency by eliminating classes, do not significantly co...