ichol 中的负支点

如何解决ichol 中的负支点

我在 Octave 中执行的实现有问题。我正在尝试使用命令 ichol(sparse(hilb(n))) 计算 Hilbert 矩阵的不完全 Cholesky 分解，当 n = 15 时，我收到以下错误： ichol:negative pivot 遇到。

我们有希尔伯特矩阵是正定对称的，所以我不知道错误来自哪里。我计算了矩阵A的特征值，第一个结果为负，我认为这是因为特征值问题是一个条件不好的问题，所以我不太关注那个结果。一个有效的解决方案是使用 diagcomp 选项，但我想知道为什么 Cholesky 的不完全分解不起作用？ A 满足获得不完全 Cholesky 分解的条件，在其他帖子中他们提到了关于支配的内容，但我不知道这会如何影响获得分解的时刻。

首先，谢谢。