如何声明一组特征稀疏矩阵

Question

我想知道以下特征稀疏矩阵的声明是否正确

#include <iomanip>
#include <iostream>
#include <sstream>
#include <fstream>
#include <cstdlib>
#include <algorithm>
#include "Eigen/Core"
#include "Eigen/LU"
#include "Eigen/Sparse"
#include "Eigen/StdVector"
#include <chrono>
#include <thread>
#include <mutex>

using namespace std;
using namespace Eigen;

#define NROW 20 
#define NODES 20000

vector<SparseMatrix<double>> Kmat(NROW, SparseMatrix<double> (NODES,NODES));

每个 NROW 的 Kmat 将被反转，我使用的是 SimplicialLLT cholesky 分解方法，这就是我声明它的方式

 static vector<SimplicialLLT<SparseMatrix<double>>> Kmat_LLT(Kmat.size());

基本上我想存储对应于 NROW 的每个逆矩阵，以便以后使用。 想知道上面Kmat_LLT的声明是否正确。 代码确实编译并运行。