24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1171  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专业求调剂 +10 月月鸟木 2026-04-01 10/500 2026-04-02 12:57 by wxiongid
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9 双马尾痞老板2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 085600,材料与化工321分调剂 +8 大馋小子 2026-03-27 10/500 2026-04-02 10:15 by kaiyanjie
[考研] 一志愿厦门大学化学工程(专硕)-数二英二406分-求调剂 +5 厦大化工 2026-04-01 5/250 2026-04-02 10:03 by jp9609
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 342求调剂 +12 Mary Keen 2026-03-28 13/650 2026-04-01 21:02 by 流情牧豪
[考研] 332求调剂 +8 Lyy930824@ 2026-03-29 8/400 2026-04-01 18:40 by 千钧澄玉宇啊
[考研] 349求调剂 +6 吃的不少 2026-04-01 6/300 2026-04-01 17:55 by JYD2011
[考研] 一志愿武理材料工程302调剂环化或化工 +15 Doleres 2026-03-31 16/800 2026-04-01 09:49 by lfj11
[考研] 一志愿 南京航空航天大学 ,080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 求调剂,一志愿北林食品与营养095500,301分,已过六级,有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
[考研] 本2一志愿C9-333分,材料科学与工程,求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +13 辰沐5211314 2026-03-26 13/650 2026-03-31 14:37 by 记事本2026
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 276求调剂 +3 赵久华 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
信息提示
请填处理意见