From 684a36f4960b1587f687e79e58fa3b12dc610b42 Mon Sep 17 00:00:00 2001
From: Claude <noreply@anthropic.com>
Date: Mon, 20 Jul 2026 13:41:48 +0000
Subject: [PATCH 1/4] refactor: use `constants/float64/nan` across
 `math/base/special`

Propagates fix from 05aaec1 ("refactor: use `constants/float64/nan`
in `math/base/special/pow`") to 100 sibling packages generating NaN
via the inline `0.0/0.0` idiom in their C implementations. Doc
example occurrences are intentionally left unchanged.

Ref: https://github.com/stdlib-js/stdlib/issues/11743

Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01EHa9r1LqZ3Tg3Gft34EGQv
---
 .../math/base/special/acos/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/acos/src/main.c |  5 ++--
 .../math/base/special/acosh/manifest.json     |  9 ++++---
 .../math/base/special/acosh/src/main.c        |  3 ++-
 .../math/base/special/asin/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/asin/src/main.c |  5 ++--
 .../math/base/special/atan2/manifest.json     |  9 ++++---
 .../math/base/special/atan2/src/main.c        |  3 ++-
 .../math/base/special/atanh/manifest.json     |  9 ++++---
 .../math/base/special/atanh/src/main.c        |  3 ++-
 .../math/base/special/bernoulli/manifest.json |  9 ++++---
 .../math/base/special/bernoulli/src/main.c    |  3 ++-
 .../math/base/special/bessely0/manifest.json  |  9 ++++---
 .../math/base/special/bessely0/src/main.c     |  3 ++-
 .../math/base/special/bessely1/manifest.json  |  9 ++++---
 .../math/base/special/bessely1/src/main.c     |  3 ++-
 .../math/base/special/beta/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/beta/src/main.c |  5 ++--
 .../math/base/special/betaln/manifest.json    |  9 ++++---
 .../math/base/special/betaln/src/main.c       |  7 +++---
 .../math/base/special/binet/manifest.json     |  9 ++++---
 .../math/base/special/binet/src/main.c        |  3 ++-
 .../math/base/special/binomcoef/manifest.json |  9 ++++---
 .../math/base/special/binomcoef/src/main.c    |  3 ++-
 .../base/special/binomcoefln/manifest.json    |  9 ++++---
 .../math/base/special/binomcoefln/src/main.c  |  3 ++-
 .../math/base/special/boxcox/manifest.json    |  9 ++++---
 .../math/base/special/boxcox/src/main.c       |  3 ++-
 .../math/base/special/boxcox1p/manifest.json  |  9 ++++---
 .../math/base/special/boxcox1p/src/main.c     |  3 ++-
 .../base/special/boxcox1pinv/manifest.json    |  9 ++++---
 .../math/base/special/boxcox1pinv/src/main.c  |  3 ++-
 .../math/base/special/boxcoxinv/manifest.json |  9 ++++---
 .../math/base/special/boxcoxinv/src/main.c    |  3 ++-
 .../math/base/special/ceilb/manifest.json     |  9 ++++---
 .../math/base/special/ceilb/src/main.c        |  3 ++-
 .../math/base/special/ceilsd/manifest.json    |  9 ++++---
 .../math/base/special/ceilsd/src/main.c       |  3 ++-
 .../math/base/special/cexp/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/cexp/src/main.c |  9 ++++---
 .../math/base/special/clamp/manifest.json     |  9 ++++---
 .../math/base/special/clamp/src/main.c        |  3 ++-
 .../math/base/special/cos/manifest.json       |  9 ++++---
 .../@stdlib/math/base/special/cos/src/main.c  |  3 ++-
 .../math/base/special/cosd/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/cosd/src/main.c |  3 ++-
 .../math/base/special/cospi/manifest.json     |  9 ++++---
 .../math/base/special/cospi/src/main.c        |  3 ++-
 .../math/base/special/digamma/manifest.json   |  9 ++++---
 .../math/base/special/digamma/src/main.c      |  5 ++--
 .../base/special/dirac-delta/manifest.json    |  9 ++++---
 .../math/base/special/dirac-delta/src/main.c  |  3 ++-
 .../base/special/dirichlet-eta/manifest.json  |  9 ++++---
 .../base/special/dirichlet-eta/src/main.c     |  3 ++-
 .../math/base/special/ellipe/manifest.json    |  9 ++++---
 .../math/base/special/ellipe/src/main.c       |  3 ++-
 .../math/base/special/ellipj/manifest.json    |  9 ++++---
 .../math/base/special/ellipj/src/main.c       | 13 +++++-----
 .../math/base/special/ellipk/manifest.json    |  9 ++++---
 .../math/base/special/ellipk/src/main.c       |  3 ++-
 .../math/base/special/erf/manifest.json       | 12 ++++++---
 .../@stdlib/math/base/special/erf/src/main.c  |  3 ++-
 .../math/base/special/erfc/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/erfc/src/main.c |  3 ++-
 .../math/base/special/erfcinv/manifest.json   |  9 ++++---
 .../math/base/special/erfcinv/src/main.c      |  5 ++--
 .../math/base/special/erfinv/manifest.json    |  9 ++++---
 .../math/base/special/erfinv/src/main.c       |  5 ++--
 .../math/base/special/exp10/manifest.json     |  9 ++++---
 .../math/base/special/exp10/src/main.c        |  3 ++-
 .../math/base/special/expit/manifest.json     |  9 ++++---
 .../math/base/special/expit/src/main.c        |  3 ++-
 .../math/base/special/factorial/manifest.json |  9 ++++---
 .../math/base/special/factorial/src/main.c    |  5 ++--
 .../base/special/factorial2/manifest.json     |  9 ++++---
 .../math/base/special/factorial2/src/main.c   |  3 ++-
 .../base/special/factorialln/manifest.json    |  9 ++++---
 .../math/base/special/factorialln/src/main.c  |  3 ++-
 .../special/falling-factorial/manifest.json   |  9 ++++---
 .../base/special/falling-factorial/src/main.c |  5 ++--
 .../base/special/fast/acosh/manifest.json     |  9 ++++---
 .../math/base/special/fast/acosh/src/main.c   |  3 ++-
 .../base/special/fast/atanh/manifest.json     |  9 ++++---
 .../math/base/special/fast/atanh/src/main.c   |  3 ++-
 .../base/special/fast/pow-int/manifest.json   |  9 ++++---
 .../math/base/special/fast/pow-int/src/main.c |  3 ++-
 .../special/fibonacci-index/manifest.json     |  9 ++++---
 .../base/special/fibonacci-index/src/main.c   |  3 ++-
 .../math/base/special/floorb/manifest.json    |  9 ++++---
 .../math/base/special/floorb/src/main.c       |  3 ++-
 .../math/base/special/floorsd/manifest.json   |  9 ++++---
 .../math/base/special/floorsd/src/main.c      |  3 ++-
 .../math/base/special/gamma/manifest.json     |  9 ++++---
 .../math/base/special/gamma/src/main.c        |  3 ++-
 .../math/base/special/gamma1pm1/manifest.json |  9 ++++---
 .../math/base/special/gamma1pm1/src/main.c    |  3 ++-
 .../math/base/special/gammainc/manifest.json  |  9 ++++---
 .../math/base/special/gammainc/src/main.c     |  3 ++-
 .../base/special/gammaincinv/manifest.json    |  9 ++++---
 .../math/base/special/gammaincinv/src/main.c  |  9 ++++---
 .../math/base/special/gcd/manifest.json       |  9 ++++---
 .../@stdlib/math/base/special/gcd/src/main.c  |  7 +++---
 .../math/base/special/heaviside/manifest.json |  9 ++++---
 .../math/base/special/heaviside/src/main.c    |  7 +++---
 .../math/base/special/hypot/manifest.json     |  9 ++++---
 .../math/base/special/hypot/src/main.c        |  3 ++-
 .../base/special/kernel-betainc/manifest.json |  9 ++++---
 .../base/special/kernel-betainc/src/main.c    | 25 ++++++++++---------
 .../special/kronecker-delta/manifest.json     |  9 ++++---
 .../base/special/kronecker-delta/src/main.c   |  3 ++-
 .../math/base/special/ln/manifest.json        | 12 ++++++---
 .../@stdlib/math/base/special/ln/src/main.c   |  3 ++-
 .../math/base/special/log10/manifest.json     |  9 ++++---
 .../math/base/special/log10/src/main.c        |  3 ++-
 .../math/base/special/log1mexp/manifest.json  |  9 ++++---
 .../math/base/special/log1mexp/src/main.c     |  3 ++-
 .../math/base/special/log1p/manifest.json     |  9 ++++---
 .../math/base/special/log1p/src/main.c        |  3 ++-
 .../math/base/special/log1pexp/manifest.json  |  9 ++++---
 .../math/base/special/log1pexp/src/main.c     |  3 ++-
 .../math/base/special/log1pmx/manifest.json   |  9 ++++---
 .../math/base/special/log1pmx/src/main.c      |  3 ++-
 .../math/base/special/log2/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/log2/src/main.c |  3 ++-
 .../math/base/special/logaddexp/manifest.json |  9 ++++---
 .../math/base/special/logaddexp/src/main.c    |  3 ++-
 .../math/base/special/logit/manifest.json     |  9 ++++---
 .../math/base/special/logit/src/main.c        |  5 ++--
 .../math/base/special/lucas/manifest.json     |  9 ++++---
 .../math/base/special/lucas/src/main.c        |  3 ++-
 .../math/base/special/max/manifest.json       |  9 ++++---
 .../@stdlib/math/base/special/max/src/main.c  |  3 ++-
 .../math/base/special/min/manifest.json       |  9 ++++---
 .../@stdlib/math/base/special/min/src/main.c  |  3 ++-
 .../math/base/special/minmax/manifest.json    |  9 ++++---
 .../math/base/special/minmax/src/main.c       |  5 ++--
 .../math/base/special/minmaxabs/manifest.json |  9 ++++---
 .../math/base/special/minmaxabs/src/main.c    |  5 ++--
 .../math/base/special/modf/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/modf/src/main.c |  5 ++--
 .../math/base/special/nanmax/manifest.json    |  9 ++++---
 .../math/base/special/nanmax/src/main.c       |  3 ++-
 .../math/base/special/nanmin/manifest.json    |  9 ++++---
 .../math/base/special/nanmin/src/main.c       |  3 ++-
 .../base/special/negafibonacci/manifest.json  |  9 ++++---
 .../base/special/negafibonacci/src/main.c     |  5 ++--
 .../math/base/special/negalucas/manifest.json |  9 ++++---
 .../math/base/special/negalucas/src/main.c    |  3 ++-
 .../math/base/special/pdiff/manifest.json     |  9 ++++---
 .../math/base/special/pdiff/src/main.c        |  3 ++-
 .../math/base/special/powm1/manifest.json     | 10 +++++---
 .../math/base/special/powm1/src/main.c        |  7 +++---
 .../math/base/special/ramp/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/ramp/src/main.c |  3 ++-
 .../math/base/special/rempio2/manifest.json   |  9 ++++---
 .../math/base/special/rempio2/src/main.c      |  5 ++--
 .../math/base/special/rempio2f/manifest.json  |  9 ++++---
 .../math/base/special/rempio2f/src/main.c     |  3 ++-
 .../base/special/riemann-zeta/manifest.json   |  9 ++++---
 .../math/base/special/riemann-zeta/src/main.c |  5 ++--
 .../special/rising-factorial/manifest.json    |  9 ++++---
 .../base/special/rising-factorial/src/main.c  |  3 ++-
 .../math/base/special/round/manifest.json     | 12 ++++++---
 .../math/base/special/round/src/main.c        |  3 ++-
 .../math/base/special/roundb/manifest.json    |  9 ++++---
 .../math/base/special/roundb/src/main.c       |  3 ++-
 .../math/base/special/roundn/manifest.json    |  9 ++++---
 .../math/base/special/roundn/src/main.c       |  3 ++-
 .../math/base/special/roundsd/manifest.json   |  9 ++++---
 .../math/base/special/roundsd/src/main.c      |  3 ++-
 .../math/base/special/sin/manifest.json       |  9 ++++---
 .../@stdlib/math/base/special/sin/src/main.c  |  3 ++-
 .../math/base/special/sinc/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/sinc/src/main.c |  3 ++-
 .../math/base/special/sincos/manifest.json    |  9 ++++---
 .../math/base/special/sincos/src/main.c       |  5 ++--
 .../math/base/special/sincospi/manifest.json  |  9 ++++---
 .../math/base/special/sincospi/src/main.c     |  5 ++--
 .../math/base/special/sind/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/sind/src/main.c |  3 ++-
 .../math/base/special/sinpi/manifest.json     |  9 ++++---
 .../math/base/special/sinpi/src/main.c        |  3 ++-
 .../math/base/special/spence/manifest.json    |  9 ++++---
 .../math/base/special/spence/src/main.c       |  3 ++-
 .../math/base/special/sqrt1pm1/manifest.json  |  9 ++++---
 .../math/base/special/sqrt1pm1/src/main.c     |  3 ++-
 .../math/base/special/tan/manifest.json       |  9 ++++---
 .../@stdlib/math/base/special/tan/src/main.c  |  3 ++-
 .../base/special/tribonacci/manifest.json     |  9 ++++---
 .../math/base/special/tribonacci/src/main.c   |  3 ++-
 .../math/base/special/trigamma/manifest.json  |  9 ++++---
 .../math/base/special/trigamma/src/main.c     |  3 ++-
 .../math/base/special/truncb/manifest.json    |  9 ++++---
 .../math/base/special/truncb/src/main.c       |  3 ++-
 .../math/base/special/truncn/manifest.json    |  9 ++++---
 .../math/base/special/truncn/src/main.c       |  3 ++-
 .../math/base/special/truncsd/manifest.json   |  9 ++++---
 .../math/base/special/truncsd/src/main.c      |  3 ++-
 .../math/base/special/wrap/manifest.json      |  9 ++++---
 .../@stdlib/math/base/special/wrap/src/main.c |  3 ++-
 200 files changed, 853 insertions(+), 451 deletions(-)

diff --git a/lib/node_modules/@stdlib/math/base/special/acos/manifest.json b/lib/node_modules/@stdlib/math/base/special/acos/manifest.json
index 3c6a30c2270f..761fd6920af4 100644
--- a/lib/node_modules/@stdlib/math/base/special/acos/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/acos/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/special/asin",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/fourth-pi"
+        "@stdlib/constants/float64/fourth-pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/special/asin",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/fourth-pi"
+        "@stdlib/constants/float64/fourth-pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/special/asin",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/fourth-pi"
+        "@stdlib/constants/float64/fourth-pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/acos/src/main.c b/lib/node_modules/@stdlib/math/base/special/acos/src/main.c
index b22d65d97279..94793e19b1e9 100644
--- a/lib/node_modules/@stdlib/math/base/special/acos/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/acos/src/main.c
@@ -35,6 +35,7 @@
 #include "stdlib/math/base/special/asin.h"
 #include "stdlib/math/base/special/sqrt.h"
 #include "stdlib/constants/float64/fourth_pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double MOREBITS = 6.123233995736765886130e-17; // pi/2 = PIO2 + MOREBITS.
@@ -81,10 +82,10 @@ static const double MOREBITS = 6.123233995736765886130e-17; // pi/2 = PIO2 + MOR
 double stdlib_base_acos( const double x ) {
 	double z;
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x < -1.0 || x > 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x > 0.5 ) {
 		return 2.0 * stdlib_base_asin( stdlib_base_sqrt( 0.5 - ( 0.5 * x ) ) );
diff --git a/lib/node_modules/@stdlib/math/base/special/acosh/manifest.json b/lib/node_modules/@stdlib/math/base/special/acosh/manifest.json
index c8420422370f..c046f462b2c5 100644
--- a/lib/node_modules/@stdlib/math/base/special/acosh/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/acosh/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/acosh/src/main.c b/lib/node_modules/@stdlib/math/base/special/acosh/src/main.c
index e998e6aeb888..d53894439ba4 100644
--- a/lib/node_modules/@stdlib/math/base/special/acosh/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/acosh/src/main.c
@@ -36,6 +36,7 @@
 #include "stdlib/math/base/special/sqrt.h"
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/ln_two.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double HUGE = 1 << 28; // 2**28
@@ -79,7 +80,7 @@ static const double HUGE = 1 << 28; // 2**28
 double stdlib_base_acosh( const double x ) {
 	double t;
 	if ( stdlib_base_is_nan( x ) || x < 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 1.0 ) {
 		return 0.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/asin/manifest.json b/lib/node_modules/@stdlib/math/base/special/asin/manifest.json
index c470fe9c8dbd..e0a2eb81c3ea 100644
--- a/lib/node_modules/@stdlib/math/base/special/asin/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/asin/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/fourth-pi"
+        "@stdlib/constants/float64/fourth-pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/fourth-pi"
+        "@stdlib/constants/float64/fourth-pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/fourth-pi"
+        "@stdlib/constants/float64/fourth-pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/asin/src/main.c b/lib/node_modules/@stdlib/math/base/special/asin/src/main.c
index 836f08fa3587..e04fef739967 100644
--- a/lib/node_modules/@stdlib/math/base/special/asin/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/asin/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/sqrt.h"
 #include "stdlib/constants/float64/fourth_pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double MOREBITS = 6.123233995736765886130e-17; // pi/2 = PIO2 + MOREBITS
@@ -165,7 +166,7 @@ double stdlib_base_asin( const double x ) {
 	double z;
 
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x > 0.0 ) {
 		sgn = 0;
@@ -175,7 +176,7 @@ double stdlib_base_asin( const double x ) {
 		a = -x;
 	}
 	if ( a > 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( a > 0.625 ) {
 		// arcsin(1-x) = pi/2 - sqrt(2x)(1+R(x))
diff --git a/lib/node_modules/@stdlib/math/base/special/atan2/manifest.json b/lib/node_modules/@stdlib/math/base/special/atan2/manifest.json
index be476ede25fb..b37fe7cd3617 100755
--- a/lib/node_modules/@stdlib/math/base/special/atan2/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/atan2/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/atan",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -63,7 +64,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/atan",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -83,7 +85,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/atan",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/atan2/src/main.c b/lib/node_modules/@stdlib/math/base/special/atan2/src/main.c
index 228c30752179..5b742c377182 100755
--- a/lib/node_modules/@stdlib/math/base/special/atan2/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/atan2/src/main.c
@@ -37,6 +37,7 @@
 #include "stdlib/math/base/special/atan.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 /**
@@ -81,7 +82,7 @@
 double stdlib_base_atan2( const double y, const double x ) {
 	double q;
 	if ( stdlib_base_is_nan( x ) ||	stdlib_base_is_nan( y ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) ) {
 		if ( x == STDLIB_CONSTANT_FLOAT64_PINF ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/atanh/manifest.json b/lib/node_modules/@stdlib/math/base/special/atanh/manifest.json
index 4678a491a177..b306081da595 100644
--- a/lib/node_modules/@stdlib/math/base/special/atanh/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/atanh/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/atanh/src/main.c b/lib/node_modules/@stdlib/math/base/special/atanh/src/main.c
index 954b9f671f17..6dc9e0ab3b9d 100644
--- a/lib/node_modules/@stdlib/math/base/special/atanh/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/atanh/src/main.c
@@ -35,6 +35,7 @@
 #include "stdlib/math/base/special/log1p.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double NEAR_ZERO = 1.0 / (1 << 28); // 2**-28
@@ -83,7 +84,7 @@ double stdlib_base_atanh( const double x ) {
 	double ax;
 	double t;
 	if ( stdlib_base_is_nan( x ) || x < -1.0 || x > 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 1.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/bernoulli/manifest.json b/lib/node_modules/@stdlib/math/base/special/bernoulli/manifest.json
index e8352697c730..3180f59fe901 100644
--- a/lib/node_modules/@stdlib/math/base/special/bernoulli/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/bernoulli/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/assert/is-odd",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/math/base/assert/is-nonnegative-integer"
+        "@stdlib/math/base/assert/is-nonnegative-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/assert/is-odd",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/math/base/assert/is-nonnegative-integer"
+        "@stdlib/math/base/assert/is-nonnegative-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/assert/is-odd",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/math/base/assert/is-nonnegative-integer"
+        "@stdlib/math/base/assert/is-nonnegative-integer",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/bernoulli/src/main.c b/lib/node_modules/@stdlib/math/base/special/bernoulli/src/main.c
index b1eb3203f5bd..c8aeaa89e989 100644
--- a/lib/node_modules/@stdlib/math/base/special/bernoulli/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/bernoulli/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_odd.h"
 #include "stdlib/constants/float64/ninf.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/bernoulli.h"
 #include <stdint.h>
 #include <stdlib.h>
@@ -172,7 +173,7 @@ int32_t MAX_BERNOULLI = 258;
 */
 double stdlib_base_bernoulli( const double n ) {
 	if ( !stdlib_base_is_nonnegative_integer( n ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( n == 1.0 ) {
 		return 0.5;
diff --git a/lib/node_modules/@stdlib/math/base/special/bessely0/manifest.json b/lib/node_modules/@stdlib/math/base/special/bessely0/manifest.json
index 0f5c4b1295c2..4795c8012d8c 100644
--- a/lib/node_modules/@stdlib/math/base/special/bessely0/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/bessely0/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pi",
-        "@stdlib/constants/float64/sqrt-pi"
+        "@stdlib/constants/float64/sqrt-pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pi",
-        "@stdlib/constants/float64/sqrt-pi"
+        "@stdlib/constants/float64/sqrt-pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pi",
-        "@stdlib/constants/float64/sqrt-pi"
+        "@stdlib/constants/float64/sqrt-pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/bessely0/src/main.c b/lib/node_modules/@stdlib/math/base/special/bessely0/src/main.c
index 9eb39650b1f1..51db37ff3fa0 100644
--- a/lib/node_modules/@stdlib/math/base/special/bessely0/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/bessely0/src/main.c
@@ -38,6 +38,7 @@
 #include "stdlib/constants/float64/ninf.h"
 #include "stdlib/constants/float64/pi.h"
 #include "stdlib/constants/float64/sqrt_pi.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double ONE_DIV_SQRT_PI = 1.0 / STDLIB_CONSTANT_FLOAT64_SQRT_PI;
 static const double TWO_DIV_PI = 2.0 / STDLIB_CONSTANT_FLOAT64_PI;
@@ -287,7 +288,7 @@ double stdlib_base_bessely0( const double x ) {
 	double z;
 
 	if ( x < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/bessely1/manifest.json b/lib/node_modules/@stdlib/math/base/special/bessely1/manifest.json
index aae8696944f8..524dcf064741 100644
--- a/lib/node_modules/@stdlib/math/base/special/bessely1/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/bessely1/manifest.json
@@ -44,7 +44,8 @@
           "@stdlib/constants/float64/pinf",
           "@stdlib/constants/float64/ninf",
           "@stdlib/constants/float64/pi",
-          "@stdlib/constants/float64/sqrt-pi"
+          "@stdlib/constants/float64/sqrt-pi",
+          "@stdlib/constants/float64/nan"
         ]
       },
       {
@@ -65,7 +66,8 @@
           "@stdlib/constants/float64/pinf",
           "@stdlib/constants/float64/ninf",
           "@stdlib/constants/float64/pi",
-          "@stdlib/constants/float64/sqrt-pi"
+          "@stdlib/constants/float64/sqrt-pi",
+          "@stdlib/constants/float64/nan"
         ]
       },
       {
@@ -86,7 +88,8 @@
           "@stdlib/constants/float64/pinf",
           "@stdlib/constants/float64/ninf",
           "@stdlib/constants/float64/pi",
-          "@stdlib/constants/float64/sqrt-pi"
+          "@stdlib/constants/float64/sqrt-pi",
+          "@stdlib/constants/float64/nan"
         ]
       }
     ]
diff --git a/lib/node_modules/@stdlib/math/base/special/bessely1/src/main.c b/lib/node_modules/@stdlib/math/base/special/bessely1/src/main.c
index 62fcc640dc1b..0560893abc62 100644
--- a/lib/node_modules/@stdlib/math/base/special/bessely1/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/bessely1/src/main.c
@@ -40,6 +40,7 @@
 #include "stdlib/constants/float64/ninf.h"
 #include "stdlib/constants/float64/pi.h"
 #include "stdlib/constants/float64/sqrt_pi.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double ONE_DIV_SQRT_PI = 1.0 / STDLIB_CONSTANT_FLOAT64_SQRT_PI;
 static const double TWO_DIV_PI = 2.0 / STDLIB_CONSTANT_FLOAT64_PI;
@@ -244,7 +245,7 @@ double stdlib_base_bessely1( double x ) {
 	double f;
 
 	if ( x < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/beta/manifest.json b/lib/node_modules/@stdlib/math/base/special/beta/manifest.json
index 64b2a1183de5..e2631db7f4aa 100644
--- a/lib/node_modules/@stdlib/math/base/special/beta/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/beta/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/pow",
         "@stdlib/constants/float64/e",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/pow",
         "@stdlib/constants/float64/e",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/pow",
         "@stdlib/constants/float64/e",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/beta/src/main.c b/lib/node_modules/@stdlib/math/base/special/beta/src/main.c
index cd7f2ae5de8f..20c10aa455fe 100644
--- a/lib/node_modules/@stdlib/math/base/special/beta/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/beta/src/main.c
@@ -38,6 +38,7 @@
 #include "stdlib/math/base/special/pow.h"
 #include "stdlib/constants/float64/e.h"
 #include "stdlib/constants/float64/eps.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double G = 10.90051099999999983936049829935654997826;
 
@@ -109,10 +110,10 @@ double stdlib_base_beta( const double a, const double b ) {
 	double c;
 
 	if ( stdlib_base_is_nan( a ) || stdlib_base_is_nan( b ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( a < 0.0 || b < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( b == 1.0 ) {
 		return 1.0 / a;
diff --git a/lib/node_modules/@stdlib/math/base/special/betaln/manifest.json b/lib/node_modules/@stdlib/math/base/special/betaln/manifest.json
index 8b449c8c22b4..85367e5c2107 100644
--- a/lib/node_modules/@stdlib/math/base/special/betaln/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/betaln/manifest.json
@@ -46,7 +46,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/constants/float64/ln-sqrt-two-pi",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -69,7 +70,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/constants/float64/ln-sqrt-two-pi",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -92,7 +94,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/constants/float64/ln-sqrt-two-pi",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/betaln/src/main.c b/lib/node_modules/@stdlib/math/base/special/betaln/src/main.c
index bd0acbe283e8..c0bc26793b46 100644
--- a/lib/node_modules/@stdlib/math/base/special/betaln/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/betaln/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/constants/float64/ln_sqrt_two_pi.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double ALGMCS[ 15 ] = {
@@ -76,7 +77,7 @@ static double dceval( const double x ) {
 	int32_t i;
 
 	if ( x < -1.1 || x > 1.1 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	b1 = 0.0;
 	b0 = 0.0;
@@ -101,7 +102,7 @@ static double dceval( const double x ) {
 */
 static double gammaCorrection( const double x ) {
 	if ( x < 10.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	// Check for underflow...
@@ -134,7 +135,7 @@ double stdlib_base_betaln( const double a, const double b ) {
 	q = stdlib_base_max( a, b );
 
 	if ( p < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( p == 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/binet/manifest.json b/lib/node_modules/@stdlib/math/base/special/binet/manifest.json
index b9595c178516..4a78423ff8cd 100644
--- a/lib/node_modules/@stdlib/math/base/special/binet/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/binet/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/math/base/special/pow",
         "@stdlib/constants/float64/phi",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/math/base/special/pow",
         "@stdlib/constants/float64/phi",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/math/base/special/pow",
         "@stdlib/constants/float64/phi",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/binet/src/main.c b/lib/node_modules/@stdlib/math/base/special/binet/src/main.c
index 07f470b67ac4..63f6acaeea0d 100644
--- a/lib/node_modules/@stdlib/math/base/special/binet/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/binet/src/main.c
@@ -23,6 +23,7 @@
 #include "stdlib/constants/float64/phi.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double SQRT_5 = 2.23606797749979;
 
@@ -49,7 +50,7 @@ double stdlib_base_binet( const double x ) {
 	double b;
 
 	if ( stdlib_base_is_nan( x ) || x == STDLIB_CONSTANT_FLOAT64_PINF || x == STDLIB_CONSTANT_FLOAT64_NINF ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	a = stdlib_base_pow( STDLIB_CONSTANT_FLOAT64_PHI, x );
 	b = stdlib_base_cospi( x ) / a;
diff --git a/lib/node_modules/@stdlib/math/base/special/binomcoef/manifest.json b/lib/node_modules/@stdlib/math/base/special/binomcoef/manifest.json
index 7b8ffed68392..fe3b44db9727 100644
--- a/lib/node_modules/@stdlib/math/base/special/binomcoef/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/binomcoef/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/max-safe-integer",
         "@stdlib/math/base/assert/is-integer",
-        "@stdlib/math/base/assert/is-odd"
+        "@stdlib/math/base/assert/is-odd",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/max-safe-integer",
         "@stdlib/math/base/assert/is-integer",
-        "@stdlib/math/base/assert/is-odd"
+        "@stdlib/math/base/assert/is-odd",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/max-safe-integer",
         "@stdlib/math/base/assert/is-integer",
-        "@stdlib/math/base/assert/is-odd"
+        "@stdlib/math/base/assert/is-odd",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/binomcoef/src/main.c b/lib/node_modules/@stdlib/math/base/special/binomcoef/src/main.c
index fd33fc6f8b0d..12326a367fda 100644
--- a/lib/node_modules/@stdlib/math/base/special/binomcoef/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/binomcoef/src/main.c
@@ -23,6 +23,7 @@
 #include "stdlib/math/base/special/gcd.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/max_safe_integer.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the binomial coefficient of two integers.
@@ -47,7 +48,7 @@ double stdlib_base_binomcoef( const double n, const double k ) {
 	double s;
 
 	if ( !stdlib_base_is_integer( n ) || !stdlib_base_is_integer( k ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	if ( k < 0 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/binomcoefln/manifest.json b/lib/node_modules/@stdlib/math/base/special/binomcoefln/manifest.json
index 6d00652ac194..6938e86c2944 100644
--- a/lib/node_modules/@stdlib/math/base/special/binomcoefln/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/binomcoefln/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/ln",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -63,7 +64,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/ln",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -83,7 +85,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/ln",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/binomcoefln/src/main.c b/lib/node_modules/@stdlib/math/base/special/binomcoefln/src/main.c
index 1383054e43b8..df42856acf90 100644
--- a/lib/node_modules/@stdlib/math/base/special/binomcoefln/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/binomcoefln/src/main.c
@@ -22,6 +22,7 @@
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the natural logarithm of the binomial coefficient of two integers.
@@ -36,7 +37,7 @@
 */
 double stdlib_base_binomcoefln( const double n, const double k ) {
 	if ( !stdlib_base_is_integer( n ) || !stdlib_base_is_integer( k ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( n < 0 ) {
 		return stdlib_base_binomcoefln( -n + k - 1, k );
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcox/manifest.json b/lib/node_modules/@stdlib/math/base/special/boxcox/manifest.json
index 464a6b0dc65d..c32050a85727 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcox/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/boxcox/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/math/base/assert/is-positive-zero",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/expm1",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/math/base/assert/is-positive-zero",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/expm1",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/math/base/assert/is-positive-zero",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/expm1",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcox/src/main.c b/lib/node_modules/@stdlib/math/base/special/boxcox/src/main.c
index 1def199b1e2c..aa7fb9b7afa4 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcox/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/boxcox/src/main.c
@@ -23,6 +23,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/math/base/special/expm1.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes a one-parameter Box-Cox transformation.
@@ -43,7 +44,7 @@
 */
 double stdlib_base_boxcox( const double x, const double lambda ) {
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( lambda ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_positive_zero( x ) && lambda < 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcox1p/manifest.json b/lib/node_modules/@stdlib/math/base/special/boxcox1p/manifest.json
index ccc37a013671..0004c8a0f759 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcox1p/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/boxcox1p/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/expm1",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/expm1",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/expm1",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcox1p/src/main.c b/lib/node_modules/@stdlib/math/base/special/boxcox1p/src/main.c
index 6aceb3a5c228..9921c0a10c25 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcox1p/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/boxcox1p/src/main.c
@@ -22,6 +22,7 @@
 #include "stdlib/math/base/special/log1p.h"
 #include "stdlib/math/base/special/expm1.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes a one-parameter Box-Cox transformation of `1+x`.
@@ -38,7 +39,7 @@ double stdlib_base_boxcox1p( const double x, const double lambda ) {
 	double lgx;
 
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( lambda ) || x < -1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == -1.0 && lambda < 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/manifest.json b/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/manifest.json
index a54b735fc6bb..d50b595987be 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/special/expm1"
+        "@stdlib/math/base/special/expm1",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/special/expm1"
+        "@stdlib/math/base/special/expm1",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/special/expm1"
+        "@stdlib/math/base/special/expm1",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/src/main.c b/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/src/main.c
index 34f8e87e6802..4b07d856aa54 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/boxcox1pinv/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/boxcox1pinv.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/expm1.h"
 #include "stdlib/math/base/special/log1p.h"
 #include "stdlib/math/base/special/abs.h"
@@ -37,7 +38,7 @@ double stdlib_base_boxcox1pinv( const double y, const double lambda ) {
 	double ly;
 
 	if ( stdlib_base_is_nan( y ) || stdlib_base_is_nan( lambda ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( lambda == 0.0 ) {
 		return stdlib_base_expm1( y );
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcoxinv/manifest.json b/lib/node_modules/@stdlib/math/base/special/boxcoxinv/manifest.json
index 0c074e7620e8..4c6fef4cf5c1 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcoxinv/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/boxcoxinv/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/boxcoxinv/src/main.c b/lib/node_modules/@stdlib/math/base/special/boxcoxinv/src/main.c
index 481e82e9edb4..32f8ba28e85a 100644
--- a/lib/node_modules/@stdlib/math/base/special/boxcoxinv/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/boxcoxinv/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/boxcoxinv.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/exp.h"
 #include "stdlib/math/base/special/log1p.h"
 #include "stdlib/math/base/assert/is_nan.h"
@@ -34,7 +35,7 @@
 */
 double stdlib_base_boxcoxinv( const double y, const double lambda ) {
 	if ( stdlib_base_is_nan( y ) || stdlib_base_is_nan( lambda ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( lambda == 0.0 ) {
 		return stdlib_base_exp( y );
diff --git a/lib/node_modules/@stdlib/math/base/special/ceilb/manifest.json b/lib/node_modules/@stdlib/math/base/special/ceilb/manifest.json
index 0ecff79da0d4..4ee3e4824ccc 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceilb/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ceilb/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/ceil",
         "@stdlib/math/base/special/ceiln",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/ceil",
         "@stdlib/math/base/special/ceiln",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/ceil",
         "@stdlib/math/base/special/ceiln",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ceilb/src/main.c b/lib/node_modules/@stdlib/math/base/special/ceilb/src/main.c
index fd3bb060f180..7b550c544f91 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceilb/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ceilb/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/ceilb.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/ceil.h"
@@ -41,7 +42,7 @@ double stdlib_base_ceilb( const double x, const int32_t n, const int32_t b ) {
 	double s;
 
 	if ( stdlib_base_is_nan( x ) || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/ceilsd/manifest.json b/lib/node_modules/@stdlib/math/base/special/ceilsd/manifest.json
index 1b5e4b62f055..160ae05e6737 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceilsd/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ceilsd/manifest.json
@@ -45,7 +45,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/ceil",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -67,7 +68,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/ceil",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -89,7 +91,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/ceil",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ceilsd/src/main.c b/lib/node_modules/@stdlib/math/base/special/ceilsd/src/main.c
index 74065dd3715c..e1ee36218813 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceilsd/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ceilsd/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/ceilsd.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/pow.h"
@@ -46,7 +47,7 @@ double stdlib_base_ceilsd( const double x, const int32_t n, const int32_t b ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || n < 1 || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/cexp/manifest.json b/lib/node_modules/@stdlib/math/base/special/cexp/manifest.json
index ddfa5646c158..759ce59d5adb 100644
--- a/lib/node_modules/@stdlib/math/base/special/cexp/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/cexp/manifest.json
@@ -45,7 +45,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/copysign",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -67,7 +68,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/copysign",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -89,7 +91,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/copysign",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/cexp/src/main.c b/lib/node_modules/@stdlib/math/base/special/cexp/src/main.c
index 44ef8126a47f..a144517293e4 100644
--- a/lib/node_modules/@stdlib/math/base/special/cexp/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/cexp/src/main.c
@@ -24,6 +24,7 @@
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/math/base/special/copysign.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/complex/float64/ctor.h"
 #include "stdlib/complex/float64/reim.h"
 
@@ -55,18 +56,18 @@ stdlib_complex128_t stdlib_base_cexp( const stdlib_complex128_t z ) {
 	stdlib_complex128_reim( z, &re, &im );
 
 	if ( stdlib_base_is_nan( re ) ) {
-		re = 0.0 / 0.0; // NaN
+		re = STDLIB_CONSTANT_FLOAT64_NAN;
 		im = ( im == 0.0 ) ? im : re;
 	} else if ( stdlib_base_is_infinite( im ) ) {
 		if ( re == STDLIB_CONSTANT_FLOAT64_PINF ) {
 			re = -re;
-			im = 0.0 / 0.0; // NaN
+			im = STDLIB_CONSTANT_FLOAT64_NAN;
 		} else if ( re == STDLIB_CONSTANT_FLOAT64_NINF ) {
 			re = -0.0;
 			im = stdlib_base_copysign( 0.0, im );
 		} else {
-			re = 0.0 / 0.0; // NaN
-			im = 0.0 / 0.0; // NaN
+			re = STDLIB_CONSTANT_FLOAT64_NAN;
+			im = STDLIB_CONSTANT_FLOAT64_NAN;
 		}
 	} else {
 		e = stdlib_base_exp( re );
diff --git a/lib/node_modules/@stdlib/math/base/special/clamp/manifest.json b/lib/node_modules/@stdlib/math/base/special/clamp/manifest.json
index a5c30237a9e4..f32c1f2e7f67 100644
--- a/lib/node_modules/@stdlib/math/base/special/clamp/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/clamp/manifest.json
@@ -40,7 +40,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/ternary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/clamp/src/main.c b/lib/node_modules/@stdlib/math/base/special/clamp/src/main.c
index defb3767570b..792246d83458 100644
--- a/lib/node_modules/@stdlib/math/base/special/clamp/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/clamp/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/clamp.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_negative_zero.h"
 
@@ -42,7 +43,7 @@ double stdlib_base_clamp( const double v, const double min, const double max ) {
 		stdlib_base_is_nan( min ) ||
 		stdlib_base_is_nan( max )
 	) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Simple cases...
 	if ( v < min ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/cos/manifest.json b/lib/node_modules/@stdlib/math/base/special/cos/manifest.json
index cf19df411559..b513eb55cdf6 100644
--- a/lib/node_modules/@stdlib/math/base/special/cos/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/cos/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/kernel-sin",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/kernel-sin",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/kernel-sin",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/cos/src/main.c b/lib/node_modules/@stdlib/math/base/special/cos/src/main.c
index ce0d6a3c74a1..e73b9d7680ea 100644
--- a/lib/node_modules/@stdlib/math/base/special/cos/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/cos/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/number/float64/base/get_high_word.h"
 #include "stdlib/constants/float64/high_word_abs_mask.h"
 #include "stdlib/constants/float64/high_word_exponent_mask.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/kernel_cos.h"
 #include "stdlib/math/base/special/kernel_sin.h"
 #include "stdlib/math/base/special/rempio2.h"
@@ -90,7 +91,7 @@ double stdlib_base_cos( const double x ) {
 	}
 	// Case: cos(Inf or NaN) is NaN */
 	if ( ix >= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_EXPONENT_MASK ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Case: Argument reduction needed...
 	n = stdlib_base_rempio2( x, &Y[ 0 ], &Y[ 1 ] );
diff --git a/lib/node_modules/@stdlib/math/base/special/cosd/manifest.json b/lib/node_modules/@stdlib/math/base/special/cosd/manifest.json
index 1c40675a5d6e..6daba4116798 100644
--- a/lib/node_modules/@stdlib/math/base/special/cosd/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/cosd/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/math/base/special/kernel-sin",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -63,7 +64,8 @@
         "@stdlib/math/base/special/kernel-sin",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -83,7 +85,8 @@
         "@stdlib/math/base/special/kernel-sin",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/cosd/src/main.c b/lib/node_modules/@stdlib/math/base/special/cosd/src/main.c
index 3e4dfd667a0d..b1fc3944300e 100644
--- a/lib/node_modules/@stdlib/math/base/special/cosd/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/cosd/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/cosd.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/kernel_sin.h"
 #include "stdlib/math/base/special/kernel_cos.h"
 #include "stdlib/math/base/special/deg2rad.h"
@@ -39,7 +40,7 @@ double stdlib_base_cosd( const double x ) {
 	double rx;
 
 	if ( stdlib_base_is_infinite( x ) || stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	rx = stdlib_base_abs( stdlib_base_fmod( x, 360.0 ) );
diff --git a/lib/node_modules/@stdlib/math/base/special/cospi/manifest.json b/lib/node_modules/@stdlib/math/base/special/cospi/manifest.json
index f96cd620858a..17a785396fde 100644
--- a/lib/node_modules/@stdlib/math/base/special/cospi/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/cospi/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/fmod",
         "@stdlib/constants/float64/pi",
-        "@stdlib/constants/float64/max-safe-integer"
+        "@stdlib/constants/float64/max-safe-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +72,8 @@
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/fmod",
         "@stdlib/constants/float64/pi",
-        "@stdlib/constants/float64/max-safe-integer"
+        "@stdlib/constants/float64/max-safe-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -95,7 +97,8 @@
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/fmod",
         "@stdlib/constants/float64/pi",
-        "@stdlib/constants/float64/max-safe-integer"
+        "@stdlib/constants/float64/max-safe-integer",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/cospi/src/main.c b/lib/node_modules/@stdlib/math/base/special/cospi/src/main.c
index b51c47377cb5..1350787be0ac 100644
--- a/lib/node_modules/@stdlib/math/base/special/cospi/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/cospi/src/main.c
@@ -26,6 +26,7 @@
 #include "stdlib/math/base/special/fmod.h"
 #include "stdlib/constants/float64/pi.h"
 #include "stdlib/constants/float64/max_safe_integer.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double MAX_INTEGER_P1 = STDLIB_CONSTANT_FLOAT64_MAX_SAFE_INTEGER + 1;
 
@@ -54,7 +55,7 @@ double stdlib_base_cospi( const double x ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_infinite( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	ax = stdlib_base_abs( x );
 	if ( ax > MAX_INTEGER_P1 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/digamma/manifest.json b/lib/node_modules/@stdlib/math/base/special/digamma/manifest.json
index bbc94721cd33..e532992995b9 100644
--- a/lib/node_modules/@stdlib/math/base/special/digamma/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/digamma/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/tan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/tan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/tan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/digamma/src/main.c b/lib/node_modules/@stdlib/math/base/special/digamma/src/main.c
index 70fdf933ccec..c6678f53a814 100644
--- a/lib/node_modules/@stdlib/math/base/special/digamma/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/digamma/src/main.c
@@ -35,6 +35,7 @@
 #include "stdlib/math/base/special/tan.h"
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pi.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double MIN_SAFE_ASYMPTOTIC = 10.0; // BIG!
 static const double root1 = 1569415565.0 / 1073741824.0;
@@ -223,7 +224,7 @@ double stdlib_base_digamma( const double x ) {
 	double xc;
 
 	if ( stdlib_base_is_nan( x ) || x == 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	// If `x` is negative, use reflection...
@@ -242,7 +243,7 @@ double stdlib_base_digamma( const double x ) {
 
 		// Check for evaluation at a negative pole:
 		if ( rem == 0.0 ) {
-			return 0.0 / 0.0; // NaN
+			return STDLIB_CONSTANT_FLOAT64_NAN;
 		}
 		tmp = STDLIB_CONSTANT_FLOAT64_PI / stdlib_base_tan( STDLIB_CONSTANT_FLOAT64_PI * rem );
 	} else {
diff --git a/lib/node_modules/@stdlib/math/base/special/dirac-delta/manifest.json b/lib/node_modules/@stdlib/math/base/special/dirac-delta/manifest.json
index 8dee620b1e2c..d113d3f6517a 100644
--- a/lib/node_modules/@stdlib/math/base/special/dirac-delta/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/dirac-delta/manifest.json
@@ -38,7 +38,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -53,7 +54,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +70,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/dirac-delta/src/main.c b/lib/node_modules/@stdlib/math/base/special/dirac-delta/src/main.c
index 7cf9b85352e6..789ecf3bdaf8 100644
--- a/lib/node_modules/@stdlib/math/base/special/dirac-delta/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/dirac-delta/src/main.c
@@ -19,6 +19,7 @@
 #include "stdlib/math/base/special/dirac_delta.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Evaluates the Dirac delta function.
@@ -32,7 +33,7 @@
 */
 double stdlib_base_dirac_delta( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/manifest.json b/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/manifest.json
index 629c9eb58123..987fb3882fd4 100644
--- a/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/special/powm1",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/riemann-zeta",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/special/powm1",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/riemann-zeta",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/special/powm1",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/riemann-zeta",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/src/main.c b/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/src/main.c
index 22fbf425edbf..be1e3a2799da 100644
--- a/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/dirichlet-eta/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/powm1.h"
 #include "stdlib/constants/float64/ln_two.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/riemann_zeta.h"
 
 /**
@@ -34,7 +35,7 @@
 */
 double stdlib_base_eta( const double s ) {
 	if ( stdlib_base_is_nan( s ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( s == 1.0 ) {
 		// Alternating harmonic series...
diff --git a/lib/node_modules/@stdlib/math/base/special/ellipe/manifest.json b/lib/node_modules/@stdlib/math/base/special/ellipe/manifest.json
index 9a8d6384fc80..6bbec8c25224 100644
--- a/lib/node_modules/@stdlib/math/base/special/ellipe/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ellipe/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/special/ellipk"
+        "@stdlib/math/base/special/ellipk",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/special/ellipk"
+        "@stdlib/math/base/special/ellipk",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/special/ellipk"
+        "@stdlib/math/base/special/ellipk",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ellipe/src/main.c b/lib/node_modules/@stdlib/math/base/special/ellipe/src/main.c
index fc6bfd9663e5..92324c26cf1a 100644
--- a/lib/node_modules/@stdlib/math/base/special/ellipe/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ellipe/src/main.c
@@ -53,6 +53,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/half_pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 /* Begin auto-generated functions. The following functions are auto-generated. Do not edit directly. */
@@ -331,7 +332,7 @@ double stdlib_base_ellipe( const double m ) {
 		return 1.0;
 	}
 	if ( x > 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x < 0.1 ) {
 		t = poly_p1( x - 0.05 );
diff --git a/lib/node_modules/@stdlib/math/base/special/ellipj/manifest.json b/lib/node_modules/@stdlib/math/base/special/ellipj/manifest.json
index 00be502cf2a4..8442d7cc1b36 100644
--- a/lib/node_modules/@stdlib/math/base/special/ellipj/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ellipj/manifest.json
@@ -55,7 +55,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/sqrt-eps",
         "@stdlib/constants/float64/eps",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -84,7 +85,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/sqrt-eps",
         "@stdlib/constants/float64/eps",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -113,7 +115,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/sqrt-eps",
         "@stdlib/constants/float64/eps",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ellipj/src/main.c b/lib/node_modules/@stdlib/math/base/special/ellipj/src/main.c
index 70393ba4e742..dcc1371f7b74 100644
--- a/lib/node_modules/@stdlib/math/base/special/ellipj/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ellipj/src/main.c
@@ -32,6 +32,7 @@
 #include "stdlib/constants/float64/sqrt_eps.h"
 #include "stdlib/constants/float64/eps.h"
 #include "stdlib/constants/float64/pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdbool.h>
 #include <stdint.h>
 
@@ -100,7 +101,7 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 		*sn = ( out[ 0 ] / out[ 2 ] ) / k1inv;
 		*cn = out[ 1 ] / out[ 2 ];
 		*dn = 1.0 / out[ 2 ];
-		*am = 0.0 / 0.0; // NaN
+		*am = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	} else if ( m > 1.0 ) {
 		// A&S 16.11.1 for reciprocal parameter, mapping m > 1 to 0 < mu < 1:
@@ -109,7 +110,7 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 		*sn = out[ 0 ] / k;
 		*cn = out[ 2 ];
 		*dn = out[ 1 ];
-		*am = 0.0 / 0.0; // NaN
+		*am = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	} else if ( m == 0.0 ) {
 		// A&S table 16.6, limiting case m = 0: circular trigonometric functions:
@@ -179,10 +180,10 @@ void stdlib_base_ellipj( const double u, const double m, double* sn, double* cn,
 			if ( N > 8 ) {
 				// Warning: Overflow encountered in iteration. Returning NaN for all output values:
 				NANFLG = true;
-				*sn = 0.0 / 0.0; // NaN
-				*cn = 0.0 / 0.0; // NaN
-				*dn = 0.0 / 0.0; // NaN
-				*am = 0.0 / 0.0; // NaN
+				*sn = STDLIB_CONSTANT_FLOAT64_NAN;
+				*cn = STDLIB_CONSTANT_FLOAT64_NAN;
+				*dn = STDLIB_CONSTANT_FLOAT64_NAN;
+				*am = STDLIB_CONSTANT_FLOAT64_NAN;
 				break;
 			}
 			atmp = ( a + b ) * 0.5;
diff --git a/lib/node_modules/@stdlib/math/base/special/ellipk/manifest.json b/lib/node_modules/@stdlib/math/base/special/ellipk/manifest.json
index 436f90cc4854..a5332ec2fcfe 100644
--- a/lib/node_modules/@stdlib/math/base/special/ellipk/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ellipk/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/constants/float64/half-pi",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/constants/float64/half-pi",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/constants/float64/half-pi",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ellipk/src/main.c b/lib/node_modules/@stdlib/math/base/special/ellipk/src/main.c
index 02df49586be0..384320054227 100644
--- a/lib/node_modules/@stdlib/math/base/special/ellipk/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ellipk/src/main.c
@@ -52,6 +52,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/half_pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double ONE_DIV_PI = 0.3183098861837907;
@@ -358,7 +359,7 @@ double stdlib_base_ellipk( const double m ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
 	}
 	if ( x > 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x < 0.1 ) {
 		t = poly_p1( x - 0.05 );
diff --git a/lib/node_modules/@stdlib/math/base/special/erf/manifest.json b/lib/node_modules/@stdlib/math/base/special/erf/manifest.json
index 461e8c261ba7..55e8af143481 100644
--- a/lib/node_modules/@stdlib/math/base/special/erf/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/erf/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -62,7 +63,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -81,7 +83,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -100,7 +103,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/erf/src/main.c b/lib/node_modules/@stdlib/math/base/special/erf/src/main.c
index e2762eab12af..4a139114eb88 100644
--- a/lib/node_modules/@stdlib/math/base/special/erf/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/erf/src/main.c
@@ -36,6 +36,7 @@
 #include "stdlib/number/float64/base/set_low_word.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double TINY = 1.0e-300;
@@ -457,7 +458,7 @@ double stdlib_base_erf( const double x ) {
 
 	// Special case: NaN
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Special case: +infinity
 	if ( x == STDLIB_CONSTANT_FLOAT64_PINF ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/erfc/manifest.json b/lib/node_modules/@stdlib/math/base/special/erfc/manifest.json
index 4c4cf755fe88..d81a9b57b5c4 100644
--- a/lib/node_modules/@stdlib/math/base/special/erfc/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/erfc/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/exp",
-        "@stdlib/number/float64/base/set-low-word"
+        "@stdlib/number/float64/base/set-low-word",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/erfc/src/main.c b/lib/node_modules/@stdlib/math/base/special/erfc/src/main.c
index cb9b4a3b4482..738d157eef6f 100644
--- a/lib/node_modules/@stdlib/math/base/special/erfc/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/erfc/src/main.c
@@ -36,6 +36,7 @@
 #include "stdlib/number/float64/base/set_low_word.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double TINY = 1.0e-300;
@@ -454,7 +455,7 @@ double stdlib_base_erfc( const double x ) {
 
 	// Special case: NaN
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Special case: +infinity
 	if ( x == STDLIB_CONSTANT_FLOAT64_PINF ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/erfcinv/manifest.json b/lib/node_modules/@stdlib/math/base/special/erfcinv/manifest.json
index 97742cc77f41..3e7183e5a7a8 100644
--- a/lib/node_modules/@stdlib/math/base/special/erfcinv/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/erfcinv/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/erfcinv/src/main.c b/lib/node_modules/@stdlib/math/base/special/erfcinv/src/main.c
index 8e4e61e66cef..ee0660485870 100644
--- a/lib/node_modules/@stdlib/math/base/special/erfcinv/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/erfcinv/src/main.c
@@ -35,6 +35,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double Y1 = 8.91314744949340820313e-2;
@@ -340,7 +341,7 @@ double stdlib_base_erfcinv( const double x ) {
 
 	// Special case: NaN
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Special case: 0
 	if ( x == 0.0 ) {
@@ -355,7 +356,7 @@ double stdlib_base_erfcinv( const double x ) {
 		return 0.0;
 	}
 	if ( x > 2.0 || x < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Argument reduction (reduce to interval [0,1]). If `x` is outside [0,1], we can take advantage of the complementary error function reflection formula: `erfc(-z) = 2 - erfc(z)`, by negating the result once finished.
 	if ( x > 1.0 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/erfinv/manifest.json b/lib/node_modules/@stdlib/math/base/special/erfinv/manifest.json
index 97742cc77f41..3e7183e5a7a8 100644
--- a/lib/node_modules/@stdlib/math/base/special/erfinv/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/erfinv/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/math/base/special/ln"
+        "@stdlib/math/base/special/ln",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/erfinv/src/main.c b/lib/node_modules/@stdlib/math/base/special/erfinv/src/main.c
index d4ec84de8e3b..a088d3d1b6fa 100644
--- a/lib/node_modules/@stdlib/math/base/special/erfinv/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/erfinv/src/main.c
@@ -35,6 +35,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double Y1 = 8.91314744949340820313e-2;
@@ -272,7 +273,7 @@ double stdlib_base_erfinv( const double x ) {
 
 	// Special case: NaN
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Special case: 1
 	if ( x == 1.0 ) {
@@ -288,7 +289,7 @@ double stdlib_base_erfinv( const double x ) {
 	}
 	// Special case: |x| > 1 (range error)
 	if ( x > 1.0 || x < -1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Argument reduction (reduce to interval [0,1]). If `x` is negative, we can safely negate the value, taking advantage of the error function being an odd function; i.e., `erf(-x) = -erf(x)`.
 	if ( x < 0.0 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/exp10/manifest.json b/lib/node_modules/@stdlib/math/base/special/exp10/manifest.json
index 22c8c12ca3bd..89e4c46efeae 100644
--- a/lib/node_modules/@stdlib/math/base/special/exp10/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/exp10/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/ldexp",
-        "@stdlib/constants/float64/min-base10-exponent"
+        "@stdlib/constants/float64/min-base10-exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/ldexp",
-        "@stdlib/constants/float64/min-base10-exponent"
+        "@stdlib/constants/float64/min-base10-exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/ldexp",
-        "@stdlib/constants/float64/min-base10-exponent"
+        "@stdlib/constants/float64/min-base10-exponent",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/exp10/src/main.c b/lib/node_modules/@stdlib/math/base/special/exp10/src/main.c
index c8eea140c600..42ff65360c46 100644
--- a/lib/node_modules/@stdlib/math/base/special/exp10/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/exp10/src/main.c
@@ -37,6 +37,7 @@
 #include "stdlib/constants/float64/max_base10_exponent.h"
 #include "stdlib/constants/float64/min_base10_exponent.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double LOG210 = 3.32192809488736234787e0;
@@ -123,7 +124,7 @@ double stdlib_base_exp10( const double x ) {
 	int32_t n;
 
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x > STDLIB_CONSTANT_FLOAT64_MAX_BASE10_EXPONENT ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/expit/manifest.json b/lib/node_modules/@stdlib/math/base/special/expit/manifest.json
index 8ca7869951dd..84884b62f51f 100644
--- a/lib/node_modules/@stdlib/math/base/special/expit/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/expit/manifest.json
@@ -38,7 +38,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/exp"
+        "@stdlib/math/base/special/exp",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -53,7 +54,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/exp"
+        "@stdlib/math/base/special/exp",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +70,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/exp"
+        "@stdlib/math/base/special/exp",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/expit/src/main.c b/lib/node_modules/@stdlib/math/base/special/expit/src/main.c
index 2ebf91c958d1..0a74ba0f7418 100644
--- a/lib/node_modules/@stdlib/math/base/special/expit/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/expit/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/expit.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/exp.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
@@ -32,7 +33,7 @@
 */
 double stdlib_base_expit( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return 1.0 / ( 1.0 + stdlib_base_exp( -x ) );
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/factorial/manifest.json b/lib/node_modules/@stdlib/math/base/special/factorial/manifest.json
index 5b4e388f472f..e33219b192a2 100644
--- a/lib/node_modules/@stdlib/math/base/special/factorial/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/factorial/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-integer",
         "@stdlib/math/base/special/gamma",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-integer",
         "@stdlib/math/base/special/gamma",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-integer",
         "@stdlib/math/base/special/gamma",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/factorial/src/main.c b/lib/node_modules/@stdlib/math/base/special/factorial/src/main.c
index c12d0da63bb0..51d438e19c8f 100644
--- a/lib/node_modules/@stdlib/math/base/special/factorial/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/factorial/src/main.c
@@ -22,6 +22,7 @@
 #include "stdlib/math/base/special/gamma.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/max_nth_factorial.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double FACTORIALS[ 171 ] = {
@@ -211,11 +212,11 @@ static const double FACTORIALS[ 171 ] = {
 */
 double stdlib_base_factorial( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_integer( x ) ) {
 		if ( x < 0.0 ) {
-			return 0.0 / 0.0; // NaN
+			return STDLIB_CONSTANT_FLOAT64_NAN;
 		}
 		if ( x <= STDLIB_CONSTANT_FLOAT64_MAX_NTH_FACTORIAL ) {
 			return FACTORIALS[ (int32_t)x ];
diff --git a/lib/node_modules/@stdlib/math/base/special/factorial2/manifest.json b/lib/node_modules/@stdlib/math/base/special/factorial2/manifest.json
index ba5469173d4c..e748db7d9200 100644
--- a/lib/node_modules/@stdlib/math/base/special/factorial2/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/factorial2/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/max-nth-double-factorial",
         "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/max-nth-double-factorial",
         "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/constants/float64/pinf",
         "@stdlib/constants/float64/max-nth-double-factorial",
         "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/factorial2/src/main.c b/lib/node_modules/@stdlib/math/base/special/factorial2/src/main.c
index e50252653a33..d440184dfa76 100644
--- a/lib/node_modules/@stdlib/math/base/special/factorial2/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/factorial2/src/main.c
@@ -22,6 +22,7 @@
 #include "stdlib/math/base/assert/is_nonnegative_integer.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/max_nth_double_factorial.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Evaluates the double factorial of `n`.
@@ -40,7 +41,7 @@ double stdlib_base_factorial2( const double n ) {
 	double i;
 
 	if ( stdlib_base_is_nan( n ) || !stdlib_base_is_nonnegative_integer( n ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( n > STDLIB_CONSTANT_FLOAT64_MAX_NTH_DOUBLE_FACTORIAL ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/factorialln/manifest.json b/lib/node_modules/@stdlib/math/base/special/factorialln/manifest.json
index 213d4fd06e9f..48f9714d9802 100644
--- a/lib/node_modules/@stdlib/math/base/special/factorialln/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/factorialln/manifest.json
@@ -38,7 +38,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-negative-integer",
-        "@stdlib/math/base/special/gammaln"
+        "@stdlib/math/base/special/gammaln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -53,7 +54,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-negative-integer",
-        "@stdlib/math/base/special/gammaln"
+        "@stdlib/math/base/special/gammaln",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +70,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-negative-integer",
-        "@stdlib/math/base/special/gammaln"
+        "@stdlib/math/base/special/gammaln",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/factorialln/src/main.c b/lib/node_modules/@stdlib/math/base/special/factorialln/src/main.c
index 596b05ad8e03..c8be4e611932 100644
--- a/lib/node_modules/@stdlib/math/base/special/factorialln/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/factorialln/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/factorialln.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_negative_integer.h"
 #include "stdlib/math/base/special/gammaln.h"
 
@@ -32,7 +33,7 @@
 */
 double stdlib_base_factorialln( const double x ) {
 	if ( stdlib_base_is_negative_integer( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return stdlib_base_gammaln( x + 1.0 );
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/falling-factorial/manifest.json b/lib/node_modules/@stdlib/math/base/special/falling-factorial/manifest.json
index a3c2f970e229..384bd26dbb0c 100644
--- a/lib/node_modules/@stdlib/math/base/special/falling-factorial/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/falling-factorial/manifest.json
@@ -45,7 +45,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/max",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -67,7 +68,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/max",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -89,7 +91,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/max",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/falling-factorial/src/main.c b/lib/node_modules/@stdlib/math/base/special/falling-factorial/src/main.c
index 1515fcb958d2..316900bfa5d9 100644
--- a/lib/node_modules/@stdlib/math/base/special/falling-factorial/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/falling-factorial/src/main.c
@@ -39,6 +39,7 @@
 #include "stdlib/constants/float64/max.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/max_nth_factorial.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 #include <stdbool.h>
 
@@ -78,7 +79,7 @@ static double risingFactorial( const double x, const int32_t n ) {
 	bool inv;
 
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	nc = n;
 	xc = x;
@@ -146,7 +147,7 @@ double stdlib_base_falling_factorial( const double x, const int32_t n ) {
 	double xc;
 
 	if ( stdlib_base_is_nan( x ) || !stdlib_base_is_nonnegative_integer( n ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		return 0.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/fast/acosh/manifest.json b/lib/node_modules/@stdlib/math/base/special/fast/acosh/manifest.json
index 26ed895f38cb..148b8bfa1e1b 100644
--- a/lib/node_modules/@stdlib/math/base/special/fast/acosh/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/fast/acosh/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/assert/is-infinite"
+        "@stdlib/math/base/assert/is-infinite",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/assert/is-infinite"
+        "@stdlib/math/base/assert/is-infinite",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/assert/is-infinite"
+        "@stdlib/math/base/assert/is-infinite",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/fast/acosh/src/main.c b/lib/node_modules/@stdlib/math/base/special/fast/acosh/src/main.c
index aadb31e8e2bc..13335f8de861 100644
--- a/lib/node_modules/@stdlib/math/base/special/fast/acosh/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/fast/acosh/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/fast/acosh.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/ln.h"
@@ -35,7 +36,7 @@
 */
 double stdlib_base_fast_acosh( const double x ) {
 	if ( x < 1.0 ) {
-		return 0.0 / 0.0 ; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_infinite( x ) ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/fast/atanh/manifest.json b/lib/node_modules/@stdlib/math/base/special/fast/atanh/manifest.json
index 3086e38c4d7e..9f25da00bea8 100644
--- a/lib/node_modules/@stdlib/math/base/special/fast/atanh/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/fast/atanh/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/assert/is-infinite"
+        "@stdlib/math/base/assert/is-infinite",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/assert/is-infinite"
+        "@stdlib/math/base/assert/is-infinite",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/math/base/assert/is-infinite"
+        "@stdlib/math/base/assert/is-infinite",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/fast/atanh/src/main.c b/lib/node_modules/@stdlib/math/base/special/fast/atanh/src/main.c
index a814f378b4b2..ef5701352fb6 100644
--- a/lib/node_modules/@stdlib/math/base/special/fast/atanh/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/fast/atanh/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/fast/atanh.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/ln.h"
@@ -37,7 +38,7 @@ double stdlib_base_fast_atanh( const double x ) {
 		return x;
 	}
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_infinite( x ) ) {
-		return 0.0 / 0.0 ; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return 0.5 * stdlib_base_ln( ( 1.0 + x ) / ( 1.0 - x ) );
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/fast/pow-int/manifest.json b/lib/node_modules/@stdlib/math/base/special/fast/pow-int/manifest.json
index 880ad9c327f1..333866c685b1 100644
--- a/lib/node_modules/@stdlib/math/base/special/fast/pow-int/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/fast/pow-int/manifest.json
@@ -40,7 +40,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/fast/pow-int/src/main.c b/lib/node_modules/@stdlib/math/base/special/fast/pow-int/src/main.c
index b5ca5a2a4c3e..3666ced52236 100644
--- a/lib/node_modules/@stdlib/math/base/special/fast/pow-int/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/fast/pow-int/src/main.c
@@ -19,6 +19,7 @@
 #include "stdlib/math/base/special/fast/pow_int.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 
 int32_t ONE = 1;
 int32_t ZERO = 0;
@@ -42,7 +43,7 @@ double stdlib_base_fast_pow( const double x, const int32_t y ) {
 	xc = x;
 	yc = y;
 	if ( stdlib_base_is_nan( xc ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( yc < ZERO ) {
 		yc = -yc;
diff --git a/lib/node_modules/@stdlib/math/base/special/fibonacci-index/manifest.json b/lib/node_modules/@stdlib/math/base/special/fibonacci-index/manifest.json
index d9f4316aeb68..31a5d80d4ce6 100644
--- a/lib/node_modules/@stdlib/math/base/special/fibonacci-index/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/fibonacci-index/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/math/base/special/round",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/assert/is-integer",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/math/base/special/round",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/assert/is-integer",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/math/base/special/round",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/assert/is-integer",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/fibonacci-index/src/main.c b/lib/node_modules/@stdlib/math/base/special/fibonacci-index/src/main.c
index 5983d6e98f30..57b4a0651300 100644
--- a/lib/node_modules/@stdlib/math/base/special/fibonacci-index/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/fibonacci-index/src/main.c
@@ -23,6 +23,7 @@
 #include "stdlib/math/base/special/round.h"
 #include "stdlib/constants/float64/phi.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdbool.h>
 
 static const double SQRT_5 = 2.23606797749979;
@@ -45,7 +46,7 @@ double stdlib_base_fibonacci_index( const double F ) {
 	double x;
 
 	if ( stdlib_base_is_nan( F ) || stdlib_base_is_integer( F ) == false || F <= 1 || F == STDLIB_CONSTANT_FLOAT64_PINF ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	x = ( F * SQRT_5 ) + 0.5;
 	return stdlib_base_round( stdlib_base_ln( x ) / stdlib_base_ln( STDLIB_CONSTANT_FLOAT64_PHI ) );
diff --git a/lib/node_modules/@stdlib/math/base/special/floorb/manifest.json b/lib/node_modules/@stdlib/math/base/special/floorb/manifest.json
index 5c67bfd871d4..80a774e7fb15 100644
--- a/lib/node_modules/@stdlib/math/base/special/floorb/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/floorb/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/floorn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/floorn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/floorn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/floorb/src/main.c b/lib/node_modules/@stdlib/math/base/special/floorb/src/main.c
index 44671bba319a..14b36aa7013e 100644
--- a/lib/node_modules/@stdlib/math/base/special/floorb/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/floorb/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/floorb.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/floor.h"
@@ -41,7 +42,7 @@ double stdlib_base_floorb( const double x, const int32_t n, const int32_t b ) {
 	double s;
 
 	if ( stdlib_base_is_nan( x ) || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/floorsd/manifest.json b/lib/node_modules/@stdlib/math/base/special/floorsd/manifest.json
index 4fce81c4238d..2a189f19e0e0 100644
--- a/lib/node_modules/@stdlib/math/base/special/floorsd/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/floorsd/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/floorsd/src/main.c b/lib/node_modules/@stdlib/math/base/special/floorsd/src/main.c
index b87b9ea9d277..455b018651bc 100644
--- a/lib/node_modules/@stdlib/math/base/special/floorsd/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/floorsd/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/floorsd.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/pow.h"
@@ -45,7 +46,7 @@ double stdlib_base_floorsd( const double x, const int32_t n, const int32_t b ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || n < 1 || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json b/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json
index 0b051015bd59..95be3c49993e 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json
@@ -49,7 +49,8 @@
         "@stdlib/constants/float64/eulergamma",
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/math/base/special/pow",
-        "@stdlib/math/base/special/exp"
+        "@stdlib/math/base/special/exp",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -75,7 +76,8 @@
         "@stdlib/constants/float64/eulergamma",
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/math/base/special/pow",
-        "@stdlib/math/base/special/exp"
+        "@stdlib/math/base/special/exp",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -101,7 +103,8 @@
         "@stdlib/constants/float64/eulergamma",
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/math/base/special/pow",
-        "@stdlib/math/base/special/exp"
+        "@stdlib/math/base/special/exp",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c b/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c
index 0fa02d54f960..cb67023dbf68 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c
@@ -42,6 +42,7 @@
 #include "stdlib/constants/float64/pi.h"
 #include "stdlib/constants/float64/eulergamma.h"
 #include "stdlib/constants/float64/sqrt_two_pi.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/pow.h"
 #include "stdlib/math/base/special/exp.h"
 #include <stdint.h>
@@ -187,7 +188,7 @@ double stdlib_base_gamma( const double x ) {
 	double z;
 
 	if ( ( stdlib_base_is_integer( x ) && x < 0 ) || x == STDLIB_CONSTANT_FLOAT64_NINF || stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		if ( stdlib_base_is_negative_zero( x ) ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma1pm1/manifest.json b/lib/node_modules/@stdlib/math/base/special/gamma1pm1/manifest.json
index 578c9d461e59..96dbc5c32258 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma1pm1/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/gamma1pm1/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/math/base/special/expm1",
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/math/base/special/expm1",
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/math/base/special/expm1",
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma1pm1/src/main.c b/lib/node_modules/@stdlib/math/base/special/gamma1pm1/src/main.c
index d582ac3ec03c..4a8c7ad04390 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma1pm1/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/gamma1pm1/src/main.c
@@ -39,6 +39,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/constants/float64/eps.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double Y1 = 0.158963680267333984375;
 static const double Y2 = 0.52815341949462890625;
@@ -282,7 +283,7 @@ static double lgammaSmallImp( const double z, const double zm1, const double zm2
 */
 double stdlib_base_gamma1pm1( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x < 0.0 ) {
 		if ( x < -0.5 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/gammainc/manifest.json b/lib/node_modules/@stdlib/math/base/special/gammainc/manifest.json
index d38ecfe1028a..ffcd69f19940 100644
--- a/lib/node_modules/@stdlib/math/base/special/gammainc/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/gammainc/manifest.json
@@ -69,7 +69,8 @@
         "@stdlib/constants/float64/e",
         "@stdlib/constants/float64/gamma-lanczos-g",
         "@stdlib/constants/float32/smallest-normal",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -111,7 +112,8 @@
         "@stdlib/constants/float64/e",
         "@stdlib/constants/float64/gamma-lanczos-g",
         "@stdlib/constants/float32/smallest-normal",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -153,7 +155,8 @@
         "@stdlib/constants/float64/e",
         "@stdlib/constants/float64/gamma-lanczos-g",
         "@stdlib/constants/float32/smallest-normal",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/gammainc/src/main.c b/lib/node_modules/@stdlib/math/base/special/gammainc/src/main.c
index b2173e5d1e2e..b61792e91d7b 100644
--- a/lib/node_modules/@stdlib/math/base/special/gammainc/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/gammainc/src/main.c
@@ -62,6 +62,7 @@
 #include "stdlib/constants/float64/max.h"
 #include "stdlib/constants/float64/pi.h"
 #include "stdlib/constants/float64/e.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdbool.h>
 #include <stdint.h>
 
@@ -1017,7 +1018,7 @@ double stdlib_base_gammainc( const double x, const double a, const bool regulari
 	bool invert;
 
 	if ( x < 0.0 || a <= 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	invert = upper;
 	if ( a >= MAX_NTH_FACTORIAL && !regularized ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/gammaincinv/manifest.json b/lib/node_modules/@stdlib/math/base/special/gammaincinv/manifest.json
index ac2aae108470..3040765fc7bb 100644
--- a/lib/node_modules/@stdlib/math/base/special/gammaincinv/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/gammaincinv/manifest.json
@@ -57,7 +57,8 @@
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/constants/float32/max",
         "@stdlib/constants/float64/two-pi",
-        "@stdlib/math/base/special/gammainc"
+        "@stdlib/math/base/special/gammainc",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -88,7 +89,8 @@
         "@stdlib/constants/float32/max",
         "@stdlib/constants/float64/two-pi",
         "@stdlib/constants/float64/eps",
-        "@stdlib/math/base/special/gammainc"
+        "@stdlib/math/base/special/gammainc",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -119,7 +121,8 @@
         "@stdlib/constants/float32/max",
         "@stdlib/constants/float64/two-pi",
         "@stdlib/constants/float64/eps",
-        "@stdlib/math/base/special/gammainc"
+        "@stdlib/math/base/special/gammainc",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/gammaincinv/src/main.c b/lib/node_modules/@stdlib/math/base/special/gammaincinv/src/main.c
index 2ef048ff8fbe..7e7ae30dc90e 100644
--- a/lib/node_modules/@stdlib/math/base/special/gammaincinv/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/gammaincinv/src/main.c
@@ -51,6 +51,7 @@
 #include "stdlib/constants/float64/two_pi.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float32/max.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdbool.h>
 #include <stdint.h>
 
@@ -1046,7 +1047,7 @@ static double compute( const double a, const double p, const double q ) {
 			invfp = 1.0 / fp;
 		} else {
 			// Overflow problems in one or more steps of the computation...
-			return 0.0 / 0.0; // NaN
+			return STDLIB_CONSTANT_FLOAT64_NAN;
 		}
 	}
 	if ( k < 2 ) {
@@ -1082,13 +1083,13 @@ static double compute( const double a, const double p, const double q ) {
 */
 double stdlib_base_gammaincinv( const double p, const double a, const bool upper ) {
 	if ( stdlib_base_is_nan( p ) || stdlib_base_is_nan( a ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( a < SMALLEST_NORMAL_FLOAT32 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( p > 1.0 || p < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Case: invert upper gamma function
 	if ( upper ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/gcd/manifest.json b/lib/node_modules/@stdlib/math/base/special/gcd/manifest.json
index 87c65f0cf256..631687723ab4 100644
--- a/lib/node_modules/@stdlib/math/base/special/gcd/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/gcd/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/special/fmod",
         "@stdlib/math/base/assert/is-integer",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/special/fmod",
         "@stdlib/math/base/assert/is-integer",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/special/fmod",
         "@stdlib/math/base/assert/is-integer",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/gcd/src/main.c b/lib/node_modules/@stdlib/math/base/special/gcd/src/main.c
index 9167bb5de935..39da1eb6161a 100644
--- a/lib/node_modules/@stdlib/math/base/special/gcd/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/gcd/src/main.c
@@ -22,6 +22,7 @@
 #include "stdlib/math/base/assert/is_integer.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 // 2^63 - 1
@@ -155,7 +156,7 @@ double stdlib_base_gcd( const double a, const double b ) {
 	double bc;
 
 	if ( stdlib_base_is_nan( a ) || stdlib_base_is_nan( b ) ) {
-		return 0.0/0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if (
 		a == STDLIB_CONSTANT_FLOAT64_PINF ||
@@ -163,10 +164,10 @@ double stdlib_base_gcd( const double a, const double b ) {
 		a == STDLIB_CONSTANT_FLOAT64_NINF ||
 		b == STDLIB_CONSTANT_FLOAT64_NINF
 	) {
-		return 0.0/0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( !( stdlib_base_is_integer( a ) && stdlib_base_is_integer( b ) ) ) {
-		return 0.0/0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	ac = a;
 	bc = b;
diff --git a/lib/node_modules/@stdlib/math/base/special/heaviside/manifest.json b/lib/node_modules/@stdlib/math/base/special/heaviside/manifest.json
index 4fcf122c4cbc..a33c9e735070 100644
--- a/lib/node_modules/@stdlib/math/base/special/heaviside/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/heaviside/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/napi/argv-int32",
         "@stdlib/napi/argv-double",
         "@stdlib/napi/create-double",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -69,7 +71,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/heaviside/src/main.c b/lib/node_modules/@stdlib/math/base/special/heaviside/src/main.c
index 4827d43530f3..6d0e9944e044 100644
--- a/lib/node_modules/@stdlib/math/base/special/heaviside/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/heaviside/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/heaviside.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
 /**
@@ -32,7 +33,7 @@
 */
 double stdlib_base_heaviside( const double x, const enum STDLIB_BASE_HEAVISIDE_CONTINUITY continuity ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x > 0.0 ) {
 		return 1.0;
@@ -47,9 +48,9 @@ double stdlib_base_heaviside( const double x, const enum STDLIB_BASE_HEAVISIDE_C
 			case STDLIB_BASE_HEAVISIDE_CONTINUITY_RIGHT_CONTINUOUS:
 				return 1.0;
 			case STDLIB_BASE_HEAVISIDE_CONTINUITY_DISCONTINUOUS:
-				return 0.0 / 0.0; // NaN
+				return STDLIB_CONSTANT_FLOAT64_NAN;
 			default:
-				return 0.0 / 0.0; // NaN
+				return STDLIB_CONSTANT_FLOAT64_NAN;
 		}
 	}
 	return 0.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/hypot/manifest.json b/lib/node_modules/@stdlib/math/base/special/hypot/manifest.json
index 755db550a9d9..bbc29643f6e7 100644
--- a/lib/node_modules/@stdlib/math/base/special/hypot/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/hypot/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/sqrt",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/hypot/src/main.c b/lib/node_modules/@stdlib/math/base/special/hypot/src/main.c
index f4c40fd94752..15c3aee873e2 100644
--- a/lib/node_modules/@stdlib/math/base/special/hypot/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/hypot/src/main.c
@@ -21,6 +21,7 @@
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/sqrt.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the hypotenuse avoiding overflow and underflow.
@@ -43,7 +44,7 @@ double stdlib_base_hypot( const double x, const double y ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
 	}
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	a = x;
 	b = y;
diff --git a/lib/node_modules/@stdlib/math/base/special/kernel-betainc/manifest.json b/lib/node_modules/@stdlib/math/base/special/kernel-betainc/manifest.json
index e9e6f5bb5914..72aa52cd35ce 100644
--- a/lib/node_modules/@stdlib/math/base/special/kernel-betainc/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/kernel-betainc/manifest.json
@@ -74,7 +74,8 @@
         "@stdlib/constants/float64/pi",
         "@stdlib/constants/float64/gamma-lanczos-g",
         "@stdlib/constants/float64/eps",
-        "@stdlib/constants/int32/max"
+        "@stdlib/constants/int32/max",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -121,7 +122,8 @@
         "@stdlib/constants/float64/pi",
         "@stdlib/constants/float64/gamma-lanczos-g",
         "@stdlib/constants/float64/eps",
-        "@stdlib/constants/int32/max"
+        "@stdlib/constants/int32/max",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -169,7 +171,8 @@
         "@stdlib/constants/float64/pi",
         "@stdlib/constants/float64/gamma-lanczos-g",
         "@stdlib/constants/float64/eps",
-        "@stdlib/constants/int32/max"
+        "@stdlib/constants/int32/max",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/kernel-betainc/src/main.c b/lib/node_modules/@stdlib/math/base/special/kernel-betainc/src/main.c
index 7731175949bd..9361c84f010c 100644
--- a/lib/node_modules/@stdlib/math/base/special/kernel-betainc/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/kernel-betainc/src/main.c
@@ -42,6 +42,7 @@
 #include "stdlib/constants/float32/smallest_normal.h"
 #include "stdlib/constants/float64/smallest_normal.h"
 #include "stdlib/constants/int32/max.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/math/base/special/asin.h"
@@ -351,7 +352,7 @@ static double ibetaPowerTerms( const double a, const double b, const double x, c
 			if ( l <= STDLIB_CONSTANT_FLOAT64_MIN_LN || l >= STDLIB_CONSTANT_FLOAT64_MAX_LN ) {
 				l += stdlib_base_ln( result );
 				if ( l >= STDLIB_CONSTANT_FLOAT64_MAX_LN ) {
-					return 0.0 / 0.0;
+					return STDLIB_CONSTANT_FLOAT64_NAN;
 				}
 				result = stdlib_base_exp( l );
 			} else {
@@ -363,7 +364,7 @@ static double ibetaPowerTerms( const double a, const double b, const double x, c
 			if ( l <= STDLIB_CONSTANT_FLOAT64_MIN_LN || l >= STDLIB_CONSTANT_FLOAT64_MAX_LN ) {
 				l += stdlib_base_ln( result );
 				if ( l >= STDLIB_CONSTANT_FLOAT64_MAX_LN ) {
-					return 0.0 / 0.0;
+					return STDLIB_CONSTANT_FLOAT64_NAN;
 				}
 				result = stdlib_base_exp( l );
 			} else {
@@ -386,7 +387,7 @@ static double ibetaPowerTerms( const double a, const double b, const double x, c
 				} else {
 					l2 += l1 + stdlib_base_ln( result );
 					if ( l2 >= STDLIB_CONSTANT_FLOAT64_MAX_LN ) {
-						return 0.0 / 0.0;
+						return STDLIB_CONSTANT_FLOAT64_NAN;
 					}
 					result = stdlib_base_exp( l2 );
 				}
@@ -398,7 +399,7 @@ static double ibetaPowerTerms( const double a, const double b, const double x, c
 				} else {
 					l2 += l1 + stdlib_base_ln( result );
 					if ( l2 >= STDLIB_CONSTANT_FLOAT64_MAX_LN ) {
-						return 0.0 / 0.0;
+						return STDLIB_CONSTANT_FLOAT64_NAN;
 					}
 					result = stdlib_base_exp( l2 );
 				}
@@ -789,21 +790,21 @@ void stdlib_base_kernel_betainc( double x, double a, double b, const bool regula
 	// Derivative not set...
 	*derivative = -1;
 	if ( stdlib_base_is_nan( x ) || x < 0.0 || x > 1.0 ) {
-		*out = 0.0 / 0.0;
-		*derivative = 0.0 / 0.0;
+		*out = STDLIB_CONSTANT_FLOAT64_NAN;
+		*derivative = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 	if ( regularized ) {
 		if ( a < 0.0 || b < 0.0 ) {
-			*out = 0.0 / 0.0;
-			*derivative = 0.0 / 0.0;
+			*out = STDLIB_CONSTANT_FLOAT64_NAN;
+			*derivative = STDLIB_CONSTANT_FLOAT64_NAN;
 			return;
 		}
 		// Extend to a few very special cases...
 		if ( a == 0.0 ) {
 			if ( b == 0.0 ) {
-				*out = 0.0 / 0.0;
-				*derivative = 0.0 / 0.0;
+				*out = STDLIB_CONSTANT_FLOAT64_NAN;
+				*derivative = STDLIB_CONSTANT_FLOAT64_NAN;
 				return;
 			}
 			if ( b > 0.0 ) {
@@ -817,8 +818,8 @@ void stdlib_base_kernel_betainc( double x, double a, double b, const bool regula
 			}
 		}
 	} else if ( a <= 0.0 || b <= 0.0 ) {
-		*out = 0.0 / 0.0;
-		*derivative = 0.0 / 0.0;
+		*out = STDLIB_CONSTANT_FLOAT64_NAN;
+		*derivative = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 	if ( x == 0.0 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/kronecker-delta/manifest.json b/lib/node_modules/@stdlib/math/base/special/kronecker-delta/manifest.json
index 229e194f54da..c0c31c9ecebf 100644
--- a/lib/node_modules/@stdlib/math/base/special/kronecker-delta/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/kronecker-delta/manifest.json
@@ -39,7 +39,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/napi/binary",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       ],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       ],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/kronecker-delta/src/main.c b/lib/node_modules/@stdlib/math/base/special/kronecker-delta/src/main.c
index 058383c53839..5e3a5d71a1dd 100644
--- a/lib/node_modules/@stdlib/math/base/special/kronecker-delta/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/kronecker-delta/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/kronecker_delta.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
 /**
@@ -36,7 +37,7 @@
 */
 double stdlib_base_kronecker_delta( const double i, const double j ) {
 	if ( stdlib_base_is_nan( i ) || stdlib_base_is_nan( j ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( i == j ) {
 		return 1.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/ln/manifest.json b/lib/node_modules/@stdlib/math/base/special/ln/manifest.json
index dcc9a2c1a87c..dcce5827f12c 100644
--- a/lib/node_modules/@stdlib/math/base/special/ln/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ln/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -62,7 +63,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -81,7 +83,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -100,7 +103,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ln/src/main.c b/lib/node_modules/@stdlib/math/base/special/ln/src/main.c
index b6a79245ca29..4e7160503c14 100644
--- a/lib/node_modules/@stdlib/math/base/special/ln/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ln/src/main.c
@@ -35,6 +35,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/constants/float64/ninf.h"
 #include "stdlib/constants/float64/exponent_bias.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double LN2_HI = 6.93147180369123816490e-01; // 3FE62E42 FEE00000
@@ -127,7 +128,7 @@ double stdlib_base_ln( const double x ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
 	}
 	if ( stdlib_base_is_nan( x ) || x < 0.0 ) {
-		return 0.0/0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	stdlib_base_float64_get_high_word( x, &hx );
 	xc = x;
diff --git a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json
index 216b8572261b..389a6f904dc3 100644
--- a/lib/node_modules/@stdlib/math/base/special/log10/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/log10/manifest.json
@@ -46,7 +46,8 @@
         "@stdlib/constants/float64/high-word-significand-mask",
         "@stdlib/number/float64/base/set-low-word",
         "@stdlib/math/base/special/kernel-log1p",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -69,7 +70,8 @@
         "@stdlib/constants/float64/high-word-significand-mask",
         "@stdlib/number/float64/base/set-low-word",
         "@stdlib/math/base/special/kernel-log1p",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -92,7 +94,8 @@
         "@stdlib/constants/float64/high-word-significand-mask",
         "@stdlib/number/float64/base/set-low-word",
         "@stdlib/math/base/special/kernel-log1p",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c
index 90eab78cc4b3..8cb0267618be 100644
--- a/lib/node_modules/@stdlib/math/base/special/log10/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/log10/src/main.c
@@ -39,6 +39,7 @@
 #include "stdlib/constants/float64/high_word_significand_mask.h"
 #include "stdlib/constants/float64/exponent_bias.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/kernel_log1p.h"
 #include <stdint.h>
@@ -87,7 +88,7 @@ double stdlib_base_log10( const double x ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || x < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	xc = x;
 	stdlib_base_float64_to_words( xc, &hx, &lx );
diff --git a/lib/node_modules/@stdlib/math/base/special/log1mexp/manifest.json b/lib/node_modules/@stdlib/math/base/special/log1mexp/manifest.json
index b3cdf711a48d..2ca746256090 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1mexp/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/log1mexp/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/ln-two",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/ln-two",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/ln-two",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/log1mexp/src/main.c b/lib/node_modules/@stdlib/math/base/special/log1mexp/src/main.c
index 73630ddeb972..92fc023557d0 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1mexp/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/log1mexp/src/main.c
@@ -25,6 +25,7 @@
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/constants/float64/ninf.h"
 #include "stdlib/constants/float64/ln_two.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the natural logarithm of \\( 1-\exp(-|x|) \\).
@@ -39,7 +40,7 @@
 double stdlib_base_log1mexp( const double x ) {
 	double ax;
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/log1p/manifest.json b/lib/node_modules/@stdlib/math/base/special/log1p/manifest.json
index 56c3a1ee7c44..00c00499553f 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1p/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/log1p/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/number/float64/base/set-high-word",
         "@stdlib/constants/float64/ninf",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/exponent-bias"
+        "@stdlib/constants/float64/exponent-bias",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/log1p/src/main.c b/lib/node_modules/@stdlib/math/base/special/log1p/src/main.c
index 52229e3256d9..628bba9e30c3 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1p/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/log1p/src/main.c
@@ -37,6 +37,7 @@
 #include "stdlib/constants/float64/ninf.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/exponent_bias.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 // High and low words of ln(2):
@@ -246,7 +247,7 @@ double stdlib_base_log1p( const double x ) {
 	double u;
 
 	if ( x < -1.0 || stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == -1.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/log1pexp/manifest.json b/lib/node_modules/@stdlib/math/base/special/log1pexp/manifest.json
index f02c16b65121..8c5d0c43ea83 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1pexp/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/log1pexp/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/exp",
-        "@stdlib/math/base/special/log1p"
+        "@stdlib/math/base/special/log1p",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/exp",
-        "@stdlib/math/base/special/log1p"
+        "@stdlib/math/base/special/log1p",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/exp",
-        "@stdlib/math/base/special/log1p"
+        "@stdlib/math/base/special/log1p",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/log1pexp/src/main.c b/lib/node_modules/@stdlib/math/base/special/log1pexp/src/main.c
index f28a49d1b48b..c2ace388ee8d 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1pexp/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/log1pexp/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/log1pexp.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/log1p.h"
 #include "stdlib/math/base/special/exp.h"
@@ -33,7 +34,7 @@
 */
 double stdlib_base_log1pexp( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x <= -37.0 ) {
 		return stdlib_base_exp( x );
diff --git a/lib/node_modules/@stdlib/math/base/special/log1pmx/manifest.json b/lib/node_modules/@stdlib/math/base/special/log1pmx/manifest.json
index f3b4915499df..9349edb6f63f 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1pmx/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/log1pmx/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/special/log1p",
         "@stdlib/math/base/special/ln",
         "@stdlib/math/base/special/abs",
-        "@stdlib/constants/float64/eps"
+        "@stdlib/constants/float64/eps",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/log1pmx/src/main.c b/lib/node_modules/@stdlib/math/base/special/log1pmx/src/main.c
index 19402e19967d..7e2c3bcc4f01 100644
--- a/lib/node_modules/@stdlib/math/base/special/log1pmx/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/log1pmx/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/constants/float64/eps.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Evaluates \\( \operatorname{log1pmx}(x) = \ln(1+x) - x \\).
@@ -48,7 +49,7 @@
 double stdlib_base_log1pmx( const double x ) {
 	double ax;
 	if ( x <= -1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	ax = stdlib_base_abs( x );
 	if ( ax > 0.95 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/log2/manifest.json b/lib/node_modules/@stdlib/math/base/special/log2/manifest.json
index 9a4bec85ced9..080336fadb79 100644
--- a/lib/node_modules/@stdlib/math/base/special/log2/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/log2/manifest.json
@@ -46,7 +46,8 @@
         "@stdlib/constants/float64/high-word-significand-mask",
         "@stdlib/constants/float64/exponent-bias",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/math/base/special/kernel-log1p"
+        "@stdlib/math/base/special/kernel-log1p",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -69,7 +70,8 @@
         "@stdlib/constants/float64/high-word-significand-mask",
         "@stdlib/constants/float64/exponent-bias",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/math/base/special/kernel-log1p"
+        "@stdlib/math/base/special/kernel-log1p",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -92,7 +94,8 @@
         "@stdlib/constants/float64/high-word-significand-mask",
         "@stdlib/constants/float64/exponent-bias",
         "@stdlib/constants/float64/ninf",
-        "@stdlib/math/base/special/kernel-log1p"
+        "@stdlib/math/base/special/kernel-log1p",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/log2/src/main.c b/lib/node_modules/@stdlib/math/base/special/log2/src/main.c
index e5d0bc22f126..94a3a18f264c 100644
--- a/lib/node_modules/@stdlib/math/base/special/log2/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/log2/src/main.c
@@ -26,6 +26,7 @@
 #include "stdlib/constants/float64/high_word_significand_mask.h"
 #include "stdlib/constants/float64/exponent_bias.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/kernel_log1p.h"
 #include <stdint.h>
 
@@ -66,7 +67,7 @@ double stdlib_base_log2( const double x ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || x < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	xc = x;
 	stdlib_base_float64_to_words( xc, &hx, &lx );
diff --git a/lib/node_modules/@stdlib/math/base/special/logaddexp/manifest.json b/lib/node_modules/@stdlib/math/base/special/logaddexp/manifest.json
index af940f4cc361..721390e55803 100644
--- a/lib/node_modules/@stdlib/math/base/special/logaddexp/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/logaddexp/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/exp",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/constants/float64/ln-two"
+        "@stdlib/constants/float64/ln-two",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/logaddexp/src/main.c b/lib/node_modules/@stdlib/math/base/special/logaddexp/src/main.c
index 26f4c74634f0..95e969884262 100644
--- a/lib/node_modules/@stdlib/math/base/special/logaddexp/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/logaddexp/src/main.c
@@ -21,6 +21,7 @@
 #include "stdlib/math/base/special/exp.h"
 #include "stdlib/math/base/special/log1p.h"
 #include "stdlib/constants/float64/ln_two.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the natural logarithm of `exp(x) + exp(y)`.
@@ -36,7 +37,7 @@
 double stdlib_base_logaddexp( const double x, const double y ) {
 	double d;
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		return 0.0/0.0;
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == y ) {
 		return x + STDLIB_CONSTANT_FLOAT64_LN2;
diff --git a/lib/node_modules/@stdlib/math/base/special/logit/manifest.json b/lib/node_modules/@stdlib/math/base/special/logit/manifest.json
index cf06f99f7f90..5ec1766c77c7 100644
--- a/lib/node_modules/@stdlib/math/base/special/logit/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/logit/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-probability",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-probability",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-probability",
         "@stdlib/constants/float64/pinf",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/logit/src/main.c b/lib/node_modules/@stdlib/math/base/special/logit/src/main.c
index 5f61a89c3a4f..965f4ae00384 100644
--- a/lib/node_modules/@stdlib/math/base/special/logit/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/logit/src/main.c
@@ -22,6 +22,7 @@
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pinf.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Evaluates the logit function.
@@ -35,10 +36,10 @@
 */
 double stdlib_base_logit( const double p ) {
 	if ( stdlib_base_is_nan( p ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( !stdlib_base_is_probability( p ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( p == 0.0 ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/lucas/manifest.json b/lib/node_modules/@stdlib/math/base/special/lucas/manifest.json
index a3c987779807..1586292f3c0a 100644
--- a/lib/node_modules/@stdlib/math/base/special/lucas/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/lucas/manifest.json
@@ -38,7 +38,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/constants/float64/max-safe-nth-lucas"
+        "@stdlib/constants/float64/max-safe-nth-lucas",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -53,7 +54,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/constants/float64/max-safe-nth-lucas"
+        "@stdlib/constants/float64/max-safe-nth-lucas",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +70,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nonnegative-integer",
-        "@stdlib/constants/float64/max-safe-nth-lucas"
+        "@stdlib/constants/float64/max-safe-nth-lucas",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/lucas/src/main.c b/lib/node_modules/@stdlib/math/base/special/lucas/src/main.c
index 6547f3ed652b..d0db07be42e2 100644
--- a/lib/node_modules/@stdlib/math/base/special/lucas/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/lucas/src/main.c
@@ -19,6 +19,7 @@
 #include "stdlib/math/base/special/lucas.h"
 #include "stdlib/math/base/assert/is_nonnegative_integer.h"
 #include "stdlib/constants/float64/max_safe_nth_lucas.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 #include <stdlib.h>
 
@@ -118,7 +119,7 @@ static const int64_t lucas_value[ 77 ] = {
 */
 double stdlib_base_lucas( const double n ) {
 	if ( !stdlib_base_is_nonnegative_integer( n ) || n > STDLIB_CONSTANT_FLOAT64_MAX_SAFE_NTH_LUCAS ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return lucas_value[ (size_t)n ];
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/max/manifest.json b/lib/node_modules/@stdlib/math/base/special/max/manifest.json
index 97a89b1debb4..2268662da987 100644
--- a/lib/node_modules/@stdlib/math/base/special/max/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/max/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-positive-zero",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-positive-zero",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-positive-zero",
-        "@stdlib/constants/float64/pinf"
+        "@stdlib/constants/float64/pinf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/max/src/main.c b/lib/node_modules/@stdlib/math/base/special/max/src/main.c
index 4170028e3031..3eea3c2a9492 100644
--- a/lib/node_modules/@stdlib/math/base/special/max/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/max/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_positive_zero.h"
 #include "stdlib/constants/float64/pinf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Returns the maximum value.
@@ -38,7 +39,7 @@
 */
 double stdlib_base_max( const double x, const double y ) {
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == STDLIB_CONSTANT_FLOAT64_PINF || y == STDLIB_CONSTANT_FLOAT64_PINF ) {
 		return STDLIB_CONSTANT_FLOAT64_PINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/min/manifest.json b/lib/node_modules/@stdlib/math/base/special/min/manifest.json
index 324b6bdf5b9e..fd591e64230d 100644
--- a/lib/node_modules/@stdlib/math/base/special/min/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/min/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-negative-zero",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-negative-zero",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-negative-zero",
-        "@stdlib/constants/float64/ninf"
+        "@stdlib/constants/float64/ninf",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/min/src/main.c b/lib/node_modules/@stdlib/math/base/special/min/src/main.c
index 392f7118dce9..4e3801208bab 100644
--- a/lib/node_modules/@stdlib/math/base/special/min/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/min/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_negative_zero.h"
 #include "stdlib/constants/float64/ninf.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Return the minimum value.
@@ -38,7 +39,7 @@
 */
 double stdlib_base_min( const double x, const double y ) {
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == STDLIB_CONSTANT_FLOAT64_NINF || y == STDLIB_CONSTANT_FLOAT64_NINF ) {
 		return STDLIB_CONSTANT_FLOAT64_NINF;
diff --git a/lib/node_modules/@stdlib/math/base/special/minmax/manifest.json b/lib/node_modules/@stdlib/math/base/special/minmax/manifest.json
index 502beab58f74..2636631021c9 100644
--- a/lib/node_modules/@stdlib/math/base/special/minmax/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/minmax/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/napi/argv-float64array",
         "@stdlib/napi/export",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -56,7 +57,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/minmax/src/main.c b/lib/node_modules/@stdlib/math/base/special/minmax/src/main.c
index d710b31d4433..3c59e30e6911 100644
--- a/lib/node_modules/@stdlib/math/base/special/minmax/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/minmax/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/minmax.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_negative_zero.h"
 
@@ -38,8 +39,8 @@
 */
 void stdlib_base_minmax( const double x, const double y, double* min, double* max ) {
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		*min = 0.0 / 0.0; // NaN
-		*max = 0.0 / 0.0; // NaN
+		*min = STDLIB_CONSTANT_FLOAT64_NAN;
+		*max = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 	if ( x == y && x == 0.0 ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/minmaxabs/manifest.json b/lib/node_modules/@stdlib/math/base/special/minmaxabs/manifest.json
index 34d444f96a5c..d9ffce8425c4 100644
--- a/lib/node_modules/@stdlib/math/base/special/minmaxabs/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/minmaxabs/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/napi/argv-float64array",
         "@stdlib/napi/export",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/abs"
+        "@stdlib/math/base/special/abs",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -56,7 +57,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/abs"
+        "@stdlib/math/base/special/abs",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/abs"
+        "@stdlib/math/base/special/abs",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/minmaxabs/src/main.c b/lib/node_modules/@stdlib/math/base/special/minmaxabs/src/main.c
index d915de16dad9..d6505689ddd2 100644
--- a/lib/node_modules/@stdlib/math/base/special/minmaxabs/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/minmaxabs/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/minmaxabs.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
@@ -41,8 +42,8 @@ void stdlib_base_minmaxabs( const double x, const double y, double* min, double*
 	double ay;
 
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		*min = 0.0 / 0.0; // NaN
-		*max = 0.0 / 0.0; // NaN
+		*min = STDLIB_CONSTANT_FLOAT64_NAN;
+		*max = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 	ax = stdlib_base_abs( x );
diff --git a/lib/node_modules/@stdlib/math/base/special/modf/manifest.json b/lib/node_modules/@stdlib/math/base/special/modf/manifest.json
index 13505de9a57a..ae9d854fab6e 100644
--- a/lib/node_modules/@stdlib/math/base/special/modf/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/modf/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/constants/float64/exponent-bias",
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/constants/float64/high-word-significand-mask",
-        "@stdlib/constants/float64/num-high-word-significand-bits"
+        "@stdlib/constants/float64/num-high-word-significand-bits",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +69,8 @@
         "@stdlib/constants/float64/exponent-bias",
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/constants/float64/high-word-significand-mask",
-        "@stdlib/constants/float64/num-high-word-significand-bits"
+        "@stdlib/constants/float64/num-high-word-significand-bits",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -89,7 +91,8 @@
         "@stdlib/constants/float64/exponent-bias",
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/constants/float64/high-word-significand-mask",
-        "@stdlib/constants/float64/num-high-word-significand-bits"
+        "@stdlib/constants/float64/num-high-word-significand-bits",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/modf/src/main.c b/lib/node_modules/@stdlib/math/base/special/modf/src/main.c
index 887504904eb3..55bc775e2780 100644
--- a/lib/node_modules/@stdlib/math/base/special/modf/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/modf/src/main.c
@@ -25,6 +25,7 @@
 #include "stdlib/constants/float64/high_word_exponent_mask.h"
 #include "stdlib/constants/float64/high_word_significand_mask.h"
 #include "stdlib/constants/float64/num_high_word_significand_bits.h"
+#include "stdlib/constants/float64/nan.h"
 
 // 4294967295 => 0xffffffff => 11111111111111111111111111111111
 static const uint32_t ALL_ONES = 4294967295;
@@ -68,8 +69,8 @@ void stdlib_base_modf( const double x, double* integral, double* frac ) {
 		return;
 	}
 	if ( stdlib_base_is_nan( x ) ) {
-		*integral = 0.0 / 0.0; // NaN
-		*frac = 0.0 / 0.0; // NaN
+		*integral = STDLIB_CONSTANT_FLOAT64_NAN;
+		*frac = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 	if ( x == STDLIB_CONSTANT_FLOAT64_PINF ) {
diff --git a/lib/node_modules/@stdlib/math/base/special/nanmax/manifest.json b/lib/node_modules/@stdlib/math/base/special/nanmax/manifest.json
index 3fdb96f82040..d4b238c916ac 100644
--- a/lib/node_modules/@stdlib/math/base/special/nanmax/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/nanmax/manifest.json
@@ -38,7 +38,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/max"
+        "@stdlib/math/base/special/max",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -53,7 +54,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/max"
+        "@stdlib/math/base/special/max",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +70,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/max"
+        "@stdlib/math/base/special/max",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/nanmax/src/main.c b/lib/node_modules/@stdlib/math/base/special/nanmax/src/main.c
index d6df587063ed..8fd0ef236f84 100644
--- a/lib/node_modules/@stdlib/math/base/special/nanmax/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/nanmax/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/nanmax.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/max.h"
 
@@ -37,7 +38,7 @@
 */
 double stdlib_base_nanmax( const double x, const double y ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return stdlib_base_is_nan( y ) ? 0.0 / 0.0 : y;
+		return stdlib_base_is_nan( y ) ? STDLIB_CONSTANT_FLOAT64_NAN : y;
 	}
 	return stdlib_base_is_nan( y ) ? x : stdlib_base_max( x, y );
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/nanmin/manifest.json b/lib/node_modules/@stdlib/math/base/special/nanmin/manifest.json
index ff2b442ff555..cb3d93b9e709 100644
--- a/lib/node_modules/@stdlib/math/base/special/nanmin/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/nanmin/manifest.json
@@ -38,7 +38,8 @@
       "dependencies": [
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/min"
+        "@stdlib/math/base/special/min",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -53,7 +54,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/min"
+        "@stdlib/math/base/special/min",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +70,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/special/min"
+        "@stdlib/math/base/special/min",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/nanmin/src/main.c b/lib/node_modules/@stdlib/math/base/special/nanmin/src/main.c
index c125e56af7c7..e12ff9bd723c 100644
--- a/lib/node_modules/@stdlib/math/base/special/nanmin/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/nanmin/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/nanmin.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/min.h"
 
@@ -37,7 +38,7 @@
 */
 double stdlib_base_nanmin( const double x, const double y ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return stdlib_base_is_nan( y ) ? 0.0 / 0.0 : y;
+		return stdlib_base_is_nan( y ) ? STDLIB_CONSTANT_FLOAT64_NAN : y;
 	}
 	return stdlib_base_is_nan( y ) ? x : stdlib_base_min( x, y );
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/negafibonacci/manifest.json b/lib/node_modules/@stdlib/math/base/special/negafibonacci/manifest.json
index 7bf1d197acf5..7780854a71ad 100644
--- a/lib/node_modules/@stdlib/math/base/special/negafibonacci/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/negafibonacci/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/constants/float64/max-safe-nth-fibonacci",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/constants/float64/max-safe-nth-fibonacci",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/constants/float64/max-safe-nth-fibonacci",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/negafibonacci/src/main.c b/lib/node_modules/@stdlib/math/base/special/negafibonacci/src/main.c
index c64339966b34..d67fe22402ea 100644
--- a/lib/node_modules/@stdlib/math/base/special/negafibonacci/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/negafibonacci/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_integer.h"
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/constants/float64/max_safe_nth_fibonacci.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdlib.h>
 
 static const double negafibonacci_value[ 79 ] = {
@@ -117,11 +118,11 @@ static const double negafibonacci_value[ 79 ] = {
 double stdlib_base_negafibonacci( const double n ) {
 	double an;
 	if ( !stdlib_base_is_integer( n ) || n > 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	an = stdlib_base_abs( n );
 	if ( an > STDLIB_CONSTANT_FLOAT64_MAX_SAFE_NTH_FIBONACCI ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return negafibonacci_value[ (size_t)an ];
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/negalucas/manifest.json b/lib/node_modules/@stdlib/math/base/special/negalucas/manifest.json
index 3476ff893e32..8a5c3ecd13ca 100644
--- a/lib/node_modules/@stdlib/math/base/special/negalucas/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/negalucas/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/max-safe-nth-lucas",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/max-safe-nth-lucas",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/special/abs",
         "@stdlib/constants/float64/max-safe-nth-lucas",
-        "@stdlib/math/base/assert/is-integer"
+        "@stdlib/math/base/assert/is-integer",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/negalucas/src/main.c b/lib/node_modules/@stdlib/math/base/special/negalucas/src/main.c
index 7099a731b6fe..4b659f2f7527 100644
--- a/lib/node_modules/@stdlib/math/base/special/negalucas/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/negalucas/src/main.c
@@ -20,6 +20,7 @@
 #include "stdlib/math/base/assert/is_integer.h"
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/constants/float64/max_safe_nth_lucas.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 #include <stdlib.h>
 
@@ -120,7 +121,7 @@ static const int64_t negalucas_value[ 77 ] = {
 double stdlib_base_negalucas( const double n ) {
 	double an = stdlib_base_abs( n );
 	if ( !stdlib_base_is_integer( n ) || n > 0.0 || an > STDLIB_CONSTANT_FLOAT64_MAX_SAFE_NTH_LUCAS ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return negalucas_value[ (size_t)an ];
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/pdiff/manifest.json b/lib/node_modules/@stdlib/math/base/special/pdiff/manifest.json
index 80656ac17472..a8cc72f1d15c 100644
--- a/lib/node_modules/@stdlib/math/base/special/pdiff/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/pdiff/manifest.json
@@ -37,7 +37,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/napi/binary",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -51,7 +52,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +67,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/pdiff/src/main.c b/lib/node_modules/@stdlib/math/base/special/pdiff/src/main.c
index e25a50eea24b..d3b6fe3b271a 100644
--- a/lib/node_modules/@stdlib/math/base/special/pdiff/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/pdiff/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/pdiff.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
 /**
@@ -36,7 +37,7 @@
 */
 double stdlib_base_pdiff( const double x, const double y ) {
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_nan( y ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x > y ) {
 		return x - y;
diff --git a/lib/node_modules/@stdlib/math/base/special/powm1/manifest.json b/lib/node_modules/@stdlib/math/base/special/powm1/manifest.json
index 4ad8fead0be1..6fc74537356f 100644
--- a/lib/node_modules/@stdlib/math/base/special/powm1/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/powm1/manifest.json
@@ -44,8 +44,8 @@
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/special/fmod"
-
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -66,7 +66,8 @@
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -87,7 +88,8 @@
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/powm1/src/main.c b/lib/node_modules/@stdlib/math/base/special/powm1/src/main.c
index eeb8ded80052..dc05d6bb306f 100644
--- a/lib/node_modules/@stdlib/math/base/special/powm1/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/powm1/src/main.c
@@ -30,6 +30,7 @@
 */
 
 #include "stdlib/math/base/special/powm1.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/abs.h"
 #include "stdlib/math/base/special/expm1.h"
 #include "stdlib/math/base/special/ln.h"
@@ -60,7 +61,7 @@ double stdlib_base_powm1( const double b, const double x ) {
 	double y;
 
 	if ( stdlib_base_is_nan( b ) || stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 0.0 ) {
 		// Any number raised to zero (including 0) is always 1 => b^0 - 1 = 0
@@ -85,11 +86,11 @@ double stdlib_base_powm1( const double b, const double x ) {
 		}
 	} else if ( stdlib_base_trunc( x ) != x ) {
 		// Exponentiation would yield a complex result...
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	result = stdlib_base_pow( bc, x ) - 1.0;
 	if ( stdlib_base_is_infinite( result ) || stdlib_base_is_nan( result ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return result;
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/ramp/manifest.json b/lib/node_modules/@stdlib/math/base/special/ramp/manifest.json
index ca59ad02f8c0..0f7e320c5fba 100644
--- a/lib/node_modules/@stdlib/math/base/special/ramp/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ramp/manifest.json
@@ -37,7 +37,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/math/base/napi/unary",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -51,7 +52,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +67,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/ramp/src/main.c b/lib/node_modules/@stdlib/math/base/special/ramp/src/main.c
index c757c3c4296d..5692eb1ad749 100644
--- a/lib/node_modules/@stdlib/math/base/special/ramp/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/ramp/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/ramp.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 
 /**
@@ -35,7 +36,7 @@
 */
 double stdlib_base_ramp( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x > 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/rempio2/manifest.json b/lib/node_modules/@stdlib/math/base/special/rempio2/manifest.json
index 69ce750a6ce5..e7b56622452f 100644
--- a/lib/node_modules/@stdlib/math/base/special/rempio2/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/rempio2/manifest.json
@@ -49,7 +49,8 @@
         "@stdlib/number/float64/base/from-words",
         "@stdlib/constants/float64/high-word-abs-mask",
         "@stdlib/constants/float64/high-word-exponent-mask",
-        "@stdlib/constants/float64/high-word-significand-mask"
+        "@stdlib/constants/float64/high-word-significand-mask",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +72,8 @@
         "@stdlib/number/float64/base/from-words",
         "@stdlib/constants/float64/high-word-abs-mask",
         "@stdlib/constants/float64/high-word-exponent-mask",
-        "@stdlib/constants/float64/high-word-significand-mask"
+        "@stdlib/constants/float64/high-word-significand-mask",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -93,7 +95,8 @@
         "@stdlib/number/float64/base/from-words",
         "@stdlib/constants/float64/high-word-abs-mask",
         "@stdlib/constants/float64/high-word-exponent-mask",
-        "@stdlib/constants/float64/high-word-significand-mask"
+        "@stdlib/constants/float64/high-word-significand-mask",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/rempio2/src/main.c b/lib/node_modules/@stdlib/math/base/special/rempio2/src/main.c
index 0384913710a3..f0a1d33b557e 100644
--- a/lib/node_modules/@stdlib/math/base/special/rempio2/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/rempio2/src/main.c
@@ -40,6 +40,7 @@
 #include "stdlib/constants/float64/high_word_abs_mask.h"
 #include "stdlib/constants/float64/high_word_exponent_mask.h"
 #include "stdlib/constants/float64/high_word_significand_mask.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double ZERO =  0.00000000000000000000e+00;       // 0x00000000, 0x00000000
@@ -595,8 +596,8 @@ int32_t stdlib_base_rempio2( const double x, double *rem1, double *rem2 ) {
 	}
 	// Case: x is NaN or infinity
 	if( ix >= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_EXPONENT_MASK ) {
-		*rem1 = 0.0 / 0.0; // NaN
-		*rem2 = 0.0 / 0.0; // NaN
+		*rem1 = STDLIB_CONSTANT_FLOAT64_NAN;
+		*rem2 = STDLIB_CONSTANT_FLOAT64_NAN;
 		return 0;
 	}
 	// Set z = scalbn(|x|, ilogb(x)-23)...
diff --git a/lib/node_modules/@stdlib/math/base/special/rempio2f/manifest.json b/lib/node_modules/@stdlib/math/base/special/rempio2f/manifest.json
index c9316fb43303..92bcea4acd07 100644
--- a/lib/node_modules/@stdlib/math/base/special/rempio2f/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/rempio2f/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/number/float32/base/from-word",
         "@stdlib/number/float32/base/to-word",
         "@stdlib/constants/float32/abs-mask",
-        "@stdlib/constants/float32/exponent-mask"
+        "@stdlib/constants/float32/exponent-mask",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -67,7 +68,8 @@
         "@stdlib/number/float32/base/from-word",
         "@stdlib/number/float32/base/to-word",
         "@stdlib/constants/float32/abs-mask",
-        "@stdlib/constants/float32/exponent-mask"
+        "@stdlib/constants/float32/exponent-mask",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -87,7 +89,8 @@
         "@stdlib/number/float32/base/from-word",
         "@stdlib/number/float32/base/to-word",
         "@stdlib/constants/float32/abs-mask",
-        "@stdlib/constants/float32/exponent-mask"
+        "@stdlib/constants/float32/exponent-mask",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/rempio2f/src/main.c b/lib/node_modules/@stdlib/math/base/special/rempio2f/src/main.c
index 6b6cd893dab5..604f9a0f5d8b 100644
--- a/lib/node_modules/@stdlib/math/base/special/rempio2f/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/rempio2f/src/main.c
@@ -38,6 +38,7 @@
 #include "stdlib/number/float32/base/to_word.h"
 #include "stdlib/constants/float32/abs_mask.h"
 #include "stdlib/constants/float32/exponent_mask.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double ZERO =  0.00000000000000000000e+00;       // 0x00000000, 0x00000000
@@ -421,7 +422,7 @@ int32_t stdlib_base_rempio2f( const float x, double *rem ) {
 	}
 	// Case: x is NaN or infinity
 	if( ix >= STDLIB_CONSTANT_FLOAT32_EXPONENT_MASK ) {
-		*rem = 0.0 / 0.0; // NaN
+		*rem = STDLIB_CONSTANT_FLOAT64_NAN;
 		return 0;
 	}
 	// Set z = scalbn(|x|, ilogb(|x|)-23)...
diff --git a/lib/node_modules/@stdlib/math/base/special/riemann-zeta/manifest.json b/lib/node_modules/@stdlib/math/base/special/riemann-zeta/manifest.json
index e9f8bb305808..bd59301693cc 100644
--- a/lib/node_modules/@stdlib/math/base/special/riemann-zeta/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/riemann-zeta/manifest.json
@@ -53,7 +53,8 @@
         "@stdlib/constants/float64/two-pi",
         "@stdlib/constants/float64/sqrt-eps",
         "@stdlib/constants/float64/ln-sqrt-two-pi",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -83,7 +84,8 @@
         "@stdlib/constants/float64/two-pi",
         "@stdlib/constants/float64/sqrt-eps",
         "@stdlib/constants/float64/ln-sqrt-two-pi",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -113,7 +115,8 @@
         "@stdlib/constants/float64/two-pi",
         "@stdlib/constants/float64/sqrt-eps",
         "@stdlib/constants/float64/ln-sqrt-two-pi",
-        "@stdlib/constants/float64/max-nth-factorial"
+        "@stdlib/constants/float64/max-nth-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/riemann-zeta/src/main.c b/lib/node_modules/@stdlib/math/base/special/riemann-zeta/src/main.c
index 4da9c0c153f1..2b03eaf1e9df 100644
--- a/lib/node_modules/@stdlib/math/base/special/riemann-zeta/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/riemann-zeta/src/main.c
@@ -47,6 +47,7 @@
 #include "stdlib/constants/float64/sqrt_eps.h"
 #include "stdlib/constants/float64/ln_sqrt_two_pi.h"
 #include "stdlib/constants/float64/max_nth_factorial.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const int32_t MAX_BERNOULLI_2N = 129;
@@ -708,12 +709,12 @@ double stdlib_base_zeta( const double s ) {
 
 	// Check for `NaN`:
 	if ( stdlib_base_is_nan( s ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	// Check for a pole:
 	if ( s == 1.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	// Check for large value:
diff --git a/lib/node_modules/@stdlib/math/base/special/rising-factorial/manifest.json b/lib/node_modules/@stdlib/math/base/special/rising-factorial/manifest.json
index 4565d1b7156a..88105a60308e 100644
--- a/lib/node_modules/@stdlib/math/base/special/rising-factorial/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/rising-factorial/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/gamma-delta-ratio",
-        "@stdlib/math/base/special/falling-factorial"
+        "@stdlib/math/base/special/falling-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/gamma-delta-ratio",
-        "@stdlib/math/base/special/falling-factorial"
+        "@stdlib/math/base/special/falling-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/gamma-delta-ratio",
-        "@stdlib/math/base/special/falling-factorial"
+        "@stdlib/math/base/special/falling-factorial",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/rising-factorial/src/main.c b/lib/node_modules/@stdlib/math/base/special/rising-factorial/src/main.c
index fccd379e2b33..89ce0fb35061 100644
--- a/lib/node_modules/@stdlib/math/base/special/rising-factorial/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/rising-factorial/src/main.c
@@ -31,6 +31,7 @@
 
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/gamma_delta_ratio.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/falling_factorial.h"
 #include <stdint.h>
 #include <stdbool.h>
@@ -67,7 +68,7 @@ double stdlib_base_rising_factorial( const double x, const int32_t n ) {
 	bool inv;
 
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	nc = n;
 	xc = x;
diff --git a/lib/node_modules/@stdlib/math/base/special/round/manifest.json b/lib/node_modules/@stdlib/math/base/special/round/manifest.json
index d52c3d1ec2b0..dca1a08f300c 100644
--- a/lib/node_modules/@stdlib/math/base/special/round/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/round/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -58,7 +59,8 @@
       "dependencies": [
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -75,7 +77,8 @@
       "dependencies": [
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -92,7 +95,8 @@
       "dependencies": [
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/round/src/main.c b/lib/node_modules/@stdlib/math/base/special/round/src/main.c
index 9137b8be571c..55978a243c59 100644
--- a/lib/node_modules/@stdlib/math/base/special/round/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/round/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/round.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/floor.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_negative_zero.h"
@@ -33,7 +34,7 @@
 */
 double stdlib_base_round( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0/0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_negative_zero( x ) || ( x >= -0.5 && x < 0.0 ) ) {
 		return -0.0; // -0
diff --git a/lib/node_modules/@stdlib/math/base/special/roundb/manifest.json b/lib/node_modules/@stdlib/math/base/special/roundb/manifest.json
index 261497985ece..8984397bd316 100644
--- a/lib/node_modules/@stdlib/math/base/special/roundb/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/roundb/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/round",
         "@stdlib/math/base/special/roundn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/round",
         "@stdlib/math/base/special/roundn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/round",
         "@stdlib/math/base/special/roundn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/roundb/src/main.c b/lib/node_modules/@stdlib/math/base/special/roundb/src/main.c
index e36234c1c3fd..5d184ca487a3 100644
--- a/lib/node_modules/@stdlib/math/base/special/roundb/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/roundb/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/roundb.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/round.h"
@@ -41,7 +42,7 @@ double stdlib_base_roundb( const double x, const int32_t n, const int32_t b ) {
 	double s;
 
 	if ( stdlib_base_is_nan( x ) || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/roundn/manifest.json b/lib/node_modules/@stdlib/math/base/special/roundn/manifest.json
index 68abf0bebf73..94ed0dd9eee7 100644
--- a/lib/node_modules/@stdlib/math/base/special/roundn/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/roundn/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/constants/float64/min-base10-exponent",
         "@stdlib/constants/float64/min-base10-exponent-subnormal",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +72,8 @@
         "@stdlib/constants/float64/min-base10-exponent",
         "@stdlib/constants/float64/min-base10-exponent-subnormal",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -95,7 +97,8 @@
         "@stdlib/constants/float64/min-base10-exponent",
         "@stdlib/constants/float64/min-base10-exponent-subnormal",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/roundn/src/main.c b/lib/node_modules/@stdlib/math/base/special/roundn/src/main.c
index 4b1ff1255a2e..ae55806f8631 100644
--- a/lib/node_modules/@stdlib/math/base/special/roundn/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/roundn/src/main.c
@@ -26,6 +26,7 @@
 #include "stdlib/constants/float64/max_base10_exponent.h"
 #include "stdlib/constants/float64/min_base10_exponent.h"
 #include "stdlib/constants/float64/min_base10_exponent_subnormal.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double MAX_INT = STDLIB_CONSTANT_FLOAT64_MAX_SAFE_INTEGER + 1.0;
@@ -117,7 +118,7 @@ double stdlib_base_roundn( const double x, const int32_t n ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	if (
diff --git a/lib/node_modules/@stdlib/math/base/special/roundsd/manifest.json b/lib/node_modules/@stdlib/math/base/special/roundsd/manifest.json
index c87297416f66..fd1aefae39b4 100644
--- a/lib/node_modules/@stdlib/math/base/special/roundsd/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/roundsd/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/number/float64/base/exponent",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/round",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +72,8 @@
         "@stdlib/number/float64/base/exponent",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/round",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -95,7 +97,8 @@
         "@stdlib/number/float64/base/exponent",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/round",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/roundsd/src/main.c b/lib/node_modules/@stdlib/math/base/special/roundsd/src/main.c
index e707479166be..d65c1f61e4d0 100644
--- a/lib/node_modules/@stdlib/math/base/special/roundsd/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/roundsd/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/roundsd.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/pow.h"
 #include "stdlib/math/base/assert/is_infinite.h"
@@ -46,7 +47,7 @@ double stdlib_base_roundsd( const double x, const int32_t n, const int32_t b ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || n < 1 ) {
-		return 0.0/0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/sin/manifest.json b/lib/node_modules/@stdlib/math/base/special/sin/manifest.json
index a5286e487f73..4fef571f3983 100644
--- a/lib/node_modules/@stdlib/math/base/special/sin/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sin/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/kernel-sin",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/kernel-sin",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/kernel-sin",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sin/src/main.c b/lib/node_modules/@stdlib/math/base/special/sin/src/main.c
index 252f144057a7..e0b1507dec04 100644
--- a/lib/node_modules/@stdlib/math/base/special/sin/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sin/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/number/float64/base/get_high_word.h"
 #include "stdlib/constants/float64/high_word_abs_mask.h"
 #include "stdlib/constants/float64/high_word_exponent_mask.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/kernel_cos.h"
 #include "stdlib/math/base/special/kernel_sin.h"
 #include "stdlib/math/base/special/rempio2.h"
@@ -90,7 +91,7 @@ double stdlib_base_sin( const double x ) {
 	}
 	// Case: x is NaN or infinity
 	if ( ix >= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_EXPONENT_MASK ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Argument reduction...
 	n = stdlib_base_rempio2( x, &Y[ 0 ], &Y[ 1 ] );
diff --git a/lib/node_modules/@stdlib/math/base/special/sinc/manifest.json b/lib/node_modules/@stdlib/math/base/special/sinc/manifest.json
index 0b6c07c7669e..90bec0abd21e 100644
--- a/lib/node_modules/@stdlib/math/base/special/sinc/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sinc/manifest.json
@@ -42,7 +42,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/constants/float64/pi",
-        "@stdlib/math/base/special/sinpi"
+        "@stdlib/math/base/special/sinpi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -61,7 +62,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/constants/float64/pi",
-        "@stdlib/math/base/special/sinpi"
+        "@stdlib/math/base/special/sinpi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -80,7 +82,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/constants/float64/pi",
-        "@stdlib/math/base/special/sinpi"
+        "@stdlib/math/base/special/sinpi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sinc/src/main.c b/lib/node_modules/@stdlib/math/base/special/sinc/src/main.c
index b1ecc5dcfdf9..1564586fe24c 100644
--- a/lib/node_modules/@stdlib/math/base/special/sinc/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sinc/src/main.c
@@ -21,6 +21,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/constants/float64/pi.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the normalized cardinal sine of a number.
@@ -53,7 +54,7 @@
 */
 double stdlib_base_sinc( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) ) {
 		return 0.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/sincos/manifest.json b/lib/node_modules/@stdlib/math/base/special/sincos/manifest.json
index 45918c065001..e3ea1d4c0be3 100644
--- a/lib/node_modules/@stdlib/math/base/special/sincos/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sincos/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/math/base/special/rempio2",
         "@stdlib/number/float64/base/get-high-word",
         "@stdlib/constants/float64/high-word-exponent-mask",
-        "@stdlib/constants/float64/high-word-abs-mask"
+        "@stdlib/constants/float64/high-word-abs-mask",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -60,7 +61,8 @@
         "@stdlib/math/base/special/rempio2",
         "@stdlib/number/float64/base/get-high-word",
         "@stdlib/constants/float64/high-word-exponent-mask",
-        "@stdlib/constants/float64/high-word-abs-mask"
+        "@stdlib/constants/float64/high-word-abs-mask",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/special/rempio2",
         "@stdlib/number/float64/base/get-high-word",
         "@stdlib/constants/float64/high-word-exponent-mask",
-        "@stdlib/constants/float64/high-word-abs-mask"
+        "@stdlib/constants/float64/high-word-abs-mask",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sincos/src/main.c b/lib/node_modules/@stdlib/math/base/special/sincos/src/main.c
index 3f38a29043dc..f8df71456e3b 100644
--- a/lib/node_modules/@stdlib/math/base/special/sincos/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sincos/src/main.c
@@ -19,6 +19,7 @@
 #include "stdlib/math/base/special/sincos.h"
 #include "stdlib/constants/float64/high_word_abs_mask.h"
 #include "stdlib/constants/float64/high_word_exponent_mask.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/number/float64/base/get_high_word.h"
 #include "stdlib/math/base/special/rempio2.h"
 #include <stdint.h>
@@ -131,8 +132,8 @@ void stdlib_base_sincos( const double x, double* sine, double* cosine ) {
 
 	// Case: x is NaN or infinity
 	if ( ix >= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_EXPONENT_MASK ) {
-		*sine = 0.0 / 0.0; // NaN
-		*cosine = 0.0 / 0.0; // NaN
+		*sine = STDLIB_CONSTANT_FLOAT64_NAN;
+		*cosine = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 
diff --git a/lib/node_modules/@stdlib/math/base/special/sincospi/manifest.json b/lib/node_modules/@stdlib/math/base/special/sincospi/manifest.json
index 04df8f3d8103..f858caf61b9a 100644
--- a/lib/node_modules/@stdlib/math/base/special/sincospi/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sincospi/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/math/base/special/fmod",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -68,7 +69,8 @@
         "@stdlib/math/base/special/fmod",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -89,7 +91,8 @@
         "@stdlib/math/base/special/fmod",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sincospi/src/main.c b/lib/node_modules/@stdlib/math/base/special/sincospi/src/main.c
index fd8be3d48632..7f0fec5ca892 100644
--- a/lib/node_modules/@stdlib/math/base/special/sincospi/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sincospi/src/main.c
@@ -25,6 +25,7 @@
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/constants/float64/pi.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Simultaneously computes the sine and cosine of a number times π.
@@ -47,8 +48,8 @@ void stdlib_base_sincospi( const double x, double* sine, double* cosine ) {
 	double r;
 
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_infinite( x ) ) {
-		*sine = 0.0 / 0.0; // NaN
-		*cosine = 0.0 / 0.0; // NaN
+		*sine = STDLIB_CONSTANT_FLOAT64_NAN;
+		*cosine = STDLIB_CONSTANT_FLOAT64_NAN;
 		return;
 	}
 	r = stdlib_base_fmod( x, 2.0 );
diff --git a/lib/node_modules/@stdlib/math/base/special/sind/manifest.json b/lib/node_modules/@stdlib/math/base/special/sind/manifest.json
index 01aba340dce1..2cfd06d9d874 100644
--- a/lib/node_modules/@stdlib/math/base/special/sind/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sind/manifest.json
@@ -44,7 +44,8 @@
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/signum",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -65,7 +66,8 @@
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/signum",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -86,7 +88,8 @@
         "@stdlib/math/base/special/kernel-cos",
         "@stdlib/math/base/special/signum",
         "@stdlib/math/base/special/abs",
-        "@stdlib/math/base/special/fmod"
+        "@stdlib/math/base/special/fmod",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sind/src/main.c b/lib/node_modules/@stdlib/math/base/special/sind/src/main.c
index d7be9229ac5a..39f62c1c8da9 100644
--- a/lib/node_modules/@stdlib/math/base/special/sind/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sind/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/sind.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/kernel_sin.h"
 #include "stdlib/math/base/special/kernel_cos.h"
 #include "stdlib/math/base/special/deg2rad.h"
@@ -41,7 +42,7 @@ double stdlib_base_sind( const double x ) {
 	double rx;
 
 	if ( stdlib_base_is_infinite( x ) || stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	rx = stdlib_base_fmod( x, 360.0 );
diff --git a/lib/node_modules/@stdlib/math/base/special/sinpi/manifest.json b/lib/node_modules/@stdlib/math/base/special/sinpi/manifest.json
index af9eaa3ea92d..7ffeb9195092 100644
--- a/lib/node_modules/@stdlib/math/base/special/sinpi/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sinpi/manifest.json
@@ -46,7 +46,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/copysign",
         "@stdlib/math/base/special/fmod",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -69,7 +70,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/copysign",
         "@stdlib/math/base/special/fmod",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -92,7 +94,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/copysign",
         "@stdlib/math/base/special/fmod",
-        "@stdlib/constants/float64/pi"
+        "@stdlib/constants/float64/pi",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sinpi/src/main.c b/lib/node_modules/@stdlib/math/base/special/sinpi/src/main.c
index 4267f2b1722e..71a04bbd7cf6 100644
--- a/lib/node_modules/@stdlib/math/base/special/sinpi/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sinpi/src/main.c
@@ -25,6 +25,7 @@
 #include "stdlib/math/base/special/copysign.h"
 #include "stdlib/math/base/special/fmod.h"
 #include "stdlib/constants/float64/pi.h"
+#include "stdlib/constants/float64/nan.h"
 
 /**
 * Computes the value of `sin(πx)`.
@@ -47,7 +48,7 @@ double stdlib_base_sinpi( const double x ) {
 	double ar;
 	double r;
 	if ( stdlib_base_is_nan( x ) || stdlib_base_is_infinite( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	// Argument reduction (reduce to [0,2))...
diff --git a/lib/node_modules/@stdlib/math/base/special/spence/manifest.json b/lib/node_modules/@stdlib/math/base/special/spence/manifest.json
index 502ee8833533..b12eb457cca7 100644
--- a/lib/node_modules/@stdlib/math/base/special/spence/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/spence/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/pi-squared"
+        "@stdlib/constants/float64/pi-squared",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/pi-squared"
+        "@stdlib/constants/float64/pi-squared",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/ln",
-        "@stdlib/constants/float64/pi-squared"
+        "@stdlib/constants/float64/pi-squared",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/spence/src/main.c b/lib/node_modules/@stdlib/math/base/special/spence/src/main.c
index 7b30a150ff28..fe1de1556b40 100644
--- a/lib/node_modules/@stdlib/math/base/special/spence/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/spence/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/ln.h"
 #include "stdlib/constants/float64/pi_squared.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 // π^2 / 6
@@ -114,7 +115,7 @@ double stdlib_base_spence( const double x ) {
 	double z;
 
 	if ( stdlib_base_is_nan( x ) || x < 0.0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( x == 1.0 ) {
 		return 0.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/manifest.json b/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/manifest.json
index ea66b59bf9ad..841491d585b9 100644
--- a/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/special/abs"
+        "@stdlib/math/base/special/abs",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/special/abs"
+        "@stdlib/math/base/special/abs",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/sqrt",
         "@stdlib/math/base/special/log1p",
-        "@stdlib/math/base/special/abs"
+        "@stdlib/math/base/special/abs",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/src/main.c b/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/src/main.c
index ee91d95abb8d..3eb8174bcbb5 100644
--- a/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/sqrt1pm1/src/main.c
@@ -30,6 +30,7 @@
 */
 
 #include "stdlib/math/base/special/sqrt1pm1.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/expm1.h"
 #include "stdlib/math/base/special/log1p.h"
@@ -48,7 +49,7 @@
 */
 double stdlib_base_sqrt1pm1( const double x ) {
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_abs( x ) > 0.75 ) {
 		return stdlib_base_sqrt( 1.0 + x ) - 1.0;
diff --git a/lib/node_modules/@stdlib/math/base/special/tan/manifest.json b/lib/node_modules/@stdlib/math/base/special/tan/manifest.json
index 047d2921f19c..035939f0712e 100644
--- a/lib/node_modules/@stdlib/math/base/special/tan/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/tan/manifest.json
@@ -43,7 +43,8 @@
         "@stdlib/constants/float64/high-word-abs-mask",
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-tan",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -63,7 +64,8 @@
         "@stdlib/constants/float64/high-word-abs-mask",
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-tan",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -83,7 +85,8 @@
         "@stdlib/constants/float64/high-word-abs-mask",
         "@stdlib/constants/float64/high-word-exponent-mask",
         "@stdlib/math/base/special/kernel-tan",
-        "@stdlib/math/base/special/rempio2"
+        "@stdlib/math/base/special/rempio2",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/tan/src/main.c b/lib/node_modules/@stdlib/math/base/special/tan/src/main.c
index 2851c48042ea..e8b8475d06de 100644
--- a/lib/node_modules/@stdlib/math/base/special/tan/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/tan/src/main.c
@@ -34,6 +34,7 @@
 #include "stdlib/number/float64/base/get_high_word.h"
 #include "stdlib/constants/float64/high_word_abs_mask.h"
 #include "stdlib/constants/float64/high_word_exponent_mask.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/kernel_tan.h"
 #include "stdlib/math/base/special/rempio2.h"
 #include <stdint.h>
@@ -89,7 +90,7 @@ double stdlib_base_tan( const double x ) {
 	}
 	// Case: tan(Inf or NaN) is NaN
 	if ( ix >= STDLIB_CONSTANT_FLOAT64_HIGH_WORD_EXPONENT_MASK ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	// Argument reduction...
 	n = stdlib_base_rempio2( x, &Y[ 0 ], &Y[ 1 ] );
diff --git a/lib/node_modules/@stdlib/math/base/special/tribonacci/manifest.json b/lib/node_modules/@stdlib/math/base/special/tribonacci/manifest.json
index f4cbee083620..a28af94964e4 100644
--- a/lib/node_modules/@stdlib/math/base/special/tribonacci/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/tribonacci/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/constants/float64/max-safe-nth-tribonacci",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-nonnegative-integer"
+        "@stdlib/math/base/assert/is-nonnegative-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/constants/float64/max-safe-nth-tribonacci",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-nonnegative-integer"
+        "@stdlib/math/base/assert/is-nonnegative-integer",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/constants/float64/max-safe-nth-tribonacci",
         "@stdlib/math/base/assert/is-nan",
-        "@stdlib/math/base/assert/is-nonnegative-integer"
+        "@stdlib/math/base/assert/is-nonnegative-integer",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/tribonacci/src/main.c b/lib/node_modules/@stdlib/math/base/special/tribonacci/src/main.c
index e62283fba283..11760080470f 100644
--- a/lib/node_modules/@stdlib/math/base/special/tribonacci/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/tribonacci/src/main.c
@@ -18,6 +18,7 @@
 
 #include "stdlib/math/base/special/tribonacci.h"
 #include "stdlib/constants/float64/max_safe_nth_tribonacci.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nonnegative_integer.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include <stdlib.h>
@@ -102,7 +103,7 @@ static const int64_t tribonacci_value[ 64 ] = {
 */
 double stdlib_base_tribonacci( const double n ) {
 	if ( stdlib_base_is_nan( n ) || !stdlib_base_is_nonnegative_integer( n ) || n > STDLIB_CONSTANT_FLOAT64_MAX_SAFE_NTH_TRIBONACCI ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	return tribonacci_value[ (size_t)n ];
 }
diff --git a/lib/node_modules/@stdlib/math/base/special/trigamma/manifest.json b/lib/node_modules/@stdlib/math/base/special/trigamma/manifest.json
index e0983f2170a0..91b106b80d49 100644
--- a/lib/node_modules/@stdlib/math/base/special/trigamma/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/trigamma/manifest.json
@@ -39,7 +39,8 @@
         "@stdlib/math/base/napi/unary",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/sinpi",
-        "@stdlib/constants/float64/pi-squared"
+        "@stdlib/constants/float64/pi-squared",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -55,7 +56,8 @@
       "dependencies": [
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/sinpi",
-        "@stdlib/constants/float64/pi-squared"
+        "@stdlib/constants/float64/pi-squared",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +73,8 @@
       "dependencies": [
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/sinpi",
-        "@stdlib/constants/float64/pi-squared"
+        "@stdlib/constants/float64/pi-squared",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/trigamma/src/main.c b/lib/node_modules/@stdlib/math/base/special/trigamma/src/main.c
index c3ca521b1ff9..d766c748062b 100644
--- a/lib/node_modules/@stdlib/math/base/special/trigamma/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/trigamma/src/main.c
@@ -33,6 +33,7 @@
 #include "stdlib/math/base/special/floor.h"
 #include "stdlib/math/base/special/sinpi.h"
 #include "stdlib/constants/float64/pi_squared.h"
+#include "stdlib/constants/float64/nan.h"
 
 static const double YOFFSET24 = 3.558437347412109375;
 
@@ -267,7 +268,7 @@ double stdlib_base_trigamma( const double x ) {
 	// Check for negative arguments and use reflection:
 	if ( x <= 0.0 ) {
 		if ( stdlib_base_floor( x ) == x ) {
-			return 0.0 / 0.0; // NaN
+			return STDLIB_CONSTANT_FLOAT64_NAN;
 		}
 		s = stdlib_base_sinpi( x );
 		z = 1.0 - x;
diff --git a/lib/node_modules/@stdlib/math/base/special/truncb/manifest.json b/lib/node_modules/@stdlib/math/base/special/truncb/manifest.json
index 28036b85938b..98eb198e3c7f 100644
--- a/lib/node_modules/@stdlib/math/base/special/truncb/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/truncb/manifest.json
@@ -41,7 +41,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/special/truncn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -59,7 +60,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/special/truncn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -77,7 +79,8 @@
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/special/truncn",
-        "@stdlib/math/base/special/pow"
+        "@stdlib/math/base/special/pow",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/truncb/src/main.c b/lib/node_modules/@stdlib/math/base/special/truncb/src/main.c
index 50b0f4885957..2d03dfef1c7d 100644
--- a/lib/node_modules/@stdlib/math/base/special/truncb/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/truncb/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/truncb.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/trunc.h"
@@ -41,7 +42,7 @@ double stdlib_base_truncb( const double x, const int32_t n, const int32_t b ) {
 	double s;
 
 	if ( stdlib_base_is_nan( x ) || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/truncn/manifest.json b/lib/node_modules/@stdlib/math/base/special/truncn/manifest.json
index 1f85b8ae0824..6824a502456d 100644
--- a/lib/node_modules/@stdlib/math/base/special/truncn/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/truncn/manifest.json
@@ -47,7 +47,8 @@
         "@stdlib/constants/float64/min-base10-exponent",
         "@stdlib/constants/float64/min-base10-exponent-subnormal",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -71,7 +72,8 @@
         "@stdlib/constants/float64/min-base10-exponent",
         "@stdlib/constants/float64/min-base10-exponent-subnormal",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -95,7 +97,8 @@
         "@stdlib/constants/float64/min-base10-exponent",
         "@stdlib/constants/float64/min-base10-exponent-subnormal",
         "@stdlib/math/base/assert/is-infinite",
-        "@stdlib/math/base/assert/is-nan"
+        "@stdlib/math/base/assert/is-nan",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/truncn/src/main.c b/lib/node_modules/@stdlib/math/base/special/truncn/src/main.c
index 8544aeb5cbc0..1ab4b9255b94 100644
--- a/lib/node_modules/@stdlib/math/base/special/truncn/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/truncn/src/main.c
@@ -26,6 +26,7 @@
 #include "stdlib/constants/float64/max_base10_exponent.h"
 #include "stdlib/constants/float64/min_base10_exponent.h"
 #include "stdlib/constants/float64/min_base10_exponent_subnormal.h"
+#include "stdlib/constants/float64/nan.h"
 #include <stdint.h>
 
 static const double MAX_INT = STDLIB_CONSTANT_FLOAT64_MAX_SAFE_INTEGER + 1.0;
@@ -96,7 +97,7 @@ double stdlib_base_truncn( const double x, const int32_t n ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 
 	if (
diff --git a/lib/node_modules/@stdlib/math/base/special/truncsd/manifest.json b/lib/node_modules/@stdlib/math/base/special/truncsd/manifest.json
index 1055163ccd52..474f91441603 100644
--- a/lib/node_modules/@stdlib/math/base/special/truncsd/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/truncsd/manifest.json
@@ -45,7 +45,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/trunc",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -67,7 +68,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/trunc",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -89,7 +91,8 @@
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/floor",
         "@stdlib/math/base/special/trunc",
-        "@stdlib/number/float64/base/exponent"
+        "@stdlib/number/float64/base/exponent",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/truncsd/src/main.c b/lib/node_modules/@stdlib/math/base/special/truncsd/src/main.c
index df49b024606b..805b612d9e6f 100644
--- a/lib/node_modules/@stdlib/math/base/special/truncsd/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/truncsd/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/truncsd.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/math/base/special/pow.h"
@@ -46,7 +47,7 @@ double stdlib_base_truncsd( const double x, const int32_t n, const int32_t b ) {
 	double y;
 
 	if ( stdlib_base_is_nan( x ) || n < 1 || b <= 0 ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	if ( stdlib_base_is_infinite( x ) || x == 0.0 ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/math/base/special/wrap/manifest.json b/lib/node_modules/@stdlib/math/base/special/wrap/manifest.json
index b2f19cbfa4f3..e7521dd7f6ff 100644
--- a/lib/node_modules/@stdlib/math/base/special/wrap/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/wrap/manifest.json
@@ -40,7 +40,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/special/fmod",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -57,7 +58,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/special/fmod",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     },
     {
@@ -74,7 +76,8 @@
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/special/trunc",
         "@stdlib/math/base/special/fmod",
-        "@stdlib/math/base/assert/is-negative-zero"
+        "@stdlib/math/base/assert/is-negative-zero",
+        "@stdlib/constants/float64/nan"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/wrap/src/main.c b/lib/node_modules/@stdlib/math/base/special/wrap/src/main.c
index 490b41f6acf3..223a64e7e5e5 100644
--- a/lib/node_modules/@stdlib/math/base/special/wrap/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/wrap/src/main.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/math/base/special/wrap.h"
+#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/special/trunc.h"
 #include "stdlib/math/base/special/fmod.h"
@@ -41,7 +42,7 @@ double stdlib_base_wrap( const double v, const double min, const double max ) {
 	double vc;
 
 	if ( stdlib_base_is_nan( v ) || stdlib_base_is_nan( min ) || stdlib_base_is_nan( max ) || max <= min ) {
-		return 0.0 / 0.0; // NaN
+		return STDLIB_CONSTANT_FLOAT64_NAN;
 	}
 	maxc = max;
 	minc = min;

From dfb27240cb38a0b6e5432a93baef7c1b92380628 Mon Sep 17 00:00:00 2001
From: Claude <noreply@anthropic.com>
Date: Mon, 20 Jul 2026 13:42:01 +0000
Subject: [PATCH 2/4] chore: flatten `arrays` parameter description in `blas`
 ndarray READMEs

Propagates fix from 2b91bfe ("chore: clean-up") to sibling packages
whose README documents a single input ndarray using the nested list
form.

Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01EHa9r1LqZ3Tg3Gft34EGQv
---
 lib/node_modules/@stdlib/blas/base/ndarray/dnrm2/README.md    | 4 +---
 lib/node_modules/@stdlib/blas/base/ndarray/dzasum/README.md   | 4 +---
 lib/node_modules/@stdlib/blas/base/ndarray/dznrm2/README.md   | 4 +---
 lib/node_modules/@stdlib/blas/base/ndarray/gnrm2/README.md    | 4 +---
 lib/node_modules/@stdlib/blas/base/ndarray/scasum/README.md   | 4 +---
 lib/node_modules/@stdlib/blas/base/ndarray/scnrm2/README.md   | 4 +---
 lib/node_modules/@stdlib/blas/base/ndarray/snrm2/README.md    | 4 +---
 .../@stdlib/blas/ext/base/ndarray/dindex-of-falsy/README.md   | 4 +---
 .../@stdlib/blas/ext/base/ndarray/gindex-of-falsy/README.md   | 4 +---
 9 files changed, 9 insertions(+), 27 deletions(-)

diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dnrm2/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/dnrm2/README.md
index 47e6b4a7b372..6db1e62b0bed 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/dnrm2/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dnrm2/README.md
@@ -66,9 +66,7 @@ var y = dnrm2( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dzasum/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/dzasum/README.md
index 728cec93cf3b..6d011a909e46 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/dzasum/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dzasum/README.md
@@ -66,9 +66,7 @@ var y = dzasum( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/dznrm2/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/dznrm2/README.md
index d0f06536c7bb..6a5f4c4436ff 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/dznrm2/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/dznrm2/README.md
@@ -66,9 +66,7 @@ var y = dznrm2( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/gnrm2/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/gnrm2/README.md
index a95e2b6cd557..393f298a04aa 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/gnrm2/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/gnrm2/README.md
@@ -66,9 +66,7 @@ var y = gnrm2( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/scasum/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/scasum/README.md
index 8b6a7500d6a8..04a1b3f0fe38 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/scasum/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/scasum/README.md
@@ -66,9 +66,7 @@ var y = scasum( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/scnrm2/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/scnrm2/README.md
index 59c4034ad6c4..df8ff5c6d4b4 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/scnrm2/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/scnrm2/README.md
@@ -66,9 +66,7 @@ var y = scnrm2( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/base/ndarray/snrm2/README.md b/lib/node_modules/@stdlib/blas/base/ndarray/snrm2/README.md
index 847e89b13662..f159043c2338 100644
--- a/lib/node_modules/@stdlib/blas/base/ndarray/snrm2/README.md
+++ b/lib/node_modules/@stdlib/blas/base/ndarray/snrm2/README.md
@@ -66,9 +66,7 @@ var y = snrm2( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 </section>
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/ndarray/dindex-of-falsy/README.md b/lib/node_modules/@stdlib/blas/ext/base/ndarray/dindex-of-falsy/README.md
index 6dd9ff3da3bb..87682c4e81d1 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/ndarray/dindex-of-falsy/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/ndarray/dindex-of-falsy/README.md
@@ -51,9 +51,7 @@ var idx = dindexOfFalsy( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 If the function is unable to find a falsy element, the function returns `-1`.
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/ndarray/gindex-of-falsy/README.md b/lib/node_modules/@stdlib/blas/ext/base/ndarray/gindex-of-falsy/README.md
index 5fdffdc55614..657227c61cf3 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/ndarray/gindex-of-falsy/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/ndarray/gindex-of-falsy/README.md
@@ -51,9 +51,7 @@ var idx = gindexOfFalsy( [ x ] );
 
 The function has the following parameters:
 
--   **arrays**: array-like object containing the following ndarrays:
-
-    -   a one-dimensional input ndarray.
+-   **arrays**: array-like object containing a one-dimensional input ndarray.
 
 If the function is unable to find a falsy element, the function returns `-1`.
 

From 504abdd0497e0a67189ac2fae21448c32560a750 Mon Sep 17 00:00:00 2001
From: Claude <noreply@anthropic.com>
Date: Mon, 20 Jul 2026 13:42:02 +0000
Subject: [PATCH 3/4] style: remove stale `eslint-disable max-len` directives
 in READMEs

Propagates fix from 2780757 ("style: re-enable lint rule") to READMEs
whose governed code blocks no longer exceed the 80-character line
limit.

Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01EHa9r1LqZ3Tg3Gft34EGQv
---
 .../@stdlib/blas/ext/base/dfirst-index-less-than/README.md    | 2 --
 lib/node_modules/@stdlib/blas/ext/base/dsome/README.md        | 2 --
 lib/node_modules/@stdlib/blas/ext/base/ssome/README.md        | 2 --
 lib/node_modules/@stdlib/lapack/base/dgttrf/README.md         | 2 --
 lib/node_modules/@stdlib/lapack/base/dlassq/README.md         | 2 --
 lib/node_modules/@stdlib/lapack/base/dpttrf/README.md         | 2 --
 lib/node_modules/@stdlib/lapack/base/spttrf/README.md         | 2 --
 lib/node_modules/@stdlib/ndarray/base/assign-scalar/README.md | 2 --
 .../@stdlib/strided/base/nullary-addon-dispatch/README.md     | 4 ----
 9 files changed, 20 deletions(-)

diff --git a/lib/node_modules/@stdlib/blas/ext/base/dfirst-index-less-than/README.md b/lib/node_modules/@stdlib/blas/ext/base/dfirst-index-less-than/README.md
index 55f07fb9a49b..6f18d1a15e71 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dfirst-index-less-than/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dfirst-index-less-than/README.md
@@ -129,8 +129,6 @@ The function has the following additional parameters:
 
 While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example, to access only the last three elements of each strided array:
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float64Array = require( '@stdlib/array/float64' );
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dsome/README.md b/lib/node_modules/@stdlib/blas/ext/base/dsome/README.md
index cd7f63b8cd87..158478bba25a 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/dsome/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/dsome/README.md
@@ -58,8 +58,6 @@ The function has the following parameters:
 
 The `N` and stride parameters determine which elements in the strided array are accessed at runtime. For example, to test every other element:
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float64Array = require( '@stdlib/array/float64' );
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/ssome/README.md b/lib/node_modules/@stdlib/blas/ext/base/ssome/README.md
index 2b9f9516dcf1..e0211f43166c 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/ssome/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/ssome/README.md
@@ -58,8 +58,6 @@ The function has the following parameters:
 
 The `N` and stride parameters determine which elements in the strided array are accessed at runtime. For example, to test every other element:
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float32Array = require( '@stdlib/array/float32' );
 
diff --git a/lib/node_modules/@stdlib/lapack/base/dgttrf/README.md b/lib/node_modules/@stdlib/lapack/base/dgttrf/README.md
index ec6ea4c7dcc1..36c986015972 100644
--- a/lib/node_modules/@stdlib/lapack/base/dgttrf/README.md
+++ b/lib/node_modules/@stdlib/lapack/base/dgttrf/README.md
@@ -241,8 +241,6 @@ The function has the following additional parameters:
 
 While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example,
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float64Array = require( '@stdlib/array/float64' );
 var Int32Array = require( '@stdlib/array/int32' );
diff --git a/lib/node_modules/@stdlib/lapack/base/dlassq/README.md b/lib/node_modules/@stdlib/lapack/base/dlassq/README.md
index 63826edd883a..426853b6c3d5 100644
--- a/lib/node_modules/@stdlib/lapack/base/dlassq/README.md
+++ b/lib/node_modules/@stdlib/lapack/base/dlassq/README.md
@@ -113,8 +113,6 @@ The function has the following additional parameters:
 
 While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example,
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float64Array = require( '@stdlib/array/float64' );
 
diff --git a/lib/node_modules/@stdlib/lapack/base/dpttrf/README.md b/lib/node_modules/@stdlib/lapack/base/dpttrf/README.md
index 89026cde78f3..86a70e04b82a 100644
--- a/lib/node_modules/@stdlib/lapack/base/dpttrf/README.md
+++ b/lib/node_modules/@stdlib/lapack/base/dpttrf/README.md
@@ -95,8 +95,6 @@ The function has the following additional parameters:
 
 While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example,
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float64Array = require( '@stdlib/array/float64' );
 
diff --git a/lib/node_modules/@stdlib/lapack/base/spttrf/README.md b/lib/node_modules/@stdlib/lapack/base/spttrf/README.md
index b267ad710e3f..47c1e7e2ac11 100644
--- a/lib/node_modules/@stdlib/lapack/base/spttrf/README.md
+++ b/lib/node_modules/@stdlib/lapack/base/spttrf/README.md
@@ -95,8 +95,6 @@ The function has the following additional parameters:
 
 While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example,
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var Float32Array = require( '@stdlib/array/float32' );
 
diff --git a/lib/node_modules/@stdlib/ndarray/base/assign-scalar/README.md b/lib/node_modules/@stdlib/ndarray/base/assign-scalar/README.md
index 1210577c042d..a359fef3100b 100644
--- a/lib/node_modules/@stdlib/ndarray/base/assign-scalar/README.md
+++ b/lib/node_modules/@stdlib/ndarray/base/assign-scalar/README.md
@@ -40,8 +40,6 @@ var assignScalar = require( '@stdlib/ndarray/base/assign-scalar' );
 
 Assigns a scalar value to every element of an output [ndarray][@stdlib/ndarray/base/ctor].
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 var scalar2ndarray = require( '@stdlib/ndarray/from-scalar' );
 var Float64Array = require( '@stdlib/array/float64' );
diff --git a/lib/node_modules/@stdlib/strided/base/nullary-addon-dispatch/README.md b/lib/node_modules/@stdlib/strided/base/nullary-addon-dispatch/README.md
index 86d3d2e7df95..1d3bc36eb546 100644
--- a/lib/node_modules/@stdlib/strided/base/nullary-addon-dispatch/README.md
+++ b/lib/node_modules/@stdlib/strided/base/nullary-addon-dispatch/README.md
@@ -105,8 +105,6 @@ where
 
 Returns a function which dispatches to a native add-on applying a nullary function using alternative indexing semantics.
 
-<!-- eslint-disable max-len -->
-
 ```javascript
 function addon( N, dtypeX, x, strideX ) {
     // Call into native add-on...
@@ -188,8 +186,6 @@ where
 
 ## Examples
 
-<!-- eslint-disable max-len -->
-
 <!-- eslint no-undef: "error" -->
 
 ```javascript

From 09327397b47d4a1c8b2aa7dd8e374ee767758339 Mon Sep 17 00:00:00 2001
From: Claude <noreply@anthropic.com>
Date: Mon, 20 Jul 2026 13:51:38 +0000
Subject: [PATCH 4/4] refactor: exclude `math/base/special/gamma` from NaN
 propagation

Reverts the `constants/float64/nan` changes for
`math/base/special/gamma`: its `src/main.c` carries pre-existing
trailing whitespace (lines 211 and 227) which fails editorconfig
linting for any PR touching the file. Dropped per the
pre-existing-lint rule rather than folding unrelated whitespace
fixes into the propagation.

Co-Authored-By: Claude Fable 5 <noreply@anthropic.com>
Claude-Session: https://claude.ai/code/session_01EHa9r1LqZ3Tg3Gft34EGQv
---
 .../@stdlib/math/base/special/gamma/manifest.json        | 9 +++------
 .../@stdlib/math/base/special/gamma/src/main.c           | 3 +--
 2 files changed, 4 insertions(+), 8 deletions(-)

diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json b/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json
index 95be3c49993e..0b051015bd59 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/manifest.json
@@ -49,8 +49,7 @@
         "@stdlib/constants/float64/eulergamma",
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/math/base/special/pow",
-        "@stdlib/math/base/special/exp",
-        "@stdlib/constants/float64/nan"
+        "@stdlib/math/base/special/exp"
       ]
     },
     {
@@ -76,8 +75,7 @@
         "@stdlib/constants/float64/eulergamma",
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/math/base/special/pow",
-        "@stdlib/math/base/special/exp",
-        "@stdlib/constants/float64/nan"
+        "@stdlib/math/base/special/exp"
       ]
     },
     {
@@ -103,8 +101,7 @@
         "@stdlib/constants/float64/eulergamma",
         "@stdlib/constants/float64/sqrt-two-pi",
         "@stdlib/math/base/special/pow",
-        "@stdlib/math/base/special/exp",
-        "@stdlib/constants/float64/nan"
+        "@stdlib/math/base/special/exp"
       ]
     }
   ]
diff --git a/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c b/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c
index cb67023dbf68..0fa02d54f960 100644
--- a/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c
+++ b/lib/node_modules/@stdlib/math/base/special/gamma/src/main.c
@@ -42,7 +42,6 @@
 #include "stdlib/constants/float64/pi.h"
 #include "stdlib/constants/float64/eulergamma.h"
 #include "stdlib/constants/float64/sqrt_two_pi.h"
-#include "stdlib/constants/float64/nan.h"
 #include "stdlib/math/base/special/pow.h"
 #include "stdlib/math/base/special/exp.h"
 #include <stdint.h>
@@ -188,7 +187,7 @@ double stdlib_base_gamma( const double x ) {
 	double z;
 
 	if ( ( stdlib_base_is_integer( x ) && x < 0 ) || x == STDLIB_CONSTANT_FLOAT64_NINF || stdlib_base_is_nan( x ) ) {
-		return STDLIB_CONSTANT_FLOAT64_NAN;
+		return 0.0 / 0.0; // NaN
 	}
 	if ( x == 0.0 ) {
 		if ( stdlib_base_is_negative_zero( x ) ) {