Fortran Coder

查看: 2002|回复: 3

[数值问题] 大数组的生成和reshape出错

[复制链接]

2

帖子

1

主题

0

精华

新人

F 币
17 元
贡献
7 点
发表于 2022-5-25 23:17:29 | 显示全部楼层 |阅读模式
本帖最后由 kyra 于 2022-5-26 08:44 编辑

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

  implicit none
  integer, parameter:: N=2000
  integer :: i,j
  real(16),  dimension(N,N):: a=0._16, b=0._16,  c0=0._16, c1=0._16
!  real(16), dimension(N*N):: aa=0, bb=0, cc=0
  character (len=20):: date,time 
  real :: tt0,tt1

   call date_and_time(date,time)
   print *, "date=", date, ";  time=", time 
   call cpu_time(tt0)
   a=reshape([(50._16-(N**2-1._16)/2.+i,i=1,N**2)],[N,N])
   b=reshape([(100._16-(N**2-1._16)/2.+i,i=1,N**2)],[N,N])
!   aa=reshape(a,[N*N])
!   bb=reshape(b,[N*N])
   call cpu_time(tt1)
   print *, "date=", date, ";  time=", time 
   print *, "reshape spend=", tt1-tt0
   call date_and_time(date,time)
   print *, "date=", date, ";  time=", time 
   call cpu_time(tt0)
   c0=matmul(a,b)
!   c0=matmul(c0,b)
!   c0=matmul(c0,b)
!   c0=matmul(c0,b)
!   c0=matmul(c0,b)
!   c0=matmul(c0,b)
!   c0=matmul(c0,b)
!   c0=matmul(c0,b)
   call date_and_time(date,time)
   call cpu_time(tt1)
   print *, "date=", date, ";  time=", time 
   print *, "multi spend=", tt1-tt0
!  print *,"complete first"

end program


gfortran 编译没问题,运行不出错。但是Flang 编译后运行会出segment fault, ifort 好像也会出。到底是啥问题呢?

2

帖子

1

主题

0

精华

新人

F 币
17 元
贡献
7 点
 楼主| 发表于 2022-5-25 23:21:30 | 显示全部楼层
貌似N=2000太大了, gfortran 编译时候加了选项-fmax-array-constructor=5000000 , flang 和ifort 要加啥选项呢

790

帖子

2

主题

0

精华

大宗师

F 币
3765 元
贡献
2255 点
发表于 2022-5-26 08:29:39 | 显示全部楼层
改为动态数组,根本上解决问题。

235

帖子

0

主题

0

精华

版主

World Analyser

F 币
630 元
贡献
464 点

新人勋章美女勋章元老勋章热心勋章规矩勋章管理勋章

QQ
发表于 2022-5-26 08:43:10 | 显示全部楼层
   a=reshape([(50._16-(N**2-1._16)/2.+i,i=1,N**2)],[N,N])
   b=reshape([(100._16-(N**2-1._16)/2.+i,i=1,N**2)],[N,N])

这两个语句非常吃堆栈,建议还是老老实实的Do循环
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-3-29 21:55

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表