A \(V\)-cycle multigrid for multilevel matrix algebras: proof of optimality (Q868670)

The convergence rate of a multigrid method for multilevel linear systems is analyzed. The coefficient matrices of these system are generated by a real and nonnegative multivariate polynomial \(f\) and belong to multilevel matrix algebras like circulant, tau, Hartley, or are of Toeplitz type. In the case of matrix algebra linear systems, the authors prove that the convergence rate is independent of the system dimension even in presence of asymptotical ill-conditioning (this happens if \(f\) takes the zero value). More precisely, if the \(d\)-level coefficient matrix has partial dimension \(n_r\) at level \(r,\) with \(r=1,\dots,d,\) then the size of the system is \(N(\mathbf{n})=\prod_{r=1}^dn_r,\;\mathbf{n}=(n_1,\dots, n_d),\) and \(O(N(\mathbf{n}))\) operations are required by the considered \(V\)-cycle multigrid in order to compute the solution within a fixed accuracy. Since the total arithmetic cost is asymptotically equivalent to the one of a matrix-vector product, the proposed method is optimal. Some numerical experiments concerning linear systems arising in 2D and 3D applications are considered and discussed.