We propose new methods for empirical point spread function characterization and implement them in Julia. We benchmark our PSF fitter against PIFF and PSFex and show that our approach bridges the accuracy and features of PIFF with the speed of PSFex.