PRECONDITIONING CHEBYSHEV SPECTRAL COLLOCATION BY FINITE-DIFFERENCE OPERATORS

摘要

In 1979 Orszag proposed a finite-difference preconditioning of the Chebyshev collocation discretization of the Poisson equation. In 1984 Haldenwang, Labrosse, Abboudi, and DeVille gave analytic formulae for the eigenvalues of this preconditioned operator in the one-dimensional case. Experimental results over many years have shown the effectiveness of this procedure and appropriate bounds on the eigenvalues in two dimensions. However, there have been no mathematical proofs describing the behavior of the eigenvalues in two or more dimensions. In this work we consider the generalized field of values (U*(A) over cap(N)U)/(U*LNU), where (A) over cap N is the matrix of the Chebyshev collocation scheme and L-N is the matrix of the finite-difference operator. For the case of the Chebyshev collocation of the Helmoltz operator, Au := -Delta u + au, a greater than or equal to 0 preconditioned by the finite-difference operator associated with the Helmholtz operator Bu := -Delta u + bu, b greater than or equal to 0, we prove that there are two constants 0 < Lambda(0) < Lambda(1) depending only on a(x, y) and b(x, y), but not on N, such that Re{(U*(A) over cap(N)U)/(U*LNU)} greater than or equal to Lambda(0) > 0 and (U*(A) over cap(N)U)/(U*LNU) less than or equal to Lambda(1). These results extend to higher dimensions and to bounds on L-N(-1)(A) over cap(N)(1,w), and A(N)(-1)L(N)(1,w), in the general case where Au := -Delta u + a(1)u(x) + a(2)u(y) + au. [References: 10]