ChebyshevGauss Subroutine

  1. SELF_Quadrature.f90
  2. SELF_Quadrature
  3. ChebyshevGauss

private subroutine ChebyshevGauss(N, nodes, weights)

Arguments

TypeIntentOptionalAttributesName
integer :: N
real(kind=real64) :: nodes(0:N)
real(kind=real64) :: weights(0:N)

Called by

proc~~chebyshevgauss~~CalledByGraph proc~chebyshevgauss ChebyshevGauss proc~chebyshevquadrature ChebyshevQuadrature proc~chebyshevquadrature->proc~chebyshevgauss proc~init_lagrange_t Init_Lagrange_t proc~init_lagrange_t->proc~chebyshevquadrature proc~init_lagrange~2 Init_Lagrange proc~init_lagrange~2->proc~chebyshevquadrature proc~init_lagrange Init_Lagrange proc~init_lagrange->proc~chebyshevquadrature

Contents

Source Code

ChebyshevGauss

Source Code

  subroutine ChebyshevGauss(N,nodes,weights)
    implicit none
    integer       :: N
    real(real64)    :: nodes(0:N)
    real(real64)    :: weights(0:N)
    ! Local
    integer    :: j

    do j = 0,N

      weights(j) = pi/(real(N,real64)+1.0_real64)
      nodes(j) = -cos(pi*(2.0_real64*real(j,real64)+1.0_real64)/(2.0_real64*real(N,real64)+2.0_real64))

    enddo

  endsubroutine ChebyshevGauss