sor迭代什么方式结束循环较好
dodo 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迭代求解一个较大的矩阵
不知道怎样结束循环
我上面用的总感觉不是很好
迭代法的收敛方式很多,一般计算的有: 残差、两次迭代结果的RMS,判断有两种,1、小于某个阈值,2、变化很小 aliouying 发表于 2014-4-18 14:12
迭代法的收敛方式很多,一般计算的有: 残差、两次迭代结果的RMS,判断有两种,1、小于某个阈值,2、变化很 ...
能不能具体点呢
我上面的就是比较了迭代之前与之后其中一个值的大小
但是矩阵中有许多值 可能恰好比较了的那个满足 就跳出了
这样结果就不精确了 可以用初始荷载的某种范数与求解得结果与失衡力同样的范数比值 作为收敛判断标准 H_S 发表于 2014-4-18 18:15
能不能具体点呢
我上面的就是比较了迭代之前与之后其中一个值的大小
但是矩阵中有许多值 可能恰好比较了 ...
你这代码有问题,为什么是判断某一个值?
我说的残差或者解的RMS都是指所有的解,即对应所有的Ek
所以判断迭代收敛与否,是在所有的k循环后计算得到Ek,然后与Ek0(即上一次的结果)做对比
PS:为什么内层循环是k=2,n-1?Ek(1)和Ek(n)不变么? 珊瑚虫 发表于 2014-4-18 18:20
可以用初始荷载的某种范数与求解得结果与失衡力同样的范数比值 作为收敛判断标准 ...
说的太专业了 看不懂:-dizzy: 本帖最后由 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)是边界条件
页:
[1]