Formal Solutions for Polarized Radiative Transfer. II. High-order Methods

When integrating the radiative transfer equation for polarized light, the necessity of high-order numerical methods is well known. In fact, well-performing high-order formal solvers enable higher accuracy and the use of coarser spatial grids. Aiming to provide a clear comparison between formal solvers, this work presents different high-order numerical schemes and applies the systematic analysis proposed by Janett et al. (2017), emphasizing their advantages and drawbacks in terms of order of accuracy, stability, and computational cost.