[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: better atan function? (sheepishly " ... nevermind ...")

<bjackel@phys.ucalgary.ca> wrote in message
> Hi Bob
>  (...snip...)
> What made you suspect problems with atan?
> Brian

Firstly, thanks for the response.
Secondly, ARGH!

I had a "float()" cast buried deep in an "argument()" function
that was lopping off my precision. To quote the immortal bard "Do'h".
There is nothing like posting a question to the entire world to make you
discover an obvious and stupid error in your code.

This error was hidden by the fact that you can do the following in IDL.
IDL> print,double(5.12345679)/double(5.12345678),format='(f40.30)'

where, the numbers are "read in" as single precision and then
cast to double precision.  I misread this as performing the action
on two double precision numbers.  One should do the following:
IDL> print,(5.12345679d)/(5.12345678d),format='(f40.30)'

Liam, thanks for the link.

bob the currently sober

PS Have you checked out www.Napster.com?  Wow eh? You can
meaningfully use the term "terrabytes of fun".