使用指针在存在特征库的情况下动态分配 c++ 中的矩阵

Using a pointer to dynamically allocate Matrices in c++ in presence of eigen library

本文关键字：动态分配 情况下 c++ 指针 存在 特征      更新时间：2023-10-16

我使用 Eigen 库来初始化矩阵并进行与对角化等相关的计算。到目前为止，我一直在使用"按值调用"方法，即：

  1. 我用MatrixXcd H(N，N(定义了一个复杂的双NxN矩阵。
  2. 然后我调用一个例程，如 SelfAdjointMatrixSolver es(H(。
  3. 最后，我使用 es.eigenvalues(( 获得特征值。

现在，我想使用指针动态分配矩阵。我写了以下代码：

#include <iostream>
#include <complex>
#include <cmath>
#include<math.h>
#include<stdlib.h>
#include<time.h>
#include<Eigen/Dense>
#include<fstream>
#include<random>
using namespace std;
using namespace Eigen;

int main()
{
int i,j,k,l;//for loops
MatrixXd *H = NULL;
H = new MatrixXd(10,10);
if (!H) 
cout << "allocation of memory failedn"; 
else
{
for(i=0;i<10;i++)
{
for(j=0;j<10;j++)
{
H[i][j]=1.0;
}
}
}
return 0;
}

我收到以下错误消息：

In file included from eigen/Eigen/Core:347:0,
from eigen/Eigen/Dense:1,
from test.cpp:7:
eigen/Eigen/src/Core/DenseCoeffsBase.h: In instantiation of ‘Eigen::DenseCoeffsBase<Derived, 1>::Scalar& Eigen::DenseCoeffsBase<Derived, 1>::operator[](Eigen::Index) [with Derived = Eigen::Matrix<double, -1, -1>; Eigen::DenseCoeffsBase<Derived, 1>::Scalar = double; Eigen::Index = long int]’:
test.cpp:32:13:   required from here
eigen/Eigen/src/Core/util/StaticAssert.h:32:40: error: static assertion failed: THE_BRACKET_OPERATOR_IS_ONLY_FOR_VECTORS__USE_THE_PARENTHESIS_OPERATOR_INSTEAD
#define EIGEN_STATIC_ASSERT(X,MSG) static_assert(X,#MSG);
^
eigen/Eigen/src/Core/DenseCoeffsBase.h:406:7: note: in expansion of macro ‘EIGEN_STATIC_ASSERT’
EIGEN_STATIC_ASSERT(Derived::IsVectorAtCompileTime,

我应该如何解决这个问题？我的错误是在定义指针还是分配值时？

正如这个断言所说THE_BRACKET_OPERATOR_IS_ONLY_FOR_VECTORS__USE_THE_PARENTHESIS_OPERATOR_INSTEAD矩阵类使用括号来索引，而不是文档中提到的方括号。所以改变

H[i][j] = 1.0;

自

H(i,j) = 1.0;

此外，矩阵已经动态分配，这里不需要new。

MatrixXd H(10, 10);
相关文章：