Openblas cblas非方阵分割错误

Question

我在运行openblas cblas测试程序时遇到了问题。我的Cblas在乘以方阵时运行得很好，但当我用非方阵尝试它时，我得到了错误“分割错误-核心转储”我检查并重新检查了维度问题，但它们似乎是正确的，所以我想知道哪里可能是错的。当我输入m=200而不是m=300时，它可以完美地工作。

例如，以下程序不起作用

    #include <iostream>
    #include <stdlib.h>

    extern "C"
    {
    #include <cblas.h>
    }



    using namespace std;

    int main()
    {
        double *a,*x, *y, *z;
        int m,k;
        m=300; k=200;

        a = (double *) malloc( m*k*sizeof( double ));
        x = (double *) malloc (k*sizeof(double));
        y = (double *) malloc (m*sizeof(double));
        z = (double *) malloc (m*sizeof(double));

        int i;

        for (i = 0; i < (m*k); ++i)
        {
            a[i] = 1;
        }

        for (i = 0; i < (k); ++i)
        {
            x[i] = 1;
        }

        for (i = 0; i < (m); ++i)
        {
            y[i] = 100 ;
        }

        cblas_dcopy(m,y,1,z,1);
        cblas_dgemv(CblasRowMajor,CblasNoTrans,m,k,1.0, a ,m ,x, 1, 1.0, z, 1);


        for (int i = 0; i<m; ++i)
        {
             cout<<z[i]<<endl;
        }

        free (a);
        free (x) ;
        free (y) ;
        free (z) ;

        return 0;
        }

提前谢谢你

Answer 1

我已经解决了这个问题: blas文档中给出的参数LDA (紧跟在矩阵后面的参数)对应于他们所说的“矩阵的前导维度”。在COLMAJOR系统(由Fortran使用)中，必须将矩阵的行数作为LDA传递，如果在C或C++中使用cblas，则必须将列数作为LDA传递，因为使用的是ROWMAJOR系统。

希望这篇文章能在未来帮助一些人