Итерационные методы решения СЛАУ: Якоби, Зейдель и релаксация
Когда размерность матрицы коэффициентов в системе линейных алгебраических уравнений достигает миллионов (что типично для задач теплопроводности, гидродинамики или квантовой химии), прямые методы решения, такие как метод Гаусса или LU-разложение, становятся бессильны. Они требуют кубического времени O(n^3) и гигантских объемов оперативной памяти. Более того, физические матрицы часто являются разреженными (состоят в основном из нулей), а прямые методы превращают эти нули в ненулевые элементы (fill-in), уничтожая структуру памяти. Решением этой проблемы являются итерационные методы. Они работают по принципу постепенного улучшения ответа: задается начальное, случайное приближение, которое шаг за шагом уточняется, пока не будет достигнута требуемая точность.