John, I think I have tracked down the problem. I had set the fit parameters to be very narrowly constrained: set_par ("egauss(1).center", 3.0325, 0, 3.03, 3.035); set_par ("egauss(1).sigma", 0.002, 0, 0.001, 0.003); because I knew that the answer would be close to those numbers and didn't want to allow the fitting algorithm to waste time on other parts of parameter space. When I relaxed the constraints: set_par ("egauss(1).center", 3.0325, 0, 2.7, 3.3); set_par ("egauss(1).sigma", 0.002, 0, 0.0, 0.01); and then fit with mpfit / chisqr:sigma=gehrels, everything works beautifully. So now my question is: suppose I really want to do a fit with the original parameter constraints. Do I need to use subplex? Is there a way to make mpfit work? Thanks, Maurice ---- You received this message because you are subscribed to the isis-users list. To unsubscribe, send a message to isis-users-request_at_email.domain.hiddenwith the first line of the message as: unsubscribeReceived on Wed Jul 07 2010 - 20:36:42 EDT
This archive was generated by hypermail 2.3.0 : Fri May 02 2014 - 08:35:46 EDT