如何使用cblas_dgemm进行向量外积？
我试图用一个行向量和一个列向量相乘。我可以用dgemm吗 换句话说，D=A*B 其中D是矩阵，a是列向量，B是行向量 我遵循了这里找到的文档。我似乎无法为cblas_ggemm获得正确的参数 这是我的尝试。在我的例子中，m=nRows，n=nCols，k=1 问题似乎是土地发展局、土地发展局和土地发展公司。我将它们分别定义为ncol、k、nRows如何使用cblas_dgemm进行向量外积？,c,matrix,vector,blas,intel-mkl,C,Matrix,Vector,Blas,Intel Mkl,我试图用一个行向量和一个列向量相乘。我可以用dgemm吗 换句话说，D=A*B 其中D是矩阵，a是列向量，B是行向量 我遵循了这里找到的文档。我似乎无法为cblas_ggemm获得正确的参数 这是我的尝试。在我的例子中，m=nRows，n=nCols，k=1 问题似乎是土地发展局、土地发展局和土地发展公司。我将它们分别定义为ncol、k、nRows #include <stdio.h> #include <time.h> #include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <mkl.h>
#define nCols 5
#define nRows 20
#define k 1
void PrintMatrix(double* pMatrix, const size_t nR, const size_t nC, const CBLAS_ORDER Order) {
unsigned int i, j;
if (Order == CblasRowMajor)
{
for (i = 0; i < nR; i++)
{
for (j = 0; j < nC; j++)
{
printf("%f \t ", pMatrix[i * nC + j]); // !!!
}
printf("\n"); // !!!
}
}
else
{
for (i = 0; i < nR; i++) {
for (j = 0; j < nC; j++) {
printf("%f \t ", pMatrix[i + j* nR ]); // !!!
}
printf("\n"); // !!!
}
}
printf("\n"); // !!!
}
int main(void) {
double A[] = { 8, 4, 7, 3, 5, 1, 1, 3, 2, 1, 2, 3, 2, 0, 1, 1 , 2, 3, 4, 1};
double B[] = { -1, 2, -1, 1, 2 };
double alpha = 1.0, beta = 0.0;
int i, lda, ldb, ldc;
double *C, *D;
D = (double*) malloc(nRows * nCols * sizeof(double));
C = (double*) malloc(nRows * nCols * sizeof(double));
for (i = 0; i < nRows*nCols; i++)
D[i] = 0.0;
for (i = 0; i < nRows*nCols; i++)
C[i] = 0.0;
lda = nCols;
ldb = k;
ldc = nRows;
cblas_dger(CblasRowMajor, nRows, nCols, alpha, A, 1, B, 1, C, nCols);
PrintMatrix(C, nRows, nCols,CblasRowMajor);
cblas_dgemm (CblasRowMajor, CblasNoTrans, CblasNoTrans, nRows, nCols, k, alpha, A, lda, B, ldb, beta, D, ldc);
PrintMatrix(D, nRows, nCols, CblasRowMajor);
free(D);
free(C);
return 0;
}
#包括
#包括
#包括
#包括
#包括
#包括
#定义NCOL5
#定义nRows 20
#定义k1
无效打印矩阵（双*P矩阵、常数大小\u t nR、常数大小\u t nC、常数CBLAS\u顺序）{
无符号整数i，j；
如果（订单==CBLASROWMARY）
{
对于（i=0；idgemmdgercblas_dgemmCblasRowMajorlda = 1;
ldb = 5;
ldc = 5;
下面的代码可以工作。我转到了major专栏，因为在Fortran BLAS文档的上下文中更容易理解主要维度问题。很抱歉，我没能准确地解决你提出的问题
#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <mkl.h>
void PrintMatrix(double* pMatrix, const size_t nR, const size_t nC, const CBLAS_ORDER Order) {
unsigned int i, j;
if (Order == CblasRowMajor) {
for (i = 0; i < nR; i++) {
for (j = 0; j < nC; j++) {
printf("%f \t ", pMatrix[i * nC + j]); // !!!
}
printf("\n"); // !!!
}
} else {
for (i = 0; i < nR; i++) {
for (j = 0; j < nC; j++) {
printf("%f \t ", pMatrix[i + j* nR ]); // !!!
}
printf("\n"); // !!!
}
}
printf("\n"); // !!!
}
int main(void)
{
const int m = 20;
const int n = 5;
const int k = 1;
double A[] = { 8, 4, 7, 3, 5, 1, 1, 3, 2, 1, 2, 3, 2, 0, 1, 1, 2, 3, 4, 1};
double B[] = { -1, 2, -1, 1, 2 };
double alpha = 1.0, beta = 0.0;
double * C = (double*) malloc(m * n * sizeof(double));
double * D = (double*) malloc(m * n * sizeof(double));
for (int i = 0; i < m*n; i++) C[i] = 0.0;
for (int i = 0; i < m*n; i++) D[i] = 0.0;
int lda = 20;
int ldb = 1;
int ldc = 20;
cblas_dger(CblasColMajor, m, n, alpha, A, 1, B, 1, C, ldc);
cblas_dgemm(CblasColMajor, CblasNoTrans, CblasNoTrans,
m, n, k,
alpha, A, lda,
B, ldb, beta,
D, ldc);
PrintMatrix(C, m, n, CblasRowMajor);
PrintMatrix(D, m, n, CblasRowMajor);
free(D);
free(C);
return 0;
}
#包括
#包括
#包括
#包括
#包括
#包括
无效打印矩阵（双*P矩阵、常数大小\u t nR、常数大小\u t nC、常数CBLAS\u顺序）{
无符号整数i，j；
如果（订单==CBLASROWMARY）{
对于（i=0；icblas.hCblasRowMajorCblasNoTrans