Skip to main content

UltrasphericalJacobiMatrix

Ultraspherical matrix for Jacobi polynomials.

Examples​

INTERFACE
  MODULE PURE SUBROUTINE UltrasphericalJacobiMatrix(n, lambda, D, E, &
    &  alphaCoeff, betaCoeff)
    INTEGER(I4B), INTENT(IN) :: n
    !! n should be greater than or equal to 1
    REAL(DFP), INTENT(IN) :: lambda
    !! lambda should be greater than -0.5
    REAL(DFP), INTENT(OUT) :: D(:)
    !! the size should be 1:n
    REAL(DFP), INTENT(OUT) :: E(:)
    !! the size should be 1:n-1
    REAL(DFP), OPTIONAL, INTENT(OUT) :: alphaCoeff(0:)
    !! recurrence coefficient of monic Ultraspherical polynomial, from 0 to n-1
    REAL(DFP), OPTIONAL, INTENT(OUT) :: betaCoeff(0:)
    !! recurrence coefficient of monic Ultraspherical polynomial, from 0 to n-1
  END SUBROUTINE UltrasphericalJacobiMatrix
END INTERFACE