adding QUAD_NAN constant

Author: SwayamInSync

quaddtype/numpy_quaddtype/src/ops.hpp

@@ -6,6 +6,7 @@
 #define QUAD_ZERO sleef_q(+0x0000000000000LL, 0x0000000000000000ULL, -16383)
 #define QUAD_ONE sleef_q(+0x1000000000000LL, 0x0000000000000000ULL, 0)
 #define QUAD_POS_INF sleef_q(+0x1000000000000LL, 0x0000000000000000ULL, 16384)
+#define QUAD_NAN sleef_q(+0x180000000000LL, 0x0000000000000000ULL, 16384)
 
 // Unary Quad Operations
 typedef Sleef_quad (*unary_op_quad_def)(const Sleef_quad *);
@@ -396,7 +397,7 @@ quad_mod(const Sleef_quad *a, const Sleef_quad *b)
 {
     // division by zero
     if (Sleef_icmpeqq1(*b, QUAD_ZERO)) {
-        return Sleef_cast_from_doubleq1(NAN);
+        return QUAD_NAN;
     }
 
     // NaN inputs
@@ -406,7 +407,7 @@ quad_mod(const Sleef_quad *a, const Sleef_quad *b)
 
     // infinity dividend -> NaN
     if (quad_isinf(a)) {
-        return Sleef_cast_from_doubleq1(NAN);
+        return QUAD_NAN;
     }
 
     // finite % inf