Singular-value (and eigenvalue) distribution and Krylov preconditioning of sequences of sampling matrices approximating integral operators

• Published in: Numerical linear algebra with applications 2014-12, Vol.21 (6), p.722-743
• Main Authors: Al-Fhaid, A.S., Serra-Capizzano, S., Sesana, D., Zaka Ullah, M.
• Format: Article
• Language: English
• Subjects: Approximation; Copyrights; eigenvalue distribution; Eigenvalues; integral equation; Linear systems; Mathematical analysis; Mathematical models; Operators; Sampling; sampling matrices; singular value distribution
• ISSN: 1070-5325; 1099-1506
• DOI: 10.1002/nla.1922

SUMMARYLet k( ⋅ , ⋅ ) be a continuous kernel defined on Ω × Ω, Ω compact subset of Rd, d ⩾1, and let us consider the integral operator K̃ from C(Ω) into C(Ω) (C(Ω) set of continuous functions on Ω) defined as the map f(x) → l(x) =∫Ωk(x,y)f(y)dy,x ∈Ω. K̃ is a compact operator and therefore its spectrum forms a bounded sequence having zero as unique accumulation point. Here, we first consider in detail the approximation of K̃ by using rectangle formula in the case where Ω = [0,1], and the step is h = 1 ∕ n. The related linear application can be represented as a matrix An of size n. In accordance with the compact character of the continuous operator, we prove that {An} ∼ σ0 and {An} ∼ λ0, that is, the considered sequence has singular values and eigenvalues clustered at zero. Moreover, the cluster is strong in perfect analogy with the compactness of K̃. Several generalizations are sketched, with special attention to the general case of pure sampling sequences, and few examples and numerical experiments are critically discussed, including the use of GMRES and preconditioned GMRES for large linear systems coming from the numerical approximation of integral equations of the form ((I −K̃)f(t))(x) = g(x),x ∈Ω, with (K̃f(t))(x) =∫Ωk(x,y)f(y)dy and datum g(x). Copyright © 2014 John Wiley & Sons, Ltd.