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

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


lapack numerical
person quant_dev    schedule 07.06.2009    source источник

Ответы (2)


arrow_upward
2
arrow_downward

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

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

person codehippo    schedule 26.06.2009
comment
Спасибо, это то, что я хотел услышать. - person quant_dev; 27.06.2009

arrow_upward
0
arrow_downward

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

  • Вероятно, это зависит от того, какую библиотеку LAPACK вы используете. Их существует несколько, оптимизированных для разных платформ. Вы используете Netlib, MKL, ACML, ??
  • Зачем вам верить на слово совершенно незнакомому человеку, когда вы можете измерить это сами?
person Die in Sente    schedule 22.06.2009
comment
Я предполагаю, что точность будет различаться гораздо меньше между реализациями, чем производительность. Что касается измерений, у меня нет опыта подготовки правильного набора тестовых матриц. - person quant_dev; 23.06.2009