ply_gauleg Subroutine

public subroutine ply_gauleg(x1, x2, x, w, nIntP)

subroutine to create gauss points and weights for one-dimensional integration on the interval [x1,x2].

Arguments

Type IntentOptional 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(inout), allocatable :: x(:)

The coordinates of the gauss points on the interval [-1,1]. The array has the length nIntP.

real(kind=rk), intent(inout), allocatable :: w(:)

The quadrature weights. The array has the length nIntP.

integer, intent(in) :: nIntP

The number of integration points.