Numerical solutions to linear transfer problems of polarized radiation II. Krylov methods and matrix-free implementation

Abstract

Context. Numerical solutions to transfer problems of polarized radiation in solar and stellar atmospheres commonly rely on stationary iterative methods, which often perform poorly when applied to large problems. In recent times, stationary iterative methods have been replaced by state-of-the-art preconditioned Krylov iterative methods for many applications. However, a general description and a convergence analysis of Krylov methods in the polarized radiative transfer context are still lacking. Aims. We describe the practical application of preconditioned Krylov methods to linear transfer problems of polarized radiation, possibly in a matrix-free context. The main aim is to clarify the advantages and drawbacks of various Krylov accelerators with respect to stationary iterative methods. Methods. We report the convergence rate and the run time of various Krylov-accelerated techniques combined with different formal solvers when applied to a 1D benchmark transfer problem of polarized radiation. In particular, we analyze the GMRES, BICGSTAB, and CGS Krylov methods, preconditioned with Jacobi, or (S)SOR. Results. Krylov methods accelerate the convergence, reduce the run time, and improve the robustness of standard stationary iterative methods. Jacobi-preconditioned Krylov methods outperform SOR-preconditioned stationary iterations in all respects. In particular, the Jacobi-GMRES method offers the best overall performance for the problem setting in use. Conclusions. Krylov methods can be more challenging to implement than stationary iterative methods. However, an algebraic formulation of the radiative transfer problem allows one to apply and study Krylov acceleration strategies with little effort. Furthermore, many available numerical libraries implement matrix-free Krylov routines, enabling an almost effortless transition to Krylov methods.

References in corpus (12)

• A New Algorithm for 2-D Transport for Astrophysical Simulations: I. General Formulation and Tests for the 1-D Spherical Case
• Formal Solutions for Polarized Radiative Transfer. I. The DELO Family
• Polarized Line Formation in Multi-Dimensional Media.III. Hanle Effect with Partial Frequency Redistribution
• Formal Solutions for Polarized Radiative Transfer. IV. Numerical Performances in Practical Problems
• Preconditioned Bi-Conjugate Gradient Method for Radiative Transfer in Spherical Media
• Formal Solutions for Polarized Radiative Transfer. II. High-order Methods
• Polarized Line Formation in Multi-dimensional Media. V. Effects of Angle-Dependent Partial Frequency Redistribution
• Numerical solutions to linear transfer problems of polarized radiation I. Algebraic formulation and stationary iterative methods
• A conjugate gradient method for the solution of the non-LTE line radiation transfer problem
• Polarized Line Formation in Arbitrary Strength Magnetic Fields: the case of a two-level atom with hyperfine structure splitting
• Multi-dimensional radiative transfer to analyze Hanle effect in Ca {\sc ii} K line at 3933 Å\,
• A new nonlocal thermodynamical equilibrium radiative transfer method for cool stars

Cited by in corpus (8)

• Modeling the scattering polarization of the solar Ca i 4227 Å line with angle-dependent partial frequency redistribution
• Novel framework for the three-dimensional NLTE inverse problem
• 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
• Numerical solutions to linear transfer problems of polarized radiation I. Algebraic formulation and stationary iterative methods
• Assessment of the CRD approximation for the observer's frame RIII redistribution matrix
• Diagnostic potential of wavelength-integrated scattering polarisation signals of the solar He II Ly-alpha line
• A numerical approach for modelling the polarisation signals of strong resonance lines with partial frequency redistribution. Numerical applications to two-term atoms and plane-parallel atmospheres
• Jacobian-Free Newton-Krylov method for multilevel NLTE radiative transfer problems