Ксения (04.06.2012 23:17, просмотров: 455) ответил USSR на Похоже, Ксения пытается вычислить QR-разложение матрицы =>
Ну да. Точнее, у меня QL-разложение, то это почти одно и тоже. Только не методом Грама-Шмидта, а сначала приведением к тридиагональной форме, а после QL-итерации по Уилкинсону. Так вот, если арифметика дурная, то большие собственные числа вычисляются правильно, а маленькие - с большой ошибкой. И всё потому, что по своей природе это процесс исчерпывания - экстракция сначала крупных вкладов, когда мелочевка остается в остатке. При короткой мантиссе одни слезы. В обычном 32-битном float вообще жить нельзя, а в 64-битном double жить можно, но пока матрица невелика.