A discrete integral transform for rapid spectral synthesis

•Novel approach for spectral synthesis resulting in typical performance increase of 2~3 orders of magnitude.•Generic approach that can be applied to any spectral code or database.•Arbitrary accuracy: Voigt function is approximated as accurately as needed.•Additional performance increase for line-of-sight spectra.•Easy to parallelize: GPU implementation yields further 2~3 orders improvement.