Abstract
This paper proposes a new preconditioning scheme for a linear system with a saddle-point structure arising from a hybrid approximation scheme on the sphere, an approximation scheme that combines (local) spherical radial basis functions and (global) spherical polynomials. In principle the resulting linear system can be preconditioned by the block-diagonal preconditioner of Murphy, Golub and Wathen. Making use of a recently derived inf–sup condition and the Brezzi stability and convergence theorem for this approximation scheme, we show that in this context the Schur complement in the above preconditioner is spectrally equivalent to a certain non-constant diagonal matrix. Numerical experiments with a non-uniform distribution of data points support the theoretically proved quality of the new preconditioner.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Aronszajn N.: Theory of reproducing kernels. Trans. Am. Math. Soc. 68, 337–404 (1950)
Benzi M., Golub G.H., Liesen J.: Numerical solution of saddle point problems. Acta Numer. 14, 1–137 (2005)
Brezzi, F.: On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers. R.A.I.R.O, R-2, 129–151 (1974)
Brezzi F.: New applications of mixed finite element methods, pp. 1335–1347. Proceedings on International Congress of Mathematicians, Berkeley (1986)
Elman H.C., Silvester D.J., Wathen A.J.: Finite Elements and Fast Iterative Solvers with Applications in Incompressible Fluid Dynamics, Numerical Mathematics and Scientific Computation. Oxford University Press, Oxford (2005)
von Golitschek M., Light W.: Interpolation by polynomials and radial basis functions on spheres. Constr. Approx. 17, 1–18 (2001)
Le Gia Q.T., Sloan I.H., Tran T.: Overlapping additive Schwarz preconditioners for elliptic PDEs on the unit sphere. Math.Comput. 78, 79–101 (2009)
Murphy M.F., Golub G.H., Wathen A.J.: A note on preconditioning for indefinite linear systems. SIAM J. Sci. Comput. 21, 1969–1972 (2000)
Narcowich F.J., Ward J.D.: Scattered data interpolation on spheres: error estimates and locally supported basis functions. SIAM J. Math. Anal. 33, 1393–1410 (2002)
Paige C.C., Saunders M.A.: Solution of sparse indefinite systems of linear equations. SIAM J. Numer. Anal. 12, 617–629 (1975)
Saff E.B., Kuijlaars A.B.J.: Distributing many points on a sphere. Math. Intell. 19, 5–11 (1997)
Sloan I.H., Sommariva A.: Approximation on the sphere using radial basis functions plus polynomials. Adv. Comput. Math. 29, 147–177 (2008)
Sloan I.H., Wendland H.: Inf–sup condition for spherical polynomials and radial basis functions on spheres. Math. Comput. 78, 1319–1331 (2009)
Verfürth R.: A combined conjugate gradient-multigrid algorithm for the numerical solution of the Stokes problem. IMA J. Numer. Anal. 4, 441–455 (1984)
Wendland H.: Scattered Data Approximation. Cambridge University Press, Cambridge (2005)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Le Gia, Q.T., Sloan, I.H. & Wathen, A.J. Stability and preconditioning for a hybrid approximation on the sphere. Numer. Math. 118, 695–711 (2011). https://doi.org/10.1007/s00211-011-0369-0
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-011-0369-0