Thats a unique idea.You should consider the following idea to enhance your idea so that it would impact uniquely: *In this note  propose a multigrid approach to the solution of (multilevel) banded circulant linear system. In particular,  discuss how to define a coarse-grid projector such that the projected matrix at lor dimension preserves the circulant structure. This approach naturally leads to an optimal algorithm having linear cost as the size N of the system and so improving the the classical one based on Fast Fourier Transforms (FFTs) that costs O(N log N) arithmetic operations (ops). It's worth mentioning that these banded circulants are used as preconditioners for elliptic and parabolic PDEs (with Dirichlet or periodic boundary conditions) and for some 2D image restoration problems where the point spread function (PSF) is numerically banded. Therefore the use of the proposed multigrid technique reduces the overall cost from O(k(e,n)N log N) to O(k(e,n)N), where k(e,n) is the number of Pre-conditioned Conjugate Gradient (PCG) iterations to reach the solution within a given accuracy of e. The full analysis of convergence and the related numerical experiments are reported in a forthcoming paper [18].
['Stefano Serra Capizzano', 'Cristina Tablino Possio']
I am writing a research paper on a idea Preliminary remarks on multigrid methods for circulant matrices
[{"from": "human", "value": "I am writing a research paper on a idea Preliminary remarks on multigrid methods for circulant matrices"}, {"from": "gpt", "value": "Thats a unique idea.You should consider the following idea to enhance your idea so that it would impact uniquely: *In this note  propose a multigrid approach to the solution of (multilevel) banded circulant linear system. In particular,  discuss how to define a coarse-grid projector such that the projected matrix at lor dimension preserves the circulant structure. This approach naturally leads to an optimal algorithm having linear cost as the size N of the system and so improving the the classical one based on Fast Fourier Transforms (FFTs) that costs O(N log N) arithmetic operations (ops). It's worth mentioning that these banded circulants are used as preconditioners for elliptic and parabolic PDEs (with Dirichlet or periodic boundary conditions) and for some 2D image restoration problems where the point spread function (PSF) is numerically banded. Therefore the use of the proposed multigrid technique reduces the overall cost from O(k(e,n)N log N) to O(k(e,n)N), where k(e,n) is the number of Pre-conditioned Conjugate Gradient (PCG) iterations to reach the solution within a given accuracy of e. The full analysis of convergence and the related numerical experiments are reported in a forthcoming paper [18]."}]