LU-разложение: оптимизация решения систем линейных уравнений
Метод Гаусса — это первый и самый известный алгоритм для решения систем линейных алгебраических уравнений (СЛАУ), который мы изучаем в университете. Однако в чистом виде (с постоянным приведением расширенной матрицы к ступенчатому виду) он не используется в серьезном программном обеспечении. С точки зрения современной вычислительной линейной алгебры, прямой ход метода Гаусса представляет собой не что иное, как факторизацию (разложение) матрицы коэффициентов в произведение двух треугольных матриц. Этот процесс называется LU-разложением (от английских слов Lower и Upper). Понимание и применение LU-разложения позволяет ускорить решение инженерных и физических задач в десятки и сотни раз, экономя колоссальные объемы вычислительных ресурсов суперкомпьютеров.