Question about the Borwein fast algorithm for certain values of Gamma

  • Context: Graduate 
  • Thread starter Thread starter mesa
  • Start date Start date
  • Tags Tags
    Algorithm Gamma
Click For Summary

Discussion Overview

The discussion revolves around Borwein's algorithm for evaluating the gamma function, particularly for small rational fractions. Participants explore the mathematical notation and components of the algorithm, including the meaning of specific symbols and the limitations of the method.

Discussion Character

  • Technical explanation
  • Conceptual clarification
  • Debate/contested

Main Points Raised

  • One participant questions the role of the curly brackets in Borwein's algorithm, wondering if they represent a multiplier or something else.
  • Another participant discusses the arithmetic geometric mean (AG) method used by Borwein, noting its effectiveness in approximating gamma function values but also its limitations in the range of values it can compute.
  • There is uncertainty regarding the meanings of "\Sigma[1]" and "AG[1]", with participants expressing confusion about these terms.
  • One participant suggests that the curly brackets may simply represent parentheses or brackets, indicating a lack of familiarity with this notation.

Areas of Agreement / Disagreement

Participants express differing views on the interpretation of the notation used in Borwein's algorithm. There is no consensus on the meaning of the curly brackets or the symbols "\Sigma[1]" and "AG[1]". The discussion remains unresolved regarding the clarity of these mathematical expressions.

Contextual Notes

Participants note that the algorithm's applicability is limited to certain values, and the process of calculating the necessary tables for the gamma function can be complex and daunting.

mesa
Gold Member
Messages
694
Reaction score
36
I am reading through Borwein's paper, "Fast evaluation of the gamma function for small rational fractions using complete elliptic integrals of the first kind" and have a question.

If we look at his algorithm's we see they are of this general form:

Gamma(1/2)=2^(-1/4)AG[1]{2^(-1/2)∑[1]}^(-1/2)

I have been able to run through everything up until the curly brackets, {}. Is this a multiplier? Or is it something else entirely?
 
Mathematics news on Phys.org
I'm sorry you are not generating any responses at the moment. Is there any additional information you can share with us? Any new findings?
 
Greg Bernhardt said:
I'm sorry you are not generating any responses at the moment. Is there any additional information you can share with us? Any new findings?

It's the first time this has happened using PF so can't really complain.

Dr. Borwein figured out a way to use the arithmetic geometric mean (which he refers to simply as AG method) to take the difference of the squares of an and bn in a summation to get some pretty astounding decimal approximations for certain values of the gamma function.

The primary issue is it seems to be extremely limited in the number of values that can be used to calculate gammas and with the exception of some extraordinary mathematical gymnastics the number of values that can be used is countable.

Regardless, it is a wonderful piece of work!
 
mesa said:
I am reading through Borwein's paper, "Fast evaluation of the gamma function for small rational fractions using complete elliptic integrals of the first kind" and have a question.

If we look at his algorithm's we see they are of this general form:

Gamma(1/2)=2^(-1/4)AG[1]{2^(-1/2)∑[1]}^(-1/2)

I have been able to run through everything up until the curly brackets, {}. Is this a multiplier? Or is it something else entirely?
I'm not clear on what "\Sigma[1]" or "AG[1]" mean. But ignoring that, because the curly brackets are to the "-1/2" power, this is
\Gamma(1/2)= \frac{1}{16}\frac{AG[1]}{\sqrt{\frac{1}{\sqrt{2}}\Sigma[1]}}
 
HallsofIvy said:
I'm not clear on what "\Sigma[1]" or "AG[1]" mean. But ignoring that, because the curly brackets are to the "-1/2" power, this is
\Gamma(1/2)= \frac{1}{16}\frac{AG[1]}{\sqrt{\frac{1}{\sqrt{2}}\Sigma[1]}}

As far as 'AG' Borwein uses it as short hand for his summation that takes the difference of the squares of an and bn from the arithmetic geometric mean operation inside a summation that is taken away from a value that is dependent on the input for gamma.

The general operation of the identity is based on calculated variable 'N' tables provided by the author to compute certain values of gamma we can 'look up' to insert to get the components of AG[N] and ∑[N]. I did not look into exactly how the author calculated these tables, the process can be daunting and does not yield much value (if someone thinks otherwise please chime in!).

Overall the scope of his identity seems very limited but for the values it can calculate the decimal approximation for accuracy is astounding and grows radically with each step.

On another note, I figured the same for the curly brackets just representing () or [], I just haven't seen them used before, is this a common thing?
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
Replies
8
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 12 ·
Replies
12
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K