Accurate Bernoulli numbers on the 41C, or "how close can you get"?

[Ángel Martin]

(03-19-2014 08:31 PM)Dieter Wrote:  "Fast and accurate" always sounds very promising. ;-) Do you have any figures here? How fast is, for instance, qf(0,01) on an actual HP41?

On V41 using the default settings: qf(0.01) = -2,326522873, in 1.5 seconds.
On a 41CL using TURBO=zero the same calculation takes 2.35 seconds.

and just 0.11002 seconds in TURBO50 ( I couldn't resist!)

(03-19-2014 08:31 PM)Dieter Wrote:  You say you used the inverse error function as implemented in the CUDA library. Do you have a link to this algorithm? I found some information on the library, but not on the actual algorithms.

cf. the paper by Mike Giles posted at:
http://people.maths.ox.ac.uk/gilesm/file...erfinv.pdf

Edited: I checked my files looking for the double-precision source code paper. It appears the link has been re-directed to another subject, so I cannot find it on-line. Unfortunately I didn't keep an electronic copy of the document so would need to go to the actual "paper vault" , I hope to find it or else would have to reverse-engineer my own MCODE

Edited: I found the document - send me your email address if you're interested.