feat: add math/base/special/ellipkf

---
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: passed
  - task: lint_license_headers
    status: passed
---

Author: anandkaranubc

lib/node_modules/@stdlib/math/base/special/ellipkf/README.md

@@ -2,7 +2,7 @@
 
 @license Apache-2.0
 
-Copyright (c) 2022 The Stdlib Authors.
+Copyright (c) 2025 The Stdlib Authors.
 
 Licensed under the Apache License, Version 2.0 (the "License");
 you may not use this file except in compliance with the License.
@@ -18,9 +18,9 @@ limitations under the License.
 
 -->
 
-# ellipk
+# ellipkf
 
-> Compute the [complete elliptic integral of the first kind][elliptic-integral].
+> Compute the [complete elliptic integral of the first kind][elliptic-integral] for a single-precision floating-point number.
 
 <section class="intro">
 
@@ -32,11 +32,6 @@ The [complete elliptic integral of the first kind][elliptic-integral] is defined
 K(m)=\int_0^\tfrac{\pi}{2} \frac{d\theta}{\sqrt{1-m \sin^2\theta}}
 ```
 
-<!-- <div class="equation" align="center" data-raw-text="K(m)=\int_0^\tfrac{\pi}{2} \frac{d\theta}{\sqrt{1-m \sin^2\theta}}" data-equation="eq:complete_elliptic_integral_first_kind">
-    <img src="https://cdn.jsdelivr.net/gh/stdlib-js/stdlib@fa2cf1d471224644251a972f9424f30bde62c2c3/lib/node_modules/@stdlib/math/base/special/ellipk/docs/img/equation_complete_elliptic_integral_first_kind.svg" alt="Complete elliptic integral of the first kind">
-    <br>
-</div> -->
-
 <!-- </equation> -->
 
 where the parameter `m` is related to the modulus `k` by `m = k^2`.
@@ -50,30 +45,30 @@ where the parameter `m` is related to the modulus `k` by `m = k^2`.
 ## Usage
 
 ```javascript
-var ellipk = require( '@stdlib/math/base/special/ellipk' );
+var ellipkf = require( '@stdlib/math/base/special/ellipkf' );
 ```
 
-#### ellipk( m )
+#### ellipkf( m )
 
-Computes the [complete elliptic integral of the first kind][elliptic-integral].
+Computes the [complete elliptic integral of the first kind][elliptic-integral] for a single-precision floating-point number.
 
 ```javascript
-var v = ellipk( 0.5 );
+var v = ellipkf( 0.5 );
 // returns ~1.854
 
-v = ellipk( -1.0 );
+v = ellipkf( -1.0 );
 // returns ~1.311
 
-v = ellipk( 2.0 );
+v = ellipkf( 2.0 );
 // returns NaN
 
-v = ellipk( Infinity );
+v = ellipkf( Infinity );
 // returns NaN
 
-v = ellipk( -Infinity );
+v = ellipkf( -Infinity );
 // returns NaN
 
-v = ellipk( NaN );
+v = ellipkf( NaN );
 // returns NaN
 ```
 
@@ -100,14 +95,14 @@ v = ellipk( NaN );
 ```javascript
 var uniform = require( '@stdlib/random/array/uniform' );
 var logEachMap = require( '@stdlib/console/log-each-map' );
-var ellipk = require( '@stdlib/math/base/special/ellipk' );
+var ellipkf = require( '@stdlib/math/base/special/ellipkf' );
 
 var opts = {
-    'dtype': 'float64'
+    'dtype': 'float32'
 };
 var m = uniform( 100, -1.0, 1.0, opts );
 
-logEachMap( 'ellipk(%0.4f) = %0.4f', m, ellipk );
+logEachMap( 'ellipkf(%0.4f) = %0.4f', m, ellipkf );
 ```
 
 </section>
@@ -137,27 +132,27 @@ logEachMap( 'ellipk(%0.4f) = %0.4f', m, ellipk );
 ### Usage
 
 ```c
