Re: Yet again, The Sky is Falling! [message #52843 is a reply to message #52842] |
Thu, 08 March 2007 20:52   |
David Fanning
Messages: 11724 Registered: August 2001
|
Senior Member |
|
|
yp writes:
> You can download the code from www.research.plymouth.ac.uk/casix/temp/test_morel_LUT.pro
>
> You will also need the LUT (morel_fq.dat) to run the code.
> www.research.plymouth.ac.uk/casix/temp/MOREL_FQ.DAT
Well, I haven't tracked the actual culprit down yet, but
the calculations in INT_LUT are doing *something* to the
variables that are being passed in by reference. If I
pass everything in by value, then there is no difference:
IDL> int_LUT,412.5D, 1D, 0.0D, 0.0D, 0.0D, 0.03D, result3
IDL> int_LUT,412.5D+0, 1D+0, nul+0, nul+0, nul+0, 0.03D+0, result4
IDL> print,'RESULT 3: ', result3[*]
IDL> print,'RESULT 4: ', result4[*]
RESULT 3: 0.090100000 0.00000000 0.00000000
0.00000000 0.00000000 0.00000000 0.00000000
RESULT 4: 0.090100000 0.00000000 0.00000000
0.00000000 0.00000000 0.00000000 0.00000000
I suspect there is a calculation in there that sets one
of the values very nearly to what you expect it to be,
but not quite.
Still weird, but more believable. :-)
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.")
|
|
|