Superfast direct inversion of the nonuniform discrete Fourier transform via hierarchically semi-separable least squares
Fuente:
arXiv
Salvato in:
| Autori principali: | , , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2024
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866908465045176320 |
|---|---|
| author | Wilber, Heather Epperly, Ethan N. Barnett, Alex H. |
| author_facet | Wilber, Heather Epperly, Ethan N. Barnett, Alex H. |
| contents | A direct solver is introduced for solving overdetermined linear systems involving nonuniform discrete Fourier transform matrices. Such matrices can be transformed into a Cauchy-like form that has hierarchical low rank structure. The rank structure of this matrix is explained, and it is shown that the ranks of the relevant submatrices grow only logarithmically with the number of columns of the matrix. A fast rank-structured hierarchical approximation method based on this analysis is developed, along with a hierarchical least-squares solver for these and related systems. This result is a direct method for inverting nonuniform discrete transforms with a complexity that is usually nearly linear with respect to the degrees of freedom in the problem.This solver is benchmarked against various iterative and direct solvers in the setting of inverting the one-dimensional type-II (or forward) transform, for a range of condition numbers and problem sizes (up to 4 x 10^6 by 2 x 10^6). These experiments demonstrate that this method is especially useful for large problems with multiple right-hand sides. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2404_13223 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | Superfast direct inversion of the nonuniform discrete Fourier transform via hierarchically semi-separable least squares Wilber, Heather Epperly, Ethan N. Barnett, Alex H. Numerical Analysis 65T50, 65F55, 65F20 A direct solver is introduced for solving overdetermined linear systems involving nonuniform discrete Fourier transform matrices. Such matrices can be transformed into a Cauchy-like form that has hierarchical low rank structure. The rank structure of this matrix is explained, and it is shown that the ranks of the relevant submatrices grow only logarithmically with the number of columns of the matrix. A fast rank-structured hierarchical approximation method based on this analysis is developed, along with a hierarchical least-squares solver for these and related systems. This result is a direct method for inverting nonuniform discrete transforms with a complexity that is usually nearly linear with respect to the degrees of freedom in the problem.This solver is benchmarked against various iterative and direct solvers in the setting of inverting the one-dimensional type-II (or forward) transform, for a range of condition numbers and problem sizes (up to 4 x 10^6 by 2 x 10^6). These experiments demonstrate that this method is especially useful for large problems with multiple right-hand sides. |
| title | Superfast direct inversion of the nonuniform discrete Fourier transform via hierarchically semi-separable least squares |
| topic | Numerical Analysis 65T50, 65F55, 65F20 |
| url | https://arxiv.org/abs/2404.13223 |