Preconditioned Bi-Conjugate Gradient Method for Radiative Transfer in Spherical Media

Abstract

A robust numerical method called the Preconditioned Bi-Conjugate Gradient (Pre-BiCG)method is proposed for the solution of radiative transfer equation in spherical geometry.A variant of this method called Stabilized Preconditioned Bi-Conjugate Gradient (Pre-BiCG-STAB) is also presented. These are iterative methods based on the construction of a set of bi-orthogonal vectors. The application of Pre-BiCG method in some benchmark tests show that the method is quite versatile, and can handle hard problems that may arise in astrophysical radiative transfer theory.

19 pages, 12 figures

Cited by in corpus (6)

• Numerical solutions to linear transfer problems of polarized radiation II. Krylov methods and matrix-free implementation
• 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
• Multi-dimensional radiative transfer to analyze Hanle effect in Ca {\sc ii} K line at 3933 Å\,
• A Nonlinear Three-Moment Model for Radiative Transfer in Spherical Symmetry
• Jacobian-Free Newton-Krylov method for multilevel NLTE radiative transfer problems