We test the accuracy of the 3D PBC solution by comparing it to the results of the 1D code for several line scattering parameters. The 1D solver uses 64 depth points, distributed logarithmically in optical depth. Figures 1-4 show the mean intensities at and the component of the emergent flux as function of wavelength for both the 1D ( symbols) and the 3D solver. The agreement is excellent for all values of from unity to , indicating that the 3D code produces an accurate solution even for extreme cases of line scattering. In the case with the continuum processes lead to earler thermalization than the classical approximation as the line strength is limited compared to the continuum. This behavior is the same as in the 1D plane-parallel comparison case. The convergence rate of the line source function (here used together with Ng acceleration) is the same as discussed in Paper II, in the case of the 3D code needed 29 iterations with the nearest-neighbor to reach a relative accuracy of using the simple starting guess . The nearest-neighbor does allow stopping the iterations earlier than a diagonal (local) due to the improved convergence rate (see paper I). This can easily cut the number of iterations by factors of two or more, even greater savings are possible if the accuracy limit is relaxed.

In addition to the mean intensities, we checked that the flux vectors have vanishing components in the and directions, typically in all voxels. We stress that this result is the result of the calculations and is not forced by the numerical scheme.