Fast solution of Sylvester-structured systems for spatial source separation of the Cosmic Microwave Background
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2022
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866910620400484352 |
|---|---|
| author | Pham, Dung Soodhalter, Kirk M. Wilson, Simon |
| author_facet | Pham, Dung Soodhalter, Kirk M. Wilson, Simon |
| contents | Implementation of many statistical methods for large, multivariate data sets requires one to solve a linear system that, depending on the method, is of the dimension of the number of observations or each individual data vector. This is often the limiting factor in scaling the method with data size and complexity. In this paper we illustrate the use of Krylov subspace methods to address this issue in a statistical solution to a source separation problem in cosmology where the data size is prohibitively large for direct solution of the required system. Two distinct approaches, adapted from techniques in the literature, are described: one that uses the method of conjugate gradients directly to the Kronecker-structured problem and another that reformulates the system as a Sylvester matrix equation. We show that both approaches produce an accurate solution within an acceptable computation time and with practical memory requirements for the data size that is currently available. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2204_08057 |
| institution | arXiv |
| publishDate | 2022 |
| record_format | arxiv |
| spellingShingle | Fast solution of Sylvester-structured systems for spatial source separation of the Cosmic Microwave Background Pham, Dung Soodhalter, Kirk M. Wilson, Simon Numerical Analysis Instrumentation and Methods for Astrophysics 65F10, 65F50, 65F08 Implementation of many statistical methods for large, multivariate data sets requires one to solve a linear system that, depending on the method, is of the dimension of the number of observations or each individual data vector. This is often the limiting factor in scaling the method with data size and complexity. In this paper we illustrate the use of Krylov subspace methods to address this issue in a statistical solution to a source separation problem in cosmology where the data size is prohibitively large for direct solution of the required system. Two distinct approaches, adapted from techniques in the literature, are described: one that uses the method of conjugate gradients directly to the Kronecker-structured problem and another that reformulates the system as a Sylvester matrix equation. We show that both approaches produce an accurate solution within an acceptable computation time and with practical memory requirements for the data size that is currently available. |
| title | Fast solution of Sylvester-structured systems for spatial source separation of the Cosmic Microwave Background |
| topic | Numerical Analysis Instrumentation and Methods for Astrophysics 65F10, 65F50, 65F08 |
| url | https://arxiv.org/abs/2204.08057 |