feat: add JS implementation for math/base/special/sicif

---
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: passed
  - task: lint_r
    status: na
  - task: lint_c_src
    status: na
  - task: lint_c_examples
    status: na
  - 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/sicif/README.md

@@ -2,7 +2,7 @@
 
 @license Apache-2.0
 
-Copyright (c) 2018 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.
 
 -->
 
-# sici
+# sicif
 
-> Compute the sine and cosine integrals.
+> Compute the sine and cosine integrals in single-precision floating-point format.
 
 <section class="intro">
 
@@ -32,11 +32,6 @@ The sine integral is defined as
 \int_0^x \frac{\sin t}{t}\; dt
 ```
 
-<!-- <div class="equation" align="center" data-raw-text="\int_0^x \frac{\sin t}{t}\; dt" data-equation="eq:sine_integral">
-    <img src="https://cdn.jsdelivr.net/gh/stdlib-js/stdlib@9d139bffba6eb2d0113c8412c1a9bf60ca59ac60/lib/node_modules/@stdlib/math/base/special/sici/docs/img/equation_sine_integral.svg" alt="Sine integral.">
-    <br>
-</div> -->
-
 <!-- </equation> -->
 
 and the cosine integral is defined as
@@ -47,11 +42,6 @@ and the cosine integral is defined as
 \gamma + \log( x ) + \int_0^x \frac{\cos t - 1}{t}\; dt
 ```
 
-<!-- <div class="equation" align="center" data-raw-text="\gamma + \log( x ) + \int_0^x \frac{\cos t - 1}{t}\; dt" data-equation="eq:cosine_integral">
-    <img src="https://cdn.jsdelivr.net/gh/stdlib-js/stdlib@9d139bffba6eb2d0113c8412c1a9bf60ca59ac60/lib/node_modules/@stdlib/math/base/special/sici/docs/img/equation_cosine_integral.svg" alt="Cosine integral.">
-    <br>
-</div> -->
-
 <!-- </equation> -->
 
 where `γ` is the [Euler-Mascheroni][eulergamma] constant.
@@ -65,38 +55,38 @@ where `γ` is the [Euler-Mascheroni][eulergamma] constant.
 ## Usage
 
 ```javascript
-var sici = require( '@stdlib/math/base/special/sici' );
+var sicif = require( '@stdlib/math/base/special/sicif' );
 ```
 
-#### sici( x )
+#### sicif( x )
 
-Computes the sine and cosine integrals.
+Computes the sine and cosine integrals in single-precision floating-point format.
 
 ```javascript
-var v = sici( 3.0 );
+var v = sicif( 3.0 );
 // returns [ ~1.849, ~0.12 ]
 
-v = sici( 0.0 );
+v = sicif( 0.0 );
 // returns [ 0.0, -Infinity ]
 
-v = sici( -9.0 );
+v = sicif( -9.0 );
 // returns [ ~-1.665, ~0.055 ]
 
-v = sici( NaN );
+v = sicif( NaN );
 // returns [ NaN, NaN ]
 ```
 
-#### sici.assign( x, out, stride, offset )
+#### sicif.assign( x, out, stride, offset )
 
-Computes the sine and cosine integrals and assigns results to a provided output array.
+Computes the sine and cosine integrals in single-precision floating-point format and assigns results to a provided output array.
 
 ```javascript
-var Float64Array = require( '@stdlib/array/float64' );
+var Float32Array = require( '@stdlib/array/float32' );
 
-var out = new Float64Array( 2 );
+var out = new Float32Array( 2 );
 
-var v = sici.assign( 3.0, out, 1, 0 );
-// returns <Float64Array>[ ~1.849, ~0.12 ]
+var v = sicif.assign( 3.0, out, 1, 0 );
+// returns <Float32Array>[ ~1.849, ~0.12 ]
 
 var bool = ( v === out );
 // returns true
@@ -114,15 +104,15 @@ var bool = ( v === out );
 
 ```javascript
 var randu = require( '@stdlib/random/base/randu' );
