cheb_dctmat2 Subroutine

public subroutine cheb_dctmat2(nx, ny, mat, rwrk)

This subroutine transforms the function values at the CGL points to the coefficients of the Chebyshev series for two variables. The matrix is not symmetric and assumes the nodal coordinates to be in ascending order.

Reference

  1. Canuto et al. (2006), 'Spectral Methods: Fundamentals in Single Domains', Chapter 2, pp. 86, eqn. (2.4.15)

  2. Gheorghiu, C. I. (2007), 'Spectral Methods for Differential Problems', Casa Cartii de Stiinta, Cluj-Napoca, ISBN 978-973-133-099-0 pp. 18, eqn. (1.40, 1.41)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: nx

Degree of polynomial along x direction
integer, intent(in) :: ny

Degree of polynomial along y direction
real(kind=rp), intent(out), dimension(:,:) :: mat

((nx+1)(ny+1),(nx+1)(ny+1)) matrix. Transform matrix.
real(kind=rp), intent(inout), dimension(:), target :: rwrk

Workspace array of minimum size (nx+1) + (ny+1)