Тонкости работы LAPACK, на c Опции

[pdk]

Имеется следующая проблема:
понадобилось мне разложение Шура, поиски привели к функции DHSEQR из LAPACK, с помощью f2c перевел с фортрана на c все что понадобилось (всего ~50 файлов), работает, но очень медленно.
Для матрицы N=38 0.32сек ( i920, MSVS без оптимизации).Сложность задачи всего 20*n^3 = 1MFLOP для N=38.
В MATLAB же даже для N=1000 функция schur выполняется мгновенно (судя по документации внутри тот же LAPACK ).

Кому приходилось сталкиваться с переносом кода из LAPACK, есть ли там тонкости?

[tulebaev]

> Справочник алгоритмов на языке АЛГОЛ.
> Вот _тут_ вот без пароля, SMS-ки, денег и совершенно без обмана
Либрусек говорит "Access denied" -- требуется абонемент?

[tulebaev]

Без проблем скачал справочник на "Библиотеке Ихтика"
http://ihtika.net/?qwe=loginfromfile&f...3256989c860bc51

Сообщение отредактировал tulebaev - Aug 27 2010, 07:10