有没有办法用BLAS,GSL或任何其他高性能库进行元素向量向量乘法？
解决方法
(以字面上的题目标题…)
是的,可以单独使用BLAS(尽管它可能不是最有效的方式).
诀窍在于将一个输入向量作为对角线矩阵：
⎡a ⎤ ⎡x⎤ ⎡ax⎤ ⎢ b ⎥ ⎢y⎥ = ⎢by⎥ ⎣ c⎦ ⎣z⎦ ⎣cz⎦
然后,您可以使用矩阵向量乘法函数中的一个,该函数可以将对角矩阵作为输入,而不需要填充,例如. SBMV
例：
void ebeMultiply(const int n,const double *a,const double *x,double *y)
{
extern void dsbmv_(const char *uplo,const int *n,const int *k,const double *alpha,const int *lda,const int *incx,const double *beta,double *y,const int *incy);
static const int k = 0; // Just the diagonal; 0 super-diagonal bands
static const double alpha = 1.0;
static const int lda = 1;
static const int incx = 1;
static const double beta = 0.0;
static const int incy = 1;
dsbmv_("L",&n,&k,&alpha,a,&lda,x,&incx,&beta,y,&incy);
}
// Test
#define N 3
static const double a[N] = {1,3,5};
static const double b[N] = {1,10,100};
static double c[N];
int main(int argc,char **argv)
{
ebeMultiply(N,b,c);
printf("Result: [%f %f %f]\n",c[0],c[1],c[2]);
return 0;
}
结果：[1.000000 30.000000 500.000000]
