Точность ЖЕЕВ и ЖЕЕВД

Question

Я использую LAPACK для диагонализации сложных эрмитовых матриц. Я могу выбирать между ЖЕЕВ и ЖЕЕВД. Какая из этих процедур является более точной для матриц размера 40 и диапазона собственных значений от 1E-2 до 1E1?

Answer 1

ZHEEVD использует метод «разделяй и властвуй» для вычисления собственных значений.

Если ваши матрицы имеют размер 40 x 40, а собственные значения находятся в диапазоне [1e-2, 1e1], тогда у вас не должно быть абсолютно никаких числовых проблем. Вы можете использовать любую рутину.

Answer 2

Я не знаю ответа, но,

• Вероятно, это зависит от того, какую библиотеку LAPACK вы используете. Есть ряд из них, оптимизированных для различных платформ. Используете ли вы Netlib, MKL, ACML, ??
• С чего бы вам взять слово незнакомца, если вы можете измерить его самостоятельно?