Fortran Coder

标题: 数组指针的堆栈问题 [打印本页]

作者: li913    时间: 2021-11-13 12:33
标题: 数组指针的堆栈问题
将已关联的数组指针赋值给另一个数组,oneAPI会对源数组进行一次copy,并默认存放于堆栈中。如果源数组过大,造成性能损失,甚至可能导致栈溢出。
1、如何将临时数据移至内存而不是堆栈中?
2、指针赋值给数组,可否有其他高效的赋值函数或手段?
[Fortran] 纯文本查看 复制代码
  program Test
  implicit none
  integer(4),allocatable,target:: a(:,:,:), b(:,:,:)
  integer(4), pointer::p(:,:,:)
  allocate(a(1024,1024,10),b(1024,1024,10)) !大小各40M
  a=0
  b=1
  p=>b !任务管理器显示使用91M
  a=p  !内存使用 131M。如果数组过大,会栈溢出
  write(*,*) a(1,1,1)
  pause
  end program Test


作者: xxdw    时间: 2021-11-14 09:49
插个TP眼。
作者: fcode    时间: 2021-11-17 08:24
没有好的思路,只能关注一下。
第一个问题,最好去intel论坛问问。
第二个问题,个人愚见,fortran的数组在内存里不能保证是连续的(可以由间距,甚至多维数组的间距也可以不一样),因此数组的赋值应该没什么提升空间。
作者: li913    时间: 2021-11-17 10:33
执行第九行a=p ( p指向b),实际会执行两次赋值: temp = b, a = temp 。这个效率低于 a=b.




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2