Eigen求解大型稀疏对称矩阵（Cholesky分解）

参考自Eigen文档
代码如下：

#include <Eigen/Sparse>typedef Eigen::SparseMatrix<double> SpMat;
typedef Eigen::Triplet<double> Trip;int main()
{//设置Triplet列表，该列表可以表示稀疏矩阵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));//用Triplet列表构造稀疏矩阵SpMat mat(3, 3);mat.setFromTriplets(tripletList.begin(), tripletList.end());//使用Cholesky分解矩阵mat(LDLT分解）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;
}

Eigen求解大型稀疏对称矩阵（Cholesky分解）相关推荐

C++实现Cholesky分解
题目: 编制程序求解矩阵 A 的 Cholesky 分解,并用程序求解方程组 Ax=b,其中代码实现: #include <iostream> #include <math.h&g ...
大型稀疏线性方程组求解技术——工业仿真的底层核心
背景在工业仿真领域,对各种现实世界的问题进行数值模拟时,如流体动力学分析.电磁场仿真.结构力学应力应变分析等,其控制方程通常是偏微分方程组,在经过不同方法的隐式离散之后最终都可转化为大型稀疏线性 ...
Eigen密集矩阵求解 1 - 线性代数及矩阵分解
简介这里介绍线性系统的解析,如何进行各种分解计算,如LU,QR,SVD,特征值分解等. 简单线性求解在一个线性系统,常如下表示,其中A,b分别是一个矩阵,需要求x: Ax=bAx \:= \: b ...
三十分钟理解：矩阵Cholesky分解，及其在求解线性方程组、矩阵逆的应用
写一篇关于Cholesky分解的文章,作为学习笔记,尽量一文看懂矩阵Cholesky分解,以及用Cholesky分解来求解对称正定线性方程组,以及求对称正定矩阵的逆的应用. 文章目录直接Choles ...
cholesky 分解加速求解线性方程组
原理求解 A x = b Ax=b Ax=b 若 A A A对称正定,对其分解: A = L L ⊤ A = LL^\top A=LL⊤ 故 L L ⊤ x = b LL^\top x = b LL ...
矩阵的Cholesky分解
原文链接:矩阵的Cholesky分解首先来复习线性代数中几个重要的概念. 1)如果一个复矩阵A = A*(共轭转置),则A称为Hermitian矩阵.(注意,矩阵A转置后仍为其本身,显然A一定是方阵 ...
几种矩阵分解算法： LU分解，Cholesky分解，QR分解，SVD分解，Jordan分解
目录 1.LU分解 2. LDLT分解法 3. Cholesky分解的形式 4. QR分解 5.SVD分解 5.1 SVD与广义逆矩阵 6. Jordan 分解参考文章: ---------我只是搬 ...
Cholesky分解—概率密度分布及累计概率分布（完整代码分享）
目录概述: 分解条件: 代码: 结果: 概述: Cholesky分解是一种分解矩阵的方法, 在线性代数中有重要的应用.Cholesky分解把矩阵分解为一个下三角矩阵以及它的共轭转置矩阵的乘积(那实数 ...
求解大型 COMSOL 模型需要多少内存？
COMSOL Multiphysics能求解多大的模型是我们最常碰到的问题之一.这个问题其实很难直接回答,因此在本篇博客中,我们将讨论内存需求.模型大小.以及用户如何预测在求解大型三维有限元问题时所需 ...

Eigen求解大型稀疏对称矩阵（Cholesky分解）

Eigen求解大型稀疏对称矩阵（Cholesky分解）相关推荐

最新文章

热门文章