A Discrete Integral Transform for Rapid Spectral Synthesis. van den Bekerom, D. C. M. & Pannier, E. Journal of Quantitative Spectroscopy and Radiative Transfer, December, 2020. ![A Discrete Integral Transform for Rapid Spectral Synthesis [link]](https://bibbase.org/img/filetypes/link.svg "link")
Paper doi abstract bibtex Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- & Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed, and a Python implementation is compared against state-of-the-art spectral code RADIS. The synthesis of a benchmark CO2 spectrum consisting of 1.8M spectral lines and 200k spectral points took only 3.1s using the proposed method (1011 lines × spectral points/s), a factor ∼300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. An experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2∼3 orders performance increase, achieving up to 5 • 1014 lines × spectral points/s.
@article{van_den_bekerom_discrete_2020,
title = {A {Discrete} {Integral} {Transform} for {Rapid} {Spectral} {Synthesis}},
issn = {0022-4073},
url = {http://www.sciencedirect.com/science/article/pii/S0022407320310049},
doi = {10.1016/j.jqsrt.2020.107476},
abstract = {Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- \& Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed, and a Python implementation is compared against state-of-the-art spectral code RADIS. The synthesis of a benchmark CO2 spectrum consisting of 1.8M spectral lines and 200k spectral points took only 3.1s using the proposed method (1011 lines × spectral points/s), a factor ∼300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. An experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2∼3 orders performance increase, achieving up to 5 • 1014 lines × spectral points/s.},
language = {en},
urldate = {2020-12-22},
journal = {Journal of Quantitative Spectroscopy and Radiative Transfer},
author = {van den Bekerom, D. C. M. and Pannier, E.},
month = dec,
year = {2020},
keywords = {mentions sympy, spectroscopy},
pages = {107476},
}
Downloads: 0
{"_id":"nRubkuaQSqSwkYbth","bibbaseid":"vandenbekerom-pannier-adiscreteintegraltransformforrapidspectralsynthesis-2020","authorIDs":[],"author_short":["van den Bekerom, D. C. M.","Pannier, E."],"bibdata":{"bibtype":"article","type":"article","title":"A Discrete Integral Transform for Rapid Spectral Synthesis","issn":"0022-4073","url":"http://www.sciencedirect.com/science/article/pii/S0022407320310049","doi":"10.1016/j.jqsrt.2020.107476","abstract":"Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- & Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed, and a Python implementation is compared against state-of-the-art spectral code RADIS. The synthesis of a benchmark CO2 spectrum consisting of 1.8M spectral lines and 200k spectral points took only 3.1s using the proposed method (1011 lines × spectral points/s), a factor ∼300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. An experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2∼3 orders performance increase, achieving up to 5 • 1014 lines × spectral points/s.","language":"en","urldate":"2020-12-22","journal":"Journal of Quantitative Spectroscopy and Radiative Transfer","author":[{"propositions":["van","den"],"lastnames":["Bekerom"],"firstnames":["D.","C.","M."],"suffixes":[]},{"propositions":[],"lastnames":["Pannier"],"firstnames":["E."],"suffixes":[]}],"month":"December","year":"2020","keywords":"mentions sympy, spectroscopy","pages":"107476","bibtex":"@article{van_den_bekerom_discrete_2020,\n\ttitle = {A {Discrete} {Integral} {Transform} for {Rapid} {Spectral} {Synthesis}},\n\tissn = {0022-4073},\n\turl = {http://www.sciencedirect.com/science/article/pii/S0022407320310049},\n\tdoi = {10.1016/j.jqsrt.2020.107476},\n\tabstract = {Accurate synthetic spectra that rely on large Line-By-Line (LBL)-databases are used in a wide range of applications such as high temperature combustion, atmospheric re-entry, planetary surveillance and laboratory plasmas. Conventionally synthetic spectra are calculated by computing a lineshape for every spectral line in the database and adding those together, which may take multiple hours for large databases. In this paper we propose a new approach for spectral synthesis based on an integral transform: the synthetic spectrum is calculated as the integral over the product of a Voigt profile and a newly proposed three-dimensional “lineshape distribution function”, which is a function of spectral position and Gaussian- \\& Lorentzian width coordinates. A fast discrete version of this transform based on the Fast Fourier Transform (FFT) is proposed, which improves performance compared to the conventional approach by several orders of magnitude while maintaining accuracy. Strategies that minimize the discretization error are discussed, and a Python implementation is compared against state-of-the-art spectral code RADIS. The synthesis of a benchmark CO2 spectrum consisting of 1.8M spectral lines and 200k spectral points took only 3.1s using the proposed method (1011 lines × spectral points/s), a factor ∼300 improvement over the state-of-the-art, with the relative improvement generally increasing for higher number of lines and/or number of spectral points. An experimental GPU-implementation of the method was also benchmarked, which demonstrated another 2∼3 orders performance increase, achieving up to 5 • 1014 lines × spectral points/s.},\n\tlanguage = {en},\n\turldate = {2020-12-22},\n\tjournal = {Journal of Quantitative Spectroscopy and Radiative Transfer},\n\tauthor = {van den Bekerom, D. C. M. and Pannier, E.},\n\tmonth = dec,\n\tyear = {2020},\n\tkeywords = {mentions sympy, spectroscopy},\n\tpages = {107476},\n}\n\n","author_short":["van den Bekerom, D. C. M.","Pannier, E."],"key":"van_den_bekerom_discrete_2020","id":"van_den_bekerom_discrete_2020","bibbaseid":"vandenbekerom-pannier-adiscreteintegraltransformforrapidspectralsynthesis-2020","role":"author","urls":{"Paper":"http://www.sciencedirect.com/science/article/pii/S0022407320310049"},"keyword":["mentions sympy","spectroscopy"],"metadata":{"authorlinks":{}},"downloads":0},"bibtype":"article","biburl":"https://bibbase.org/zotero-group/nicoguaro/525293","creationDate":"2020-12-29T13:44:13.902Z","downloads":0,"keywords":["mentions sympy","spectroscopy"],"search_terms":["discrete","integral","transform","rapid","spectral","synthesis","van den bekerom","pannier"],"title":"A Discrete Integral Transform for Rapid Spectral Synthesis","year":2020,"dataSources":["YtBDXPDiQEyhyEDZC","fhHfrQgj3AaGp7e9E","qzbMjEJf5d9Lk78vE","45tA9RFoXA9XeH4MM","MeSgs2KDKZo3bEbxH","nSXCrcahhCNfzvXEY","ecatNAsyr4f2iQyGq","tpWeaaCgFjPTYCjg3"]}