| 查看: 1071 | 回复: 8 | ||
[求助]
初学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,LINTEGER::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是现成的一个子程序。。。 求大侠教教我。 |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有13人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求fortran教学视频
已经有5人回复- 小弟初学abaqus二次开发,有个小问题，请大家不要笑话
已经有1人回复
- 求助一个问题，初学Fortran
已经有6人回复
- 跪求初学Fortran的方法及资料。。谢谢呀。。。
已经有1人回复
- abaqus焊接模拟子程序中 SUBROUTINE DFLUX(FLUX,……,JLTYP,是什么意思？
已经有4人回复
- Fortran初学者经验分享
已经有3人回复
- 求visual fortran 6.5编译器
已经有7人回复
- 求助fortran大神，我是初学者，帮我分析一下下面的代码
已经有5人回复
- exp中有循环变量，结果出现infinity，怎么处理
已经有9人回复
- 编一个fortran程序
已经有17人回复
- 关于fortran语言，初学不懂程序，请教各位高手解答。
已经有7人回复
- 初学vasp，编译vasplib出错------dlexlib.f
已经有7人回复
- fortran执行出错，可是看不懂错误的问题
已经有9人回复
- 求助学习fortran的书籍
已经有5人回复
- 在文本和屏幕不能同时输出，这是为何啊，求指教
已经有2人回复
- 初学MPI，以下Fortran程序请教
已经有2人回复
- 百度文库里的，转给大家看看，Python的学习
已经有12人回复
- 【求助】初学fortran，
已经有15人回复
- 链接失效，被关闭
已经有9人回复
2楼2015-11-08 23:03:55
3楼2015-11-09 16:29:30
4楼2015-11-09 23:21:46
5楼2015-11-10 08:12:46
6楼2015-11-10 09:53:56
7楼2015-11-10 15:59:06
8楼2015-11-10 16:21:37
9楼2015-11-10 17:29:25