diff --git a/lib/node_modules/@stdlib/math/special/data/unary.json b/lib/node_modules/@stdlib/math/special/data/unary.json
index 40ed9f32ca27..f2ca8fcfa491 100644
--- a/lib/node_modules/@stdlib/math/special/data/unary.json
+++ b/lib/node_modules/@stdlib/math/special/data/unary.json
@@ -1008,7 +1008,81 @@
     ],
     "extra_keywords": []
   },
-  "@stdlib/math/base/special/gammaln": {},
+  "@stdlib/math/base/special/gammaln": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "gammaln",
+    "alias": "gammaln",
+    "pkg_desc": "evaluate the natural logarithm of the gamma function",
+    "desc": "evaluates the natural logarithm of the gamma function",
+    "short_desc": "natural logarithm of the gamma function",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "double",
+          "dtype": "float64"
+        },
+        "domain": [
+          {
+            "min": "-infinity",
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            -10,
+            10
+          ]
+        },
+        "example_values": [
+          1,
+          3.5,
+          4.5,
+          -0.5,
+          2,
+          3,
+          -3.5,
+          0.1,
+          4,
+          1.5,
+          5,
+          0.5,
+          2.5,
+          -1.5,
+          -2.5,
+          50,
+          100,
+          -50.5,
+          -100.5,
+          2.2
+        ]
+      }
+    ],
+    "output_policy": "real_floating_point_and_generic",
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "gamma",
+      "factorial",
+      "natural",
+      "logarithm",
+      "log",
+      "ln",
+      "lgamma"
+    ],
+    "extra_keywords": []
+  },
   "@stdlib/math/base/special/acosf": {
     "$schema": "math/base@v1.0",
     "base_alias": "acos",
@@ -8271,7 +8345,136 @@
       "math.exp"
     ]
   },
-  "@stdlib/math/base/special/cexp": {},
+  "@stdlib/math/base/special/cexp": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "cexp",
+    "alias": "cexp",
+    "pkg_desc": "evaluate the exponential function for a double-precision complex floating-point number",
+    "desc": "evaluates the exponential function for a double-precision complex floating-point number",
+    "short_desc": "exponential function",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "Complex128",
+          "jsdoc": "Complex128",
+          "c": "stdlib_complex128_t",
+          "dtype": "complex128"
+        },
+        "domain": null,
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            [
+              -10,
+              10
+            ],
+            [
+              -10,
+              10
+            ]
+          ]
+        },
+        "example_values": [
+          {
+            "re": -3.14,
+            "im": -1.5
+          },
+          {
+            "re": 0,
+            "im": 0
+          },
+          {
+            "re": -1.5,
+            "im": 2.5
+          },
+          {
+            "re": 2.5,
+            "im": -1.5
+          },
+          {
+            "re": 0,
+            "im": -3.7
+          },
+          {
+            "re": 4.2,
+            "im": 0
+          },
+          {
+            "re": 21.2,
+            "im": 3
+          },
+          {
+            "re": 11,
+            "im": -5
+          },
+          {
+            "re": 33,
+            "im": -14.67
+          },
+          {
+            "re": -42,
+            "im": 9.3
+          },
+          {
+            "re": -3,
+            "im": 3
+          },
+          {
+            "re": 73,
+            "im": 31
+          },
+          {
+            "re": -2.45,
+            "im": 1.23
+          },
+          {
+            "re": 2.45,
+            "im": -1.23
+          },
+          {
+            "re": 1.77,
+            "im": -3.14
+          },
+          {
+            "re": -7.5,
+            "im": 8.2
+          },
+          {
+            "re": 5.5,
+            "im": -12.3
+          },
+          {
+            "re": -15.8,
+            "im": 0.4
+          },
+          {
+            "re": 0.99,
+            "im": -0.99
+          }
+        ]
+      }
+    ],
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "Complex128",
+        "jsdoc": "Complex128",
+        "c": "stdlib_complex128_t",
+        "dtype": "complex128"
+      }
+    },
+    "keywords": [
+      "complex",
+      "cmplx",
+      "exponential",
+      "exp",
+      "pow",
+      "cexp"
+    ],
+    "extra_keywords": []
+  },
   "@stdlib/math/base/special/exp10": {
     "$schema": "math/base@v1.0",
     "base_alias": "exp10",
@@ -8502,15 +8705,13 @@
       "math.exp"
     ]
   },
