Стаття, яку ви цитуєте, захищає Гауссова елімінація, кажучи, що, хоча вона є чисельною нестабільністю, вона, як правило, добре працює на випадкових матрицях, і оскільки більшість матриць, про які можна думати, схожі на випадкові матриці, у нас повинно бути гаразд. Це ж твердження можна сказати про безліч нестабільних методів.
Розглянемо простір усіх матриць. Ці методи прекрасно працюють майже скрізь. Тобто 99,999 ...% усіх матриць, які можна створити, не матиме проблем із нестабільними методами. Існує лише дуже мала частка матриць, з якою GE та інші матимуть труднощі.
Проблеми, які хвилюють дослідників, як правило, полягають у тій невеликій частці.
Ми не будуємо матриці випадковим чином. Ми будуємо матриці з дуже особливими властивостями, які відповідають дуже спеціальним, невипадковим системам. Ці матриці часто є умовними.
Геометрично ви можете розглянути лінійний простір усіх матриць. Через цей простір існує підпростір нуля об'єм / міра сингулярних матриць. Багато проблем, які ми будуємо, згруповані навколо цього підпростору. Вони не розподіляються випадковим чином.
Як приклад розглянемо рівняння або дисперсію тепла. Ці системи, як правило, видаляють інформацію із системи (усі початкові стани тяжіють до єдиного кінцевого стану), і в результаті матриці, що описують ці рівняння, надзвичайно одиничні. Цей процес є малоймовірним у випадковій ситуації, але є всюдисущим у фізичних системах.