-#include "stdlib/math/base/special/ellipk.h"
+#include "stdlib/math/base/special/ellipkf.h"
 ```
 
-#### stdlib_base_ellipk( m )
+#### stdlib_base_ellipkf( m )
 
-Computes the [complete elliptic integral of the first kind][elliptic-integral].
+Computes the [complete elliptic integral of the first kind][elliptic-integral] for a single-precision floating-point number.
 
 ```c
-double out = stdlib_base_ellipk( 0.5 );
-// returns ~1.854
+float out = stdlib_base_ellipkf( 0.5f );
+// returns ~1.854f
 
-out = stdlib_base_ellipk( -1.0 );
-// returns ~1.311
+out = stdlib_base_ellipkf( -1.0f );
+// returns ~1.311f
 ```
 
 The function accepts the following arguments:
 
--   **x**: `[in] double` input value.
+-   **x**: `[in] float` input value.
 
 ```c
-double stdlib_base_ellipk( const double m );
+float stdlib_base_ellipkf( const float m );
 ```
 
 </section>
@@ -179,19 +174,19 @@ double stdlib_base_ellipk( const double m );
 ### Examples
 
 ```c
-#include "stdlib/math/base/special/ellipk.h"
+#include "stdlib/math/base/special/ellipkf.h"
 #include <stdlib.h>
 #include <stdio.h>
 
 int main( void ) {
-    double m;
-    double v;
+    float m;
+    float v;
     int i;
 
     for ( i = 0; i < 100; i++ ) {
-        m = -1.0 + ( ( (double)rand() / (double)RAND_MAX ) * 2.0 );
-        v = stdlib_base_ellipk( m );
-        printf( "ellipk(%lf) = %lf\n", m, v );
+        m = -1.0f + ( ( (float)rand() / (float)RAND_MAX ) * 2.0f );
+        v = stdlib_base_ellipkf( m );
+        printf( "ellipkf(%f) = %f\n", m, v );
     }
 }
 ```
@@ -208,11 +203,6 @@ int main( void ) {
 
 <section class="references">
 
-## References
-
--   Fukushima, Toshio. 2009. "Fast computation of complete elliptic integrals and Jacobian elliptic functions." _Celestial Mechanics and Dynamical Astronomy_ 105 (4): 305. doi:[10.1007/s10569-009-9228-z][@fukushima:2009a].
--   Fukushima, Toshio. 2015. "Precise and fast computation of complete elliptic integrals by piecewise minimax rational function approximation." _Journal of Computational and Applied Mathematics_ 282 (July): 71–76. doi:[10.1016/j.cam.2014.12.038][@fukushima:2015a].
-
 </section>
 
 <!-- /.references -->
@@ -221,13 +211,6 @@ int main( void ) {
 
 <section class="related">
 
-* * *
-
-## See Also
-
--   <span class="package-name">[`@stdlib/math/base/special/ellipe`][@stdlib/math/base/special/ellipe]</span><span class="delimiter">: </span><span class="description">compute the complete elliptic integral of the second kind.</span>
--   <span class="package-name">[`@stdlib/math/base/special/ellipj`][@stdlib/math/base/special/ellipj]</span><span class="delimiter">: </span><span class="description">compute the Jacobi elliptic functions sn, cn, and dn.</span>
-
 </section>
 
 <!-- /.related -->
@@ -238,16 +221,8 @@ int main( void ) {
 
 [elliptic-integral]: https://en.wikipedia.org/wiki/Elliptic_integral
 
-[@fukushima:2009a]: https://doi.org/10.1007/s10569-009-9228-z
-
-[@fukushima:2015a]: https://doi.org/10.1016/j.cam.2014.12.038
-
 <!-- <related-links> -->
 
-[@stdlib/math/base/special/ellipe]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/math/base/special/ellipe
-
-[@stdlib/math/base/special/ellipj]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/math/base/special/ellipj
-
 <!-- </related-links> -->
 
 </section>

lib/node_modules/@stdlib/math/base/special/ellipkf/benchmark/benchmark.js

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2019 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -21,27 +21,32 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
-var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var uniform = require( '@stdlib/random/array/uniform' );
+var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var pkg = require( './../package.json' ).name;
-var ellipk = require( './../lib' );
+var ellipkf = require( './../lib' );
 
 
 // MAIN //
 
 bench( pkg, function benchmark( b ) {
+	var x;
 	var y;
 	var i;
 
+	x = uniform( 100, -1.0, 1.0, {
+		'dtype': 'float32'
+	});
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = ellipk( randu() );
-		if ( isnan( y ) ) {
+		y = ellipkf( x[ i%x.length ] );
+		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/ellipkf/benchmark/benchmark.native.js

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2022 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -22,35 +22,40 @@
 
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
-var randu = require( '@stdlib/random/base/randu' );
-var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var uniform = require( '@stdlib/random/array/uniform' );
+var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var tryRequire = require( '@stdlib/utils/try-require' );
 var pkg = require( './../package.json' ).name;
 
 
 // VARIABLES //
 
-var ellipk = tryRequire( resolve( __dirname, './../lib/native.js' ) );
+var ellipkf = tryRequire( resolve( __dirname, './../lib/native.js' ) );
 var opts = {
-	'skip': ( ellipk instanceof Error )
+	'skip': ( ellipkf instanceof Error )
 };
 
 
 // MAIN //
 
 bench( pkg+'::native', opts, function benchmark( b ) {
+	var x;
 	var y;
 	var i;
 
+	x = uniform( 100, -1.0, 1.0, {
+		'dtype': 'float32'
+	});
+
 	b.tic();
 	for ( i = 0; i < b.iterations; i++ ) {
-		y = ellipk( randu() );
-		if ( isnan( y ) ) {
+		y = ellipkf( x[ i%x.length ] );
+		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/ellipkf/benchmark/c/native/Makefile

@@ -1,7 +1,7 @@
 #/
 # @license Apache-2.0
 #
-# Copyright (c) 2022 The Stdlib Authors.
+# Copyright (c) 2025 The Stdlib Authors.
 #
 # Licensed under the Apache License, Version 2.0 (the "License");
 # you may not use this file except in compliance with the License.

lib/node_modules/@stdlib/math/base/special/ellipkf/benchmark/c/native/benchmark.c

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2022 The Stdlib Authors.
+* Copyright (c) 2025 The Stdlib Authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
@@ -16,14 +16,14 @@
 * limitations under the License.
 */
 
-#include "stdlib/math/base/special/ellipk.h"
+#include "stdlib/math/base/special/ellipkf.h"
 #include <stdlib.h>
 #include <stdio.h>
 #include <math.h>
 #include <time.h>
 #include <sys/time.h>
 
-#define NAME "ellipk"
+#define NAME "ellipkf"
 #define ITERATIONS 1000000
 #define REPEATS 3
 
@@ -75,13 +75,15 @@ static double tic( void ) {
 }
 
 /**
-* Generates a random number on the interval [0,1).
+* Generates a random number on the interval [min,max).
 *
-* @return random number
+* @param min    minimum value (inclusive)
+* @param max    maximum value (exclusive)
+* @return       random number
 */
-static double rand_double( void ) {
-	int r = rand();
-	return (double)r / ( (double)RAND_MAX + 1.0 );
+static float random_uniform( const float min, const float max ) {
+	float v = (float)rand() / ( (float)RAND_MAX + 1.0f );
+	return min + ( v*(max-min) );
 }
 
 /**
@@ -90,16 +92,19 @@ static double rand_double( void ) {
 * @return elapsed time in seconds
 */
 static double benchmark( void ) {
+	float x[ 100 ];
 	double elapsed;
-	double x;
-	double y;
 	double t;
+	float y;
 	int i;
 
+	for ( i = 0; i < 100; i++ ) {
+		x[ i ] = random_uniform( -1.0f, 1.0f );
+	}
+
 	t = tic();
 	for ( i = 0; i < ITERATIONS; i++ ) {
-		x = rand_double();
-		y = stdlib_base_ellipk( x );
+		y = stdlib_base_ellipkf( x[ i%100 ] );
 		if ( y != y ) {
 			printf( "should not return NaN\n" );
 			break;