A PDE-Level Iterative Method for the Numerical Solution of Helmholtz Equation in Underwater Acoustics