有谁用过ORHTOMIN解线性方程？请进

moujianye

我用ORTHOMIN解压力泊松方程，很不稳定，当更新速度后，有时不能收敛，请教是什么原因？谢谢！

dadeng

大概是程序编错了，或者是没有用好预条件。

caiqd

orthomin方法中有一个参数，即保存的搜索方向的数目，调整一下这个参数，换不同的数目试一试，有时这个数太小，收敛可能很慢。

moujianye

在ORTHOMIN中，为什么不能把前一次的结果用做下一次调用ORTHOMIN的初值呢？在算法中有这么一条语句：
for (IB = 0; IB <= IJKM - 1; IB++)
{
　DP[IB] = 0;
　Rsq = Rsq + RN[IB] * RN[IB];
}
即每次调用ORTHOMIN时，初值为０。请问怎样能把前一次的结果用做下一次调用ORTHOMIN的初值。谢谢！

moujianye

我用其他方法解过此方程组，结果是正确的，所以问题出在ORTHOMIN中。