Circumvent persistent FPE for vectorized sqrt

jppelteret · Jun 12, 2022 · 351d069 · 351d069
1 parent 391dd74
commit 351d069
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/include/weak_forms/unary_operators.h b/include/weak_forms/unary_operators.h
@@ -1607,6 +1607,10 @@ private:                                                                 \
       {
         using namespace std;
 
+        // Circumvent persistent FPE
+        if (value == vectorized_value_type<ScalarType, width>{})
+          return vectorized_value_type<ScalarType, width>{};
+
 #ifdef WEAK_FORMS_VECTORIZATION_FPE_SQRT_OF_ZERO
         // We need to check if the vector is a zero vector. If so, simply
         // avoid the operation that would result in an FPE as we know the