Method

In the following discussion we use notation of Hauschildt (1992). First, we will describe the process for the formal solution, then we will describe how we construct the non-local approximate $\Lambda$ operator, $\ifmmode{\Lambda^*}\else\hbox{$\Lambda^*$}\fi$, and methods to solve the necessary linear equations in the operator splitting step.