Particular features of implementation of an unsaturated numerical method for the exterior axisymmetric Neumann problem (Q2436093)

The author continues the study from [Dokl. Math. 76, No. 3, 882--885 (2007); translation from Dokl. Akad. Nauk, Ross. Akad. Nauk 417, No. 4, 442--445 (2007; Zbl 1188.65159)] and [Sib. Math. J. 52, No. 6, 980--994 (2011); translation from Sib. Mat. Zh. 52, No. 6, 1234--1252 (2011; Zbl 1248.65125)] of a numerical method for solving the exterior axisymmetric Neumann problem for the Laplace equation. The problem is transformed to a one-dimensional weakly singular integral equation which is solved with even trigonometric polynomials. It was shown, that this method is unsaturated, i.e., the principal error term is determined by the smoothness of the solution, resulting in exponential convergence for smooth solutions. In the present paper, the author applies the method for solving the spectral problem for the corresponding operator. The deviation of the first eigenvalue of the discretized problem from the eigenvalue of the Neumann operator is estimated with efficiently computable parameters. This allows to choose the parameters of the proposed iteration method for solving the linear system with a non-symmetric dense matrix such that the accuracy is preserved.