循环问题求教

kif117

计算要求是这样的:
有两种文件:
第一种文件数据为 a1,a2,a3....a100
第二种文件数据为 b1,b2,b3...b100
; 其中每种文件各有1000个，文件名分别是x1,x2,....x100(x对应第一种文件)和y1,y2,...y100(y对应第二种文件)，分别放入名为x和y的两个文件夹内。

计算过程是:
以各自文件夹的第一个文件为例子，对应文件x1和y1，
b1-a(1，100)
b2-a(1，100)
...
b100-a(1，100)
分别用b1减去全部数据a, b2减去全部数据a,.....直到第a100; 这样得到共10000数据个差。
对它进行运算，对最后的运算结果求和，为输出的第一个数......循环操作全部1000个文件，最后得到共1000个输出结果。

这是经过帮助优化后的计算程序:

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

查看代码

打印?

01	PROGRAM test

02	IMPLICIT NONE

03	REAL, ALLOCATABLE :: a(:),b(:),hhs(:)

04	INTEGER :: i, z, j, k, tn, num, numm

05	REAL    :: summ, cos2, r1, q1, w1

06	character(160):: filename

07	Real(8) time

08	Integer(4) time0, time1, dtime

09

10

num=100

11	numm=num*num

12

tn=1000

13	ALLOCATE(a(num),b(num),hhs(numm))

14	Print*, 'Running...'

15	Call system_clock(time0)

16	OPEN(10, FILE='output.dat',status='unknown' )

17

18

Do i=1,tn

19	write(filename,'(i6)') i

20	Open(200,File="D:\x\"//trim(adjustl(filename))//'.dat',status='old')

21	READ(200,*) a

22	close(200)

23	Open(300,File="D:\y\"//trim(adjustl(filename))//'.dat',status='old')

24	READ(300,*) b

25	close(300)

26

z=0

27

summ=0.0

28

r1 = 20

29	do  j = 1, num

30	do k = 1, num

31

z=z+1

32	q1 = b(j)-a(k)

33	cos2 = (r1*q1+r1*r1+q1*q1)**2/3

34	w1=sqrt(q1*q1+r1*r1)

35	if (w1.GT.0.and.w1.LE.20) then

36	hhs(z)=cos2/(w1*w1*w1)

37

endif

38

enddo

39

enddo

40	summ=sum(hhs)

41	WRITE(10,'(1X,F12.6)') summ

42

enddo

43	DEALLOCATE(a,b)

44

close(10)

45	Call system_clock(time1, dtime)

46	time = 1D0*(time1-time0)/dtime

47	Write (*, '(a7,f16.7)') 'Time = ', time

48	END PROGRAM

现在我想改变一下运算，就是将其中的 r1 = 20由定量换成变量，让它等于以各自文件夹的第一个文件操作所得到的数据差，
对应文件x1和y1，
b1-a(1，100)
b2-a(1，100)
...
b100-a(1，100)
分别用b1减去全部数据a, b2减去全部数据a,.....直到第a100; 这样得到共10000数据。以最初的10000个数据作为r1，然后循环其他文件和它进行计算。
我想问能不能不再另外写循环实现这个目的? 我现在的写法是另外独立写一个循环读取，这样又降低了计算速度。
感谢指导!

fcode

这逻辑有点复杂，还是没理清。

r1 等于变量没有问题，你可以用任意变量对其赋值。其他逻辑嘛，有点乱，帮不上你。

kif117

fcode 发表于 2015-6-2 10:17
这逻辑有点复杂，还是没理清。

r1 等于变量没有问题，你可以用任意变量对其赋值。其他逻辑嘛，有点乱，帮 ...

之前我说得不好，现在我努力分成部分地说一说问题:

首先，假设有两种文件，分别放在两个文件夹里，分成是A和B。对应文件是x1,x2,...x1000，以及y1, y2,...y1000，假设每种各1000个文件。

其中，每个单独的文件，无论种类，分别包含10000个数据，数据共有三列(假设)。x对应a数据，y对应b数据。以第一列数据为例，取x1和y1:

b1-a(1,10000)
b2-a(1,10000)
...
b10000-a(1,10000)

对应一组(a1,b1)一共有10000*10000个数据，这是已经写好的程序循环，之后操作全部的x和y，一共有1000组10000*10000个数据。

我想用2-1000组10000*10000个数据和第1组10000*10000个数据进行减法计算，为了提高运算速度，想请教能不能不再另外写循环就实现目的。现在我的办法是把第1组10000*10000个数据另外写到一个文件里，然后去循环和read.


感谢您的帮助!

fcode

开一个 10000*10000 的数组，把它存储起来。就不需要写到文件里，循环read了。
但是2-1000组数据与第一组的减法，还是要循环
do i = 2 , 1000
的