c – BLAS中的元素向量向量乘法？

是的,可以单独使用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]