Hi,
I am using OneAPI 2025.0 and 2015.1 to compile the following Fortran code:
solv.f:
program main
implicit none
integer(kind=4) :: info,lwork,i
real(kind=8), dimension(12) :: work
real(kind=8), dimension(3) :: er,ei
real(kind=8), dimension(3,3) :: mat,vl,vvr
lwork=12
mat=0.d0
do i=1,3
mat(i,i)=1.d0
enddo
call dgeev('N','V',3,mat,3,er,ei,vl,3,vvr,3,work,12,info)
end program main
My Makefile is:
solv:
ifx -O0 -g -fpe0 -gdwarf-2 -c solv.f
ifx -O0 -g -fpe0 -gdwarf-2 -o solv solv.o -lmkl_lapack95_lp64 -lmkl_intel_lp64 -lmkl_core -lmkl_sequential -lpthread
clean:
rm -f solv solv.o
Running the program produces 'floating divide by zero' error. The program runs well if compiled by OneAPI 2024.2, or gfortran, or other Intel compiler versions dated 2020,2018, etc. Also, the error does not show up if the code is compiled by other options like "-O2".
Similar error is also found for zhegvx. But I do not have a testing program to reproduce it. The dependencies on the versions of the Intel compilers & MKL libraries, as well as the compiler options, are all the same as zgeev.
I would be grateful if you can look into it. Thanks in advance.
Dr. Lin-Hui Ye
Peking University.
