Simultaneous analysis of large INTEGRAL/SPI11Based on observations with INTEGRAL, an ESA project with instruments and science data center funded by ESA member states (especially the PI countries: Denmark, France, Germany, Italy, Spain, and Switzerland), Czech Republic and Poland with participation of Russia and the USA. datasets: Optimizing the computation of the solution and its variance using sparse matrix algorithms

Nowadays, analyzing and reducing the ever larger astronomical datasets is becoming a crucial challenge, especially for long cumulated observation times. The INTEGRAL/SPI  X/γ-ray spectrometer is an instrument for which it is essential to process many exposures at the same time in order to increase t...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Astronomy and computing Ročník 1; s. 59 - 69
Hlavní autoři: Bouchet, L., Amestoy, P., Buttari, A., Rouet, F.-H., Chauvin, M.
Médium: Journal Article
Jazyk:angličtina
Vydáno: Elsevier B.V 01.02.2013
Témata:
Gamma-rays: general
Methods: data analysis
Methods: numerical
Techniques: imaging spectroscopy
Techniques: miscellaneous
Techniques: miscellaneous
Methods: data analysis
Techniques: imaging spectroscopy
Gamma-rays: general
Methods: numerical
ISSN:2213-1337, 2213-1345
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí:Nowadays, analyzing and reducing the ever larger astronomical datasets is becoming a crucial challenge, especially for long cumulated observation times. The INTEGRAL/SPI  X/γ-ray spectrometer is an instrument for which it is essential to process many exposures at the same time in order to increase the low signal-to-noise ratio of the weakest sources. In this context, the conventional methods for data reduction are inefficient and sometimes not feasible at all. Processing several years of data simultaneously requires computing not only the solution of a large system of equations, but also the associated uncertainties. We aim at reducing the computation time and the memory usage. Since the SPI transfer function is sparse, we have used some popular methods for the solution of large sparse linear systems; we briefly review these methods. We use the Multifrontal Massively Parallel Solver (MUMPS) to compute the solution of the system of equations. We also need to compute the variance of the solution, which amounts to computing selected entries of the inverse of the sparse matrix corresponding to our linear system. This can be achieved through one of the latest features of the MUMPS software that has been partly motivated by this work. In this paper we provide a brief presentation of this feature and evaluate its effectiveness on astrophysical problems requiring the processing of large datasets simultaneously, such as the study of the entire emission of the Galaxy. We used these algorithms to solve the large sparse systems arising from SPI data processing and to obtain both their solutions and the associated variances. In conclusion, thanks to these newly developed tools, processing large datasets arising from SPI is now feasible with both a reasonable execution time and a low memory usage. •INTEGRAL/SPI X/γ-ray spectrometer data analysis.•Large astronomical datasets arising from the simultaneous analysis of years of data.•Resolution of a large sparse system of equations; the solution and its variance.•The Multifrontal Massively Parallel Solver (MUMPS) to solve the equations.•MUMPS A−1 feature to compute selected inverse entries (variance of the solution etc.).
ISSN:2213-1337
2213-1345
DOI:10.1016/j.ascom.2013.03.002