My model crashes at this single line of code, with an erroneous arithmetic operation statement (preceded by the statement Floating-point exception) : gcc1 = exp(-84.165_DP / t9a13) I have trapped the value of t9a13 at 0.271320760 Now, this does produce an arithmetic result of 1.90445382861070831…

This code works on my side with GFortran, with no compiler option. It also gives a correct result with -184.165_DP, and with -284.165_DP it finally prints 0.0000000000000000

Thank you. I am doing another run with the code; I think I may have spotted the error. Basically, the actual variable definition was real :: gcc1, t9a13 and not real(DP) :: gcc1, t9a13; as I had it in my example above and in the code model I was looking at at the time. I realised my error (looking …

Again default kinds malediction…

EXP failure with very small negative value

Help

tyranids June 11, 2023, 6:48pm 5

Yeah take a look at this thread for similar issue/explanation https://fortran-lang.discourse.group/t/can-floating-point-literals-be-adapted-by-the-compiler-to-double-precision-variable/

Topic		Replies	Views
Resolving floating point exceptions signalled after program termination Help	9	2173	April 25, 2023
Flang-new 19.0: "warning: overflow on power with INTEGER exponent" Flang	2	299	April 9, 2024
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO Uncategorized	16	8086	April 4, 2023
Simple Mathematics Help	2	527	February 17, 2022
Exponent(inf) trouble Help	21	518	April 19, 2026

EXP failure with very small negative value

Related topics