24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1248  |  回复: 8

444874095

新虫 (初入文坛)

[求助] 初学fortran,运行后出现this name does not match the unit name

我尝试运行了一个高斯赛德尔迭代
代码如下:
PROGRAM MAIN
IMPLICIT NONE
DIMENSION A(4,4),B(4),X(4)
DOUBLE PRECISION A,B,X
REAL:PS,EPS,L
INTEGER::I
DATA A/7.0,9.0,-2.0,1.0,2.0,15.0,-2.0,3.0,1.0,3.0,11.0,2.0,-2.0,-2.0,5.0,13.0/  
DATA B/4.0,7.0,-1.0,0.0/
DPS=1.0E-06
CALL AGSDL(A,B,4,X,EPS,L)
IF (L.NE.0) THEN
  WRITE(*,10) (X(I),I=1,4)
END IF
10 FORMAT(1X,D15.6)
end
CONTAINS
SUBROUTINE AGSDL(A,B,N,X,EPS,L)
DIMENSION A(N,N),B(N),X(N)
DOUBLE PRECISION A,B,X,T,S,P,Q
INTEGER::N
DO 5 I=1,N
  IF (ABS(A(I,I))+1.0.EQ.1.0) THEN
L=0
WRITE(*,100)
RETURN
  END IF
5  CONTINUE
100  FORMAT(1X,'FAIL')
L=100
DO 10 I=1,N
10   X(I)=0.0
20   P=0.0
L=L-1
DO 50 I=1,N
  T=X(I)
  S=0.0
  DO 30 J=1,N
IF(J.NE.I) S=S+A(I,J)*X(J)
30  CONTINUE
      X(I)=(B(I)-S)/A(I,I)
      Q=ABS(X(I)-T)/(1+ABS(X(I)))
      IF (Q.GT.P) P=Q
50  CONTINUE
      IF((P.GE.EPS).AND.(L.NE.0)) GOTO 20
      IF(L.EQ.0) WRITE(*,100)
      RETURN
END SUBROUTINE AGSDL
END PROGRAM MAIN
其中的引用的子例子程序agsdl是现成的一个子程序。。。
求大侠教教我。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444874095

新虫 (初入文坛)

那个笑脸那里是一个D。。。是声明DPS是real型的。不知怎么变笑脸了。。。大侠请忽略这个细节。。
2楼2015-11-08 23:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444874095

新虫 (初入文坛)

3楼2015-11-09 16:29:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444874095

新虫 (初入文坛)

4楼2015-11-09 23:21:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

指环王jyj

新虫 (小有名气)

还在用77语言……那个double precision命令啥意思?

发自小木虫Android客户端
5楼2015-11-10 08:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444874095

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by 指环王jyj at 2015-11-10 08:12:46
还在用77语言……那个double precision命令啥意思?

双精度数值类型

发自小木虫Android客户端
6楼2015-11-10 09:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kimjungil

铁杆木虫 (小有名气)

将其中一个循环的unit=5改为55试试。

发自小木虫Android客户端
7楼2015-11-10 15:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444874095

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by kimjungil at 2015-11-10 15:59:06
将其中一个循环的unit=5改为55试试。

还是一样的错误。。这样改有什么依据吗?

发自小木虫Android客户端
8楼2015-11-10 16:21:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

444874095

新虫 (初入文坛)

9楼2015-11-10 17:29:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 444874095 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 上海大学实验技术岗位非升即走 +8 嘻嘻哈哈乐呵呵 2026-05-15 8/400 2026-05-18 11:53 by ning312
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[考博] 现在不知道怎么办,感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
信息提示
请填处理意见