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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
15815/16返回列表首页上一页111213141516下一页
查看: 12637  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] FileCode能看出啥？ +8 要乐观耀哥 2026-08-10 21/1050 2026-08-12 16:08 by 云上清扬
[基金申请] 长年满屏的广告，版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
关闭
请填处理意见
关闭