Как вычислить сумму одномерного массива с BLAS?

Question

На уровне BLAS 1 есть * ASUM и * NRM2, которые вычисляют нормы векторов L1 и L2, но как вычислить (векторную) сумму вектора?Должно быть что-то лучше, чем заполнить другой вектор, полный единиц, и сделать * DOT ...

Answer 1

BLAS не обеспечивает операцию горизонтальной суммы, как вы ищете, потому что это не операция, которая часто требуется библиотекам линейной алгебры.

Многие библиотеки DSP делают обеспечивают эту операцию; например, в OS X и iOS вы должны использовать функцию vDSP_sve( ), предоставляемую платформой Accelerate. К сожалению, доступные библиотеки DSP, как правило, сильно различаются от платформы к платформе, поэтому нам нужно знать больше о том, на какую платформу вы ориентируетесь.

Answer 2

Вы можете создать точечное произведение, в котором второй вектор имеет приращение нуля.С использованием C это будет выглядеть так:

int n;
int ix = 1;
int iy = 0;
double y = 1.0;

ddot_(&n, x, &ix, &y, &iy);

Answer 3

Одним из способов является использование точечного произведения с вектором единиц, более конкретно, использование функции cblas_caxpy.

Answer 4

Как видно из http://www.netlib.org/blas/blasqr.pdf, xAXPY поддерживает векторное суммирование.