AllStarFit
AllStarFit (ascl:2201.005) is (or was) an experimental R package for source detection, PSF and multi-component galaxy fitting. Specifically, I used it to perform simultaneous fits of multiple stars in a single image to model the PSF, multi-band fits of galaxies, and to fit multi-component discs with offsets along the minor axis to mimic an inclined disc with a finite thickness/vertical scale height. I also experimented with simultaneously fitting the PSF and sources (stars and galaxies), which was an abject failure and is not recommended.
Relatedly, I also provided cmaeshpc, an extension of the https://cran.r-project.org/web/packages/cmaes/index.html that allows for parallel optimization with runtime limits in an HPC setting (where jobs typically are allocated multiple CPU cores and are terminated after a fixed maximum runtime).