Coarse space (numerical analysis)
- This article deals with a component of numerical methods. For coarse space in topology, see coarse structure.
In numerical analysis, coarse problem is an auxiliary system of equations used in an iterative method for the solution of a given larger system of equations. A coarse problem is basically a version of the same problem at a lower resolution, retaining its essential characteristics, but with fewer variables. The purpose of the coarse problem is to propagate information throughout the whole problem globally.
In multigrid methods for partial differential equations, the coarse problem is typically obtained as a discretization of the same equation on a coarser grid (usually, in finite difference methods) or by a Galerkin approximation on a subspace, called a coarse space. In finite element methods, the Galerkin approximation is typically used, with the coarse space generated by larger elements on the same domain. Typically, the coarse problem corresponds to a grid that is twice or three times coarser.
In domain decomposition methods, the construction of a coarse problem follows the same principles as in multigrid methods, but the coarser problem has much fewer unknowns, generally only one or just a few unknowns per subdomain or substructure, and the coarse space can be of a quite different type that the original finite element space, e.g. piecewise constants with averaging in balancing domain decomposition or built from energy minimal functions in BDDC. The construction of the coarse problem in FETI is unusual in that it is not obtained as a Galerkin approximation of the original problem, however.
In Algebraic Multigrid Methods and in iterative aggregation methods in mathematical economics and Markov chains, the coarse problem is generally obtained by the Galerkin approximation on a subspace. In mathematical economics, the coarse problem may be obtained by the aggregation of products or industries into a coarse description with fewer variables. In Markov chains, a coarse Markov chain may be obtained by aggregating states.
The speed of convergence of multigrid and domain decomposition methods for elliptic partial differential equations without a coarse problem deteriorates with decreasing mesh step (or decreasing element size, or increasing number of subdomains or substructures), thus making a coarse problem necessary for a scalable algorithm.
References
- Jan Mandel and Bedrich Sousedik, Coarse space over the ages, Nineteenth International Conference on Domain Decomposition, Springer-Verlag, submitted, 2009. arXiv:0911.5725
- Olof B. Widlund, The Development of Coarse Spaces for Domain Decomposition Algorithms, in: Domain Decomposition Methods in Science and Engineering XVIII, Bercovier, M. and Gander, M.J. and Kornhuber, R. and Widlund, O. (eds.), Lecture Notes in Computational Science and Engineering 70, Springer-Verlag, 2009, Proceedings of 18th International Conference on Domain Decomposition, Jerusalem, Israel, January 2008. article