-  "@stdlib/math/base/special/expm1": {},
-  "@stdlib/math/base/special/expm1rel": {},
-  "@stdlib/math/base/special/factorial": {
+  "@stdlib/math/base/special/expm1": {
     "$schema": "math/base@v1.0",
-    "base_alias": "factorial",
-    "alias": "factorial",
-    "pkg_desc": "compute the factorial",
-    "desc": "computes the factorial",
-    "short_desc": "factorial",
+    "base_alias": "expm1",
+    "alias": "expm1",
+    "pkg_desc": "compute `exp(x)-1`, where `exp(x)` is the natural exponential function",
+    "desc": "computes `exp(x)-1`, where `exp(x)` is the natural exponential function",
+    "short_desc": "exp(x)-1",
     "parameters": [
       {
         "name": "x",
@@ -8523,44 +8724,43 @@
         },
         "domain": [
           {
-            "min": 0,
+            "min": "-infinity",
             "max": "infinity"
           }
         ],
         "rand": {
-          "prng": "random/base/discrete-uniform",
+          "prng": "random/base/uniform",
           "parameters": [
-            0,
-            100
+            -10,
+            10
           ]
         },
         "example_values": [
-          1,
+          -1.2,
           2,
-          3,
-          4,
-          5,
-          6,
-          7,
-          8,
-          9,
-          10,
-          11,
-          12,
-          13,
-          14,
-          15,
-          16,
-          17,
-          21,
-          34,
-          99
+          -3.1,
+          -4.7,
+          5.5,
+          6.7,
+          8.9,
+          -10.2,
+          11.3,
+          -12.4,
+          13.5,
+          14.6,
+          -15.7,
+          16.8,
+          -17.9,
+          18.1,
+          -19.11,
+          20.12,
+          -21.15,
+          23.78
         ]
       }
     ],
-    "output_policy": "real_floating_point_and_generic",
     "returns": {
-      "desc": "factorial",
+      "desc": "function value",
       "type": {
         "javascript": "number",
         "jsdoc": "number",
@@ -8569,12 +8769,167 @@
       }
     },
     "keywords": [
-      "factorial",
-      "fact"
+      "expm1",
+      "power",
+      "natural",
+      "exp",
+      "exponential",
+      "function",
+      "euler"
     ],
-    "extra_keywords": []
+    "extra_keywords": [
+      "math.expm1"
+    ]
   },
-  "@stdlib/math/base/special/factorial2f": {
+  "@stdlib/math/base/special/expm1rel": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "expm1rel",
+    "alias": "expm1rel",
+    "pkg_desc": "compute the relative error exponential",
+    "desc": "computes the relative error exponential",
+    "short_desc": "relative error exponential",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "double",
+          "dtype": "float64"
+        },
+        "domain": [
+          {
+            "min": "-infinity",
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            -10,
+            10
+          ]
+        },
+        "example_values": [
+          -1.2,
+          2,
+          -3.1,
+          -4.7,
+          5.5,
+          6.7,
+          8.9,
+          -10.2,
+          11.3,
+          -12.4,
+          13.5,
+          14.6,
+          -15.7,
+          16.8,
+          -17.9,
+          18.1,
+          -19.11,
+          20.12,
+          -21.15,
+          23.78
+        ]
+      }
+    ],
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "expm1",
+      "power",
+      "natural",
+      "exp",
+      "relative",
+      "error",
+      "err",
+      "exponential",
+      "euler"
+    ],
+    "extra_keywords": [
+      "math.expm1"
+    ]
+  },
+  "@stdlib/math/base/special/factorial": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "factorial",
+    "alias": "factorial",
+    "pkg_desc": "compute the factorial",
+    "desc": "computes the factorial",
+    "short_desc": "factorial",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "double",
+          "dtype": "float64"
+        },
+        "domain": [
+          {
+            "min": 0,
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/discrete-uniform",
+          "parameters": [
+            0,
+            100
+          ]
+        },
+        "example_values": [
+          1,
+          2,
+          3,
+          4,
+          5,
+          6,
+          7,
+          8,
+          9,
+          10,
+          11,
+          12,
+          13,
+          14,
+          15,
+          16,
+          17,
+          21,
+          34,
+          99
+        ]
+      }
+    ],
+    "output_policy": "real_floating_point_and_generic",
+    "returns": {
+      "desc": "factorial",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "factorial",
+      "fact"
+    ],
+    "extra_keywords": []
+  },
+  "@stdlib/math/base/special/factorial2f": {
     "$schema": "math/base@v1.0",
     "base_alias": "factorial2",
     "alias": "factorial2f",
@@ -9650,10 +10005,305 @@
     ],
     "extra_keywords": []
   },
