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系统。

希望这可以帮助将来的某个人

非方阵

1 个答案: