[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Documentation Suggestion
From: |
m elmer |
Subject: |
Documentation Suggestion |
Date: |
Wed, 18 Jan 2023 15:01:43 -0600 (CST) |
Howdy GSL Team!
First, I just wanted to thank y'all for maintaining this library. Its C
approach makes its many functionalities incredibly convenient. I only ran into
one small issue.
In the Gauss-Legendre integration
https://www.gnu.org/software/gsl/doc/html/integration.html#gauss-legendre-integration
section, I couldn't find any information on how the
gsl_integration_glfixed_table stored the nodes and weights. Since it wasn't
documented, I assumed that all n nodes and n weights were stored in the x and w
arrays, respectively, and got some very confusing results. After printing n
values from each of the x and w arrays, I discovered that it doesn't store the
negative values, presumably to save memory. I realize that I should have been
using gsl_integration_glfixed_point in the first place, but there was nothing
that said so or warned me against directly accessing
gsl_integration_glfixed_table.
So that others don't make the same mistake I did, would anyone be willing to
add a few sentences to the description of gsl_integration_glfixed_table? My
suggestion is as follows:
"To save memory, this table does not store the negative nodes or weights for
the quadrature, since they are the same absolute value as the positive nodes
and weights, respectively. The non-negative values stored are ordered by
increasing value, and a desired point is easily retrieved for any integration
interval using the function gsl_integration_glfixed_point."
Thank you for considering my suggestion.
-Matt
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- Documentation Suggestion,
m elmer <=