-var sici = require( '@stdlib/math/base/special/sici' );
+var sicif = require( '@stdlib/math/base/special/sicif' );
 
 var x;
 var y;
 var i;
 
 for ( i = 0; i < 100; i++ ) {
     x = randu() * 100.0;
-    y = sici( x );
+    y = sicif( x );
     console.log( 'si(%d) = %d, ci(%d) = %d', x, y[ 0 ], x, y[ 1 ] );
 }
 ```

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

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2018 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.
@@ -23,9 +23,9 @@
 var bench = require( '@stdlib/bench' );
 var uniform = require( '@stdlib/random/array/uniform' );
 var isArray = require( '@stdlib/assert/is-array' );
-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 sici = require( './../lib' );
+var sicif = require( './../lib' );
 
 
 // MAIN //
@@ -35,20 +35,22 @@ 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 = sici( x[ i%x.length ] );
-		if ( isnan( y[ 0 ] ) ) {
+		y = sicif( x[ i%x.length ] );
+		if ( isnanf( y[ 0 ] ) ) {
 			b.fail( 'should not return NaN' );
 		}
 	}
 	b.toc();
 	if ( !isArray( y ) ) {
 		b.fail( 'should return an array' );
 	}
-	if ( isnan( y[ 0 ] ) || isnan( y[ 1 ] ) ) {
+	if ( isnanf( y[ 0 ] ) || isnanf( y[ 1 ] ) ) {
 		b.fail( 'should not return NaN' );
 	}
 	b.pass( 'benchmark finished' );
@@ -62,20 +64,22 @@ bench( pkg+':assign', function benchmark( b ) {
 	var i;
 
 	out = [ 0.0, 0.0 ];
-	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 = sici.assign( x[ i%x.length ], out, 1, 0 );
-		if ( isnan( y[ 0 ] ) ) {
+		y = sicif.assign( x[ i%x.length ], out, 1, 0 );
+		if ( isnanf( y[ 0 ] ) ) {
 			b.fail( 'should not return NaN' );
 		}
 	}
 	b.toc();
 	if ( !isArray( y ) || y !== out ) {
 		b.fail( 'should return the output array' );
 	}
-	if ( isnan( y[ 0 ] ) || isnan( y[ 1 ] ) ) {
+	if ( isnanf( y[ 0 ] ) || isnanf( y[ 1 ] ) ) {
 		b.fail( 'should not return NaN' );
 	}
 	b.pass( 'benchmark finished' );

lib/node_modules/@stdlib/math/base/special/sicif/benchmark/c/cephes/Makefile

@@ -1,7 +1,7 @@
 #/
 # @license Apache-2.0
 #
-# Copyright (c) 2018 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/sicif/benchmark/c/cephes/benchmark.c

@@ -1,7 +1,7 @@
 /**
 * @license Apache-2.0
 *
-* Copyright (c) 2018 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,14 +22,14 @@
 #include <time.h>
 #include <sys/time.h>
 
-#define NAME "sici"
+#define NAME "sicif"
 #define ITERATIONS 1000000
 #define REPEATS 3
 
 /**
 * Define prototypes for external functions.
 */
-extern int sici( double x, double *si, double *ci );
+extern int sici( float x, float *si, float *ci );
 
 /**
 * Prints the TAP version.
@@ -79,13 +79,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) );
 }
 
 /**
@@ -94,15 +96,15 @@ static double rand_double( void ) {
 * @return elapsed time in seconds
 */
 static double benchmark( void ) {
-	double x[ 100 ];
 	double elapsed;
-	double y;
-	double z;
+	float x[ 100 ];
 	double t;
+	float y;
+	float z;
 	int i;
 
 	for ( i = 0; i < 100; i++ ) {
-		x[ i ] = ( 100.0*rand_double() ) - 50.0;
+		x[ i ] = random_uniform( -50.0f, 50.0f );
 	}
 
 	t = tic();