Floating divide by zero error in dgeev and possibly also others

yelh2015 — Thu, 08 May 2025 02:57:39 GMT

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.

Re: Floating divide by zero error in dgeev and possibly also others

yelh2015 — Thu, 08 May 2025 06:46:40 GMT

By the way, the 'floating divide by zero' error of dgeev and zhegvx happened under xubuntu 22.04 operating system.

Dr. Lin-Hui Ye

Peking University

Re: Floating divide by zero error in dgeev and possibly also others

Steve_Lionel — Thu, 08 May 2025 15:14:28 GMT

Please ask this in Intel® oneAPI Math Kernel Library - Intel Community as it is not a Fortran compiler issue.

topic Floating divide by zero error in dgeev and possibly also others in Intel® oneAPI Math Kernel Library

Floating divide by zero error in dgeev and possibly also others

Re: Floating divide by zero error in dgeev and possibly also others

Re: Floating divide by zero error in dgeev and possibly also others