Многопоточный LAPACK

Для атласа надо линковать либу ptf77blas. Прирост будет только на алгоритмах, использубщих blas level 3.

Нахождение обратной матрице к симметричной неопределенной вроде использует level3. Или нужно напрямую использовать blas, чтобы профит появился?

Нет, нужно только правильно слинковать. Проверь через ldd свой бинарь.

Показывает liblapacke, libblas, liblapack, даже когда делаю -L/path/ -lpt77blas

Вот это плохо, скорее всего придется пересобирать lapack.

Еще есть вариант статикой слинковать liblapack.a libatlas.a libpt77blas.a

Мб можно перестать юзать lapacke?
Я использую lapack (atlas), который сам собрал.

многопоточный lapack

Intel MKL

проприетарно, организация не позволяет:-)

gotoBLAS, MKL, Sun perf lib, ACML и еще куча альтернативных реализаций

gotoBLAS практически рипнулась, заменена на openblas, который даже медленне нетюненного ATLAS, sub и acml проприетарно.

Я использую lapack (atlas), который сам собрал.

Ничего не понял. Раз собрал сам lapack, то и используй его и не линкуй системный. Вообще, atlas при компиляции может и lapack собрать, но делает он это исключительно для линковки lapack'овских so'шников с atlas'ом. lapack весь построен на блочных алгоритмах и если blas level3 у тебя распараллелена, то и lapack оказывается автоматически распараллеленным.

Не можешь схематично расписать как все это друг в друга соединяется и что еще подключаешь?

Да, использую собранный и аккуратно упакованный в пакет atlas, но он использует только одно ядро. Сейчас таки собрал plasma, она работает.

Для атласа надо линковать либу ptf77blas

где она в дебиане?

lapack работает с матрицами _только_ с помощью вызовов blas, при этом многие алгоритмы из его состава работают с матрицами блочно и используют функции blas level3, которые эффективно параллелятся. Соответственно, если нижележащая реализация blas распараллелена, то будет использоваться несколько ядер. Ничего дополнительно подключать не надо.

Откуда мне знать, в mageia есть. Вообще, atlas ооочень рекомендуется собирать самому под свой проц.

Я бы тогда проверил наличие параллельности на нижележащем слое, запустив dgemm.

пока я ограничился опенбласом :) ладно не судьба увидеть атлас в паралельном виде.

Работает в параллельном виде он отлично. Я выжимал из него для своего проца gflops близкие к тем, что указаны на сайте интеля :)

проприетарно, организация не позволяет:-)

Однакож. Религиозная организация?

о, увидел на гитхабе опенбласа как установить многопоточность. спс.

мы чё исходный код будем выкладывать что ли?

мы чё исходный код будем выкладывать что ли?

Ну а ежели не будете, то какая разница? Под линукс интели дают за так, но для некоммерческого пользования.

Ну а организации коммерческое использование, так что вот ТСу и не подходит

а openblas ты собирал или использовал пакет?

я про атлас, openblas не использовал

ок, ясно.

А использовал только blas3 функции? PLASMA использует не многопоточный BLAS, распараллеливает на верхнем уровне, и соответствующее ускорение есть. Интересно, что будет, если использовать LAPACKE с параллельным BLAS.

Собрал многопоточный OpenBLAS, PLASMA получается лучше.

А использовал только blas3 функции?

Да, blas1 и 2 не параллелятся почти никак.

Кастую Reset, ты сравнивал OpenBLAS и ATLAS по скорости?

OpenBLAS вообще не щупал. Сравнивал mkl, acml, atlas на функции dgemm (умножение произвольных матриц). На моих q8200 и q9550 получилось mkl<atlas<acml. Удивительно, но родная интелевская библиотека показала худшие результаты. Сравнивал где-то полтора года назад на текущих на тот момент версиях.

Ясно, просто мне нужно под виндой проект собрать, использующий atlas, пока не получается, думал на другое перейти. Спасибо!

acml под виндой работает

Лицензия не та:-(