参考自Eigen文档
代码如下:
#include <Eigen/Sparse>
typedef Eigen::SparseMatrix<double> SpMat;
typedef Eigen::Triplet<double> Trip;
int main()
{
std::vector<Trip> tripletList(9);
tripletList.push_back(Trip(0, 0, 10));
tripletList.push_back(Trip(0, 1, -2));
tripletList.push_back(Trip(0, 2, -2));
tripletList.push_back(Trip(1, 0, -2));
tripletList.push_back(Trip(1, 1, 10));
tripletList.push_back(Trip(1, 2, -1));
tripletList.push_back(Trip(2, 0, -2));
tripletList.push_back(Trip(2, 1, -1));
tripletList.push_back(Trip(2, 2, 3));
SpMat mat(3, 3);
mat.setFromTriplets(tripletList.begin(), tripletList.end());
Eigen::SimplicialLDLT<SpMat> chol(mat);
Eigen::VectorXd b(3);
b.coeffRef(0) = 1.0;
b.coeffRef(1) = 0.5;
b.coeffRef(2) = 1.0;
Eigen::VectorXd x = chol.solve(b);
for (int i = 0; i < 3; i++)
{
std::cout << x.coeffRef(i) << std::endl;
}
return 0;
}
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)