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

Submitted by:  Susi Lehtola <jussilehtola>
Submitted on:  Mon 24 Jan 2011 04:51:39 PM UTC  
 
Category: NoneSeverity: 3 - Normal
Operating System: Status: Fixed
Assigned to: Rhys Ulerich <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>
Project AdministratorIn charge of this item.
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>

 

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

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

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 4 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Tue 29 Mar 2011 11:04:29 PM UTCbjgOpen/ClosedOpen=>Closed
    Mon 14 Feb 2011 07:25:05 PM UTCrhysuStatusIn Progress=>Fixed
    Mon 14 Feb 2011 03:15:38 PM UTCrhysuStatusNone=>In Progress
      Assigned toNone=>rhysu

    Back to the top


    Powered by Savane 3.1-cleanup1