-  "@stdlib/math/base/special/gamma-lanczos-sum": {},
-  "@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf": {},
-  "@stdlib/math/base/special/gamma-lanczos-sum-expg-scaled": {},
-  "@stdlib/math/base/special/gamma1pm1": {},
+  "@stdlib/math/base/special/gamma-lanczos-sum": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "gammaLanczosSum",
+    "alias": "gammaLanczosSum",
+    "pkg_desc": "calculate the Lanczos sum for the approximation of the gamma function",
+    "desc": "calculates the Lanczos sum for the approximation of the gamma function",
+    "short_desc": "Lanczos sum for the approximation of the gamma function",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "double",
+          "dtype": "float64"
+        },
+        "domain": [
+          {
+            "min": "-infinity",
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            -10,
+            10
+          ]
+        },
+        "example_values": [
+          1,
+          3.5,
+          4.5,
+          -0.5,
+          2,
+          3,
+          -3.5,
+          0.1,
+          4,
+          1.5,
+          5,
+          0.5,
+          2.5,
+          -1.5,
+          -2.5,
+          50,
+          100,
+          -50.5,
+          -100.5,
+          2.2
+        ]
+      }
+    ],
+    "output_policy": "real_floating_point_and_generic",
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "gamma",
+      "lanczos",
+      "approx",
+      "approximation",
+      "polynomial",
+      "rational",
+      "factorial"
+    ],
+    "extra_keywords": []
+  },
+  "@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "gammaLanczosSumExpGScaled",
+    "alias": "gammaLanczosSumExpGScaledf",
+    "pkg_desc": "calculate a scaled Lanczos sum for the approximation of the gamma function as a single precision floating-point number",
+    "desc": "calculates a scaled Lanczos sum for the approximation of the gamma function as a single precision floating-point number",
+    "short_desc": "scaled Lanczos sum for the approximation of the gamma function",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "float",
+          "dtype": "float32"
+        },
+        "domain": [
+          {
+            "min": "-infinity",
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            -10,
+            10
+          ]
+        },
+        "example_values": [
+          1,
+          3.5,
+          4.5,
+          -0.5,
+          2,
+          3,
+          -3.5,
+          0.1,
+          4,
+          1.5,
+          5,
+          0.5,
+          2.5,
+          -1.5,
+          -2.5,
+          50,
+          100,
+          -50.5,
+          -100.5,
+          2.2
+        ]
+      }
+    ],
+    "output_policy": "real_floating_point_and_generic",
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "float",
+        "dtype": "float32"
+      }
+    },
+    "keywords": [
+      "gamma",
+      "lanczos",
+      "scaled",
+      "approx",
+      "approximation",
+      "polynomial",
+      "rational",
+      "factorial"
+    ],
+    "extra_keywords": []
+  },
+  "@stdlib/math/base/special/gamma-lanczos-sum-expg-scaled": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "gammaLanczosSumExpGScaled",
+    "alias": "gammaLanczosSumExpGScaled",
+    "pkg_desc": "calculate a scaled Lanczos sum for the approximation of the gamma function",
+    "desc": "calculates a scaled Lanczos sum for the approximation of the gamma function",
+    "short_desc": "scaled Lanczos sum for the approximation of the gamma function",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "double",
+          "dtype": "float64"
+        },
+        "domain": [
+          {
+            "min": "-infinity",
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            -10,
+            10
+          ]
+        },
+        "example_values": [
+          1,
+          3.5,
+          4.5,
+          -0.5,
+          2,
+          3,
+          -3.5,
+          0.1,
+          4,
+          1.5,
+          5,
+          0.5,
+          2.5,
+          -1.5,
+          -2.5,
+          50,
+          100,
+          -50.5,
+          -100.5,
+          2.2
+        ]
+      }
+    ],
+    "output_policy": "real_floating_point_and_generic",
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "gamma",
+      "lanczos",
+      "scaled",
+      "approx",
+      "approximation",
+      "polynomial",
+      "rational",
+      "factorial"
+    ],
+    "extra_keywords": []
+  },
+  "@stdlib/math/base/special/gamma1pm1": {
+    "$schema": "math/base@v1.0",
+    "base_alias": "gamma1pm1",
+    "alias": "gamma1pm1",
+    "pkg_desc": "compute gamma(x+1) - 1",
+    "desc": "computes gamma(x+1) - 1",
+    "short_desc": "gamma(x+1) - 1",
+    "parameters": [
+      {
+        "name": "x",
+        "desc": "input value",
+        "type": {
+          "javascript": "number",
+          "jsdoc": "number",
+          "c": "double",
+          "dtype": "float64"
+        },
+        "domain": [
+          {
+            "min": "-infinity",
+            "max": "infinity"
+          }
+        ],
+        "rand": {
+          "prng": "random/base/uniform",
+          "parameters": [
+            -10,
+            10
+          ]
+        },
+        "example_values": [
+          1,
+          3.5,
+          4.5,
+          -0.5,
+          2,
+          3,
+          -3.5,
+          0.1,
+          4,
+          1.5,
+          5,
+          0.5,
+          2.5,
+          -1.5,
+          -2.5,
+          50,
+          100,
+          -50.5,
+          -100.5,
+          2.2
+        ]
+      }
+    ],
+    "output_policy": "real_floating_point_and_generic",
+    "returns": {
+      "desc": "function value",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "gamma",
+      "gammap1m1",
+      "gamma1pm1",
+      "factorial"
+    ],
+    "extra_keywords": []
+  },
   "@stdlib/math/base/special/gammasgnf": {
     "$schema": "math/base@v1.0",
     "base_alias": "gammasgn",