Асимптотически оптимальное решение модельной задачи для экранированного уравнения Пуассона

Abstract

Экранированное уравнение рассматривается на прямоугольной области со смешанными краевыми условиями. При численном решении этой задачи предлагается использовать итерационную факторизацию после фиктивного продолжения дискретной задачи аппроксимирующей решаемую задачу. В итоге решение основывается на решении систем линейных алгебраических уравнений с матрицами треугольного вида, в которых ненулевых элементов не боле трех в каждой строке. При достаточно малой погрешности аппроксимации рассматриваемой задачи требуемая относительная погрешность предлагаемого итерационного процесса достигается за количество итераций, независящее от параметров дискретизации. Итерационный процесс оказывается методом, дающим оптимальную асимптотику по количеству операций в арифметических действиях. Разработанный итерационный процесс основывается на характерных особенностях указанной модельной задачи. Эта задача может быть получена в методах фиктивных компонент, пространств, когда решают краевые задачи для эллиптических уравнений в областях сложной формы. Приводится алгоритм реализации итерационного метода с выбором итерационных параметров в автоматическом режиме, с применением метода минимальных невязок, поправок. Это дает критерий для остановки итерационного процесса при получении указанной предварительно относительной погрешности. Приведен простейший тестовый пример для вычислительных экспериментов, подтверждающих асимптотическую оптимальность для итерационного метода в количестве вычислительных затрат. Реализация метода существенно основывается на использовании комплексного анализа. The screened equation is considered on rectangular area, with the mixed boundary conditions. At the numerical solution of this task it is suggested to use iterative factorization after fictitious continuation of the discrete task approximating the solved task. As a result, the solution is based on the solution of systems of the linear algebraic equations with matrixes of the triangular type, which contain no more than three nonzero elements in every line. At rather small error of approximation of the considered task, the demanded relative error of the suggested iterative process is reached in over the number of iterations independent of the discretization parameters. The iterative process turns out to be the method giving an optimal asymptotics as per the number of operations in arithmetic actions. The developed iterative process is based on the characteristic specifics of the stated model task. This task can be obtained in the methods of fictitious components, spaces, when boundary problems for the elliptic equations in areas of a complex shape are being solved. The algorithm is given for the fulfillment of the iterative method with choosing of the iterative parameters in the automatic mode, by applying the minimal residual, corrections method. This gives a criterion to stop the iterative process when the specified preliminarily relative error is obtained. An elementary test example is given on the computing experiments confirming the asymptotic optimality for the iterative method in the number of computing expenses. The fulfillment of the method is substantially based on the use of complex analysis.