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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
15815/16返回列表首页上一页111213141516下一页
查看: 13160  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗？ 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[职场人生] 40+中年人失业后，还能找到好工作吗？ +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
信息提示
关闭
请填处理意见
关闭