Fortran Coder

标题: eoshift函数gfortran运行正常，intel Fortran运行却出错 [打印本页]

---

作者: 愤怒的三炮    时间: 2021-2-5 21:23
标题: eoshift函数gfortran运行正常，intel Fortran运行却出错
看到Fortran内置的一个数组操作函数“eoshift”，功能是将数组进行挪位。
就想测试一下它和普通方法的速度差异，代码如下：

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

program eoshift_test

    implicit none

    integer:: i

    integer,dimension(:,:,:),allocatable:: a 

    integer,dimension(:,:,:),allocatable:: temp !临时数组

    integer:: grid=1000,dim1=30 ! 数组大小

    real:: time_start, time_end ! 计时



    allocate(a(dim1,grid,grid))

    allocate(temp(dim1,grid,grid))



    a=0

    temp=0

    !===== 用 eoshift函数 ==================

    call cpu_time(time_start)

        a=eoshift(a,shift=1,boundary=0,dim=1)

    call cpu_time(time_end)

    write(*,*) 'function eoshift:'

    write(*,*) time_end-time_start,'seconds elapsed'



    !====== 用常规方法 =======================

    call cpu_time(time_start)

        do i =1, dim1-1

            temp(i,:,:)=a(i+1,:,:)

        end do

        a=temp

        a(dim1,:,:)=0

    call cpu_time(time_end)

    write(*,*) 'normal method:'

    write(*,*) time_end-time_start,'seconds elapsed'

end

用gfortran编译运行没有问题；
但是用Intel Fortran编译通过，运行时却报错“StackOverflow”；
在Linux下也是同样报错；

既然是栈溢出的问题，于是我将数组大小调整得很小，果然不报错了，但问题是：
1. 调用内置函数eoshift为什么gfortran没有StackOverflow，而ifort却StackOverflow了？
2. 我只知道a是我分配到内存里的，而内置函数eoshift为什么会使用栈？
希望有前辈能够指导一下，感激不尽！

{31522AA7-DFF2-402A-A5FF-7E202B389BF3}.png.jpg (99.48 KB, 下载次数: 370)

下载附件

2021-2-5 21:22 上传

Linux下也报错

{63225753-1D20-44ED-B8E2-D1A1AD76F591}.png.jpg (78.36 KB, 下载次数: 344)

下载附件

2021-2-5 21:22 上传

visual studio StackOverflow

---

作者: fcode    时间: 2021-2-10 12:26
建议不要用intrinsic的函数处理数据量较大的数组。容易堆栈溢出，不同编译器处理方式不同。

---

作者: 愤怒的三炮    时间: 2021-2-14 15:56
好的，谢谢雪球，祝新春愉快！

---

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)

Powered by Discuz! X3.2