Optimized parameterized Uzawa methods for solving complex Helmholtz equations

Abstract

As we know that, the parameterized Uzawa (PU) method can be very efficient when used to solve the standard saddle point problem, especially, when we have good and accurate estimation of preconditioned Schur complement matrix. In this paper, by taking full use of the special structure of coefficient matrix arising from complex Helmholtz equations, two types of optimized PU (OPU) methods are discussed theoretically and experimentally. Specifically, the convergence factors of these two OPU methods are less than 0.172, and the optimal result of first OPU method can reach 0.0396 with σ1≥σ2>0, which is currently the best theoretical result in the literatures. Moreover, the second OPU method has better computational advantages compared with the first OPU method as it avoids the inverse computation of W at each step, indicating the less CPU time will be costed by the second OPU method. In addition, the optimal parameters involved in our algorithms consist of constants, reducing the computational complexity associated with parameter selection. Finally, numerical results are given, not only show the effectiveness of OPU methods, but also confirm the rationality of theoretical analysis.