| 查看: 1381 | 回复: 13 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】50金币求大家帮忙改一段程序
|
|||
|
有一组数据如下所示，第一列代表istep变化范围（1-5000），第二列ia（1-30000），第三列ib（1-800），现在我想用一个三维数组H（istep,ia,ib)来标记这些数据，同一行数字定义的H=1，其他的H均为0.也就是说当istep=1时，ia=781，ib=455时H=1。我的主要目的就是能够将istep，ia，ib用一个函数联系起来，因为在程序后面要用这个函数作为进一步的判断条件。我写了一个程序： program exam implicit none integer ia,ib,istep integer ::h(5000,30000,800) integer i h=0 open (10,file='test.out') do i=1,99999999 read(10,*,end=2,err=2) istep,ia,ib h(istep,ia,ib)=1 enddo 2 continue close(10) do ia=1,30000 do ib=1,800 do istep=1,5000 h(istep)=h(istep,ia,ib) enddo do ntel=1,nsteps if(h(ntel)/=0) then 进行后面的程序 ----------------- endif enddo enddo enddo end 但是现在的问题是，定义的三维数组h太大了，程序不能运行， 请大家帮忙看看应该怎么办？先谢了 数据： 1 758 22 1 820 245 1 2050 399 1 2587 232 1 15733 34 1 18193 65 1 20718 256 2 820 245 2 2050 399 2 2587 232 2 20251 125 3 758 22 3 820 245 3 1817 625 3 2050 399 3 16723 555 3 18193 65 3 20718 256 --- ---- ---- |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有22人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【急】请大家帮忙看一下这个offer是否合乎要求，新手散光所有金币求助~~
已经有12人回复
- 帮忙修改一段我的翻译，好难翻啊555
已经有4人回复
- 也请高人指点，帮忙看看每位评委的意见是什么？（50金币，回帖1金币）
已经有84人回复
- 帮忙修改一段翻译！谢谢！
已经有3人回复
- 求高手帮忙改一下英文摘要
已经有12人回复
- 【求助】撒点金币请大家帮忙出主意，回国都带什么礼物啊
已经有20人回复
- 急求，30分钟内翻译一段文字（药学），30分钟内完成6金币再追加50金币
已经有1人回复
- 再求一段英译中，在线等，急，50金币
已经有5人回复
- 【求助】10金币请帮忙修改一下MATLAB程序
已经有9人回复
» 抢金币啦！回帖就可以得到:
-
坐标济南，来碰碰运气
+1/474
-
成都理工大学全国重点实验室公开诚聘绿色有机合成方向联培生及科研助理
+1/78
-
内蒙古大学能源材料化学研究院招聘2026年博士生
+1/76
-
中科院长春光机所 招收计算材料学博士/硕士研究生（含机器学习辅助材料设计方向）
+1/75
-
希望你在这里
+1/64
-
香港理工大学-应用生物与化学科技学系 招收2025年博士研究生
+2/50
-
中国科学院深海所 招收2026秋入学博士生1名 申请-考核制
+1/43
-
急招碳材料相关特任研究人员/博士后/科研助理/26级博士和硕士
+1/43
-
华中科技大学龚江研究员课题组诚招博士研究生、科研助理和博士后
+2/34
-
厦门大学航空航天学院智能制造课题组招2026年申请审核制博士生1-2名
+1/33
-
SCI，计算机相关可以写
+1/28
-
征女友 @长安
+1/24
-
长江学者团队招聘药学/生物信息学等方向高校教师7名(地点杭州、有事业编)+博后5名
+1/10
-
求博导收留
+1/5
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/4
-
中科院动物所招收2026年博士生（优先少干专项计划、化学或生命科学背景）
+1/4
-
电子科技大学激光与光子制造团队招硕士博士
+1/4
-
中国矿业大学黄赳课题组联合中国科学院南京土壤研究所朱晓芳研究员诚聘博士后
+1/2
-
海南大学化学院—功能分子器件团队2026博士/研究助理招生
+1/2
-
求博导收留
+1/1
|
可能是我的表达有问题。 我的数据： 1 758 22 1 820 245 1 2050 399 1 2587 232 1 15733 34 1 18193 65 1 20718 256 2 820 245 2 2050 399 2 2587 232 2 20251 125 3 758 22 3 820 245 3 1817 625 3 2050 399 3 16723 555 3 18193 65 3 20718 256 现有一数据文件test.out,文件结构：第一列istep，第二列ia，第三列ib。istep的变化范围是从1到5000，ia的变化范围是从1到30000，ib的变化范围是从1到800。定义一个三维数组h（5000,30000,800），当h的三维取文件中一行数据时，h=1，否者h=0。我的想法是：从test.out文件读入数据，每读一行（即istep，ia，ib），就将相应的h（istep，ia，ib）赋值为1.相应的程序我是这样写的： program exam implicit none integer ia,ib,istep integer ::h(5000,30000,800) integer i h=0 open (10,file='test.out') do i=1,99999999 read(10,*,end=2,err=2) istep,ia,ib h(istep,ia,ib)=1 enddo 2 continue close(10) 当读完文件，就有了三维数组h中每个元素的值（0或1），然后用h是否等于0来做判断条件，进行下面的计算。相应的程序为 do istep=1,5000 do ia=1,30000 do ib=1,800 if (h(istep,ia,ib)/=0) then 进行计算--- endif enddo enddo enddo 但是如果像我这样定义三位数组的话，需要的内存太大，程序不运行，所以想请大家帮忙看看，如何修改程序，才能使程序执行判断条件后面的计算。举个例子，当istep=1,ia,ib只有取文件中的值也就是下面这些数据对应h才等于1，可以进行后面的计算。取其他的值时h均为0，不进行后面的计算。 1 758 22 1 820 245 1 2050 399 1 2587 232 1 15733 34 1 18193 65 1 20718 256 |
9楼2011-03-31 01:08:16
4楼2011-03-28 23:54:35
|
以我的第一行数据 1 758 22 为例，1代表istep，758代表ia，22代表ib。 我想用函数将istep，ia，ib联系起来，只有当取我数据中一行上的数字时，定义的函数h=1，其他的h均为0.以第一行为例， istep=1，ia= 758，ib=22时，h=1.或者用多个函数也行，只要能完成后面循环中判断（ if(h(ntel)/=0) then）就行。 do ia=1,30000 do ib=1,800 do istep=1,5000 h(istep)=h(istep,ia,ib) enddo do ntel=1,nsteps if(h(ntel)/=0) then 进行后面的程序 ----------------- endif enddo enddo enddo |
5楼2011-03-29 03:05:47
6楼2011-03-29 15:16:15