Create Gauss-Legendre integration points and weights for one-dimensional integration on the interval [x1,x2].
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | x1 | lower limit of integration interval |
||
real(kind=rk), | intent(in) | :: | x2 | upper limit of integration interval |
||
real(kind=rk), | intent(out) | :: | x(:) | The coordinates of the gauss points on the interval [x1,x2]. The array has the length nIntP. |
||
real(kind=rk), | intent(out) | :: | w(:) | The quadrature weights. The array has the length nIntP. |
||
integer, | intent(in) | :: | nIntP | The number of integration points. |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
real(kind=rk), | private | :: | z1 | some working variables |
|||
real(kind=rk), | private | :: | z | some working variables |
|||
real(kind=rk), | private | :: | xm | some working variables |
|||
real(kind=rk), | private | :: | xl | some working variables |
|||
real(kind=rk), | private | :: | pp | some working variables |
|||
real(kind=rk), | private | :: | p3 | some working variables |
|||
real(kind=rk), | private | :: | p2 | some working variables |
|||
real(kind=rk), | private | :: | p1 | some working variables |
|||
real(kind=rk), | private | :: | EPS | the relative precision of the points |
|||
integer, | private | :: | m | ||||
integer, | private | :: | i | ||||
integer, | private | :: | j |