Как лучше всего диагонализировать треугольную матрицу в Numpy/Scipy? ⇐ Python

[Anonymous]

У меня есть матрица, которую я знаю, треугольная, и я хочу его диагонализировать. Поскольку собственные значения треугольной матрицы легко найти (только диагональные значения), эта проблема теоретически точно решается, поэтому я решил, что использование numpy.linalg.eig может быть излишним и может быть медленнее/не настолько точно, как это может быть. То есть, конечно, если eig достаточно умный, чтобы эффективно работать с треугольной матрицей, но я не смог найти никакой документации, которая уверяет это. Матрицы (~ 200) и приличное количество собственных значений (также ~ 200), поэтому я чувствую, что это может быть не оптимальным.

Подробнее здесь: https://stackoverflow.com/questions/797 ... umpy-scipy