sor迭代什么方式结束循环较好

H_S

[Fortran] 纯文本查看 复制代码

do
            do k=2,n-1
                Ek0=E(k)
                P0=(P(k)*Dz(k)+P(k-1)*Dz(k-1))/(Dz(k)+Dz(k-1))
                Dz0=(Dz(k)+Dz(k-1))/2D0
                E(k)=E(k)*(1D0-W0)+W0*(E(k+1)/(Dz0*Dz(k))+E(k-1)/(Dz0*Dz(k-1)))/(1.0D0/(Dz(k)*Dz0)+1.0D0/(Dz(k-1)*Dz0)+i*W*U*P0)
                T=(abs(dble(Ek0)-dble(E(k)))+abs(aimag(Ek0)-aimag(E(k))))/2D0
                if(T<1e-45)goto 100
            enddo
        enddo
100     H2=(E(3)-E(2))/(Dz(2)*i*W*U)
        Hs=P(2)*Dz(2)*(E(2)+(E(3)-E(2))/4)/2+H2
        Pa=(E(2)/Hs)**2/(W*U)

如题 作业中
用到sor迭代求解一个较大的矩阵
不知道怎样结束循环
我上面用的总感觉不是很好

aliouying

迭代法的收敛方式很多，一般计算的有： 残差、两次迭代结果的RMS，判断有两种，1、小于某个阈值，2、变化很小

H_S

aliouying 发表于 2014-4-18 14:12
迭代法的收敛方式很多，一般计算的有： 残差、两次迭代结果的RMS，判断有两种，1、小于某个阈值，2、变化很 ...

能不能具体点呢
我上面的就是比较了迭代之前与之后其中一个值的大小
但是矩阵中有许多值 可能恰好比较了的那个满足 就跳出了
这样结果就不精确了

珊瑚虫

可以用初始荷载的某种范数与求解得结果与失衡力同样的范数比值 作为收敛判断标准

aliouying

H_S 发表于 2014-4-18 18:15
能不能具体点呢
我上面的就是比较了迭代之前与之后其中一个值的大小
但是矩阵中有许多值 可能恰好比较了 ...

你这代码有问题，为什么是判断某一个值？

我说的残差或者解的RMS都是指所有的解，即对应所有的Ek
所以判断迭代收敛与否，是在所有的k循环后计算得到Ek，然后与Ek0（即上一次的结果）做对比

PS:为什么内层循环是k=2,n-1？Ek(1)和Ek(n)不变么？

H_S

珊瑚虫 发表于 2014-4-18 18:20
可以用初始荷载的某种范数与求解得结果与失衡力同样的范数比值 作为收敛判断标准 ...

说的太专业了 看不懂

H_S

aliouying 发表于 2014-4-19 10:29
你这代码有问题，为什么是判断某一个值？

我说的残差或者解的RMS都是指所有的解，即对应所有的Ek

是就判断了一个值 觉得不是很好 就来问了

看我同学的是
DO
   E0=E
   SUM=0
   DO k=2,n-1
        E(k)=E(k)+....
   ENDDO
   DO i=1,n
        SUM=SUM+abs(E(i)-E0(k))
   ENDDO
   IF(SUM<1E-15)EXIT
ENDDO

大概就是这样的 凭记忆打的

我就觉得又多了一个数组 好麻烦

E(1)和E(n)是边界条件