diff --git a/lib/node_modules/@stdlib/_tools/changelog/parse-commits/scripts/commits.sh b/lib/node_modules/@stdlib/_tools/changelog/parse-commits/scripts/commits.sh
index ec2094057b68..f6424dbeac10 100755
--- a/lib/node_modules/@stdlib/_tools/changelog/parse-commits/scripts/commits.sh
+++ b/lib/node_modules/@stdlib/_tools/changelog/parse-commits/scripts/commits.sh
@@ -35,7 +35,7 @@ if [ -z "$1" ]; then
 fi
 
 if [ -z "$GIT_COMMIT_SEP" ]; then
-    GIT_COMMIT_SEP="^---^"; # Default separator
+	GIT_COMMIT_SEP="^---^"; # Default separator
 fi
 
 git log --name-only --no-merges --notes --pretty=format:"%H|%ad|%aN <%aE>|%B|%N|" "$1" | awk -v sep="$GIT_COMMIT_SEP" '/^$/{p=1;next} /^[0-9a-f]{40}\|/{if (p==1) print sep; p=0} {print}';
diff --git a/lib/node_modules/@stdlib/math/base/special/ceiln/manifest.json b/lib/node_modules/@stdlib/math/base/special/ceiln/manifest.json
index 2a9e59a90ffa..a85af2e9172b 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceiln/manifest.json
+++ b/lib/node_modules/@stdlib/math/base/special/ceiln/manifest.json
@@ -41,6 +41,7 @@
         "@stdlib/math/base/napi/binary",
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/ceil",
+        "@stdlib/math/base/special/pow",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/constants/float64/max-base10-exponent",
@@ -65,6 +66,7 @@
       "dependencies": [
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/ceil",
+        "@stdlib/math/base/special/pow",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/constants/float64/max-base10-exponent",
@@ -89,6 +91,7 @@
       "dependencies": [
         "@stdlib/math/base/special/abs",
         "@stdlib/math/base/special/ceil",
+        "@stdlib/math/base/special/pow",
         "@stdlib/math/base/assert/is-nan",
         "@stdlib/math/base/assert/is-infinite",
         "@stdlib/constants/float64/max-base10-exponent",
diff --git a/lib/node_modules/@stdlib/math/base/special/ceiln/src/ceiln.c b/lib/node_modules/@stdlib/math/base/special/ceiln/src/ceiln.c
index 4e3449324892..02b1d1ed6d99 100644
--- a/lib/node_modules/@stdlib/math/base/special/ceiln/src/ceiln.c
+++ b/lib/node_modules/@stdlib/math/base/special/ceiln/src/ceiln.c
@@ -19,6 +19,7 @@
 #include "stdlib/math/base/special/ceiln.h"
 #include "stdlib/math/base/special/ceil.h"
 #include "stdlib/math/base/special/abs.h"
+#include "stdlib/math/base/special/pow.h"
 #include "stdlib/math/base/assert/is_nan.h"
 #include "stdlib/math/base/assert/is_infinite.h"
 #include "stdlib/constants/float64/max_base10_exponent.h"
@@ -27,7 +28,6 @@
 #include "stdlib/constants/float64/max_safe_integer.h"
 #include "stdlib/constants/float64/pinf.h"
 #include <stdint.h>
-#include <math.h>
 
 
 // VARIABLES //
@@ -90,14 +90,14 @@ double stdlib_base_ceiln( const double x, const int32_t n ) {
 	}
 	// If we overflow, return `x`, as the number of digits to the right of the decimal is too small (i.e., `x` is too large / lacks sufficient fractional precision) for there to be any effect when rounding...
 	if ( n < STDLIB_CONSTANT_FLOAT64_MIN_BASE10_EXPONENT ) {
-		s = pow( 10.0, -( n + STDLIB_CONSTANT_FLOAT64_MAX_BASE10_EXPONENT ) ); // TODO: replace use of `pow` once have stdlib equivalent
+		s = stdlib_base_pow( 10.0, -( n + STDLIB_CONSTANT_FLOAT64_MAX_BASE10_EXPONENT ) );
 		y = ( x * HUGE_VALUE ) * s; // order of operation matters!
 		if ( stdlib_base_is_infinite( y ) ) {
 			return x;
 		}
 		return ( stdlib_base_ceil( y ) / HUGE_VALUE ) / s;
 	}
-	s = pow( 10.0, -n ); // TODO: replace use of `pow` once have stdlib equivalent
+	s = stdlib_base_pow( 10.0, -n );
 	y = x * s;
 	if ( stdlib_base_is_infinite( y ) ) {
 		return x;
diff --git a/lib/node_modules/@stdlib/stats/base/dists/poisson/median/README.md b/lib/node_modules/@stdlib/stats/base/dists/poisson/median/README.md
index b4f9f8f90807..4234b6d995dd 100644
--- a/lib/node_modules/@stdlib/stats/base/dists/poisson/median/README.md
+++ b/lib/node_modules/@stdlib/stats/base/dists/poisson/median/README.md
@@ -31,7 +31,7 @@ The [median][median] for a [Poisson][poisson-distribution] random variable is
 <!-- <equation class="equation" label="eq:poisson_median" align="center" raw="\operatorname{Median}\left( X \right) = \lfloor \lambda+1/3-0.02/\lambda \rfloor " alt="Median for a Poisson distribution."> -->
 
 ```math
-\mathop{\mathrm{Median}}\left( X \right) = \lfloor \lambda+1/3-0.02/\lambda \rfloor 
+\mathop{\mathrm{Median}}\left( X \right) = \lfloor \lambda+1/3-0.02/\lambda \rfloor
 ```
 
 <!-- <div class="equation" align="center" data-raw-text="\operatorname{Median}\left( X \right) = \lfloor \lambda+1/3-0.02/\lambda \rfloor " data-equation="eq:poisson_median">