comp.lang.idl-pvwave archive: archive » mpfit: what's different

Home » Public Forums » archive » mpfit: what's different

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: mpfit: what's different [message #58703 is a reply to message #58600]

Tue, 12 February 2008 09:40

Wox
Messages: 184
Registered: August 2006

Senior Member

Btw, can it be that mpcurvefit expects a different analytical
derivative (dResid/dP) than curvefit (dFunc/dP)? I think this should
be added in mpcurvefit_eval:

;; Compute the deviates, applying the weights
result = (y-f)*wts
if n_params() GT 1 then $
dp*=rebin(wts,n_elements(x),n_elements(p),/sample)

Actually, it should be
dp*=rebin(-wts,n_elements(x),n_elements(p),/sample)
but since mpfit_fdjac2 flips the sign (why?), we have to do it here to
correct for that.

Report message to a moderator

[Message index]

		mpfit: what's different By: Wox on Wed, 06 February 2008 08:35
		Re: mpfit: what's different By: Wox on Mon, 11 February 2008 09:13
		Re: mpfit: what's different By: Craig Markwardt on Tue, 12 February 2008 02:19
		Re: mpfit: what's different By: Craig Markwardt on Mon, 11 February 2008 05:39
		Re: mpfit: what's different By: Wox on Tue, 12 February 2008 09:40

Previous Topic:	Re: One ellipse to rule them all
Next Topic:	Array has a corrupted descriptor

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 14:55:00 PDT 2025

Total time taken to generate the page: 2.27483 seconds