[Top][All Lists]

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

[Bug-gsl] [bug #32233] Gauss-Legendre routine does not work in gsl-1.14

From: Rhys Ulerich
Subject: [Bug-gsl] [bug #32233] Gauss-Legendre routine does not work in gsl-1.14
Date: Mon, 24 Jan 2011 15:32:40 +0000
User-agent: Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv: Gecko/20101206 Ubuntu/10.04 (lucid) Firefox/3.6.13

Update of bug #32233 (project gsl):

                  Status:                    None => Need Info              
             Assigned to:                    None => rhysu                  


Follow-up Comment #1:

Good morning Jussi,

Three comments:

1) You've filed a bug report on an internal implementation detail-- you'll
note that the GSL manual says nothing definite about the contents of
gsl_integration_glfixed_table.  That was deliberate.

I believe that if you check the code against a 2nd order polynomial it will
behave correctly.  The GSL source code contains exactly such a test within
integration/test.c near a comment saying "Test the fixed-order Gauss-Legendre
rules with a monomial".

2) The weights and points are symmetric about zero leading to only about half
of the points needing to be stored.  That 0.57735 value you show is used
twice, once at 0.57735 and once at -0.57735 each time with weight 1.  This
sums to 2.  The point x = 1 has w = 0, as your table shows, and does not
contribute to such a sum.  But again, this is an implementation detail.

3) If you have a need to access the points and weights, please consider
filing a feature request.  It should be easy to provide an API call that will
return the values that MathWorld lists.

Does this address your concerns about the correctness?  May I close this

- Rhys


Reply to this item at:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]