Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aape_v1_0.gz(10 Kbytes)
Manuscript Title: GISP: Gauss-integration subroutine pack.
Authors: M. Rysavy
Program title: GISP
Catalogue identifier: AAPE_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 32(1984)191
Programming language: Fortran.
Computer: ES 1040.
Operating system: OS/21.8F - MVT.
RAM: 48K words
Word size: 8
Keywords: General purpose, Numerical, Quadrature, Orthogonal polynomials, Gauss method.
Classification: 4.11.

Nature of problem:
A calculation of definite integrals appears in almost all physical problems. In many cases the integrals have to be evaluated numerically. The present pack enables one to employ the Gauss integration method with an optional number of nodes for four different weight functions, i.e. for four different types of integrals.

Solution method:
The integration may be based on the Legendre, Jacobi, Hermite, or Laguerre polynomials. The nodes of the required polynomial are roughly found by a direct scan of a suitable interval (which is automatically chosen by the program) and then corrected by the interval-halving method. The corresponding Gaussian weights are determined and the integral is evaluated. Since the determination of the nodes and weights is time-consuming, the program provides the possibility to skip their evaluation if they are already known (e.g. from the previous call).

Restrictions:
In principle, there are no restrictions, but the program was tested up to 250 nodes only.

Running time:
The running time depends on several factors and more detailed information is provided in the sect. 7 of the long write-up. The test run took about 5.5 min.