Sweeping Preconditioner for the 3 D 1 Helmholtz Equation

This paper introduces the recursive sweeping preconditioner for the numerical solu4 tion of the Helmholtz equation in 3D. This is based on the earlier work of the sweeping preconditioner 5 with the moving perfectly matched layers (PMLs). The key idea is to apply the sweeping precondi6 tioner recursively to the quasi-2D auxiliary problems introduced in the 3D sweeping preconditioner. 7 Compared to the non-recursive 3D sweeping preconditioner, the setup cost of this new approach 8 drops from O(N4/3) to O(N), the application cost per iteration drops from O(N logN) to O(N), 9 and the iteration count only increases mildly when combined with the standard GMRES solver. 10 Several numerical examples are tested and the results are compared with the non-recursive sweeping 11 preconditioner to demonstrate the efficiency of the new approach. 12