vector - BLASを使用して1D配列の合計をどのように計算しますか？

Question

BLASレベル1には、ベクトルのL1およびL2ノルムを計算する*ASUMおよび*NRM2がありますが、ベクトルの（符号付き）合計をどのように計算するのでしょうか。別のベクトルを1つでいっぱいにして、*DOT...を実行するよりも優れたものが必要です。

Answer 1

BLASは、線形代数ライブラリで頻繁に必要とされる演算ではないため、求めているような水平和演算を提供しません。

多くのDSPライブラリはこの操作を提供します。たとえば、OS XおよびiOSではvDSP_sve( )、Accelerateフレームワークによって提供される機能を使用します。残念ながら、利用可能なDSPライブラリはプラットフォームごとに大きく異なる傾向があるため、ターゲットにしているプラ​​ットフォームについて詳しく知る必要があります。

Answer 2

2番目のベクトルの増分がゼロである内積を行うことができます。Cを使用すると、次のようになります。

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

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

Answer 3

1つの方法は、1のベクトルを持つ内積を使用することです。具体的には、cblas_caxpy関数を使用します。

Answer 4

http://www.netlib.org/blas/blasqr.pdfに見られるように、xAXPYベクトルの合計をサポートします。