Matrix 带矩阵的Blas daxpy例程 我正在研究C++中的一些矩阵相关问题。我想解决这个问题：Y=aX+Y，其中X和Y是矩阵，a是常数。我曾考虑过使用daxpy BLAS例程，但是，根据文档，daxpy是一个向量例程，我没有得到与在matlab中解决相同问题时相同的结果

我目前正在运行以下程序：

F77NAME(daxpy)(N, a, X, 1, Y, 1);

---

当需要执行操作Y=a*X+Y时，X'，Y`是1D矩阵还是2D矩阵并不重要，因为操作是按元素进行的

因此，如果将矩阵分配给单指针double A[]=new[]M*N；，然后可以通过将向量的维数定义为M*N来使用daxpy

堆栈二维矩阵双A[3][2]也是如此；因为这个内存是按顺序分配的

否则，您需要使用for循环并分别添加每一行

int MN = M*N;
int one = 1;
F77NAME(daxpy)(&MN, &a, &X, &one, &Y, &one);