src/gray_tables.f90: avoid floating point exception

Due to Fortran non-shortcircuiting expressions, if info /= 0 and psi_n = NaN or is undefined, it will still raise an exception.
2024-11-04 10:53:45 +01:00 · 2024-11-04 10:53:45 +01:00 · 2368a0ba33
commit 2368a0ba33
parent 72eb224568
1 changed files with 2 additions and 1 deletions
--- a/src/gray_tables.f90
+++ b/src/gray_tables.f90
@ -410,7 +410,8 @@ contains
        psi_n = sol(1)

        ! Handle spurious or no solutions
-        if (info /= 1 .or. psi_n < 0 .or. psi_n > 1) cycle
+        if (info /= 1) cycle
+        if (psi_n < 0 .or. psi_n > 1) cycle
      end block

      ! Compute and print the ψ_n(R,z) = ψ_n₀ contour