comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » "infinite" nested for
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: "infinite" nested for [message #86645 is a reply to message #86642] Mon, 25 November 2013 10:25 Go to previous messageGo to previous message
Mats Löfdahl is currently offline  Mats Löfdahl
Messages: 263
Registered: January 2012
Senior Member
On 2013-11-25 17:18, Chrisss wrote:
> Hi Mats Löfdahl, thanks for your interest!
> my problem is that my "Im" is a 3D matrix in which every pixel is a curve with 60 elements.

That much I got. :o)

> Inner loop is for these ones: it does the fitting after the exponential curve inversion. I must follow a scholastic approach that tests all of 60 inversion time T0 within inner loop.

In the loop you do

T0=t
yCurr=y
if (t ne 0) then yCurr[0:T0]=-y[0:To]

So you are changing the sign of the first t elements in y. (I'm assuming
you don't have To in your program and just happened to write it instead
of T0.)

What I don't really see is why you'd want to try fitting an exponential
to negative numbers, hence the question about just using abs(). Is there
some noise that can make the measured (I assume) numbers be negative?

I guess I don't know what an "exponential curve inversion" is. Maybe you
should make t0 a parameter that you can use in myfunc instead, so
instead of testing all possible points where you can flip the sign of
your data, you fit your unchanged data to a function that is an
exponential where you flip the sign at some point t0, that is part of
the fit.
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: lamda default value in CURVEFIT function
Next Topic: Secret of Failed IDL Business Plan

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 17:06:51 PDT 2025

Total time taken to generate the page: 0.00481 seconds