Salamin-Brent algorithm



Net researching algorithms for calculating pi, I found the wikipedia
page of that title and learned of the Salamin-Brent algorithm (I'm new
to this particular mathematical esoterica). Therein it states that
the algorithm is O(N log(N) log(log(N))) (i.e., time to calculate N
digits of pi is proportional to N log(N) log(log(N))). My question
is, is the log base the number system base, e.g., the time to
calculate the binary expression of pi is O(N log2(N) log2(log2(N))),
whereas the time to compute the hexadecimal expression is O(N log16(N)
log16(log16(N))), etc. (neglecting any benefits the binary calculation
might gain from taking place on a binary digital computer)?

Oh, and while I have someone's attention, what is done about the
"floating point epsilon" problem, i.e., when the precision of the
calculation reaches the point where it would exceed the representable
precision of the implemented floating point standard? Is it as simple
as just adequately scaling all the numbers involved? Or is there more
to it than that?

Thanks!

DG

.



Relevant Pages

  • Re: Floating point arithmetic question
    ... If I do this calculation in floating point arithmetic (of some specific ... The precison of the 'side with respect to a given plane' test result ... precision calculation will give you wrong answer. ...
    (comp.graphics.algorithms)
  • Re: Microcontroller choice for research project
    ... effort when you're fiddling with the algorithm). ... getting tight for simple PID if you add on the usual saturation, ... I found that the greatest boon conferred by floating point is freedom ... number of bits you need for calculation and the number of bits ...
    (sci.engr.control)
  • Re: General Opinion on a how to?
    ... So different prize ... not have to change the database as you add new algorithm to the system. ... The way i see my calculation working with the aid of Jameys tips is this: ... Load into a prize pool object instance for that game ...
    (microsoft.public.dotnet.languages.csharp)
  • Re: Bug in % (Float)?
    ... floating point number and it ends up slightly larger. ... calculation is taken into account, ... with the algorithm used. ... The discontinuity of % is certainly at the root of the problem, ...
    (comp.lang.ruby)
  • Feature suggestion: sum() ought to use a compensated summation algorithm
    ... I did the following calculation: Generated a list of a million random numbers between 0 and 1, constructed a new list by subtracting the mean value from each number, and then calculated the mean again. ... However, I noticed that the simple Python program below gives a result of ~ 10^-14, while an equivalent Mathematica program (also using double precision) gives a result of ~ 10^-17, i.e. three orders of magnitude more precise. ... A little research shows that Mathematica uses a "compensated summation" algorithm. ...
    (comp.lang.python)

Loading