Тонкости работы 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, есть ли там тонкости?

[Oldring]

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

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

Классическая книжка Голуба по матричным вычислениям вас спасет.
Там проще написать с нуля, чем переносить. Там десяток строк всё разложение Шура, а не полсотни файлов.

[pdk]

В книгу обязательно загляну, но фраза из документации LAPACK намекает что все может быть сложнее:

Although we started from a working experimental code and
we reused some of the older Lapack 3.0 code (in xLAHQR), nevertheless it
took many weeks of full time eort to write, debug and tune a nal production
version of the new QR algorithm