Fortran Coder

标题: sor迭代什么方式结束循环较好 [打印本页]

作者: H_S 时间: 2014-4-18 13:19
标题: sor迭代什么方式结束循环较好

[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 时间: 2014-4-18 14:12
迭代法的收敛方式很多，一般计算的有：残差、两次迭代结果的RMS，判断有两种，1、小于某个阈值，2、变化很小

作者: H_S 时间: 2014-4-18 18:15

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

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

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

作者: aliouying 时间: 2014-4-19 10:29

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-19 12:17

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

说的太专业了看不懂

作者: H_S 时间: 2014-4-19 12:29
本帖最后由 H_S 于 2014-4-19 19:19 编辑

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)是边界条件

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)