高斯牛顿法，迭代无法终止（死循环）

jz_cheng · 发表于 2014-10-17 11:12:45

在用高斯牛顿法计算数据时，当输入的数据是连续时（比如2、3、4、5、6）程序能跳出循环，但是当输入数据不连续时，会出现死循环情况。由于本人能力有限，无法找出其中的原因，希望那位大侠能够帮我指出。由于代码比较长，所以我就打包上传。
本人使用广义逆修正。

jz_cheng · 发表于 2014-10-17 12:46:10

这是我运行的结果

jz_cheng · 发表于 2014-10-17 13:31:21

楚香饭发表于 2014-10-17 13:09
你用的什么编译器？我试过 ivf，cvf，gfortran，都无法跳出循环。

谢谢，我用的是ivf

jz_cheng · 发表于 2014-10-17 13:38:46

楚香饭发表于 2014-10-17 13:29
OK ，改了以下代码就可以了

[mw_shl_code=fortran,true]REAL*8::G(row,col),Gt(col,row),GtG(col,col)=0. ...

谢谢你，但是初始化之后，输入数据还是不能

[其他行业算法] 高斯牛顿法，迭代无法终止（死循环）