ifort bug example

run_me.f90

! ifort run_me.f90 -o run_me
! gfortran run_me.f90 -o run_me
program run_me

    implicit none
    
    double precision, allocatable, dimension(:,:) :: M
    integer :: icol, irow
    

    ! This will work with ifort and gfortran
    irow = 512
    icol = 512
    
    allocate(M(irow,icol))
    
    M(:,:) = 1.0d0
    
    write (*,*) "M", shape(M)
    write (*,*) shape(matmul(transpose(M), M))
    
    deallocate(M)

    ! This will fail with ifort but not gfortran
    irow = 768
    icol = 768
    
    allocate(M(irow,icol))

    M(:,:) = 1.0d0
    
    write (*,*) "M", shape(M)
    write (*,*) shape(matmul(transpose(M), M))

    deallocate(M)
    
    
end program run_me