Chebyshev1JacobiRadauMatrix
Radau matrix for Chebyshev1 polynomials.
INTERFACE
MODULE PURE SUBROUTINE Chebyshev1JacobiRadauMatrix(a, n, D, E, alphaCoeff, &
& betaCoeff)
REAL(DFP), INTENT(IN) :: a
!! one of the end of the domain
INTEGER(I4B), INTENT(IN) :: n
!! n should be greater than or equal to 1
REAL(DFP), INTENT(OUT) :: D(:)
!! the size should be 1:n+1
REAL(DFP), INTENT(OUT) :: E(:)
!! the size should be 1:n
REAL(DFP), OPTIONAL, INTENT(OUT) :: alphaCoeff(0:)
REAL(DFP), OPTIONAL, INTENT(OUT) :: betaCoeff(0:)
END SUBROUTINE Chebyshev1JacobiRadauMatrix
END INTERFACE
Examples
- ️܀ See example
- ↢
program main
use easifemBase
implicit none
INTEGER( I4B ), parameter :: n = 5
REAL( DFP ) :: D(n+1), E(n+1), alphaCoeff(n+1), betaCoeff(n+1)
alphaCoeff=0.0_DFP
betaCoeff = 0.0_DFP
CALL Chebyshev1JacobiRadauMatrix(-1.0_DFP, n, 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.50000
E:
-------
0.70711
0.50000
0.50000
0.50000
0.50000
0.00000
alphaCoeff:
------------
0.00000E+00
0.00000E+00
0.00000E+00
0.00000E+00
0.00000E+00
betaCoeff:
-----------
3.14159
0.50000
0.25000
0.25000
0.25000
0.00000