GNU Scientific Library - Bugs: bug #32237, RFE: Add function for extracting...
You are not allowed to post comments on this tracker with your current authentication level.
bug #32237: RFE: Add function for extracting nodes and weights for Gauss-Legendre integration
Submitter: | Susi Lehtola <jussilehtola> | ||
Submitted: | Mon 24 Jan 2011 04:51:39 PM UTC | ||
Category: | None | Severity: | 3 - Normal |
Operating System: | Status: | Fixed | |
Assigned to: | rhysu | Open/Closed: | Closed |
Release: |
Mon 14 Feb 2011 07:25:05 PM UTC, comment #1: |
Rhys Ulerich <rhysu> |
Mon 24 Jan 2011 04:51:39 PM UTC, original submission:
It would be useful if there was some function in gsl_integration_glfixed that would return the nodes and weights for the wanted order of GL quadrature.
|
Susi Lehtola <jussilehtola> |
No files currently attached
Depends on the following items: None found
Items that depend on this one: None found
Carbon-Copy List
Follow 4 latest changes.
Date | Changed by | Updated Field | Previous Value | => | Replaced by |
---|---|---|---|---|---|
2011-03-29 | bjg | Open/Closed | Open | Closed | |
2011-02-14 | rhysu | Status | In Progress | Fixed | |
2011-02-14 | rhysu | Status | None | In Progress | |
Assigned to | None | rhysu |
Powered by Savane 3.13-758e.
Corresponding source code
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.