ply_nodeset_legendre Function

public function ply_nodeset_legendre(nPoints) result(x)

Compute Gauss-Legendre integration points on the interval [-1,1].

Arguments

TypeIntentOptionalAttributesName
integer, intent(in) :: nPoints

The number of integration points.

Return Value real(kind=rk)(nPoints)

The coordinates of the Legendre points on the interval [-1,1]. The array has to have the length nPoints.


Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: z1
real(kind=rk), private :: z
real(kind=rk), private :: pp
real(kind=rk), private :: p3
real(kind=rk), private :: p2
real(kind=rk), private :: p1
real(kind=rk), private :: EPS
integer, private :: m
integer, private :: i
integer, private :: j