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

Home » Public Forums » archive » Bug in NR_QROMB?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Bug in NR_QROMB? [message #3522] Sun, 05 February 1995 12:49 Go to previous message
mallozzi is currently offline  mallozzi
Messages: 60
Registered: August 1994
Member
Hi all,
I needed to compute an integral, and the integrand itself contained a
second integral, so I wrote a simple program to test this case using
the IDL-supplied NR_QROMB. In order to use NR_QROMB, you must write a
function that contains the integrand, and then you call NR_QROMB with this
function name as a parameter. In my case, within my first function, I
placed another call to QROMB to evaluate the inner integral (this is not a
double integral), the result of which becomes part of the first integrand.

If I make BOTH calls to NR_QROMB in single or double precision, I get the
wrong result. If I make one call single and the other call double, I get
the correct result. Does anyone know what is happening?
[Message index]
 
Read Message
Read Message
Previous Topic: Re: IDL routine for compress raster-8 image? (lossless)
Next Topic: Wanted: IDL code to read DAOPHOT2 PSFs

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

Current Time: Sat Oct 11 07:10:45 PDT 2025

Total time taken to generate the page: 0.96312 seconds