Indexed by:
Abstract:
Chebyshev polynomial spectral methods are very accurate, but are plagued by the cost and ill-conditioning of dense discretization matrices. Modified schemes, collectively known as “integration sparsification”, have mollified these problems by discretizing the highest derivative as a diagonal matrix. Here, we examine five case studies where the highest derivative diagonalization fails. Nevertheless, we show that Galerkin discretizations do yield banded matrices that retain most of the advantages of “integration sparsification”. Symbolic computer algebra greatly extends the reach of spectral methods. When spectral methods are implemented using exact rational arithmetic, as is possible for small truncation N in Maple, Mathematica and their ilk, roundoff error is irrelevant, and sparsification failure is not worrisome. When the discretization contains a parameter L, symbolic algebra spectral methods return, as answer to an eigenproblem, not discrete numbers but rather a plane algebraic curve defined as the zero set of a bivariate polynomial P(λ,L); the optimal approximations to the eigenvalues λj are in the middle of the straight portions of the zero contours of P(λ;L) where the isolines are parallel to the L axis. © 2018 International Association for Mathematics and Computers in Simulation (IMACS)
Keyword:
Reprint Author's Address:
Email:
Source :
Mathematics and Computers in Simulation
ISSN: 0378-4754
Year: 2019
Volume: 160
Page: 82-102
1 . 6 2
JCR@2019
2 . 4 6 3
JCR@2020
ESI Discipline: ENGINEERING;
ESI HC Threshold:83
JCR Journal Grade:2
CAS Journal Grade:3
Cited Count:
WoS CC Cited Count: 0
SCOPUS Cited Count:
ESI Highly Cited Papers on the List: 0 Unfold All
WanFang Cited Count:
Chinese Cited Count:
30 Days PV: 8
Affiliated Colleges: