问题描述
我在 Octave 中执行的实现有问题。我正在尝试使用命令 ichol(sparse(hilb(n))) 计算 Hilbert 矩阵的不完全 Cholesky 分解，当 n = 15 时，我收到以下错误： ichol:negative pivot 遇到。
我们有希尔伯特矩阵是正定对称的，所以我不知道错误来自哪里。我计算了矩阵A的特征值，第一个结果为负，我认为这是因为特征值问题是一个条件不好的问题，所以我不太关注那个结果。一个有效的解决方案是使用 diagcomp 选项，但我想知道为什么 Cholesky 的不完全分解不起作用？ A 满足获得不完全 Cholesky 分解的条件，在其他帖子中他们提到了关于支配的内容，但我不知道这会如何影响获得分解的时刻。
首先，谢谢。
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）