comp.lang.idl-pvwave archive: archive » Re: mpfit: multivariate fit

Home » Public Forums » archive » Re: mpfit: multivariate fit

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

Re: mpfit: multivariate fit [message #53917]

Tue, 08 May 2007 07:45

David Fanning
Messages: 11724
Registered: August 2001

Senior Member

Dave writes:

> Thanks hradilv, that does help. I also found that I could get it to
> work if I expand the dimensions and make everything a vector prior to
> fitting. I'd be interested in knowing from Craig if your fix is ok as
> I greatly prefer it to the fiddling below.

I don't speak for Craig, but I'm going to guess that Vince's
"fix" is not going to go over well. :-)

Here is the line that breaks:

fjac(0,j) = (fp-fvec)/h(j)

When it breaks, fjac is a [300,9] array. And what
you are trying to stuff into it is a [100,3] array.
But, j is 7, which means you are trying to do something
like:

fjac[0:99, 7:9] = thisThing

Of course, that "9" is what is causing the grief.

I think Vince's fix just makes it possible to stuff
the data into fjac *somewhere*. This seems a dubious
proposition when you are looking for good results. :-)

Cheers,

David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")

Report message to a moderator

[Message index]

		Re: mpfit: multivariate fit By: David Fanning on Tue, 08 May 2007 07:45
		Re: mpfit: multivariate fit By: Dave[3] on Tue, 08 May 2007 07:23
		Re: mpfit: multivariate fit By: Vince Hradil on Tue, 08 May 2007 07:11
		Re: mpfit: multivariate fit By: Craig Markwardt on Fri, 11 May 2007 01:05

Previous Topic:	How to print 100 in exponential format as '10^2'
Next Topic:	Re: How to print 100 in exponential format as '10^2'

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 09:00:51 PDT 2025

Total time taken to generate the page: 1.03996 seconds