для любителей копаться с кодом квантовохимических программ

[Yurii]

Я понимаю, что таких любителей - единицы, но как человек. относящийся к этой вымерающей когорте, даю ссылку на модификацию кода одной функции, которая используется не только в диагонализации:
ru: http://software.intel.com/ru-ru/forums/ ... p?p=128709
en: http://software.intel.com/en-us/forums/ ... 1&o=d&s=lr
может кому и сгодится.

[VTur]

Юрий, старый вопрос. Не встретились алгоритмы диагонализации эрмитовых матриц? Не лапаковские или линпаковские.

[Yurii]

Юрий, старый вопрос. Не встретились алгоритмы диагонализации эрмитовых матриц? Не лапаковские или линпаковские.

Может это поможет: http://www.netlib.org/utk/people/JackDo ... la-sw.html

[VTur]

Спасибо!

[Yurii]

Спасибо!

У меня для эрмитовых матриц еще более впечатляющие результаты, чем для вещественных матриц, но только для трехдиагонализации (а она "жрет" основную долю времени) и при этом для упакованных матриц, что очень существенно, т.к. эрмитовы матрицы "потребляют" в два раза больше памяти.

[VTur]

Как Вы сводите эрмитовы матрицы общего вида (не спарз) к трёхдиагональным или треугольным, а потом к трёхдиагональным?

[Yurii]

Как Вы сводите эрмитовы матрицы общего вида (не спарз) к трёхдиагональным или треугольным, а потом к трёхдиагональным?

к трёхдиагональным

[ivan_chem]

Юрий, а есть уже что-нибудь готовое, что можно быстро попробовать (скомпилировать/прилинковать)? Я пишу для gpu, и там настолько все быстро считается что для больших систем основная часть времени уходит на blas и диагонализацию. Получается что скорость начинает определяться эйгенсолвером. Сейчас пользую dsyevd из mkl -- самое быстрое что смог найти

[Yurii]

Юрий, а есть уже что-нибудь готовое, что можно быстро попробовать (скомпилировать/прилинковать)? Я пишу для gpu, и там настолько все быстро считается что для больших систем основная часть времени уходит на blas и диагонализацию. Получается что скорость начинает определяться эйгенсолвером. Сейчас пользую dsyevd из mkl -- самое быстрое что смог найти

dsyevd и dsyev - на настоящий момент единственнные функции в пакете Intel MKL, которым можно доверять. dsyev требует в три раза меньше памяти, чем dsyevd, но значительно медленней. Можно попробовать еще dstemr, в которой реализован RRR алгоритм, разработанный в конце прошлого века. Раньше он работал очень неустойчиво, но сейчас они его подтянули, хотя со 100% точностью за него не ручаюсь. Да и точность собственных векторов в нем не на высоте. Что касается исходников, то можно обратиться к Intel: http://software.intel.com/en-us/forums/ ... p?p=130621

[Yurii]

Вот опубликовал выдержки из своей статьи: http://software.intel.com/ru-ru/forums/ ... p?p=133175