Matrix approximations and solvers using tensor products and non-standard wavelet transforms related to irregular grids

J. M. Ford, I. V. Oseledets, E. E. Tyrtyshnikov

Результат исследований: Вклад в журналСтатьярецензирование

11 Цитирования (Scopus)

Аннотация

Dense large-scale matrices coming from integral equations and tensor-product grids can be approximated by a sum of Kronecker products with further sparsification of the factors via discrete wavelet transforms, which results in reduced storage and computational costs and also in good preconditioners in the case of uniform one-dimensional grids. However, irregular grids lead to a loss of approximation quality and, more significantly, to a severe deterioration in efficiency of the preconditioners that have been considered previously (using a sparsification of the inverse to one Kronecker product or an incomplete factorization approach). In this paper we propose to use non-standard wavelet transforms related to the irregular grids involved and, using numerical examples, we show that the new transforms provide better compression than the Daubechies wavelets. A further innovation is a scaled two-level circulant preconditioner that performs well on irregular grids. The proposed approximation and preconditioning techniques have been applied to a hypersingular integral equation modelling flow around a thin aerofoil and made it possible to solve linear systems with more than 1 million unknowns in 15-20 minutes even on a personal computer.

Язык оригиналаАнглийский
Страницы (с-по)185-204
Число страниц20
ЖурналRussian Journal of Numerical Analysis and Mathematical Modelling
Том19
Номер выпуска2
DOI
СостояниеОпубликовано - 2004
Опубликовано для внешнего пользованияДа

Fingerprint

Подробные сведения о темах исследования «Matrix approximations and solvers using tensor products and non-standard wavelet transforms related to irregular grids». Вместе они формируют уникальный семантический отпечаток (fingerprint).

Цитировать