The incomplete Cholesky decomposition is known as an excellent smoother in a multigrid iteration and as a preconditioner for the conjugate gradient method. However, the existence of the decomposition is only ensured if the system matrix is an M-matrix. It is well-known that finite element methods usually do not lead to M-matrices. In contrast to this restricting fact, numerical experiments show that, even in cases where the system matrix is not an M-matrix the behaviour of the incomplete Cholesky decomposition apparently does not depend on the structure of the grid. In this paper the behaviour of the method is investigated theoretically for a model problem, where the M-matrix condition is violated systematically by a suitable perturbation. It is shown that in this example the stability of the incomplete Cholesky decomposition is independent of the perturbation and that the analysis of the smoothing property can be carried through. This can be considered as a generalization of the results for the so called square-grid triangulation, as has been established by Wittum in [12] and [11].

Sauter, S (1995). *On the stability of the incomplete Cholesky decomposition for a singular perturbed problem, where the coefficient matrix is not an M-matrix.* Numerical Linear Algebra with Applications, 2(1):17-28.

## Abstract

The incomplete Cholesky decomposition is known as an excellent smoother in a multigrid iteration and as a preconditioner for the conjugate gradient method. However, the existence of the decomposition is only ensured if the system matrix is an M-matrix. It is well-known that finite element methods usually do not lead to M-matrices. In contrast to this restricting fact, numerical experiments show that, even in cases where the system matrix is not an M-matrix the behaviour of the incomplete Cholesky decomposition apparently does not depend on the structure of the grid. In this paper the behaviour of the method is investigated theoretically for a model problem, where the M-matrix condition is violated systematically by a suitable perturbation. It is shown that in this example the stability of the incomplete Cholesky decomposition is independent of the perturbation and that the analysis of the smoothing property can be carried through. This can be considered as a generalization of the results for the so called square-grid triangulation, as has been established by Wittum in [12] and [11].

## Citations

## Altmetrics

## Additional indexing

Item Type: | Journal Article, refereed, original work |
---|---|

Communities & Collections: | 07 Faculty of Science > Institute of Mathematics |

Dewey Decimal Classification: | 510 Mathematics |

Language: | English |

Date: | 1995 |

Deposited On: | 29 Nov 2010 16:28 |

Last Modified: | 05 Apr 2016 13:27 |

Publisher: | Wiley-Blackwell |

ISSN: | 1070-5325 |

Publisher DOI: | https://doi.org/10.1002/nla.1680020103 |

Related URLs: | http://www.zentralblatt-math.org/zbmath/search/?q=an%3A0826.65101 |

## Download

Full text not available from this repository.View at publisher

TrendTerms displays relevant terms of the abstract of this publication and related documents on a map. The terms and their relations were extracted from ZORA using word statistics. Their timelines are taken from ZORA as well. The bubble size of a term is proportional to the number of documents where the term occurs. Red, orange, yellow and green colors are used for terms that occur in the current document; red indicates high interlinkedness of a term with other terms, orange, yellow and green decreasing interlinkedness. Blue is used for terms that have a relation with the terms in this document, but occur in other documents.

You can navigate and zoom the map. Mouse-hovering a term displays its timeline, clicking it yields the associated documents.