BLAS и поэлементые действия и срезы

Зачем для такой элементарщины брать BLAS, если можно просто в цикле умножать. А для распараллеливания взять openmp.

А для распараллеливания взять openmp.

На данной задаче он будет давать эффект только с -O0.

А вообще blas это не про это, а cblas это вообще не про эффективность. Если нужна эффективность, то надо брать реализацию от вендора железа или atlas, собранный под твое железо.

Так вот над этим и размышляю ). Будет ли BLAS быстрее циклов? Ведь в нем есть зачем-то alpha x + y - эту операцию тоже можно без проблем сделать в цикле.

P.S. хорошо бы учебник какой-нить основательный найти, и чтоб для чайников

atlas

ну это уже второе дело, как я понимаю atlas или mkl используют CBLAS/CLAPACK для Си.

первое дело --- ускорит ли он поэлементые действия для срезов массивов? и есть ли возможности работать со срезами ... ?

Будет ли BLAS быстрее циклов?

Сильно сомневаюсь.

Почему решил подумать над усовершенствованием кода:

существуют же octave, numpy (может быть matlab)... которые используют BLAS, LAPACK и утверждают, что поэлементые и матричные операции очень эффективны.

octave

Там вообще ни о какой производительности говорить не приходится! Тот еще тормоз.

octave, numpy (может быть matlab)... которые используют BLAS, LAPACK

Они используют любую совместимую библиотеку. OpenBLAS и ATLAS совместимы со стандартными реализациями.

Тот еще тормоз.

Только blas и lapack тут совершенно не при чём.

atlas, собранный под твое железо.

А оно не умеет собираться стразу под все оптимизации и динамически выбирать математическое ядро, как openblas?

А не проще посмотреть на что-то вроде http://arma.sourceforge.net ?

как я понимаю atlas или mkl используют CBLAS/CLAPACK для Си.

Не используют. От них там только api.

первое дело --- ускорит ли он поэлементые действия для срезов массивов?

Не ускорит.

Оно почти так и делает, только итоговый код выбирает на основе запуска тестов на твоем железе, поэтому в итоговый бинарь попадает не всё.

Поэлементное умножение не входит в blas api.

хорошо бы учебник какой-нить основательный найти

Деммель Вычислительная линейная алгебра.

и чтоб для чайников

Такого не бывает :)

первое дело --- ускорит ли он поэлементые действия для срезов массивов?

Не ускорит

тогда зачем в BLAS функция alpha * x + y, ведь это же поэлементное сложение двух векторов (с поэлементным умножение на скаляр)? если, в циклы работают не медлее блазовской функции, то такого рода операцию лучше делать просто в цикле.

Поэлементное умножение не входит в blas api.

А со срезами как дела в BLAS-подобных и LAPACK-подобных библиотеках?

получится у меня подматрицу извлечь

например

int x[10];

x = y[10:20] // здесь какая-то функция из BLAS/LAPACK вместо оператора []

За книжку спасибо, устанавливает связи между численными методами и LAPACK!

memcpy тебя чем не устраивает?

// правда, для многомерных будет долго и нудно

BLAS/LAPACK умеет суммировать элементы вектора?

Откуда ж мне знать? Читай мануал. Я не использую внешние библиотеки для такой элементарщины. Я вообще напрямую BLAS/LAPACK никогда не использовал, только через интерфейсы GSL. Да и нет у меня практически таких задач. Разве только аппроксимация наименьшими квадратами. Фурье я при помощи libfftw3 делаю, через Фурье же и свертку реализую; остальное все либо велосипежу с openmp, либо на куде (если есть возможность сильно распараллелить, а вычисления достаточно нудные).

1. поддерживает ли CBLAS/CLAPACK поэлементное умножение векторов ??? другие операции (+-/) тоже интересуют.

Нет, но есть MKL интеловская, там есть.

2. можно ли делать срезы ???
3. имеется ли в этих библиотеках какой-нить удобный интерфейс для многомерных массивов, например

Если устроит на одном ядре, см Eigen.