【急】ICCG算法1万阶迭代次数为160次正常吗？
10000阶，非零元素28万。精度10e-9。迭代160次。我看别人的论文
350000阶，非零元素46万，迭代80次。
20000阶，7次（不会吧？！）
360000阶，非零元素400万，迭代63次。
所以感觉我的程序有误，不过最终得到的解是正确的。也就是说可能IC分解错了，没达到降低条件数的目的。
代码查了几遍也没找到错误，1万160次到底正常吗？ 找到两个bug，一个是精度，这个bug导致我的程序是梯度的精度10e-9，而不是b-Ax。另一个是IC分解时稀疏矩阵的遍历上的bug，即同一行的非零元素，列数大的不一定被记录在列数小的元素的前面。现在1万阶迭代40多次，经matlab验证是正常的。100万阶迭代了200多次，这个矩阵太大了就没用matlab验证。
幸好有matlab这东西，要不都摸不着头。这两个bug都不会导致最后的结果出错。尤其第二个bug，含有bug的IC分解后迭代次数确实比直接用CG少了，让人觉得是有效果的。
页:
[1]