subroutine ChebyshevQuadrature(N,nodes,weights,quadType)
implicit none
integer,intent(in) :: N
real(prec),intent(out) :: nodes(0:N)
real(prec),intent(out) :: weights(0:N)
integer,intent(in) :: QuadType
! Local
real(real64) :: nodesLocal(0:N)
real(real64) :: weightsLocal(0:N)
integer :: i
if(QuadType == CHEBYSHEV_GAUSS_LOBATTO) then
call ChebyshevGaussLobatto(N,nodesLocal,weightsLocal)
elseif(QuadType == CHEBYSHEV_GAUSS) then
call ChebyshevGauss(N,nodesLocal,weightsLocal)
endif
do i = 0,N
nodes(i) = real(nodesLocal(i),prec)
weights(i) = real(weightsLocal(i),prec)
enddo
endsubroutine ChebyshevQuadrature