Hi friends,
I am using the OneAPI 2025.0 and 2015.1 HPC bundles 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 (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 dgeev.
I am using Xubuntu 22.04 operating system.
I would be grateful if you can look into it. Thank you very much.
Dr. Lin-Hui Ye
Peking University.
