大数组的生成和reshape出错

XiaOaiX · 发表于 2022-5-25 23:17:29

本帖最后由 kyra 于 2022-5-26 08:44 编辑

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

01program arraymult
02 
03  implicit none
04  integer, parameter:: N=2000
05  integer :: i,j
06  real(16),  dimension(N,N):: a=0._16, b=0._16,  c0=0._16, c1=0._16
07!  real(16), dimension(N*N):: aa=0, bb=0, cc=0
08  character (len=20):: date,time
09  real :: tt0,tt1
10 
11   call date_and_time(date,time)
12   print *, "date=", date, ";  time=", time
13   call cpu_time(tt0)
14   a=reshape([(50._16-(N**2-1._16)/2.+i,i=1,N**2)],[N,N])
15   b=reshape([(100._16-(N**2-1._16)/2.+i,i=1,N**2)],[N,N])
16!   aa=reshape(a,[N*N])
17!   bb=reshape(b,[N*N])
18   call cpu_time(tt1)
19   print *, "date=", date, ";  time=", time
20   print *, "reshape spend=", tt1-tt0
21   call date_and_time(date,time)
22   print *, "date=", date, ";  time=", time
23   call cpu_time(tt0)
24   c0=matmul(a,b)
25!   c0=matmul(c0,b)
26!   c0=matmul(c0,b)
27!   c0=matmul(c0,b)
28!   c0=matmul(c0,b)
29!   c0=matmul(c0,b)
30!   c0=matmul(c0,b)
31!   c0=matmul(c0,b)
32   call date_and_time(date,time)
33   call cpu_time(tt1)
34   print *, "date=", date, ";  time=", time
35   print *, "multi spend=", tt1-tt0
36!  print *,"complete first"
37 
38end program