bugGNU 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: 

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.

Rhys Ulerich <rhysu>
Group administrator
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.

Although gsl_integration_glfixed can be used to calculate simple one-dimensional integrals, in more dimensions (e.g. integration on the surface of the sphere) it's often a lot easier and more efficient not to use the GSL API.

Susi Lehtola <jussilehtola>

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by rhysu (Updated the item)
  • -email is unavailable- added by jussilehtola (Submitted the item)
  •  

    Follow 4 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2011-03-29 bjg Open/ClosedOpen Closed
    2011-02-14 rhysu StatusIn Progress Fixed
    2011-02-14 rhysu StatusNone In Progress
        Assigned toNone rhysu

    Back to the top

    Powered by Savane 3.13-758e.
    Corresponding source code