UltrasphericalJacobiLobattoMatrix
Ultraspherical Lobatto matrix for Jacobi polynomial.
Examples
- ܀ Interface
- ️܀ See example
- ↢
INTERFACE
MODULE PURE SUBROUTINE UltrasphericalJacobiLobattoMatrix(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+2
REAL(DFP), INTENT(OUT) :: E(:)
!! the size should be 1:n+1
REAL(DFP), OPTIONAL, INTENT(OUT) :: alphaCoeff(0:)
REAL(DFP), OPTIONAL, INTENT(OUT) :: betaCoeff(0:)
END SUBROUTINE UltrasphericalJacobiLobattoMatrix
END INTERFACE
program main
use easifemBase
implicit none
INTEGER( I4B ), parameter :: n = 5
REAL( DFP ) :: D(n+2), E(n+2), alphaCoeff(n+2), betaCoeff(n+2)
CALL UltrasphericalJacobiLobattoMatrix(n, 0.5_DFP, D, E, alphaCoeff, betaCoeff)
CALL Display(D, "D: ")
CALL Display(E, "E: ")
CALL Display(alphaCoeff, "alphaCoeff: ")
CALL Display(betaCoeff, "betaCoeff: ")
end program main
See results
results
D:
-------
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
E:
-------
0.57735
0.51640
0.50709
0.50395
0.50252
0.73855
0.00000
alphaCoeff:
------------
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
betaCoeff:
-----------
2.00000
0.33333
0.26667
0.25714
0.25397
0.25253
0.00000