返回
Multigrid methods for multilevel circulant matrices
DOI:10.1137/S1064827501388509.png)
摘要
En 中文
We introduce a multigrid technique for the solution of multilevel circulant linear systems whose coefficient matrix has eigenvalues of the form f(x(j)([n])), where f is continuous and independent of n = (n(1),..., n(d)), and x(j)([n]) = 2pij/n = (2pij(1)/n(1),..., 2pij(d)/n(d)), 0 less than or equal to j(r) less than or equal to n(r) - 1. The interest of the proposed technique pertains to the multilevel banded case, where the total cost is optimal, i.e., O(N) arithmetic operations (ops), N = Pi(r=1)(d) n(r), instead of O(N log N) ops arising from the use of FFTs. In fact, multilevel banded circulants are used as preconditioners for elliptic and parabolic PDEs (with Dirichlet or periodic boundary conditions) and for some two-dimensional image restoration problems where the point spread function (PSF) is numerically banded, so that the overall cost is reduced from O(k(epsilon, n)N log N) to O(k(epsilon, n)N), where k(epsilon, n) is the number of PCG iterations to reach the solution within an accuracy of epsilon. Several numerical experiments concerning one-rank regularized circulant discretization of elliptic 2q-differential operators over one-dimensional and two-dimensional square domains with mixed boundary conditions are performed and discussed.
Keyword:
circulant algebra
two-grid and multigrid iterations
multilevel matrices
AI总结
对已上传原文的论文进行重点信息的提取,主要内容包括:简要概述、研究摘要、背景介绍、关键亮点、图文解析、展望与总结。
期刊
IF:
2.6
论文数:
5.1K
被引数:
1.8W
机构
暂无机构信息

