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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
15811/16返回列表首页上一页91011121314下一页末页
查看: 12519  |  回复: 157
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

jihe01

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
微尘、梦想(金币+2): 谢谢回贴交流…… 2011-04-19 15:13:08
看了maomao1210开的这个帖子，很受启发，对Fortran学习帮助很大。有个程序问题请教一下，这是网上一个产生随机数的程序，用95编写的，当我在调用FUNCTION时（用6.6版本），出现Access volidation错误，问题出在11 continue后面一行  IDUM=1.0上，但我不明白什么原因，请帮忙看一下，谢谢！


FUNCTION RAN1(IDUM)


      DIMENSION R(97)                                         
      PARAMETER(M1=259200,IA1=7141,IC1=54773,RM1=1./M1)
      PARAMETER(M2=134456,IA2=8121,IC2=28411,RM2=1./M2)
      PARAMETER(M3=243000,IA3=4561,IC3=51349)
      
          DATA IFF/0/

      IF(IDUM.LT.0.OR.IFF.EQ.0) THEN
        IFF=1
        IX1=MOD(IC1-IDUM,M1)
        IX1=MOD(IA1*IX1+IC1,M1)
        IX2=MOD(IX1,M2)
        IX1=MOD(IA1*IX1+IC1,M1)
        IX3=MOD(IX1,M3)

        DO 11 J=1,97
        IX1=MOD(IA1*IX1+IC1,M1)
        IX2=MOD(IA2*IX2+IC2,M2)
        R(J)=(real(IX1)+real(IX2)*RM2)*RM1
   11   CONTINUE
        IDUM=1.0
      ENDIF

      IX1=MOD(IA1*IX1+IC1,M1)
      IX2=MOD(IA2*IX2+IC2,M2)
      IX3=MOD(IA3*IX3+IC3,M3)
      J=1+(97*IX3)/M3
      IF(J.GT.97.OR.J.LT.1)PAUSE
      RAN1=R(J)
      R(J)=(real(IX1)+real(IX2)*RM2)*RM1

      RETURN
      END function
赞一下(2人) 回复此楼
101楼2011-04-11 15:24:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★ ★
微尘、梦想(金币+2): 谢谢回复！ 2011-04-19 15:13:35
引用回帖:
Originally posted by jihe01 at 2011-04-11 15:24:00:
看了maomao1210开的这个帖子，很受启发，对Fortran学习帮助很大。有个程序问题请教一下，这是网上一个产生随机数的程序，用95编写的，当我在调用FUNCTION时（用6.6版本），出现Access volidation错误，问题出在11 ...

我刚测试了一下（用的是你说的vf6.6），结果一切都不错，估计程序不是出现这个函数中吧，是你操作有问题，或者其他模块的程序有问题。打开资源管理器，把相应的exe关掉，再运行一下，看是否还会出现类似的错误。
赞一下(1人) 回复此楼
102楼2011-04-11 15:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jihe01

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
dubo(金币+1): 欢迎讨论 2011-08-27 16:39:20
谢谢maomao的测试，可能确实不是这个function的问题，可能是其他程序调用的问题，这个程序有点复杂，我是从网上下的，是关于MCMC方法计算的，当我调试时，总是到这个地方出现措施（ＩＤＵＭ＝１）。附件是原程序和错误提示，请看一下，谢谢。
赞一下(2人) 回复此楼
103楼2011-04-11 19:57:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qwer_2009

铜虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 欢迎常来程序语言版，鼓励新虫多发话题！ 2011-04-19 15:15:09
program main
implicit none
type node
real x   !设置随机数
type(node),pointer::next
end type
type(node),pointer::head,q,p
real rand
integer i,j
integer::m=0,n=0


call random_seed( )
nullify(head)


do i=1,10
call random_number(rand)
allocate(q)
q%x=rand
q%next=>head
head=>q
end do

q=>head
do while(associated(q))
print*,q%x
q=>q%next
end do

!判断是否大于0.50
q=>head

do while(associated(q))


if(q%x<0.50.and.associated(q,head))then
n=n+1  !判断小于0.5的个数
head=>q%next
deallocate(q)
q=>head
else if(q%x<0.50.and..not.associated(q,head))then
n=n+1
p=>q%next
deallocate(q)
q=>p
else
m=m+1
q=>q%next
end if

end do


write(*,*) "大于0.5的个数m=", m
write(*,*) "小于0.5的个数n=", n

!输出大于0.5的数
q=>head
do while(associated(q))

print*,q%x
q=>q%next
end do

end

请大家帮忙看看这个程序问题出在哪里了，在输出m和n的值时还是正确的，为什么加上最后几句 显示大于0.5的数时就会无限循环下去，而且得到的结果还是错误的呢，谢谢
赞一下(2人) 回复此楼
104楼2011-04-19 11:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abc__1203

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
可否提供一段Fortran调用MATLAB函数的实例程序
例如想调用MATLAB里边的besselj函数
谢谢！
赞一下(1人) 回复此楼
105楼2011-05-16 21:29:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suxingtao

铁杆木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-08-27 16:40:22
我有个Fortran插值问题，http://muchong.com/bbs/viewthread.php?tid=3305448&fpage=1，向您请教！
赞一下(2人) 回复此楼
106楼2011-06-16 14:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿水哥王

铜虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-08-27 16:40:11
楼主你好  我目前在算一个球星量子点矩阵的本征值和本征矢  我想问下这个矩阵元的积分里怎样编译这个沿径向函数的积分元
赞一下(2人) 回复此楼
沉默越久 爆发越大
107楼2011-06-20 15:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shaofain

铁虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
jjdg(金币+1): 感谢参与 2011-08-01 17:02:39
楼主你好：能不能推荐一个关于fortran工程算法的网站或是书？谢谢
赞一下(2人) 回复此楼
108楼2011-08-01 15:24:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anshenyu

铜虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 欢迎讨论 2011-08-27 16:39:57
请教下fortran能不能处理这样的矩阵代数方程?
P=(I-PH^TH/r)(APA^T+BB^Tq)
其中大写字母均代表矩阵
赞一下(2人) 回复此楼
109楼2011-08-25 16:19:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+1): 鼓励发帖交流~ 2011-09-07 11:54:33
本帖仅楼主可见
赞一下
110楼2011-09-05 09:23:29
已阅   申请程序强帖   回复此楼   编辑   查看我的主页
相关版块跳转 我要订阅楼主 maomao1210 的主题更新
15811/16返回列表首页上一页91011121314下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 bjm2e4yior 2026-07-18 3/150 2026-07-19 16:43 by 0u3sh5hhfi
[考博] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 16:10 by 0u3sh5hhfi
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:03 by 0u3sh5hhfi
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:03 by 0u3sh5hhfi
[基金申请] 国社科通讯评审 +5 潜水冠军123 2026-07-18 6/300 2026-07-19 12:28 by lfy8008
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +14 botar 2026-07-14 15/750 2026-07-19 00:15 by beefly
[基金申请] 5个%2F，非固定段22个字母 +8 jaasxulu 2026-07-15 9/450 2026-07-18 17:40 by dxcgbm
[基金申请] 生命口C13面上会评时间 +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by 虫卜1988
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 投英文期刊如何查重 10+5 不知道取什么名� 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
信息提示
关闭
请填处理意见
关闭