QUADRATURE_RULES_HERMITE_PHYSICIST is a dataset directory which contains examples of quadrature rules of Gauss-Hermite type using the physicist weight.
Gauss-Hermite quadrature rules are designed to approximate integrals on the infinite interval (-oo,+oo).
The Gauss Hermite quadrature assumes that the integrand we are considering has a form like:
Integral ( -oo < x < +oo ) w(x) * f(x) dxwhere the factor w(x) is regarded as a weight factor.
There are three common variations of the rule, depending on the form of the weight factor w(x):
Integral ( -oo < x < +oo ) f(x) dx
Integral ( -oo < x < +oo ) exp(-x*x) f(x) dx
Integral ( -oo < x < +oo ) exp(-x*x/2) f(x) dx
The corresponding Gauss-Hermite rule that uses order points will approximate the integral by
sum ( 1 <= i <= order ) w(i) * f(x(i))where, confusingly, w(i) is a vector of quadrature weights, which has no connection with the w(x) weight function.
For this directory, a quadrature rule is stored as three files, containing the weights, the points, and a file containing two points defining the endpoints of the region. Since the Hermite rules are defined on an infinite region, we set the endpoints to a very large negative and positive values respectively, and hope the program will understand what we mean.
If a Gauss-Hermite rule of order n is computed and evaluated with no numerical error, then it should produce the exact integral of any (weighted) polynomial of degree 2*n-1 or less.
As a simple check, for the physicist rule, the weights should add to sqrt(pi) = 1.7724538509055160273
We consider a physicist weighted Gauss-Hermite quadrature rule of order 3.
Here is the text of the "W" file storing the weights of such a rule:
0.2954089751509195 1.181635900603677 0.2954089751509196
Here is the text of the "X" file storing the abscissas of such a rule:
-1.224744871391589 -0.9862844991098402E-16 1.224744871391589
Here is the text of the "R" file storing the lower and upper limits of the region. These are formally -oo and +oo, but here we simply give them as large values.
The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.
HERMITE_POLYNOMIAL, a C++ library which evaluates the physicist's Hermite polynomial, the probabilist's Hermite polynomial, the Hermite function, and related functions.
HERMITE_RULE, a C++ program which can compute and print a Gauss-Hermite quadrature rule.
INT_EXACTNESS_HERMITE, a C++ program which tests the polynomial exactness of Gauss-Hermite quadrature rules.
QUADRATURE_RULES_GEN_HERMITE, a dataset directory which contains quadrature rules for integration on an infinite interval, using a generalized Gauss-Hermite rule.
QUADRATURE_RULES_HERMITE_PROBABILIST, a dataset directory which contains Gauss-Hermite quadrature rules, for integration on the interval (-oo,+oo), with weight function exp(-x^2/2).
QUADRATURE_RULES_HERMITE_UNWEIGHTED, a dataset directory which contains Gauss-Hermite quadrature rules, for integration on the interval (-oo,+oo), with weight function 1.
SPARSE_GRID_HERMITE, a dataset directory which contains M-dimensional Smolyak sparse grids based on the 1D Gauss-Hermite rule;
TEST_INT_HERMITE, a C++ library which defines test integrands for integration over (-oo,+oo).
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 1:
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 3:
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 7:
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 15:
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 31:
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 63:
Gauss-Hermite Quadrature Rule, Physicist Weight, Order 127:
You can go up one level to the DATASETS page.