Update mfit for noparams branch of lmfit
Add add constrain p1_center < p2_center in factory_two_aym_gaussians