topic Floating divide by zero error in dfeast_scsrgv in IntelĀ® Fortran Compiler
https://community.intel.com/t5/Intel-Fortran-Compiler/Floating-divide-by-zero-error-in-dfeast-scsrgv/m-p/1219813#M152314
<P>Hi,</P>
<P>I am using the FEAST eigenvalue solver to solve a generalized eigenvalue problem with symmetric sparse matrices in Fortran (Compiler: ifort 2020). I store the upper triangular parts of the matrices, i.e. set UPLO = 'U'.</P>
<P>Currently, I receive an error from ifort "forrtl: error (73): floating divide by zero", which is traced back to the dfeast_scsrgv call. I know that the matrices in the generalized eigenvalue problem may contain small values (i.e. values close to zero). So the sparse patterns might contain zero values. Could this be the cause of the error?</P>
<P>I have tried setting fpm(27) = 1 and fpm(28) = 1 to check input matrices, but this is not recognized by the dfeast_scsrgv call (it isn't listed in the output: List of input parameters fpm(1:64)-- if different from default). Is this a known issue?</P>
<P>Thank you,</P>
<P>Anna</P>Wed, 21 Oct 2020 14:42:51 GMTDalklint__Anna2020-10-21T14:42:51ZFloating divide by zero error in dfeast_scsrgv
https://community.intel.com/t5/Intel-Fortran-Compiler/Floating-divide-by-zero-error-in-dfeast-scsrgv/m-p/1219813#M152314
<P>Hi,</P>
<P>I am using the FEAST eigenvalue solver to solve a generalized eigenvalue problem with symmetric sparse matrices in Fortran (Compiler: ifort 2020). I store the upper triangular parts of the matrices, i.e. set UPLO = 'U'.</P>
<P>Currently, I receive an error from ifort "forrtl: error (73): floating divide by zero", which is traced back to the dfeast_scsrgv call. I know that the matrices in the generalized eigenvalue problem may contain small values (i.e. values close to zero). So the sparse patterns might contain zero values. Could this be the cause of the error?</P>
<P>I have tried setting fpm(27) = 1 and fpm(28) = 1 to check input matrices, but this is not recognized by the dfeast_scsrgv call (it isn't listed in the output: List of input parameters fpm(1:64)-- if different from default). Is this a known issue?</P>
<P>Thank you,</P>
<P>Anna</P>Wed, 21 Oct 2020 14:42:51 GMThttps://community.intel.com/t5/Intel-Fortran-Compiler/Floating-divide-by-zero-error-in-dfeast-scsrgv/m-p/1219813#M152314Dalklint__Anna2020-10-21T14:42:51ZRe: Floating divide by zero error in dfeast_scsrgv
https://community.intel.com/t5/Intel-Fortran-Compiler/Floating-divide-by-zero-error-in-dfeast-scsrgv/m-p/1219906#M152321
<P>You may want to try the option -no-ftz</P>
<P><A href="https://software.intel.com/content/www/us/en/develop/documentation/fortran-compiler-developer-guide-and-reference/top/compiler-reference/compiler-options/compiler-option-details/floating-point-options/ftz-qftz.html" target="_blank">https://software.intel.com/content/www/us/en/develop/documentation/fortran-compiler-developer-guide-and-reference/top/compiler-reference/compiler-options/compiler-option-details/floating-point-options/ftz-qftz.html</A></P>
<P>this may or may not help. If you think you are in the realm of denormals FTZ could be causing the zeroes.</P>Wed, 21 Oct 2020 19:50:10 GMThttps://community.intel.com/t5/Intel-Fortran-Compiler/Floating-divide-by-zero-error-in-dfeast-scsrgv/m-p/1219906#M152321Ron_Green2020-10-21T19:50:10Z