Итерационная реализация спектрального метода Чебышева для двумерных эллиптических уравнений с переменными коэффициентами

Построены и исследованы два варианта итерационных алгоритмов реализации спектрального метода Чебышева для двумерных эллиптических уравнений с переменными коэффициентами. Рассмотренные алгоритмы основаны на использовании стабилизированной версии итерационного метода бисопряженных градиентов с комбинированным переобусловливателем в виде диагональной матрицы коэффициентов уравнения и дискретного аналога оператора Лапласа, представленного конечно-разностной или спектральной аппроксимацией. Для обработки дискретного аналога оператора Лапласа в первом случае реализован итерационный метод переменных направлений с оптимальным набором итерационных параметров, а во втором случае – алгоритм Бартельса-Стюарта. На основе численных экспериментов показана высокая эффективность предлагаемых алгоритмов. В обоих случаях количество итераций практически не зависит от размерности сетки и умеренно возрастает при увеличении степени неоднородности коэффициентов задачи. Вычислительная сложность алгоритмов характеризуется величиной $O(N\sqrt{N})$, где $N$ – количество узлов сетки. Несмотря на существенную субоптимальность вычислительной сложности, при размерности сетки $N=n\times n, n\leq 300$, временные затраты на реализацию алгоритмов демонстрируют значения не хуже, чем у алгоритмов оптимальной вычислительной сложности $O(N)$.