A conjugate gradient method for the solution of the non-LTE line radiation transfer problem

Abstract

This study concerns the fast and accurate solution of the line radiation transfer problem, under non-LTE conditions. We propose and evaluate an alternative iterative scheme to the classical ALI-Jacobi method, and to the more recently proposed Gauss-Seidel and Successive Over-Relaxation (GS/SOR) schemes. Our study is indeed based on the application of a preconditioned bi-conjugate gradient method (BiCG-P). Standard tests, in 1D plane parallel geometry and in the frame of the two-level atom model, with monochromatic scattering, are discussed. Rates of convergence between the previously mentioned iterative schemes are compared, as well as their respective timing properties. The smoothing capability of the BiCG-P method is also demonstrated.

Research note: 4 pages, 5 figures, accepted to A&A

References in corpus (1)

• Fast 2D non-LTE radiative modelling of prominences I. Numerical methods and benchmark results

Cited by in corpus (6)

• Numerical solutions to linear transfer problems of polarized radiation II. Krylov methods and matrix-free implementation
• Numerical non-LTE 3D radiative transfer using a multigrid method
• The transfer of polarized radiation in resonance lines with partial frequency redistribution, J-state interference, and arbitrary magnetic fields. A radiative transfer code and useful approximations
• Gauss-Seidel and Successive Overrelaxation Methods for Radiative Transfer with Partial Frequency Redistribution
• Numerical radiative transfer with state-of-the-art iterative methods made easy
• A new nonlocal thermodynamical equilibrium radiative transfer method for cool stars