| NaN & mpfit question.. [message #51114] |
Sun, 29 October 2006 16:41 |
D.Kochman@gmail.com
Messages: 7
Registered: October 2006
|
Junior Member |
|
|
I understand this is fairly well covered territory, however, I couldn't
find a solution by searching through the archives.
Anyways, as I call mpfit.pro, it goes into an infinite loop where I can
only break it using the ctrl+break. I was curious as to what was
causing this and it just keeps cycling with the same number values for
the parameters, and Chi-Squared=NaN
So this got me thinking that somehow I was giving mpfit an array with
NaN, but I can't figure out where. I even tried a hackish type thing
to force NaNs to be zero, but still no luck.
Heres the input of the data:
if(sFile ne "")then $
begin
z=dblarr(2,1024)
dg.path=outpath
header=read_binary(sFile,DATA_TYPE=1,DATA_DIMS=32,DATA_START =0)
result_load=read_binary(sFile,DATA_TYPE=3,DATA_START=32)
result_truncated=result_load[dg.starttimeval:dg.starttimeval +1023]
result_reversed=REVERSE(result_truncated)
z(0,0:*)=result_reversed
result=total(z,1,/nan)
chn_data=result
dg.file1=chn_data
and some more stuff
I was hoping to force 0's onto the situation by totalling an array with
0's. But in this situation, I still get Chi-Square = NaN? So I'm
guessing its not with the data, but with something else.
So is there something else that I'm missing? dg.file1 is the one that
ends up getting returned to mpfit.pro at the (y-f)/err
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
