We optimise the non-linear parameters of the shapelet basis (scale and basis center) using the Particle Swarm Optimiser [1]} implemented in Lenstronomy, which is particularly suited for finding global minima in large parameter spacesWe also used the Nelder-Mead simplex algorithm (from the scipy library), but found that it quickly fails at exploring the parameter space once more complex light profiles such as shapelets+Sérsic or large shapelet basis are being optimised, given the constraints of HST images..
