UMFPACK学习（一）

最新推荐文章于 2022-05-18 17:26:44 发布

weixin_34137799

最新推荐文章于 2022-05-18 17:26:44 发布

阅读量1.5k

点赞数 1

文章标签： python matlab

原文链接：http://www.cnblogs.com/silver1116/archive/2012/08/04/2623083.html

版权

本文介绍了UMFPACK库，用于求解稀疏矩阵方程，特别是非对称矩阵。UMFPACK基于非对称多波前算法，支持LU分解。文章详细讲解了如何在C程序中使用UMFPACK，包括symbolic和numeric分解，以及solve、free_symbolic和free_numeric五个步骤，并给出了稀疏矩阵的compressed column form表示方法和一个小示例。

摘要由CSDN通过智能技术生成

引子

最近毕设论文中需要求解如下的稀疏矩阵方程：Ax=b，

现在采用Eigen中的SparseLDLT进行求解，无奈速度有点慢。所以开始寻找其他的求解库，找到了UMFPACK。

UMFPACK就是求解类似于Ax=b这样问题的一个库，来自佛罗里达州立大学。可以直接到http://www.cise.ufl.edu/research/sparse/umfpack/去下载对应的包然后编译得到Windows下的lib，包含到自己的工程就可以了。今天下午看了看UMFPACK的guide ，稍微总结一下。

什么是UMFPACK

UMFPACK是专门求解类似于Ax=b这样的稀疏矩阵方程的一个库，一般情况下A 是稀疏非对称的矩阵。它是基于非对称的多波前算法（ Unsymmetric-pattern MultiFrontal method ）对稀疏矩阵方程进行求解的。UMFPACk可以对PAQ，PRAQ，PR-1AQ 这样的矩阵进行 LU分解（L 和 U 分别是下三角矩阵和上三角矩阵，P和Q是可置换矩阵，R是对角阵）。这里有一个概念：reduce fill-in (减少注入元)，稀疏矩阵的注入元是指执行算法后从初始的零值变为非零值的那些元素。