在 Matlab mex 文件中使用特征向量

Question

我正在尝试将 Eigen 库用于我想在 Matlab 中制作成 mex 文件的 c++。我将代码分解为下面的基本要素。我需要创建一个返回向量，其长度对应于输入矩阵的行数。

使用下面的代码，我在编译时遇到以下错误：

double y_OUT[nrow] = {};

• 错误 C2057：预期的常量表达式错误
• 错误 C2466：无法分配常量大小为 0 的数组

我不知道为什么我不能分配长度为 nrow 的向量。 通过取消注释特定行并打印 nrow 我验证它实际上包含正确的数字。有没有人可以给出任何指示（双关语）？

     void mexFunction(
             int          nlhs,
             mxArray      *plhs[],
             int          nrhs,
             const mxArray *prhs[]
             )
    {



      double *x_IN;
      int nrow,ncols;

      /* Check for proper number of arguments */
      //...

      x_IN = mxGetPr(prhs[0]);
      nrow = (int)mxGetM(prhs[0]);
      ncols = (int)mxGetN(prhs[0]);

      double y_OUT[nrow] = {};

      MatrixXd x=Map<MatrixXd>(x_IN,nrow,ncols);
      VectorXd Respons=VectorXd::Zero(nrow);

      Map<VectorXd>(y_OUT,nrow)=Respons.array();
      return;
    }

Answer 1

你没有在里面分配任何东西。 double y_OUT[nrow] = {}; 行包含两个错误。如错误消息所述。

1. nrow 不是 constexpr。它不能在编译时进行评估。
2. 通过大括号初始化程序，您建议向量具有0 元素，这也没有意义。

动态分配的正确语法是：double* y_OUT = new double[nrow];。我想这就是你想要的。您使用的语法double y_OUT[nrow] 用于静态分配，换句话说，在编译时必须知道数组的大小。

旁注：当然你需要记住释放你保留的内存。而且问题与IMO主题完全无关=）