|
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].
|
|
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 ) |
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 1851 of file pesg.f90.
The documentation for this interface was generated from the following file:
- /scratch4/NCEPDEV/nems/Brian.Curtis/git/BrianCurtis-NOAA/UFS_UTILS/v1_14_0/sorc/grid_tools.fd/regional_esg_grid.fd/pesg.f90