|
subroutine | gaulegh (m, x, w) |
| This Gauss-Legendre quadrature integrates exactly any even polynomial up to degree m*4-2 in the half-interval [0,1]. More...
|
|
Definition at line 58 of file pesg.f90.
◆ gaulegh()
subroutine pesg::gaulegh::gaulegh |
( |
integer(spi), intent(in) |
m, |
|
|
real(dp), dimension(m), intent(out) |
x, |
|
|
real(dp), dimension(m), intent(out) |
w |
|
) |
| |
|
private |
This Gauss-Legendre quadrature integrates exactly any even polynomial up to degree m*4-2 in the half-interval [0,1].
This code is liberally adapted from the algorithm given in Press et al., Numerical Recipes.
- Parameters
-
m | number of nodes in half-interval |
x | nodes and weights |
w | nodes and weights |
- Author
- R. J. Purser
Definition at line 1852 of file pesg.f90.
The documentation for this interface was generated from the following file:
- /scratch1/NCEPDEV/da/George.Gayno/ufs_utils.git/UFS_UTILS.upstream/sorc/grid_tools.fd/regional_esg_grid.fd/pesg.f90