应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
北京石油化工学院2026年研究生招生接收调剂公告
15815/16返回列表首页上一页111213141516下一页
查看: 11518  |  回复: 157
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

568286346

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我想请问在centos6.3linux系统中如何安装ifort和icc啊？
赞一下 回复此楼
坚持才能实现自己的理想！
141楼2012-11-22 17:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

未央八月

铁虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送鲜花一朵
您好，问一下关于编译fortran77程序的问题：
f77 ***.f  $lplt -g -o ***
楼主能普及一下-g -o及$lplt的含义吗？
谢啦！
赞一下 回复此楼
142楼2012-12-02 02:05:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhg19

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by maomao1210 at 2009-06-01 12:41:49
哦。我给你举个例子吧，比如对称矩阵 A
                              A(1,1)     A(1,2)     A(1,3)
                              A(2,1)     A(2,2)     A(2,3)
                              A(3,1)      ...

是否应该II=MAX（I,J）*（MAX（I,J）-1）/2+MIN（I,J）

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
143楼2013-07-09 13:34:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我有个疑问想请教下，
不是说单精度的表示范围：10e-37到10e37吗，那我为什么在fortran里给单精度变量赋值 a=2.1*10**(-15）写出来的a却是0呢，或者当a当分母的时候结果是NAN。这两个情况是不是矛盾了呀？
赞一下 回复此楼
144楼2013-07-09 16:23:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
引用回帖:
144楼: Originally posted by btx97 at 2013-07-09 16:23:19
我有个疑问想请教下，
不是说单精度的表示范围：10e-37到10e37吗，那我为什么在fortran里给单精度变量赋值 a=2.1*10**(-15）写出来的a却是0呢，或者当a当分母的时候结果是NAN。这两个情况是不是矛盾了呀？

10**(-1) = 0
10**(-2) =0
.....
....
....
...
10**(-15) = 0
-15是整数，这个肯定是零的呀。
可以修改为 10**（-15.0）
赞一下 回复此楼
145楼2013-07-09 19:01:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
引用回帖:
144楼: Originally posted by btx97 at 2013-07-09 16:23:19
我有个疑问想请教下，
不是说单精度的表示范围：10e-37到10e37吗，那我为什么在fortran里给单精度变量赋值 a=2.1*10**(-15）写出来的a却是0呢，或者当a当分母的时候结果是NAN。这两个情况是不是矛盾了呀？

10**(-1) = 0
10**(-2) =0
.....
....
....
...
10**(-15) = 0
-15是整数，这个肯定是零的呀。
可以修改为 10**（-15.0）
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

btx97
146楼2013-07-09 19:01:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

btx97

金虫 (小有名气)
送红花一朵
引用回帖:
146楼: Originally posted by maomao1210 at 2013-07-09 12:01:38
10**(-1) = 0
10**(-2) =0
.....
....
....
...
10**(-15) = 0
-15是整数，这个肯定是零的呀。
可以修改为 10**（-15.0）...

谢谢！
回复此楼
147楼2013-07-11 15:22:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yydwhu

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
请问下,3维数组.10801*5401*66,好象编译不了,怎么处理?
赞一下 回复此楼
148楼2013-07-28 23:06:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

血汗泪

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
前辈您好，我写了一个计算复数矩阵行列式的小程序，出现错误了，望前辈指点！
        SUBROUTINE CDET(A,N,DET)
        USE MSIMSL
      IMPLICIT INTEGER*4 (I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
        COMPLEX*16 A(N,N),DET
        DIMENSION RAR(N,N),RAI(N,N)
        REAL*8 SS
!C=====================================================================
        RAR=DREAL(A)
        RAI=DIMAG(A)
        F  =1.0
        DET=1.0
!C=====================================================================
        DO L=1,N-1
                RD=0.0
                DO I=L,N
!C=====================================================================
                        DO J=L,N
!C---------------------------------------------------------------------
                                RM=RAR(I,J)*RAR(I,J)+RAI(I,J)*RAI(I,J)
!C---------------------------------------------------------------------
                                IF (RM.GT.RD) THEN
                                        RD=RM
                                        IS=I
                                        JS=J
                                END IF
!C---------------------------------------------------------------------
                        ENDDO
                ENDDO
!C=====================================================================
                IF (RD+1.0.EQ.1.0) THEN
                        DET=0.0
                        RETURN
                END IF
!C=====================================================================
                IF (IS.NE.L) THEN
                        F=-F
!C---------------------------------------------------------------------
                        DO J=L,N
                                P=RAR(L,J)
                                RAR(L,J)=RAR(IS,J)
                                RAR(IS,J)=P
                                P=RAI(L,J)
                                RAI(L,J)=RAI(IS,J)
                                RAI(IS,J)=P
                        ENDDO
!C---------------------------------------------------------------------
                ENDIF
!C=====================================================================
                IF (JS.NE.L) THEN
                        F=-F
!C---------------------------------------------------------------------
                        DO I=L,N
                                P=RAR(I,L)
                                RAR(I,L)=RAR(I,JS)
                                RAR(I,JS)=P
                                P=RAI(I,L)
                                RAI(I,L)=RAI(I,JS)
                                RAI(I,JS)=P
                        ENDDO
!C---------------------------------------------------------------------
                END IF
                DET=DET*DCMPLX(RAR(L,L),RAI(L,L))
!C=====================================================================
                DO J=L+1,N
                        P=RAR(L,J)*RAR(L,L)
                        Q=-RAI(L,J)*RAI(L,L)
                        SS=(RAR(L,L)-RAI(L,L))*(RAR(L,J)+RAI(L,J))
                        RAR(L,J)=(P-Q)/RD
                        RAI(L,J)=(SS-P-Q)/RD
                ENDDO
                DO I=L+1,N
                        DO J=L+1,N
                                P=RAR(I,L)*RAR(L,J)
                                Q=RAI(I,L)*RAI(L,J)
                                SS=(RAR(I,L)+RAI(I,L))*(RAR(L,J)+RAI(L,J))
                                RAR(I,J)=RAR(I,J)-P+Q
                                RAI(I,J)=RAI(I,J)-SS+P+Q
                        ENDDO
                ENDDO
!C---------------------------------------------------------------------
        ENDDO
!C=====================================================================
        DET=F*DET*DCMPLX(RAR(N,N),RAI(N,N))
        RETURN
        END SUBROUTINE
!C=====================================================================
错误为Error: The shapes of the array expressions do not conform.   [DET]
请您看一下，万分感谢！
赞一下 回复此楼
149楼2013-07-29 17:04:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leixie

铜虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by snoopyzhao at 2009-06-16 09:42:00
个人认为，如果选择现成的子程序的话，一定要优先考虑 blas  或 lapack 等经过历史检验的，其次才是徐士良、NC 等……...

Fortran菜鸟长见识了。哈哈。
赞一下 回复此楼
真正的平静，不是避开车马喧嚣，而是在心中修篱种菊.
150楼2013-11-07 11:37:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maomao1210 的主题更新
15815/16返回列表首页上一页111213141516下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336求调剂，一志愿中科大 +4 墨彧 yuyu 2026-04-06 4/200 2026-04-06 20:03 by lijunpoly
[考研] 286求调剂 +12 Faune 2026-04-06 12/600 2026-04-06 17:40 by dongzh2009
[考研] 288环境专硕,求调材料方向 +16 lllllos 2026-04-04 18/900 2026-04-06 16:41 by 蓝云思雨
[考研] 求调剂 +5 wos666 2026-04-03 5/250 2026-04-06 10:13 by 蓝云思雨
[考研] 求调剂 +5 chenxrlkx 2026-04-05 7/350 2026-04-06 07:54 by houyaoxu
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 341求调剂 +3 洛多罗 2026-04-02 4/200 2026-04-04 21:36 by 智能智慧
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 化工调剂303分，过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
[考研] 材料工程322分 +8 哈哈哈吼吼吼哈 2026-04-01 8/400 2026-04-02 11:53 by 3041
[考研] 一志愿9初试366 本双非求调剂 +4 运气来得若有似� 2026-04-02 4/200 2026-04-02 09:56 by guanxin1001
[考研] 285求调剂 +11 AZMK 2026-04-01 11/550 2026-04-01 22:40 by peike
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
信息提示
关闭
请填处理意见
关闭