特征稀疏基质决定因素为零

Eigen sparse matrix determinant is zero

本文关键字：决定因素特征更新时间：2023-10-16

我正在尝试计算我使用的稀疏矩阵是正定的。为此，我试图使用Sylvester标准，这意味着领先的未成年人是积极的。

要计算矩阵的决定因素，我正在构造矩阵每个块的sparselu求解器，然后可以为我提供矩阵的决定因素。但是从一定维度开始（左右130*130），我得到的结果是所有决定因素。对于某种截面算法，本特定算法的决定因素简单低于某些阈值。

我对这种机制的搜索没有带来不错的结果。我的矩阵的尺寸约为16k*16k，所有非零元素都位于对角线附近的96个元素上。

eigen中是否实现了任何截断机制，我可以以某种方式控制其阈值吗？

这很可能是由于底流，即，确定因素被计算为小于1.0的大量数字的产物。如果计算0.5周围130个值的乘积，您就靠近可以用单个精密浮子代表的边界。

您可以使用方法logAbsDeterminant和signDeterminant获得有意义的结果。