Issue |
A&A
Volume 561, January 2014
|
|
---|---|---|
Article Number | A88 | |
Number of page(s) | 7 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/201322177 | |
Published online | 07 January 2014 |
Compressed convolution
1
Department of Physics and AstronomyUniversity College London,
London
WC1E 6BT,
UK
e-mail:
f.elsner@ucl.ac.uk
2
Institut d’Astrophysique de Paris, UMR 7095, CNRS - Université
Pierre et Marie Curie (Univ Paris 06), 98 bis blvd Arago, 75014
Paris,
France
3
Departments of Physics and Astronomy, University of Illinois at
Urbana-Champaign, Urbana
IL
61801,
USA
Received:
29
June
2013
Accepted:
11
December
2013
We introduce the concept of compressed convolution, a technique to convolve a given data set with a large number of non-orthogonal kernels. In typical applications our technique drastically reduces the effective number of computations. The new method is applicable to convolutions with symmetric and asymmetric kernels and can be easily controlled for an optimal trade-off between speed and accuracy. It is based on linear compression of the collection of kernels into a small number of coefficients in an optimal eigenbasis. The final result can then be decompressed in constant time for each desired convolved output. The method is fully general and suitable for a wide variety of problems. We give explicit examples in the context of simulation challenges for upcoming multi-kilo-detector cosmic microwave background (CMB) missions. For a CMB experiment with detectors with similar beam properties, we demonstrate that the algorithm can decrease the costs of beam convolution by two to three orders of magnitude with negligible loss of accuracy. Likewise, it has the potential to allow the reduction of disk space required to store signal simulations by a similar amount. Applications in other areas of astrophysics and beyond are optimal searches for a large number of templates in noisy data, e.g. from a parametrized family of gravitational wave templates; or calculating convolutions with highly overcomplete wavelet dictionaries, e.g. in methods designed to uncover sparse signal representations.
Key words: methods: data analysis / methods: statistical / methods: numerical / cosmic background radiation
© ESO, 2014
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.