Accelerated higher order radiative perturbation computation: an application of GDOM

In a previous paper by the authors, the computational expressions for the higher order terms of the radiative perturbation series have been developed, and numerical results have been obtained. In this paper, new computational expressions are developed that use the analytic Green's function and the GDOM code, which were developed recently by Qin and Box. Our analysis and numerical computations indicate that the new scheme dramatically improves the efficiency of the computation (by reducing the CPU-time to less than 2% of that required by the previous scheme), and the huge demand on memory usage has been completely removed. This new scheme allows us to expand the high order perturbation computations to more general cases, for example, to include azimuth dependence, and to use the perturbation theory in more potential applications.