[Bug-gsl] [bug #32237] RFE: Add function for extracting nodes and weights for Gauss-Legendre integration

[Rhys Ulerich]

Update of bug #32237 (project gsl):

                  Status:             In Progress => Fixed                  

    _______________________________________________________

Follow-up Comment #1:

I just ran into this need too...

In revision 4711 I implemented gsl_integration_glfixed_point to allow
retrieving Gauss-Legendre points and weights from a
gsl_integration_glfixed_table.  The points are strictly increasing which
should allow you to compare them against other implementations.  Test cases in
integration/test.c check their validity for n = 1, 2, 3, 4, and 5 against some
analytic results as well as ensuring they're scaled correctly for non-[-1,1]
intervals.  Both pass on GCC's -O3.

I'm closing this feature.  Please re-open if you run into any trouble.

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?32237>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/