向量乘BLAS'catlas_saxpby不能正常工作

Question

我想要有兩個任意長度的向量（典型長度將是2048）並且乘以元素。所以對於所有n，Z [n] = X [n] * Y [n]。

我所建立的測試代碼是相當簡單：

float inputX[4] = { 2, 4, 8, 16 }; 
float inputY[4] = { 2, 4, 8, 16 }; 

catlas_saxpby(4, 1, inputX, 1, 1, inputY, 1); 

結果進入inputY，其結果是

4.000000, 8.000000, 16.000000, 32.000000 

其中，如果他們被乘以它應該是4，16 ，64,256。但它看起來像添加。

所以這不是我所期望的，文檔沒有給我足夠的信息來確定它在做什麼。

任何想法？

Apple's documentation for BLAS says this: 

Computes the product of two vectors, scaling each one separately (single-precision). 

void catlas_saxpby (
    const int N, 
    const float alpha, 
    const float *X, 
    const int incX, 
    const float beta, 
    float *Y, 
    const int incY 
); 
Parameters 
N 
Number of elements in the vector. 
alpha 
Scaling factor for X. 
X 
Input vector X. 
incX 
Stride within X. For example, if incX is 7, every 7th element is used. 
beta 
Scaling factor for Y. 
Y 
Input vector Y. 
incY 
Stride within Y. For example, if incY is 7, every 7th element is used. 
Discussion 
On return, the contents of vector Y are replaced with the result. 

Answer 1

Apple文檔有誤。 saxpby函數計算標量爲alpha和beta的表達式alpha*X + beta*Y和矢量X和Y。

我不認爲有一個函數可用於計算兩個向量的元素明智的乘積，因爲這不是線性代數中的常見操作。您可以採用外部產品的對角線，但由於它計算整個外部產品（而不是N），所以這是一種嚴重的工作量浪費。

Answer 2

As Adam Rosenfield表示，文檔不正確。請提交一個錯誤。

除此之外，他的答案還有一些更正。首先，saxpby計算alpha * X + beta * Y。其次，對你更有用：在BLAS中沒有你想做的功能，但在vDSP中也有這樣的功能，它也是Accelerate.framework的一部分：vDSP_vmul。