Workaround rocm 5.7 compiler bug in debug mode when evaluating an if condition

diff --git a/source/base/polynomials_piecewise.cc b/source/base/polynomials_piecewise.cc
index 0514f334b22586bc1d728f1e95410fc1c8a333e4..67e96a0fabf16e7c9c936496f9c0d6a3da00b904 100644 (file)
--- a/source/base/polynomials_piecewise.cc
+++ b/source/base/polynomials_piecewise.cc
@@ -139,7 +139,16 @@ namespace Polynomials
         else
           {
             const double offset = step * interval;
-            if (x < offset || x > offset + step)
+            // ROCm 5.7 throw a floating point exception in debug when trying to
+            // evaluate (x < offset && x > offset + step). Separating the
+            // conditions fixes the issue.
+            if (x < offset)
+              {
+                for (unsigned int k = 0; k <= n_derivatives; ++k)
+                  values[k] = 0;
+                return;
+              }
+            else if (x > offset + step)
               {
                 for (unsigned int k = 0; k <= n_derivatives; ++k)
                   values[k] = 0;