Разложение A=LU представляет квадратную матрицу в виде нижнетреугольной L и верхнетреугольной U. При единичной диагонали L прямой ход особенно прост. Необратимая матрица может не иметь такого разложения без перестановок строк: нулевой ведущий элемент может потребовать перестановки. Для PA=LU решают Ly=Pb, затем Ux=y. Перестановка действует как на правую часть, так и на матрицу коэффициентов. Ненулевые главные миноры оправдывают обычное исключение без перестановок, а не только сам факт обратимости.
В LU=b вычислите y сверху вниз, вычитая уже известные члены; делите на текущий диагональный элемент, если он не равен единице. В Ux=y вычисляйте x снизу вверх. Для каждого уравнения подставьте значения обратно в исходную строку для проверки. Повторное использование LU для множества правых частей позволяет избежать повторения процесса исключения: плотная факторизация требует работы кубического порядка по размерности, тогда как каждое треугольное решение — квадратичного. Точная арифметика может проверять малые примеры; числовое ведение пивота снижает некоторые проблемы округления, но не устраняет присущую плохую обусловленность.
Для B:V→W и A:W→Z ядро B содержится в ядре (A∘B), но последнее может быть больше. Дополнительные векторы — это те, которые отображаются B в ядро A. Ограничьте B на ядре (A∘B): его образ равен im B∩ker A, а ядро равно ker B. Формула ранга-нуля для этого ограниченного отображения дает dim ker(A∘B)=dim ker B+dim(im B∩ker A). Именно пересечение, а не всё ядро A, определяет добавление к нулю.
Для эндоморфизмов R⁶ с нуль-рангами A=2 и B=3, ранг B=3, а размерность пересечения может варьироваться от 0 до 2. Следовательно, нуль-ранг A∘B варьируется от 3 до 5. Границы зависят от общего промежуточного пространства: для подпространств размерностей r и s в m-мерном пространстве их пересечение имеет размерность между max(0,r+s−m) и min(r,s). Обратное умножение может изменить нуль-ранг, хотя AB и BA оба определены. Выбирайте совместимые пространства перед применением этих формул.