摘要
In a previous work, the author and D.C. Dobson proposed a numerical method for solving the complex Helmholtz equation based on the minimization variational principles developed by Milton, Seppecher, and Bouchitte. This method results in a system of equations with a symmetric positive definite coefficient matrix, but at the same time requires solving simultaneously for the solution and its gradient. Herein is presented a method based on the saddle point variational principles of Milton, Seppecher, and Bouch- itte, which produces symmetric positive definite systems of equations, but eliminates the necessity of solving for the gradient of the solution. The result is a method for a wide class of Helmholtz problems based completely on the Conjugate Gradient algorithm.
In a previous work, the author and D.C. Dobson proposed a numerical method for solving the complex Helmholtz equation based on the minimization variational principles developed by Milton, Seppecher, and Bouchitte. This method results in a system of equations with a symmetric positive definite coefficient matrix, but at the same time requires solving simultaneously for the solution and its gradient. Herein is presented a method based on the saddle point variational principles of Milton, Seppecher, and Bouch- itte, which produces symmetric positive definite systems of equations, but eliminates the necessity of solving for the gradient of the solution. The result is a method for a wide class of Helmholtz problems based completely on the Conjugate Gradient algorithm.
