一个批量运算的相关问题

kif117 · 发表于 2015-3-30 18:47:38

问题是这样的:

数据样本为:

1    5.402000    -4.100000    6.486000
2    5.106000    10.426000    14.167000
3    8.123000    20.117000    -1.005000
4 -11.142000    17.711000    3.485000
5 -17.692000    8.982000 -10.812000
.....
160 ...

每个文件包含160个ID，一共有几万个那么多。但是不需要全部使用，我单挑选那些文件号为10的倍数的来用，这样需要处理的文件个数为几千个。我希望结果是以ID号作为RESULT_(ID号)文件名，然后里面包含提取出几千个的数据。

我参考论坛里另外一个帖子(数据分类问题求助)，仿照着fcode同学给的例子写了一个:

[Fortran] 纯文本查看 复制代码

01Program main
02       Implicit None
03       Character(160) :: Filename
04       Character(160) :: fn1,fn2,fn3,tmp
05       Integer :: I, J,num
06       real :: A(4), B(4), C(4)
07       integer ios
08 
09       fn1='c6c8_ID_'
10         fn2='c6_'
11         fn3='c8_'
12 
13         num=160  ! total number of result files
14 
15       print*, 'Running'
16       do J = 1, num
17         write (tmp,*)J
18 C        print*, J
19       open (j,File="D:\c6c8\"//trim(adjustl(fn1))//trim(adjustl(tmp))//'.dat',status='unknown')   ! save result here
20         write (j,'(3f8.3)') 0.0,0.0,1.0
21         enddo
22 
23       open (20, File='IDlist.dat',status='old') ! 里面存放的是10的倍数，共计9200个。
24      do I = 1, 9200
25       read (20, *,iostat=ios) Filename
26       if (ios /=0) then
27       exit
28       endif
29 c        print*, filename
30       open (200,File="D:\C6\"//trim(adjustl(fn2)) //trim(adjustl(filename))//'.dat')
31       open (300,File="D:\C8\"//trim(adjustl(fn3)) //trim(adjustl(filename))//'.dat')
32 
33       Do J = 1, num
34       Read (200, *) A(:)  ! first row=ID number， then x,y,z
35         Read (300, *) B(:)  ! first row=ID number
36         C(:)=B(:)-A(:)
37 c        print*, c(:)
38       Write (c(1), '(3f8.3)') C(2:4) ! first row=ID number 
39       Close (300)
40         Close (200)
41       enddo
42       enddo
43       Close (20)
44       Do J = 1, num
45       Close (J)
46       End Do
47       End Program

请您指导如何修改，谢谢!

fcode · 发表于 2015-3-30 19:15:32

[Fortran] 纯文本查看 复制代码

01Program main
02  Implicit None
03  Character (160) :: filename
04  Character (160) :: fn1, fn2, fn3, tmp
05  Integer :: i, j, num
06  Real :: a(4), b(4), c(4)
07  Integer ios
08 
09  fn1 = 'c6c8_ID_'
10  fn2 = 'c6_'
11  fn3 = 'c8_'
12 
13  num = 160 ! total number of result files
14 
15  Print *, 'Running'
16  Do j = 1, num
17    Write (tmp, *) j
18    !c print*, j
19    Open (j, File='D:\c6c8\'//trim(adjustl(fn1))//trim(adjustl(tmp))//'.dat', Status='unknown') ! save result here
20    Write (j, '(3f8.3)') 0.0, 0.0, 1.0
21  End Do
22 
23  Open (20, File='IDlist.dat', Status='old') ! 里面存放的是10的倍数，共计9200个。
24  Do i = 1, 9200
25    Read (20, *, Iostat=ios) filename
26    If (ios/=0) Exit
27    ! print*, filename
28    Open (200, File='D:\C6\'//trim(adjustl(fn2))//trim(adjustl(filename))//'.dat')
29    Open (300, File='D:\C8\'//trim(adjustl(fn3))//trim(adjustl(filename))//'.dat')
30 
31    Do j = 1, num
32      Read (200, *) a(:) ! first row=ID number， then x,y,z
33      Read (300, *) b(:) ! first row=ID number
34      c(:) = b(:) - a(:)
35      !c print*, c(:)
36      Write ( j , '(3f8.3)') c(2:4) ! write( j ,
37    End Do
38    Close (300)!// Close 在循环外面
39    Close (200)
40  End Do
41  Close (20)
42  Do j = 1, num
43    Close (j)
44  End Do
45End Program main

kif117 · 发表于 2015-3-30 19:38:55

本帖最后由 kif117 于 2015-3-30 19:41 编辑

fcode 发表于 2015-3-30 19:15
[mw_shl_code=fortran,true]Program main
Implicit None
Character (160) :: filename

您好! 我这样改过以后，目前的result是160个ID文件里只有(0，0，1)和一行不知道具体从哪里提取的(X,Y,Z)。

又新发现: 原来存为9200个挑选的文件号的IDList 被重写成这样:
10
0.415 -4.906 -3.266

fcode · 发表于 2015-3-30 19:56:28

你把 Open(20 , Read(20 , Close(20 都改为 201

kif117 · 发表于 2015-3-30 20:39:08

本帖最后由 kif117 于 2015-3-31 02:49 编辑

fcode 发表于 2015-3-30 19:56
你把 Open(20 , Read(20 , Close(20 都改为 201

应该是文件号的问题，这种compile时不出现的错误实在很难察觉。我接下来把问题做得‘更难’了一些，是这样的:

现在我想这样求: 比照之前统计各个ID的情况，现在想做这样的运算:

针对一种类数据(c6)的ID1，计算它和另一种类数据(c8)其他159个ID的c(:)=B(:)-A(:)，再求平均。(比如说C6_1中ID1的数值，和C8_1的ID(2至160)共159个其他ID的数值做差后求平均，作为result_1的第一行数据，...以此类推，一直做完全部9200个文件。)

是不是只能通过进行比较 if (id...) 之类的判断语句来实现?

vvt · 发表于 2015-4-2 08:13:35

sum( c8_1_ID( 2:160 ) - c6_1_ID(1) ) / 159
这样就行了

[求助] 一个批量运算的相关问题

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章