feat: add math/base/special/gamma-lanczos-sum-expg-scaledf

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: passed
  - task: lint_package_json
    status: passed
  - task: lint_repl_help
    status: passed
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: passed
  - task: lint_javascript_tests
    status: passed
  - task: lint_javascript_benchmarks
    status: passed
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: passed
  - task: lint_c_examples
    status: passed
  - task: lint_c_benchmarks
    status: passed
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: passed
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

Author: anandkaranubc

lib/node_modules/@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf/README.md

@@ -20,7 +20,7 @@ limitations under the License.
 
 # Gamma Scaled Lanczos Sum
 
-> Calculate a scaled Lanczos sum for the approximation of the [gamma function][gamma-function].
+> Calculate a scaled Lanczos sum for the approximation of the [gamma function][gamma-function] as a single precision floating-point number.
 
 <section class="intro">
 
@@ -32,26 +32,16 @@ The [Lanczos approximation][lanczos-approximation] for the [gamma function][gamm
 \Gamma ( n ) = \frac{(n+g-0.5)^{n-0.5}}{e^{n+g-0.5}} L_g(n)
 ```
 
-<!-- <div class="equation" align="center" data-raw-text="\Gamma ( n ) = \frac{(n+g-0.5)^{n-0.5}}{e^{n+g-0.5}} L_g(n)" data-equation="eq:lanczos_approximation">
-    <img src="https://cdn.jsdelivr.net/gh/stdlib-js/stdlib@bb29798906e119fcb2af99e94b60407a270c9b32/lib/node_modules/@stdlib/math/base/special/gamma-lanczos-sum-expg-scaled/docs/img/equation_lanczos_approximation.svg" alt="Lanczos approximation for gamma function.">
-    <br>
-</div> -->
-
 <!-- </equation> -->
 
-where `g` is an [arbitrary constant][@stdlib/constants/float64/gamma-lanczos-g] and `L_g(n)` is the Lanczos sum. The scaled Lanczos sum is given by
+where `g` is an [arbitrary constant][@stdlib/constants/float32/gamma-lanczos-g] and `L_g(n)` is the Lanczos sum. The scaled Lanczos sum is given by
 
 <!-- <equation class="equation" label="eq:scaled_lanczos_sum" align="center" raw="L_g(n) \cdot \exp(-g)" alt="Scaled Lanczos sum."> -->
 
 ```math
 L_g(n) \cdot \exp(-g)
 ```
 
-<!-- <div class="equation" align="center" data-raw-text="L_g(n) \cdot \exp(-g)" data-equation="eq:scaled_lanczos_sum">
-    <img src="https://cdn.jsdelivr.net/gh/stdlib-js/stdlib@bb29798906e119fcb2af99e94b60407a270c9b32/lib/node_modules/@stdlib/math/base/special/gamma-lanczos-sum-expg-scaled/docs/img/equation_scaled_lanczos_sum.svg" alt="Scaled Lanczos sum.">
-    <br>
-</div> -->
-
 <!-- </equation> -->
 
 </section>
@@ -62,23 +52,27 @@ L_g(n) \cdot \exp(-g)
 
 ## Usage
 
+<!-- eslint-disable id-length -->
+
 ```javascript
-var gammaLanczosSumExpGScaledf = require( '@stdlib/math/base/special/gamma-lanczos-sum-expg-scaled' );
+var gammaLanczosSumExpGScaledf = require( '@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf' );
 ```
 
 #### gammaLanczosSumExpGScaledf( x )
 
-Calculates the Lanczos sum for the approximation of the [gamma function][gamma-function] (scaled by `exp(-g)`, where `g = 10.900511`).
+Calculates the Lanczos sum for the approximation of the [gamma function][gamma-function] (scaled by `exp(-g)`, where `g = 1.42845618724823`) as a single precision floating-point number.
+
+<!-- eslint-disable id-length -->
 
 ```javascript
 var v = gammaLanczosSumExpGScaledf( 4.0 );
-// returns ~0.018
+// returns ~0.748
 
 v = gammaLanczosSumExpGScaledf( -1.5 );
-// returns ~25.337
+// returns ~0.193
 
 v = gammaLanczosSumExpGScaledf( -0.5 );
-// returns ~-12.911
+// returns ~-0.558
 
 v = gammaLanczosSumExpGScaledf( 0.5 );
 // returns ~1.772
@@ -100,13 +94,15 @@ v = gammaLanczosSumExpGScaledf( NaN );
 
 <!-- eslint no-undef: "error" -->
 
+<!-- eslint-disable id-length -->
+
 ```javascript
 var uniform = require( '@stdlib/random/array/uniform' );
 var logEachMap = require( '@stdlib/console/log-each-map' );
-var gammaLanczosSumExpGScaledf = require( '@stdlib/math/base/special/gamma-lanczos-sum-expg-scaled' );
+var gammaLanczosSumExpGScaledf = require( '@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf' );
 
 var opts = {
-    'dtype': 'float64'
+    'dtype': 'float32'
 };
 var x = uniform( 100, -10.0, 10.0, opts );
 
@@ -140,27 +136,27 @@ logEachMap( 'x: %0.4f, f(x): %0.4f', x, gammaLanczosSumExpGScaledf );
 ### Usage
 
 ```c
-#include "stdlib/math/base/special/gamma_lanczos_sum_expg_scaled.h"
+#include "stdlib/math/base/special/gamma_lanczos_sum_expg_scaledf.h"
 ```
 
 #### stdlib_base_gamma_lanczos_sum_expg_scaledf( x )
 
-Calculates the Lanczos sum for the approximation of the [gamma function][gamma-function] (scaled by `exp(-g)`, where `g = 10.900511`).
+Calculates the Lanczos sum for the approximation of the [gamma function][gamma-function] (scaled by `exp(-g)`, where `g = 1.42845618724823`) as a single precision floating-point number.
 
 ```c
-double out = stdlib_base_gamma_lanczos_sum_expg_scaledf( 4.0 );
-// returns ~0.018
+float out = stdlib_base_gamma_lanczos_sum_expg_scaledf( 4.0f );
+// returns ~0.018f
 
-out = stdlib_base_gamma_lanczos_sum_expg_scaledf( -1.5 );
-// returns ~25.337
+out = stdlib_base_gamma_lanczos_sum_expg_scaledf( -1.5f );
+// returns ~25.337f
 ```
 
 The function accepts the following arguments:
 
--   **x**: `[in] double` input value.
+-   **x**: `[in] float` input value.
 
 ```c
-double stdlib_base_gamma_lanczos_sum_expg_scaledf( const double x );
+float stdlib_base_gamma_lanczos_sum_expg_scaledf( const float x );
 ```
 
 </section>
@@ -182,18 +178,17 @@ double stdlib_base_gamma_lanczos_sum_expg_scaledf( const double x );
 ### Examples
 
 ```c
-#include "stdlib/math/base/special/gamma_lanczos_sum_expg_scaled.h"
-#include <stdlib.h>
+#include "stdlib/math/base/special/gamma_lanczos_sum_expg_scaledf.h"
 #include <stdio.h>
 
 int main( void ) {
-    const double x[] = { 4.0, -1.5, -0.5, 0.5 };
+    const float x[] = { 4.0f, -1.5f, -0.5f, 0.5f };
 
-    double y;
+    float y;
     int i;
     for ( i = 0; i < 4; i++ ) {
         y = stdlib_base_gamma_lanczos_sum_expg_scaledf( x[ i ] );
-        printf( "gamma_lanczos_sum_expg_scaled(%lf) = %lf\n", x[ i ], y );
+        printf( "gamma_lanczos_sum_expg_scaledf(%f) = %f\n", x[ i ], y );
     }
 }
 ```
@@ -210,13 +205,6 @@ int main( void ) {
 
 <section class="related">
 
-* * *
-
-## See Also
-
--   <span class="package-name">[`@stdlib/math/base/special/gamma`][@stdlib/math/base/special/gamma]</span><span class="delimiter">: </span><span class="description">gamma function.</span>
--   <span class="package-name">[`@stdlib/math/base/special/gamma-lanczos-sum`][@stdlib/math/base/special/gamma-lanczos-sum]</span><span class="delimiter">: </span><span class="description">calculate the Lanczos sum for the approximation of the gamma function.</span>
-
 </section>
 
 <!-- /.related -->
@@ -225,18 +213,14 @@ int main( void ) {
 
 <section class="links">
 
-[@stdlib/constants/float64/gamma-lanczos-g]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/constants/float64/gamma-lanczos-g
+[@stdlib/constants/float32/gamma-lanczos-g]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/constants/float32/gamma-lanczos-g
 
 [gamma-function]: https://en.wikipedia.org/wiki/Gamma_function
 
 [lanczos-approximation]: https://en.wikipedia.org/wiki/Lanczos_approximation
 
 <!-- <related-links> -->
 
-[@stdlib/math/base/special/gamma]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/math/base/special/gamma
-
-[@stdlib/math/base/special/gamma-lanczos-sum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/math/base/special/gamma-lanczos-sum
-
 <!-- </related-links> -->
 
 </section>

lib/node_modules/@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf/benchmark/benchmark.js

@@ -22,7 +22,7 @@
 
 var bench = require( '@stdlib/bench' );
 var uniform = require( '@stdlib/random/array/uniform' );
-var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var pkg = require( './../package.json' ).name;
 var gammaLanczosSumExpGScaledf = require( './../lib' );
 
@@ -34,17 +34,19 @@ bench( pkg, function benchmark( b ) {
 	var y;
 	var i;
 
-	x = uniform( 100, -50.0, 50.0 );
+	x = uniform( 100, -50.0, 50.0, {
+		'dtype': 'float32'
+	});
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
 		y = gammaLanczosSumExpGScaledf( x[ i%x.length ] );
-		if ( isnan( y ) ) {
+		if ( isnanf( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
 	}
 	b.toc();
-	if ( isnan( y ) ) {
+	if ( isnanf( y ) ) {
 		b.fail( 'should not return NaN' );
 	}
 	b.pass( 'benchmark finished' );

lib/node_modules/@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf/benchmark/benchmark.native.js

@@ -23,7 +23,7 @@
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
 var uniform = require( '@stdlib/random/array/uniform' );
-var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var pkg = require( './../package.json' ).name;
 
@@ -43,17 +43,19 @@ bench( pkg+'::native', opts, function benchmark( b ) {
 	var y;
 	var i;
 
-	x = uniform( 100, -50.0, 50.0 );
+	x = uniform( 100, -50.0, 50.0, {
+		'dtype': 'float32'
+	});
 
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
 		y = gammaLanczosSumExpGScaledf( x[ i%x.length ] );
-		if ( isnan( y ) ) {
+		if ( isnanf( y ) ) {
 			b.fail( 'should not return NaN' );
 		}
 	}
 	b.toc();
-	if ( isnan( y ) ) {
+	if ( isnanf( y ) ) {
 		b.fail( 'should not return NaN' );
 	}
 	b.pass( 'benchmark finished' );

lib/node_modules/@stdlib/math/base/special/gamma-lanczos-sum-expg-scaledf/benchmark/c/native/benchmark.c

@@ -16,14 +16,14 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/special/gamma_lanczos_sum_expg_scaled.h"
+#include "stdlib/math/base/special/gamma_lanczos_sum_expg_scaledf.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
 #include <time.h>
 #include <sys/time.h>
 
-#define NAME "gamma_lanczos_sum_expg_scaled"
+#define NAME "gamma_lanczos_sum_expg_scaledf"
 #define ITERATIONS 1000000
 #define REPEATS 3
 
@@ -79,9 +79,9 @@ static double tic( void ) {
 *
 * @return random number
 */
-static double rand_double( void ) {
+static float rand_float( void ) {
 	int r = rand();
-	return (double)r / ( (double)RAND_MAX + 1.0 );
+	return (float)r / ( (float)RAND_MAX + 1.0f );
 }
 
 /**
@@ -91,13 +91,13 @@ static double rand_double( void ) {
 */
 static double benchmark( void ) {
 	double elapsed;
-	double x[ 100 ];
-	double y;
+	float x[ 100 ];
 	double t;
+	float y;
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = ( 100.0 * rand_double() ) - 50.0;
+		x[ i ] = ( 100.0f*rand_float() ) - 50.0f;
 	}
 
 	t = tic();