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

Question

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

Answer 1

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

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

Answer 2

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

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