From e85f3947aaca3cbbd6ab4f84e41577928e0a10a4 Mon Sep 17 00:00:00 2001
From: Aman Bhansali <92033532+aman-095@users.noreply.github.com>
Date: Thu, 19 Sep 2024 03:01:58 +0530
Subject: [PATCH 1/5] feat: add C `ndarray` implementation for
 `blas/base/saxpy`

PR-URL: https://github.com/stdlib-js/stdlib/pull/2918
Ref: https://github.com/stdlib-js/stdlib/issues/2039
Reviewed-by: Athan Reines <kgryte@gmail.com>
---
 .../@stdlib/blas/base/daxpy/manifest.json     |   6 +
 .../@stdlib/blas/base/saxpy/README.md         | 136 ++++++++++++++++++
 .../base/saxpy/benchmark/c/benchmark.length.c |  44 +++++-
 .../blas/base/saxpy/examples/c/example.c      |   8 ++
 .../saxpy/include/stdlib/blas/base/saxpy.h    |   9 +-
 .../include/stdlib/blas/base/saxpy_cblas.h    |   4 +-
 .../@stdlib/blas/base/saxpy/lib/ndarray.js    |   4 +-
 .../blas/base/saxpy/lib/ndarray.native.js     |  19 +--
 .../@stdlib/blas/base/saxpy/manifest.json     | 101 ++++++++++---
 .../@stdlib/blas/base/saxpy/src/addon.c       |  26 +++-
 .../@stdlib/blas/base/saxpy/src/saxpy.c       |  56 +-------
 .../@stdlib/blas/base/saxpy/src/saxpy_cblas.c |  24 +++-
 .../@stdlib/blas/base/saxpy/src/saxpy_f.c     |  22 ++-
 .../blas/base/saxpy/src/saxpy_ndarray.c       |  83 +++++++++++
 14 files changed, 445 insertions(+), 97 deletions(-)
 create mode 100644 lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_ndarray.c

diff --git a/lib/node_modules/@stdlib/blas/base/daxpy/manifest.json b/lib/node_modules/@stdlib/blas/base/daxpy/manifest.json
index 17dc2243c60..87718cfd396 100644
--- a/lib/node_modules/@stdlib/blas/base/daxpy/manifest.json
+++ b/lib/node_modules/@stdlib/blas/base/daxpy/manifest.json
@@ -134,6 +134,7 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
         "@stdlib/strided/base/min-view-buffer-index"
       ]
     },
@@ -154,6 +155,7 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
         "@stdlib/strided/base/min-view-buffer-index"
       ]
     },
@@ -262,6 +264,7 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
         "@stdlib/strided/base/min-view-buffer-index"
       ]
     },
@@ -281,6 +284,7 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
         "@stdlib/strided/base/min-view-buffer-index"
       ]
     },
@@ -328,6 +332,7 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
         "@stdlib/strided/base/min-view-buffer-index"
       ]
     },
@@ -348,6 +353,7 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
         "@stdlib/strided/base/min-view-buffer-index"
       ]
     },
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/README.md b/lib/node_modules/@stdlib/blas/base/saxpy/README.md
index 71bb966ad8e..efac1e63552 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/README.md
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/README.md
@@ -163,6 +163,142 @@ console.log( y );
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/blas/base/saxpy.h"
+```
+
+#### c_saxpy( N, alpha, \*X, strideX, \*Y, strideY )
+
+Multiplies a vector `X` by a constant and adds the result to `Y`.
+
+```c
+const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f };
+float y[] = { 0.0f, 0.0f, 0.0f, 0.0f };
+
+c_saxpy( 4, 5.0f, x, 1, y, 1 );
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **alpha**: `[in] float` scalar constant.
+-   **X**: `[in] float*` input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **Y**: `[inout] float*` output array.
+-   **strideY**: `[in CBLAS_INT` index increment for `Y`.
+
+```c
+void c_saxpy( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, float *Y, const CBLAS_INT strideY );
+```
+
+#### c_saxpy_ndarray( N, alpha, \*X, strideX, offsetX, \*Y, strideY, offsetY )
+
+Multiplies a vector `X` by a constant and adds the result to `Y` using alternative indexing semantics.
+
+```c
+const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f };
+float y[] = { 0.0f, 0.0f, 0.0f, 0.0f };
+
+c_saxpy_ndarray( 4, 5.0f, x, 1, 0, y, 1, 0 );
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **alpha**: `[in] float` scalar constant.
+-   **X**: `[in] float*` input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **offsetX**: `[in] CBLAS_INT` starting index for `X`.
+-   **Y**: `[inout] float*` output array.
+-   **strideY**: `[in CBLAS_INT` index increment for `Y`.
+-   **offsetY**: `[in] CBLAS_INT` starting index for `Y`.
+
+```c
+void c_saxpy_ndarray( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/blas/base/saxpy.h"
+#include <stdio.h>
+
+int main( void ) {
+    // Create strided arrays:
+    const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f, 5.0f, 6.0f, 7.0f, 8.0f };
+    float y[] = { 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f, 0.0f };
+
+    // Specify the number of elements:
+    const int N = 4;
+
+    // Specify stride lengths:
+    const int strideX = 2;
+    const int strideY = -2;
+
+    // Compute `a*x + y`:
+    c_saxpy( N, 5.0f, x, strideX, y, strideY );
+
+    // Print the result:
+    for ( int i = 0; i < 8; i++ ) {
+        printf( "y[ %i ] = %f\n", i, y[ i ] );
+    }
+
+    // Compute `a*x + y`:
+    c_saxpy_ndarray( N, 5.0f, x, strideX, 1, y, strideY, 7 );
+
+    // Print the result:
+    for ( int i = 0; i < 8; i++ ) {
+        printf( "y[ %i ] = %f\n", i, y[ i ] );
+    }
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/base/saxpy/benchmark/c/benchmark.length.c
index 23f658eafee..ed96c9959f9 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/benchmark/c/benchmark.length.c
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/benchmark/c/benchmark.length.c
@@ -94,7 +94,7 @@ static float rand_float( void ) {
 * @param len          array length
 * @return elapsed time in seconds
 */
-static double benchmark( int iterations, int len ) {
+static double benchmark1( int iterations, int len ) {
 	double elapsed;
 	float x[ len ];
 	float y[ len ];
@@ -120,6 +120,39 @@ static double benchmark( int iterations, int len ) {
 	return elapsed;
 }
 
+/**
+* Runs a benchmark.
+*
+* @param iterations   number of iterations
+* @param len          array length
+* @return elapsed time in seconds
+*/
+static double benchmark2( int iterations, int len ) {
+	double elapsed;
+	float x[ len ];
+	float y[ len ];
+	double t;
+	int i;
+
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( rand_float()*100.0f ) - 200.0f;
+		y[ i ] = ( rand_float()*10000.0f ) - 20000.0f;
+	}
+	t = tic();
+	for ( i = 0; i < iterations; i++ ) {
+		c_saxpy_ndarray( len, 5.0f, x, 1, 0, y, 1, 0 );
+		if ( y[ 0 ] != y[ 0 ] ) {
+			printf( "should not return NaN\n" );
+			break;
+		}
+	}
+	elapsed = tic() - t;
+	if ( y[ 0 ] != y[ 0 ] ) {
+		printf( "should not return NaN\n" );
+	}
+	return elapsed;
+}
+
 /**
 * Main execution sequence.
 */
@@ -142,7 +175,14 @@ int main( void ) {
 		for ( j = 0; j < REPEATS; j++ ) {
 			count += 1;
 			printf( "# c::%s:len=%d\n", NAME, len );
-			elapsed = benchmark( iter, len );
+			elapsed = benchmark1( iter, len );
+			print_results( iter, elapsed );
+			printf( "ok %d benchmark finished\n", count );
+		}
+		for ( j = 0; j < REPEATS; j++ ) {
+			count += 1;
+			printf( "# c::%s:ndarray:len=%d\n", NAME, len );
+			elapsed = benchmark2( iter, len );
 			print_results( iter, elapsed );
 			printf( "ok %d benchmark finished\n", count );
 		}
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/saxpy/examples/c/example.c
index f60b7ad5c39..efbbd38a0a0 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/examples/c/example.c
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/examples/c/example.c
@@ -38,4 +38,12 @@ int main( void ) {
 	for ( int i = 0; i < 8; i++ ) {
 		printf( "y[ %i ] = %f\n", i, y[ i ] );
 	}
+
+	// Compute `a*x + y`:
+	c_saxpy_ndarray( N, 5.0f, x, strideX, 1, y, strideY, 7 );
+
+	// Print the result:
+	for ( int i = 0; i < 8; i++ ) {
+		printf( "y[ %i ] = %f\n", i, y[ i ] );
+	}
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy.h b/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy.h
index 8c85c830b01..250510d2372 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy.h
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy.h
@@ -22,6 +22,8 @@
 #ifndef SAXPY_H
 #define SAXPY_H
 
+#include "stdlib/blas/base/shared.h"
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -32,7 +34,12 @@ extern "C" {
 /**
 * Multiplies a vector `X` by a constant and adds the result to `Y`.
 */
-void c_saxpy( const int N, const float alpha, const float *X, const int strideX, float *Y, const int strideY );
+void API_SUFFIX(c_saxpy)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, float *Y, const CBLAS_INT strideY );
+
+/**
+* Multiplies a vector `X` by a constant and adds the result to `Y` using alternative indexing semantics.
+*/
+void API_SUFFIX(c_saxpy_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy_cblas.h b/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy_cblas.h
index 1cd7e454741..048ce2edaa4 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy_cblas.h
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/include/stdlib/blas/base/saxpy_cblas.h
@@ -22,6 +22,8 @@
 #ifndef SAXPY_CBLAS_H
 #define SAXPY_CBLAS_H
 
+#include "stdlib/blas/base/shared.h"
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -32,7 +34,7 @@ extern "C" {
 /**
 * Multiplies a vector `X` by a constant and adds the result to `Y`.
 */
-void cblas_saxpy( const int N, const float alpha, const float *X, const int strideX, float *Y, const int strideY );
+void API_SUFFIX(cblas_saxpy)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, float *Y, const CBLAS_INT strideY );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.js
index 2dd79fa9dd7..90972ba9d3d 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.js
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.js
@@ -37,10 +37,10 @@ var M = 4;
 * @param {number} alpha - scalar
 * @param {Float32Array} x - input array
 * @param {integer} strideX - `x` stride length
-* @param {NonNegativeInteger} offsetX - starting `x` index
+* @param {NonNegativeInteger} offsetX - starting index for `x`
 * @param {Float32Array} y - output array
 * @param {integer} strideY - `y` stride length
-* @param {NonNegativeInteger} offsetY - starting `y` index
+* @param {NonNegativeInteger} offsetY - starting index for `y`
 * @returns {Float32Array} output array
 *
 * @example
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.native.js
index b373e6fdfeb..ef785b4d69f 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.native.js
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/lib/ndarray.native.js
@@ -20,9 +20,7 @@
 
 // MODULES //
 
-var minViewBufferIndex = require( '@stdlib/strided/base/min-view-buffer-index' );
-var offsetView = require( '@stdlib/strided/base/offset-view' );
-var addon = require( './saxpy.native.js' );
+var addon = require( './../src/addon.node' );
 
 
 // MAIN //
@@ -34,10 +32,10 @@ var addon = require( './saxpy.native.js' );
 * @param {number} alpha - scalar
 * @param {Float32Array} x - input array
 * @param {integer} strideX - `x` stride length
-* @param {NonNegativeInteger} offsetX - starting `x` index
+* @param {NonNegativeInteger} offsetX - starting index for `x`
 * @param {Float32Array} y - output array
 * @param {integer} strideY - `y` stride length
-* @param {NonNegativeInteger} offsetY - starting `y` index
+* @param {NonNegativeInteger} offsetY - starting index for `y`
 * @returns {Float32Array} output array
 *
 * @example
@@ -51,16 +49,7 @@ var addon = require( './saxpy.native.js' );
 * // y => <Float32Array>[ 6.0, 11.0, 16.0, 21.0, 26.0 ]
 */
 function saxpy( N, alpha, x, strideX, offsetX, y, strideY, offsetY ) {
-	var viewX;
-	var viewY;
-
-	offsetX = minViewBufferIndex( N, strideX, offsetX );
-	offsetY = minViewBufferIndex( N, strideY, offsetY );
-
-	viewX = offsetView( x, offsetX );
-	viewY = offsetView( y, offsetY );
-
-	addon( N, alpha, viewX, strideX, viewY, strideY );
+	addon.ndarray( N, alpha, x, strideX, offsetX, y, strideY, offsetY );
 	return y;
 }
 
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/manifest.json b/lib/node_modules/@stdlib/blas/base/saxpy/manifest.json
index 86bce79a68f..ce201234aee 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/manifest.json
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/manifest.json
@@ -43,6 +43,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -56,14 +58,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -71,14 +77,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -98,6 +108,8 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -121,7 +133,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -139,7 +154,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -157,6 +175,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -170,14 +190,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -185,14 +209,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -211,6 +239,8 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -233,7 +263,10 @@
         "-lblas"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -250,7 +283,10 @@
         "-lblas"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -270,6 +306,8 @@
       ],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -293,7 +331,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -311,7 +352,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -320,7 +364,8 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
@@ -328,6 +373,8 @@
       "libraries": [],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -341,14 +388,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -356,14 +407,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -372,14 +427,18 @@
       "blas": "",
       "wasm": true,
       "src": [
-        "./src/saxpy.c"
+        "./src/saxpy.c",
+        "./src/saxpy_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     }
   ]
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/src/addon.c b/lib/node_modules/@stdlib/blas/base/saxpy/src/addon.c
index da508a76138..6178332b966 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/src/addon.c
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/src/addon.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/blas/base/saxpy.h"
+#include "stdlib/blas/base/shared.h"
 #include "stdlib/napi/export.h"
 #include "stdlib/napi/argv.h"
 #include "stdlib/napi/argv_float.h"
@@ -39,8 +40,29 @@ static napi_value addon( napi_env env, napi_callback_info info ) {
 	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 5 );
 	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 );
 	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, Y, N, strideY, argv, 4 );
-	c_saxpy( N, alpha, X, strideX, Y, strideY );
+	API_SUFFIX(c_saxpy)( N, alpha, X, strideX, Y, strideY );
 	return NULL;
 }
 
-STDLIB_NAPI_MODULE_EXPORT_FCN( addon )
+/**
+* Receives JavaScript callback invocation data.
+*
+* @param env    environment under which the function is invoked
+* @param info   callback data
+* @return       Node-API value
+*/
+static napi_value addon_method( napi_env env, napi_callback_info info ) {
+	STDLIB_NAPI_ARGV( env, info, argv, argc, 8 );
+	STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 );
+	STDLIB_NAPI_ARGV_FLOAT( env, alpha, argv, 1 );
+	STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 );
+	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 6 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetY, argv, 7 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, Y, N, strideY, argv, 5 );
+	API_SUFFIX(c_saxpy_ndarray)( N, alpha, X, strideX, offsetX, Y, strideY, offsetY );
+	return NULL;
+}
+
+STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method )
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy.c b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy.c
index 2868c9f1141..384cd84639f 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy.c
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy.c
@@ -17,6 +17,8 @@
 */
 
 #include "stdlib/blas/base/saxpy.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/stride2offset.h"
 
 /**
 * Multiplies a vector `X` by a constant and adds the result to `Y`.
@@ -28,54 +30,8 @@
 * @param Y        output array
 * @param strideY  Y stride length
 */
-void c_saxpy( const int N, const float alpha, const float *X, const int strideX, float *Y, const int strideY ) {
-	int ix;
-	int iy;
-	int i;
-	int m;
-
-	if ( N <= 0 ) {
-		return;
-	}
-	// If `alpha` is `0`, then `y` is unchanged...
-	if ( alpha == 0.0f ) {
-		return;
-	}
-	// If both strides are equal to `1`, use unrolled loops...
-	if ( strideX == 1 && strideY == 1 ) {
-		m = N % 4;
-
-		// If we have a remainder, do a clean-up loop...
-		if ( m > 0 ) {
-			for ( i = 0; i < m; i++ ) {
-				Y[ i ] += alpha * X[ i ];
-			}
-			if ( N < 4 ) {
-				return;
-			}
-		}
-		for ( i = m; i < N; i += 4 ) {
-			Y[ i ] += alpha * X[ i ];
-			Y[ i+1 ] += alpha * X[ i+1 ];
-			Y[ i+2 ] += alpha * X[ i+2 ];
-			Y[ i+3 ] += alpha * X[ i+3 ];
-		}
-		return;
-	}
-	if ( strideX < 0 ) {
-		ix = (1-N) * strideX;
-	} else {
-		ix = 0;
-	}
-	if ( strideY < 0 ) {
-		iy = (1-N) * strideY;
-	} else {
-		iy = 0;
-	}
-	for ( i = 0; i < N; i++ ) {
-		Y[ iy ] += alpha * X[ ix ];
-		ix += strideX;
-		iy += strideY;
-	}
-	return;
+void API_SUFFIX(c_saxpy)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, float *Y, const CBLAS_INT strideY ) {
+	CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX );
+	CBLAS_INT oy = stdlib_strided_stride2offset( N, strideY );
+	API_SUFFIX(c_saxpy_ndarray)( N, alpha, X, strideX, ox, Y, strideY, oy );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_cblas.c b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_cblas.c
index 31469a0306b..d57d9095121 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_cblas.c
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_cblas.c
@@ -18,6 +18,8 @@
 
 #include "stdlib/blas/base/saxpy.h"
 #include "stdlib/blas/base/saxpy_cblas.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Multiplies a vector `X` by a constant and adds the result to `Y`.
@@ -29,6 +31,24 @@
 * @param Y        output array
 * @param strideY  Y stride length
 */
-void c_saxpy( const int N, const float alpha, const float *X, const int strideX, float *Y, const int strideY ) {
-	cblas_saxpy( N, alpha, X, strideX, Y, strideY );
+void API_SUFFIX(c_saxpy)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, float *Y, const CBLAS_INT strideY ) {
+	API_SUFFIX(cblas_saxpy)( N, alpha, X, strideX, Y, strideY );
+}
+
+/**
+* Multiplies a vector `X` by a constant and adds the result to `Y` using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param alpha    scalar
+* @param X        input array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        output array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+*/
+void API_SUFFIX(c_saxpy_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	X += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	Y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
+	API_SUFFIX(cblas_saxpy)( N, alpha, X, strideX, Y, strideY );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_f.c b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_f.c
index 29dbe869b56..64f49315bb8 100644
--- a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_f.c
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_f.c
@@ -18,6 +18,8 @@
 
 #include "stdlib/blas/base/saxpy.h"
 #include "stdlib/blas/base/saxpy_fortran.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Multiplies a vector `X` by a constant and adds the result to `Y`.
@@ -29,6 +31,24 @@
 * @param Y        output array
 * @param strideY  Y stride length
 */
-void c_saxpy( const int N, const float alpha, const float *X, const int strideX, float *Y, const int strideY ) {
+void API_SUFFIX(c_saxpy)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, float *Y, const CBLAS_INT strideY ) {
+	saxpy( &N, &alpha, X, &strideX, Y, &strideY );
+}
+
+/**
+* Multiplies a vector `X` by a constant and adds the result to `Y`.
+*
+* @param N        number of indexed elements
+* @param alpha    scalar
+* @param X        input array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        output array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+*/
+void API_SUFFIX(c_saxpy_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	X += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	Y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
 	saxpy( &N, &alpha, X, &strideX, Y, &strideY );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_ndarray.c b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_ndarray.c
new file mode 100644
index 00000000000..0b4922baec2
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/saxpy/src/saxpy_ndarray.c
@@ -0,0 +1,83 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2018 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.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+#include "stdlib/blas/base/saxpy.h"
+#include "stdlib/blas/base/shared.h"
+
+static const CBLAS_INT M = 4;
+
+/**
+* Multiplies a vector `X` by a constant and adds the result to `Y` using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param alpha    scalar
+* @param X        input array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        output array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+*/
+void API_SUFFIX(c_saxpy_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	CBLAS_INT ix;
+	CBLAS_INT iy;
+	CBLAS_INT i;
+	CBLAS_INT m;
+
+	if ( N <= 0 ) {
+		return;
+	}
+	// If `alpha` is `0`, then `y` is unchanged...
+	if ( alpha == 0.0f ) {
+		return;
+	}
+	ix = offsetX;
+	iy = offsetY;
+
+	// If both strides are equal to `1`, use unrolled loops...
+	if ( strideX == 1 && strideY == 1 ) {
+		m = N % M;
+
+		// If we have a remainder, do a clean-up loop...
+		if ( m > 0 ) {
+			for ( i = 0; i < m; i++ ) {
+				Y[ iy ] += alpha * X[ ix ];
+				ix += strideX;
+				iy += strideY;
+			}
+			if ( N < M ) {
+				return;
+			}
+		}
+		for ( i = m; i < N; i += M ) {
+			Y[ iy ] += alpha * X[ ix ];
+			Y[ iy+1 ] += alpha * X[ ix+1 ];
+			Y[ iy+2 ] += alpha * X[ ix+2 ];
+			Y[ iy+3 ] += alpha * X[ ix+3 ];
+			ix += M;
+			iy += M;
+		}
+		return;
+	}
+	for ( i = 0; i < N; i++ ) {
+		Y[ iy ] += alpha * X[ ix ];
+		ix += strideX;
+		iy += strideY;
+	}
+	return;
+}

From 2bedae91dbdd0080b38ae0413f5e8f26b88db95b Mon Sep 17 00:00:00 2001
From: Aman Bhansali <92033532+aman-095@users.noreply.github.com>
Date: Thu, 19 Sep 2024 03:14:18 +0530
Subject: [PATCH 2/5] feat: add C `ndarray` implementation for `blas/base/sdot`

PR-URL: https://github.com/stdlib-js/stdlib/pull/2919
Ref: https://github.com/stdlib-js/stdlib/issues/2039
Co-authored-by: Athan Reines <kgryte@gmail.com>
Reviewed-by: Athan Reines <kgryte@gmail.com>
Signed-off-by: Athan Reines <kgryte@gmail.com>
---
 .../@stdlib/blas/base/sdot/README.md          | 132 ++++++++++++++++++
 .../base/sdot/benchmark/c/benchmark.length.c  |  46 +++++-
 .../blas/base/sdot/examples/c/example.c       |   6 +
 .../base/sdot/include/stdlib/blas/base/sdot.h |   9 +-
 .../include/stdlib/blas/base/sdot_cblas.h     |   4 +-
 .../blas/base/sdot/lib/ndarray.native.js      |  15 +-
 .../@stdlib/blas/base/sdot/manifest.json      | 101 +++++++++++---
 .../@stdlib/blas/base/sdot/src/addon.c        |  25 +++-
 .../@stdlib/blas/base/sdot/src/sdot.c         |  51 +------
 .../@stdlib/blas/base/sdot/src/sdot_cblas.c   |  24 +++-
 .../@stdlib/blas/base/sdot/src/sdot_f.c       |  27 +++-
 .../@stdlib/blas/base/sdot/src/sdot_ndarray.c |  84 +++++++++++
 12 files changed, 436 insertions(+), 88 deletions(-)
 create mode 100644 lib/node_modules/@stdlib/blas/base/sdot/src/sdot_ndarray.c

diff --git a/lib/node_modules/@stdlib/blas/base/sdot/README.md b/lib/node_modules/@stdlib/blas/base/sdot/README.md
index 9812103ade9..63bece13641 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/README.md
+++ b/lib/node_modules/@stdlib/blas/base/sdot/README.md
@@ -177,6 +177,138 @@ console.log( out );
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/blas/base/sdot.h"
+```
+
+#### c_sdot( N, \*X, strideX, \*Y, strideY )
+
+Computes the dot product of two single-precision floating-point vectors.
+
+```c
+const float x[] = { 4.0f, 2.0f, -3.0f, 5.0f, -1.0f };
+const float y[] = { 2.0f, 6.0f, -1.0f, -4.0f, 8.0f };
+
+float d = c_sdot( 5, x, 1, y, 1 );
+// returns -5.0f
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **X**: `[in] float*` first input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **Y**: `[in] float*` second input array.
+-   **strideY**: `[in] CBLAS_INT` index increment for `Y`.
+
+```c
+float c_sdot( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY );
+```
+
+#### c_sdot_ndarray( N, \*X, strideX, offsetX, \*Y, strideY, offsetY )
+
+Computes the dot product of two single-precision floating-point vectors using alternative indexing semantics.
+
+```c
+const float x[] = { 4.0f, 2.0f, -3.0f, 5.0f, -1.0f };
+const float y[] = { 2.0f, 6.0f, -1.0f, -4.0f, 8.0f };
+
+float d = c_sdot_ndarray( 3, x, 1, 2, y, 1, 2 );
+// returns -25.0f
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **X**: `[in] float*` first input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **offsetX**: `[in] CBLAS_INT` starting index for `X`.
+-   **Y**: `[in] float*` second input array.
+-   **strideY**: `[in] CBLAS_INT` index increment for `Y`.
+-   **offsetY**: `[in] CBLAS_INT` starting index for `Y`.
+
+```c
+float c_sdot_ndarray( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/blas/base/sdot.h"
+#include <stdio.h>
+
+int main( void ) {
+    // Create strided arrays:
+    const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+    const float y[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+
+    // Specify the number of indexed elements:
+    const int N = 8;
+
+    // Specify strides:
+    const int strideX = 1;
+    const int strideY = -1;
+
+    // Compute the dot product:
+    float d = c_sdot( N, x, strideX, y, strideY );
+
+    // Print the result:
+    printf( "dot product: %f\n", d );
+
+    // Compute the dot product:
+    d = c_sdot_ndarray( N, x, strideX, 0, y, strideY, 7 );
+
+    // Print the result:
+    printf( "dot product: %f\n", d );
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/base/sdot/benchmark/c/benchmark.length.c
index dd1970a400d..726b6d6bad8 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/benchmark/c/benchmark.length.c
+++ b/lib/node_modules/@stdlib/blas/base/sdot/benchmark/c/benchmark.length.c
@@ -94,7 +94,7 @@ static float rand_float( void ) {
 * @param len          array length
 * @return elapsed time in seconds
 */
-static double benchmark( int iterations, int len ) {
+static double benchmark1( int iterations, int len ) {
 	double elapsed;
 	float x[ len ];
 	float y[ len ];
@@ -122,6 +122,41 @@ static double benchmark( int iterations, int len ) {
 	return elapsed;
 }
 
+/**
+* Runs a benchmark.
+*
+* @param iterations   number of iterations
+* @param len          array length
+* @return elapsed time in seconds
+*/
+static double benchmark2( int iterations, int len ) {
+	double elapsed;
+	float x[ len ];
+	float y[ len ];
+	float z;
+	double t;
+	int i;
+
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( rand_float()*20000.0f ) - 10000.0f;
+		y[ i ] = ( rand_float()*20000.0f ) - 10000.0f;
+	}
+	z = 0.0f;
+	t = tic();
+	for ( i = 0; i < iterations; i++ ) {
+		z = c_sdot_ndarray( len, x, 1, 0, y, 1, 0 );
+		if ( z != z ) {
+			printf( "should not return NaN\n" );
+			break;
+		}
+	}
+	elapsed = tic() - t;
+	if ( z != z ) {
+		printf( "should not return NaN\n" );
+	}
+	return elapsed;
+}
+
 /**
 * Main execution sequence.
 */
@@ -144,7 +179,14 @@ int main( void ) {
 		for ( j = 0; j < REPEATS; j++ ) {
 			count += 1;
 			printf( "# c::%s:len=%d\n", NAME, len );
-			elapsed = benchmark( iter, len );
+			elapsed = benchmark1( iter, len );
+			print_results( iter, elapsed );
+			printf( "ok %d benchmark finished\n", count );
+		}
+		for ( j = 0; j < REPEATS; j++ ) {
+			count += 1;
+			printf( "# c::%s:ndarray:len=%d\n", NAME, len );
+			elapsed = benchmark2( iter, len );
 			print_results( iter, elapsed );
 			printf( "ok %d benchmark finished\n", count );
 		}
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/sdot/examples/c/example.c
index dd07f29a32a..500617bb4ef 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/examples/c/example.c
+++ b/lib/node_modules/@stdlib/blas/base/sdot/examples/c/example.c
@@ -36,4 +36,10 @@ int main( void ) {
 
 	// Print the result:
 	printf( "dot product: %f\n", d );
+
+	// Compute the dot product:
+	d = c_sdot_ndarray( N, x, strideX, 0, y, strideY, 7 );
+
+	// Print the result:
+	printf( "dot product: %f\n", d );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot.h b/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot.h
index ca3453ff658..8a08025cba0 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot.h
+++ b/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot.h
@@ -22,6 +22,8 @@
 #ifndef SDOT_H
 #define SDOT_H
 
+#include "stdlib/blas/base/shared.h"
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -32,7 +34,12 @@ extern "C" {
 /**
 * Computes the dot product of two single-precision floating-point vectors.
 */
-float c_sdot( const int N, const float *X, const int strideX, const float *Y, const int strideY );
+float API_SUFFIX(c_sdot)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY );
+
+/**
+* Computes the dot product of two single-precision floating-point vectors using alternative indexing semantics.
+*/
+float API_SUFFIX(c_sdot_ndarray)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot_cblas.h b/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot_cblas.h
index d02eaf1817f..5bae0337f38 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot_cblas.h
+++ b/lib/node_modules/@stdlib/blas/base/sdot/include/stdlib/blas/base/sdot_cblas.h
@@ -22,6 +22,8 @@
 #ifndef SDOT_CBLAS_H
 #define SDOT_CBLAS_H
 
+#include "stdlib/blas/base/shared.h"
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -32,7 +34,7 @@ extern "C" {
 /**
 * Computes the dot product of two single-precision floating-point vectors.
 */
-float cblas_sdot( const int N, const float *X, const int strideX, const float *Y, const int strideY );
+float API_SUFFIX(cblas_sdot)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/sdot/lib/ndarray.native.js
index b032bdcc820..0aa3602a73c 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/lib/ndarray.native.js
+++ b/lib/node_modules/@stdlib/blas/base/sdot/lib/ndarray.native.js
@@ -20,9 +20,7 @@
 
 // MODULES //
 
-var minViewBufferIndex = require( '@stdlib/strided/base/min-view-buffer-index' );
-var offsetView = require( '@stdlib/strided/base/offset-view' );
-var addon = require( './sdot.native.js' );
+var addon = require( './../src/addon.node' );
 
 
 // MAIN //
@@ -49,16 +47,7 @@ var addon = require( './sdot.native.js' );
 * // returns -5.0
 */
 function sdot( N, x, strideX, offsetX, y, strideY, offsetY ) {
-	var viewX;
-	var viewY;
-
-	offsetX = minViewBufferIndex( N, strideX, offsetX );
-	offsetY = minViewBufferIndex( N, strideY, offsetY );
-
-	viewX = offsetView( x, offsetX );
-	viewY = offsetView( y, offsetY );
-
-	return addon( N, viewX, strideX, viewY, strideY );
+	return addon.ndarray( N, x, strideX, offsetX, y, strideY, offsetY );
 }
 
 
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/manifest.json b/lib/node_modules/@stdlib/blas/base/sdot/manifest.json
index 7d0e5b19ee2..2f2fd17c43d 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/manifest.json
+++ b/lib/node_modules/@stdlib/blas/base/sdot/manifest.json
@@ -45,6 +45,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
@@ -57,14 +59,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -72,14 +78,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -100,6 +110,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
@@ -122,7 +134,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -140,7 +155,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -160,6 +178,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
@@ -172,14 +192,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -187,14 +211,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -214,6 +242,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
@@ -235,7 +265,10 @@
         "-lblas"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -252,7 +285,10 @@
         "-lblas"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -273,6 +309,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
@@ -295,7 +333,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -313,7 +354,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -322,7 +366,8 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
@@ -331,6 +376,8 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
@@ -343,14 +390,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -358,14 +409,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -374,14 +429,18 @@
       "blas": "",
       "wasm": true,
       "src": [
-        "./src/sdot.c"
+        "./src/sdot.c",
+        "./src/sdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     }
   ]
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/src/addon.c b/lib/node_modules/@stdlib/blas/base/sdot/src/addon.c
index 5fa37f0dd74..712e45afb38 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/src/addon.c
+++ b/lib/node_modules/@stdlib/blas/base/sdot/src/addon.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/blas/base/sdot.h"
+#include "stdlib/blas/base/shared.h"
 #include "stdlib/napi/export.h"
 #include "stdlib/napi/argv.h"
 #include "stdlib/napi/argv_int64.h"
@@ -38,8 +39,28 @@ static napi_value addon( napi_env env, napi_callback_info info ) {
 	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 4 );
 	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 1 );
 	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, Y, N, strideY, argv, 3 );
-	STDLIB_NAPI_CREATE_DOUBLE( env, (double)c_sdot( N, X, strideX, Y, strideY ), v );
+	STDLIB_NAPI_CREATE_DOUBLE( env, (double)API_SUFFIX(c_sdot)( N, X, strideX, Y, strideY ), v );
 	return v;
 }
 
-STDLIB_NAPI_MODULE_EXPORT_FCN( addon )
+/**
+* Receives JavaScript callback invocation data.
+*
+* @param env    environment under which the function is invoked
+* @param info   callback data
+* @return       Node-API value
+*/
+static napi_value addon_method( napi_env env, napi_callback_info info ) {
+	STDLIB_NAPI_ARGV( env, info, argv, argc, 7 );
+	STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 );
+	STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 2 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 3 );
+	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 5 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetY, argv, 6 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 1 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, Y, N, strideY, argv, 4 );
+	STDLIB_NAPI_CREATE_DOUBLE( env, (double)API_SUFFIX(c_sdot_ndarray)( N, X, strideX, offsetX, Y, strideY, offsetY ), v );
+	return v;
+}
+
+STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method )
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot.c b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot.c
index 8ed48cfb8c0..40782617c30 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot.c
+++ b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot.c
@@ -22,6 +22,8 @@
 * @see <a href="http://www.netlib.org/lapack/expolore-html/df/d28/group__single__blas__level1.html">sdot</a>
 */
 #include "stdlib/blas/base/sdot.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/stride2offset.h"
 
 /**
 * Computes the dot product of two single-precision floating-point vectors.
@@ -33,50 +35,9 @@
 * @param strideY  Y stride length
 * @return         the dot product
 */
-float c_sdot( const int N, const float *X, const int strideX, const float *Y, const int strideY ) {
-	float dot;
-	int ix;
-	int iy;
-	int m;
-	int i;
-
-	dot = 0.0f;
-	if ( N <= 0 ) {
-		return dot;
-	}
-	// If both strides are equal to `1`, use unrolled loops...
-	if ( strideX == 1 && strideY == 1 ) {
-		m = N % 5;
-
-		// If we have a remainder, do a clean-up loop...
-		if ( m > 0 ) {
-			for ( i = 0; i < m; i++ ) {
-				dot += X[ i ] * Y[ i ];
-			}
-		}
-		if ( N < 5 ) {
-			return dot;
-		}
-		for ( i = m; i < N; i += 5 ) {
-			dot += ( X[i]*Y[i] ) + ( X[i+1]*Y[i+1] ) + ( X[i+2]*Y[i+2] ) + ( X[i+3]*Y[i+3] ) + ( X[i+4]*Y[i+4] );
-		}
-		return dot;
-	}
-	if ( strideX < 0 ) {
-		ix = (1-N) * strideX;
-	} else {
-		ix = 0;
-	}
-	if ( strideY < 0 ) {
-		iy = (1-N) * strideY;
-	} else {
-		iy = 0;
-	}
-	for ( i = 0; i < N; i++ ) {
-		dot += X[ ix ] * Y[ iy];
-		ix += strideX;
-		iy += strideY;
-	}
-	return dot;
+float API_SUFFIX(c_sdot)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY ) {
+	CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX );
+	CBLAS_INT oy = stdlib_strided_stride2offset( N, strideY );
+	return API_SUFFIX(c_sdot_ndarray)( N, X, strideX, ox, Y, strideY, oy );
 }
 
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_cblas.c b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_cblas.c
index a725236c56c..37af408a28b 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_cblas.c
+++ b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_cblas.c
@@ -18,6 +18,8 @@
 
 #include "stdlib/blas/base/sdot.h"
 #include "stdlib/blas/base/sdot_cblas.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Computes the dot product of two single-precision floating-point vectors.
@@ -29,6 +31,24 @@
 * @param strideY  Y stride length
 * @return         the dot product
 */
-float c_sdot( const int N, const float *X, const int strideX, const float *Y, const int strideY ) {
-	return cblas_sdot( N, X, strideX, Y, strideY );
+float API_SUFFIX(c_sdot)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY ) {
+	return API_SUFFIX(cblas_sdot)( N, X, strideX, Y, strideY );
+}
+
+/**
+* Computes the dot product of two single-precision floating-point vectors using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param X        first array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+* @return         the dot product
+*/
+float API_SUFFIX(c_sdot_ndarray)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	X += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	Y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
+	return API_SUFFIX(cblas_sdot)( N, X, strideX, Y, strideY );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_f.c b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_f.c
index 743baf7d52b..df545e7782c 100644
--- a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_f.c
+++ b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_f.c
@@ -23,6 +23,8 @@
 */
 #include "stdlib/blas/base/sdot.h"
 #include "stdlib/blas/base/sdot_fortran.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Computes the dot product of two single-precision floating-point vectors.
@@ -36,8 +38,31 @@
 * @param strideY  Y stride length
 * @return         the dot product
 */
-float c_sdot( const int N, const float *X, const int strideX, const float *Y, const int strideY ) {
+float API_SUFFIX(c_sdot)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY ) {
 	float dot;
 	sdotsub( &N, X, &strideX, Y, &strideY, &dot );
 	return dot;
 }
+
+/**
+* Computes the dot product of two single-precision floating-point vectors using alternative indexing semantics.
+*
+* Arguments are passed by reference to a Fortran subroutine implementing `sdot`.
+*
+* @param N        number of indexed elements
+* @param X        first array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second array
+* @param strideY  Y stride length
+* @param strideY  starting index for Y
+* @return         the dot product
+*/
+float API_SUFFIX(c_sdot_ndarray)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	float dot;
+
+	X += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	Y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
+	sdotsub( &N, X, &strideX, Y, &strideY, &dot );
+	return dot;
+}
diff --git a/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_ndarray.c b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_ndarray.c
new file mode 100644
index 00000000000..5829d725994
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/sdot/src/sdot_ndarray.c
@@ -0,0 +1,84 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2023 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.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+/**
+* Compute the dot product of two single-precision floating-point vectors.
+*
+* @see <a href="http://www.netlib.org/lapack/expolore-html/df/d28/group__single__blas__level1.html">sdot</a>
+*/
+#include "stdlib/blas/base/sdot.h"
+#include "stdlib/blas/base/shared.h"
+
+static const CBLAS_INT M = 5;
+
+/**
+* Computes the dot product of two single-precision floating-point vectors using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param X        first array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+* @return         the dot product
+*/
+float API_SUFFIX(c_sdot_ndarray)( const CBLAS_INT N, const float *X, const CBLAS_INT strideX, CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, CBLAS_INT offsetY ) {
+	CBLAS_INT ix;
+	CBLAS_INT iy;
+	CBLAS_INT m;
+	CBLAS_INT i;
+	float dot;
+
+	dot = 0.0f;
+	if ( N <= 0 ) {
+		return dot;
+	}
+	ix = offsetX;
+	iy = offsetY;
+	
+	// If both strides are equal to `1`, use unrolled loops...
+	if ( strideX == 1 && strideY == 1 ) {
+		m = N % M;
+
+		// If we have a remainder, do a clean-up loop...
+		if ( m > 0 ) {
+			for ( i = 0; i < m; i++ ) {
+				dot += X[ ix ] * Y[ iy ];
+				ix += strideX;
+				iy += strideY;
+			}
+		}
+		if ( N < M ) {
+			return dot;
+		}
+		for ( i = m; i < N; i += M ) {
+			dot += ( X[ ix ]*Y[ iy ] ) + ( X[ ix+1 ]*Y[ iy+1 ] ) + ( X[ ix+2 ]*Y[ iy+2 ] ) + ( X[ ix+3 ]*Y[ iy+3 ] ) + ( X[ ix+4 ]*Y[ iy+4 ] );
+			ix += M;
+			iy += M;
+		}
+		return dot;
+	}
+	for ( i = 0; i < N; i++ ) {
+		dot += X[ ix ] * Y[ iy];
+		ix += strideX;
+		iy += strideY;
+	}
+	return dot;
+}
+

From f7edbe044289b877586d0140467cb43cdfd18c19 Mon Sep 17 00:00:00 2001
From: Aman Bhansali <92033532+aman-095@users.noreply.github.com>
Date: Thu, 19 Sep 2024 03:14:53 +0530
Subject: [PATCH 3/5] feat: add C `ndarray` implementation for
 `blas/base/sdsdot`

PR-URL: https://github.com/stdlib-js/stdlib/pull/2921
Ref: https://github.com/stdlib-js/stdlib/issues/2039
Co-authored-by: Athan Reines <kgryte@gmail.com>
Reviewed-by: Athan Reines <kgryte@gmail.com>
Signed-off-by: Athan Reines <kgryte@gmail.com>
---
 .../@stdlib/blas/base/sdsdot/README.md        | 134 ++++++++++++++++++
 .../sdsdot/benchmark/c/benchmark.length.c     |  46 +++++-
 .../blas/base/sdsdot/examples/c/example.c     |   6 +
 .../sdsdot/include/stdlib/blas/base/sdsdot.h  |   9 +-
 .../include/stdlib/blas/base/sdsdot_cblas.h   |   4 +-
 .../blas/base/sdsdot/lib/ndarray.native.js    |  15 +-
 .../@stdlib/blas/base/sdsdot/manifest.json    | 113 +++++++++++----
 .../@stdlib/blas/base/sdsdot/src/addon.c      |  26 +++-
 .../@stdlib/blas/base/sdsdot/src/sdsdot.c     |  51 +------
 .../blas/base/sdsdot/src/sdsdot_cblas.c       |  25 +++-
 .../@stdlib/blas/base/sdsdot/src/sdsdot_f.c   |  28 +++-
 .../blas/base/sdsdot/src/sdsdot_ndarray.c     |  85 +++++++++++
 12 files changed, 448 insertions(+), 94 deletions(-)
 create mode 100644 lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_ndarray.c

diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/README.md b/lib/node_modules/@stdlib/blas/base/sdsdot/README.md
index 4f041688095..d6b8f9131d5 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/README.md
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/README.md
@@ -178,6 +178,140 @@ console.log( out );
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/blas/base/sdsdot.h"
+```
+
+#### c_sdsdot( N, scalar, \*X, strideX, \*Y, strideY )
+
+Calculates the dot product of vectors `x` and `y` with extended accumulation.
+
+```c
+const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+const float y[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+
+float v = c_sdsdot( 5, 0.0f, x, 1, y, -1 );
+// returns -120.0f
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **scalar**: `[in] float` scalar constant to add to dot product.
+-   **X**: `[in] float*` first input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **Y**: `[in] float*` second input array.
+-   **strideY**: `[in] CBLAS_INT` index increment for `Y`.
+
+```c
+float c_sdsdot( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY );
+```
+
+#### c_sdsdot_ndarray( N, scalar, \*X, strideX, offsetX, \*Y, strideY, offsetY )
+
+Calculates the dot product of vectors `x` and `y` with extended accumulation using alternative indexing semantics.
+
+```c
+const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+const float y[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+
+float v = c_sdsdot_ndarray( 5, 0.0f, x, 1, 0, y, -1, 7 );
+// returns -80.0f
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **scalar**: `[in] float` scalar constant to add to dot product.
+-   **X**: `[in] float*` first input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **offsetX**: `[in] CBLAS_INT` starting index for `X`.
+-   **Y**: `[in] float*` second input array.
+-   **strideY**: `[in] CBLAS_INT` index increment for `Y`.
+-   **offsetY**: `[in] CBLAS_INT` starting index for `Y`.
+
+```c
+float c_sdsdot_ndarray( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/blas/base/sdsdot.h"
+#include <stdio.h>
+
+int main( void ) {
+    // Create strided arrays:
+    const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+    const float y[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+
+    // Specify the number of indexed elements:
+    const int N = 8;
+
+    // Specify strides:
+    const int strideX = 1;
+    const int strideY = -1;
+
+    // Compute the dot product:
+    float d = c_sdsdot( N, 0.0f, x, strideX, y, strideY );
+
+    // Print the result:
+    printf( "dot product: %f\n", d );
+
+    // Compute the dot product:
+    d = c_sdsdot_ndarray( N, 0.0f, x, strideX, 0, y, strideY, 7 );
+
+    // Print the result:
+    printf( "dot product: %f\n", d );
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 * * *
 
 <section class="references">
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/base/sdsdot/benchmark/c/benchmark.length.c
index 82512090bde..565da1c5833 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/benchmark/c/benchmark.length.c
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/benchmark/c/benchmark.length.c
@@ -94,7 +94,7 @@ static float rand_float( void ) {
 * @param len          array length
 * @return elapsed time in seconds
 */
-static double benchmark( int iterations, int len ) {
+static double benchmark1( int iterations, int len ) {
 	double elapsed;
 	float x[ len ];
 	float y[ len ];
@@ -122,6 +122,41 @@ static double benchmark( int iterations, int len ) {
 	return elapsed;
 }
 
+/**
+* Runs a benchmark.
+*
+* @param iterations   number of iterations
+* @param len          array length
+* @return elapsed time in seconds
+*/
+static double benchmark2( int iterations, int len ) {
+	double elapsed;
+	float x[ len ];
+	float y[ len ];
+	float z;
+	double t;
+	int i;
+
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( rand_float()*20000.0f ) - 10000.0f;
+		y[ i ] = ( rand_float()*20000.0f ) - 10000.0f;
+	}
+	z = 0.0f;
+	t = tic();
+	for ( i = 0; i < iterations; i++ ) {
+		z = c_sdsdot_ndarray( len, 0.0f, x, 1, 0, y, 1, 0 );
+		if ( z != z ) {
+			printf( "should not return NaN\n" );
+			break;
+		}
+	}
+	elapsed = tic() - t;
+	if ( z != z ) {
+		printf( "should not return NaN\n" );
+	}
+	return elapsed;
+}
+
 /**
 * Main execution sequence.
 */
@@ -144,7 +179,14 @@ int main( void ) {
 		for ( j = 0; j < REPEATS; j++ ) {
 			count += 1;
 			printf( "# c::%s:len=%d\n", NAME, len );
-			elapsed = benchmark( iter, len );
+			elapsed = benchmark1( iter, len );
+			print_results( iter, elapsed );
+			printf( "ok %d benchmark finished\n", count );
+		}
+		for ( j = 0; j < REPEATS; j++ ) {
+			count += 1;
+			printf( "# c::%s:ndarray:len=%d\n", NAME, len );
+			elapsed = benchmark2( iter, len );
 			print_results( iter, elapsed );
 			printf( "ok %d benchmark finished\n", count );
 		}
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/sdsdot/examples/c/example.c
index 0238b0ff825..c46ff68ceb7 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/examples/c/example.c
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/examples/c/example.c
@@ -36,4 +36,10 @@ int main( void ) {
 
 	// Print the result:
 	printf( "dot product: %f\n", d );
+
+	// Compute the dot product:
+	d = c_sdsdot_ndarray( N, 0.0f, x, strideX, 0, y, strideY, 7 );
+
+	// Print the result:
+	printf( "dot product: %f\n", d );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot.h b/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot.h
index d4832b5a3b7..96c7690223d 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot.h
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot.h
@@ -22,6 +22,8 @@
 #ifndef SDSDOT_H
 #define SDSDOT_H
 
+#include "stdlib/blas/base/shared.h"
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -32,7 +34,12 @@ extern "C" {
 /**
 * Computes the dot product of two single-precision floating-point vectors with extended accumulation.
 */
-float c_sdsdot( const int N, const float scalar, const float *X, const int strideX, const float *Y, const int strideY );
+float API_SUFFIX(c_sdsdot)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY );
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation using alternative indexing semantics.
+*/
+float API_SUFFIX(c_sdsdot_ndarray)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot_cblas.h b/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot_cblas.h
index b2b0274eef4..31a15464d1c 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot_cblas.h
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/include/stdlib/blas/base/sdsdot_cblas.h
@@ -22,6 +22,8 @@
 #ifndef SDSDOT_CBLAS_H
 #define SDSDOT_CBLAS_H
 
+#include "stdlib/blas/base/shared.h"
+
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
 */
@@ -32,7 +34,7 @@ extern "C" {
 /**
 * Computes the dot product of two single-precision floating-point vectors with extended accumulation.
 */
-float cblas_sdsdot( const int N, const float scalar, const float *X, const int strideX, const float *Y, const int strideY );
+float API_SUFFIX(cblas_sdsdot)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/sdsdot/lib/ndarray.native.js
index 5a2d55a1a1e..8ab1aa3dd7f 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/lib/ndarray.native.js
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/lib/ndarray.native.js
@@ -20,9 +20,7 @@
 
 // MODULES //
 
-var minViewBufferIndex = require( '@stdlib/strided/base/min-view-buffer-index' );
-var offsetView = require( '@stdlib/strided/base/offset-view' );
-var addon = require( './sdsdot.native.js' );
+var addon = require( './../src/addon.node' );
 
 
 // MAIN //
@@ -50,16 +48,7 @@ var addon = require( './sdsdot.native.js' );
 * // returns -5.0
 */
 function sdsdot( N, scalar, x, strideX, offsetX, y, strideY, offsetY ) {
-	var viewX;
-	var viewY;
-
-	offsetX = minViewBufferIndex( N, strideX, offsetX );
-	offsetY = minViewBufferIndex( N, strideY, offsetY );
-
-	viewX = offsetView( x, offsetX );
-	viewY = offsetView( y, offsetY );
-
-	return addon( N, scalar, viewX, strideX, viewY, strideY );
+	return addon.ndarray( N, scalar, x, strideX, offsetX, y, strideY, offsetY );
 }
 
 
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/manifest.json b/lib/node_modules/@stdlib/blas/base/sdsdot/manifest.json
index 417ab8fdfff..5dd3391921a 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/manifest.json
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/manifest.json
@@ -45,9 +45,11 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
-        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-float",
+        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
         "@stdlib/napi/create-double"
       ]
@@ -58,14 +60,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -73,14 +79,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -101,9 +111,11 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
-        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-float",
+        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
         "@stdlib/napi/create-double"
       ]
@@ -124,7 +136,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -142,7 +157,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -162,9 +180,11 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
-        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-float",
+        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
         "@stdlib/napi/create-double"
       ]
@@ -175,14 +195,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -190,14 +214,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -217,9 +245,11 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
-        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-float",
+        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
         "@stdlib/napi/create-double"
       ]
@@ -239,7 +269,10 @@
         "-lblas"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -256,7 +289,10 @@
         "-lblas"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -277,9 +313,11 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index",
         "@stdlib/napi/argv",
-        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-float",
+        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
         "@stdlib/napi/create-double"
       ]
@@ -300,7 +338,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
     {
       "task": "examples",
@@ -318,7 +359,10 @@
         "-lpthread"
       ],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/min-view-buffer-index"
+      ]
     },
 
     {
@@ -327,7 +371,8 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
@@ -336,9 +381,11 @@
       "libpath": [],
       "dependencies": [
         "@stdlib/napi/export",
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset",
         "@stdlib/napi/argv",
-        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-float",
+        "@stdlib/napi/argv-int64",
         "@stdlib/napi/argv-strided-float32array",
         "@stdlib/napi/create-double"
       ]
@@ -349,14 +396,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
@@ -364,14 +415,18 @@
       "blas": "",
       "wasm": false,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
 
     {
@@ -380,14 +435,18 @@
       "blas": "",
       "wasm": true,
       "src": [
-        "./src/sdsdot.c"
+        "./src/sdsdot.c",
+        "./src/sdsdot_ndarray.c"
       ],
       "include": [
         "./include"
       ],
       "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     }
   ]
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/src/addon.c b/lib/node_modules/@stdlib/blas/base/sdsdot/src/addon.c
index 637e4371357..597ffa2c161 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/src/addon.c
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/src/addon.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/blas/base/sdsdot.h"
+#include "stdlib/blas/base/shared.h"
 #include "stdlib/napi/export.h"
 #include "stdlib/napi/argv.h"
 #include "stdlib/napi/argv_int64.h"
@@ -40,8 +41,29 @@ static napi_value addon( napi_env env, napi_callback_info info ) {
 	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 5 );
 	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 );
 	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, Y, N, strideY, argv, 4 );
-	STDLIB_NAPI_CREATE_DOUBLE( env, (double)c_sdsdot( N, scalar, X, strideX, Y, strideY ), v );
+	STDLIB_NAPI_CREATE_DOUBLE( env, (double)API_SUFFIX(c_sdsdot)( N, scalar, X, strideX, Y, strideY ), v );
 	return v;
 }
 
-STDLIB_NAPI_MODULE_EXPORT_FCN( addon )
+/**
+* Receives JavaScript callback invocation data.
+*
+* @param env    environment under which the function is invoked
+* @param info   callback data
+* @return       Node-API value
+*/
+static napi_value addon_method( napi_env env, napi_callback_info info ) {
+	STDLIB_NAPI_ARGV( env, info, argv, argc, 8 );
+	STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 );
+	STDLIB_NAPI_ARGV_FLOAT( env, scalar, argv, 1 );
+	STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 );
+	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 6 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetY, argv, 7 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, Y, N, strideY, argv, 5 );
+	STDLIB_NAPI_CREATE_DOUBLE( env, (double)API_SUFFIX(c_sdsdot_ndarray)( N, scalar, X, strideX, offsetX, Y, strideY, offsetY ), v );
+	return v;
+}
+
+STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method )
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot.c b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot.c
index ba3d601c02d..1d745b1cc2e 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot.c
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot.c
@@ -22,6 +22,8 @@
 * @see <a href="http://www.netlib.org/lapack/expolore-html/df/d28/group__single__blas__level1.html">sdsdot</a>
 */
 #include "stdlib/blas/base/sdsdot.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/stride2offset.h"
 
 /**
 * Computes the dot product of two single-precision floating-point vectors with extended accumulation.
@@ -34,50 +36,9 @@
 * @param strideY  Y stride length
 * @return         dot product
 */
-float c_sdsdot( const int N, const float scalar, const float *X, const int strideX, const float *Y, const int strideY ) {
-	double dot;
-	int ix;
-	int iy;
-	int m;
-	int i;
-
-	dot = (double)scalar;
-	if ( N <= 0 ) {
-		return dot;
-	}
-	// If both strides are equal to `1`, use unrolled loops...
-	if ( strideX == 1 && strideY == 1 ) {
-		m = N % 5;
-
-		// If we have a remainder, do a clean-up loop...
-		if ( m > 0 ) {
-			for ( i = 0; i < m; i++ ) {
-				dot += (double)X[ i ] * (double)Y[ i ];
-			}
-		}
-		if ( N < 5 ) {
-			return dot;
-		}
-		for ( i = m; i < N; i += 5 ) {
-			dot += ( (double)X[i]*(double)Y[i] ) + ( (double)X[i+1]*(double)Y[i+1]) + ( (double)X[i+2]*(double)Y[i+2] ) + ( (double)X[i+3]*(double)Y[i+3] ) + ( (double)X[i+4]*(double)Y[i+4] );
-		}
-		return dot;
-	}
-	if ( strideX < 0 ) {
-		ix = (1-N) * strideX;
-	} else {
-		ix = 0;
-	}
-	if ( strideY < 0 ) {
-		iy = (1-N) * strideY;
-	} else {
-		iy = 0;
-	}
-	for ( i = 0; i < N; i++ ) {
-		dot += (double)X[ ix ] * (double)Y[ iy ];
-		ix += strideX;
-		iy += strideY;
-	}
-	return dot;
+float API_SUFFIX(c_sdsdot)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY ) {
+	CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX );
+	CBLAS_INT oy = stdlib_strided_stride2offset( N, strideY );
+	return API_SUFFIX(c_sdsdot_ndarray)( N, scalar, X, strideX, ox, Y, strideY, oy );
 }
 
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_cblas.c b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_cblas.c
index bd47331e133..8a8fdeca4fe 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_cblas.c
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_cblas.c
@@ -18,6 +18,8 @@
 
 #include "stdlib/blas/base/sdsdot.h"
 #include "stdlib/blas/base/sdsdot_cblas.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Computes the dot product of two single-precision floating-point vectors with extended accumulation.
@@ -30,6 +32,25 @@
 * @param strideY  Y stride length
 * @return         dot product
 */
-float c_sdsdot( const int N, const float scalar, const float *X, const int strideX, const float *Y, const int strideY ) {
-	return cblas_sdsdot( N, scalar, X, strideX, Y, strideY );
+float API_SUFFIX(c_sdsdot)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY ) {
+	return API_SUFFIX(cblas_sdsdot)( N, scalar, X, strideX, Y, strideY );
+}
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param scalar   scalar constant added to the dot product
+* @param X        first array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+* @return         dot product
+*/
+float API_SUFFIX(c_sdsdot_ndarray)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	X += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	Y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
+	return API_SUFFIX(cblas_sdsdot)( N, scalar, X, strideX, Y, strideY );
 }
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_f.c b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_f.c
index 8e5de1ca0cd..e99160eaa69 100644
--- a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_f.c
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_f.c
@@ -23,6 +23,8 @@
  */
 #include "stdlib/blas/base/sdsdot.h"
 #include "stdlib/blas/base/sdsdot_fortran.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Computes the dot product of two single-precision floating-point vectors with extended accumulation.
@@ -37,8 +39,32 @@
 * @param strideY  Y stride length
 * @return         dot product
 */
-float c_sdsdot( const int N, const float scalar, const float *X, const int strideX, const float *Y, const int strideY ) {
+float API_SUFFIX(c_sdsdot)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const float *Y, const CBLAS_INT strideY ) {
 	float dot;
 	sdsdotsub( &N, &scalar, X, &strideX, Y, &strideY, &dot );
 	return dot;
 }
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation using alternative indexing semantics.
+*
+* Arguments are passed by reference to a Fortran subroutine implementing `sdsdot`.
+*
+* @param N        number of indexed elements
+* @param scalar   scalar constant added to the dot product
+* @param X        first array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+* @return         dot product
+*/
+float API_SUFFIX(c_sdsdot_ndarray)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	float dot;
+
+	X += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); // adjust array pointer
+	Y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); // adjust array pointer
+	sdsdotsub( &N, &scalar, X, &strideX, Y, &strideY, &dot );
+	return dot;
+}
diff --git a/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_ndarray.c b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_ndarray.c
new file mode 100644
index 00000000000..6aed77a944c
--- /dev/null
+++ b/lib/node_modules/@stdlib/blas/base/sdsdot/src/sdsdot_ndarray.c
@@ -0,0 +1,85 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2023 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.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+/**
+* Compute the dot product of two single-precision floating-point vectors with extended accumulation.
+*
+* @see <a href="http://www.netlib.org/lapack/expolore-html/df/d28/group__single__blas__level1.html">sdsdot</a>
+*/
+#include "stdlib/blas/base/sdsdot.h"
+#include "stdlib/blas/base/shared.h"
+
+static const CBLAS_INT M = 5;
+
+/**
+* Computes the dot product of two single-precision floating-point vectors with extended accumulation using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param scalar   scalar constant added to the dot product
+* @param X        first array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+* @return         dot product
+*/
+float API_SUFFIX(c_sdsdot_ndarray)( const CBLAS_INT N, const float scalar, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const float *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	CBLAS_INT ix;
+	CBLAS_INT iy;
+	CBLAS_INT m;
+	CBLAS_INT i;
+	double dot;
+
+	dot = (double)scalar;
+	if ( N <= 0 ) {
+		return dot;
+	}
+	ix = offsetX;
+	iy = offsetY;
+
+	// If both strides are equal to `1`, use unrolled loops...
+	if ( strideX == 1 && strideY == 1 ) {
+		m = N % M;
+
+		// If we have a remainder, do a clean-up loop...
+		if ( m > 0 ) {
+			for ( i = 0; i < m; i++ ) {
+				dot += (double)X[ ix ] * (double)Y[ iy ];
+				ix += strideX;
+				iy += strideY;
+			}
+		}
+		if ( N < M ) {
+			return dot;
+		}
+		for ( i = m; i < N; i += M ) {
+			dot += ( (double)X[ ix ]*(double)Y[ iy ] ) + ( (double)X[ ix+1 ]*(double)Y[ iy+1 ]) + ( (double)X[ ix+2 ]*(double)Y[ iy+2 ] ) + ( (double)X[ ix+3 ]*(double)Y[ iy+3 ] ) + ( (double)X[ ix+4 ]*(double)Y[ iy+4 ] );
+			ix += M;
+			iy += M;
+		}
+		return dot;
+	}
+	for ( i = 0; i < N; i++ ) {
+		dot += (double)X[ ix ] * (double)Y[ iy ];
+		ix += strideX;
+		iy += strideY;
+	}
+	return dot;
+}
+

From fd41e1bd1dc6b7c351702bdbf95978b343215b0b Mon Sep 17 00:00:00 2001
From: Muhammad Haris <101793258+headlessNode@users.noreply.github.com>
Date: Thu, 19 Sep 2024 03:05:58 +0500
Subject: [PATCH 4/5] feat: add C `ndarray` API to `blas/ext/base/sfill`

PR-URL: https://github.com/stdlib-js/stdlib/pull/2923
Co-authored-by: Athan Reines <kgryte@gmail.com>
Reviewed-by: Athan Reines <kgryte@gmail.com>
---
 .../@stdlib/blas/ext/base/sfill/README.md     | 139 ++++++++++++++++--
 .../ext/base/sfill/benchmark/benchmark.js     |   9 +-
 .../base/sfill/benchmark/benchmark.native.js  |   9 +-
 .../base/sfill/benchmark/benchmark.ndarray.js |   9 +-
 .../benchmark/benchmark.ndarray.native.js     |   9 +-
 .../base/sfill/benchmark/c/benchmark.length.c |  42 +++++-
 .../@stdlib/blas/ext/base/sfill/docs/repl.txt |  20 +--
 .../blas/ext/base/sfill/docs/types/index.d.ts |  12 +-
 .../blas/ext/base/sfill/examples/c/example.c  |   8 +-
 .../blas/ext/base/sfill/examples/index.js     |   9 +-
 .../include/stdlib/blas/ext/base/sfill.h      |   9 +-
 .../blas/ext/base/sfill/lib/ndarray.js        |  14 +-
 .../blas/ext/base/sfill/lib/ndarray.native.js |  15 +-
 .../@stdlib/blas/ext/base/sfill/lib/sfill.js  |   6 +-
 .../blas/ext/base/sfill/lib/sfill.native.js   |   6 +-
 .../@stdlib/blas/ext/base/sfill/manifest.json |  30 ++--
 .../@stdlib/blas/ext/base/sfill/src/addon.c   |  27 +++-
 .../ext/base/sfill/src/{sfill.c => main.c}    |  59 +++++---
 18 files changed, 314 insertions(+), 118 deletions(-)
 rename lib/node_modules/@stdlib/blas/ext/base/sfill/src/{sfill.c => main.c} (52%)

diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/README.md b/lib/node_modules/@stdlib/blas/ext/base/sfill/README.md
index a58f14ca374..a35e5d3dcf1 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/README.md
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/README.md
@@ -30,7 +30,7 @@ limitations under the License.
 var sfill = require( '@stdlib/blas/ext/base/sfill' );
 ```
 
-#### sfill( N, alpha, x, stride )
+#### sfill( N, alpha, x, strideX )
 
 Fills a single-precision floating-point strided array `x` with a specified scalar constant `alpha`.
 
@@ -48,7 +48,7 @@ The function has the following parameters:
 -   **N**: number of indexed elements.
 -   **alpha**: scalar constant.
 -   **x**: input [`Float32Array`][@stdlib/array/float32].
--   **stride**: index increment.
+-   **strideX**: index increment.
 
 The `N` and stride parameters determine which elements in the strided array are accessed at runtime. For example, to fill every other element
 
@@ -77,7 +77,7 @@ sfill( 3, 5.0, x1, 2 );
 // x0 => <Float32Array>[ 1.0, 5.0, 3.0, 5.0, 5.0, 5.0 ]
 ```
 
-#### sfill.ndarray( N, alpha, x, stride, offset )
+#### sfill.ndarray( N, alpha, x, strideX, offsetX )
 
 Fills a single-precision floating-point strided array `x` with a specified scalar constant `alpha` using alternative indexing semantics.
 
@@ -92,9 +92,9 @@ sfill.ndarray( x.length, 5.0, x, 1, 0 );
 
 The function has the following additional parameters:
 
--   **offset**: starting index.
+-   **offsetX**: starting index.
 
-While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying `buffer`, the offset parameter supports indexing semantics based on a starting index. For example, to access only the last three elements of the strided array
+While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameter supports indexing semantics based on a starting index. For example, to access only the last three elements of the strided array
 
 ```javascript
 var Float32Array = require( '@stdlib/array/float32' );
@@ -126,13 +126,12 @@ sfill.ndarray( 3, 5.0, x, 1, x.length-3 );
 <!-- eslint no-undef: "error" -->
 
 ```javascript
-var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ).factory;
-var filledarrayBy = require( '@stdlib/array/filled-by' );
+var discreteUniform = require( '@stdlib/random/array/discrete-uniform' );
 var sfill = require( '@stdlib/blas/ext/base/sfill' );
 
-var rand = discreteUniform( -100, 100 );
-
-var x = filledarrayBy( 10, 'float32', rand );
+var x = discreteUniform( 10, -100, 100, {
+    'dtype': 'float32'
+});
 console.log( x );
 
 sfill( x.length, 5.0, x, 1 );
@@ -143,6 +142,126 @@ console.log( x );
 
 <!-- /.examples -->
 
+<!-- C interface documentation. -->
+
+* * *
+
+<section class="c">
+
+## C APIs
+
+<!-- Section to include introductory text. Make sure to keep an empty line after the intro `section` element and another before the `/section` close. -->
+
+<section class="intro">
+
+</section>
+
+<!-- /.intro -->
+
+<!-- C usage documentation. -->
+
+<section class="usage">
+
+### Usage
+
+```c
+#include "stdlib/blas/ext/base/sfill.h"
+```
+
+#### c_sfill( N, alpha, \*X, strideX )
+
+Fills a single-precision floating-point strided array `X` with a specified scalar constant `alpha`.
+
+```c
+float x[] = { 1.0f, 2.0f, 3.0f, 4.0f };
+
+c_sfill( 4, 5.0f, x, 1 );
+
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **alpha**: `[in] float` scalar constant.
+-   **X**: `[out] float*` input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+
+```c
+void c_sfill( const CBLAS_INT N, const float alpha, float *X, const CBLAS_INT strideX );
+```
+
+#### c_sfill_ndarray( N, alpha, \*X, strideX, offsetX )
+
+Fills a single-precision floating-point strided array `X` with a specified scalar constant `alpha` using alternative indexing semantics.
+
+```c
+float x[] = { 1.0f, 2.0f, 3.0f, 4.0f };
+
+c_sfill_ndarray( 4, 5.0f, x, 1, 0 );
+```
+
+The function accepts the following arguments:
+
+-   **N**: `[in] CBLAS_INT` number of indexed elements.
+-   **alpha**: `[in] float` scalar constant.
+-   **X**: `[out] float*` input array.
+-   **strideX**: `[in] CBLAS_INT` index increment for `X`.
+-   **offsetX**: `[in] CBLAS_INT` starting index for `X`.
+
+```c
+void c_sfill_ndarray( const CBLAS_INT N, const float alpha, float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX );
+```
+
+</section>
+
+<!-- /.usage -->
+
+<!-- C API usage notes. Make sure to keep an empty line after the `section` element and another before the `/section` close. -->
+
+<section class="notes">
+
+</section>
+
+<!-- /.notes -->
+
+<!-- C API usage examples. -->
+
+<section class="examples">
+
+### Examples
+
+```c
+#include "stdlib/blas/ext/base/sfill.h"
+#include <stdio.h>
+
+int main( void ) {
+    // Create a strided array:
+    float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
+
+    // Specify the number of indexed elements:
+    const int N = 8;
+
+    // Specify a stride:
+    const int strideX = 1;
+
+    // Fill the array:
+    c_sfill( N, 5.0f, x, strideX );
+
+    // Print the result:
+    for ( int i = 0; i < 8; i++ ) {
+        printf( "x[ %i ] = %f\n", i, x[ i ] );
+    }
+}
+```
+
+</section>
+
+<!-- /.examples -->
+
+</section>
+
+<!-- /.c -->
+
 <!-- Section for related `stdlib` packages. Do not manually edit this section, as it is automatically populated. -->
 
 <section class="related">
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.js
index 5a94f4d5bc9..436cdad2cc7 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.js
@@ -21,8 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var uniform = require( '@stdlib/random/base/uniform' ).factory;
-var filledarrayBy = require( '@stdlib/array/filled-by' );
+var uniform = require( '@stdlib/random/array/uniform' );
 var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var pow = require( '@stdlib/math/base/special/pow' );
 var pkg = require( './../package.json' ).name;
@@ -31,7 +30,9 @@ var sfill = require( './../lib/sfill.js' );
 
 // VARIABLES //
 
-var rand = uniform( -10, 10 );
+var options = {
+	'dtype': 'float32'
+};
 
 
 // FUNCTIONS //
@@ -44,7 +45,7 @@ var rand = uniform( -10, 10 );
 * @returns {Function} benchmark function
 */
 function createBenchmark( len ) {
-	var x = filledarrayBy( len, 'float32', rand );
+	var x = uniform( len, -10.0, 10.0, options );
 	return benchmark;
 
 	function benchmark( b ) {
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.native.js
index 8e103f5af93..f9e933e061b 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.native.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.native.js
@@ -22,8 +22,7 @@
 
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
-var uniform = require( '@stdlib/random/base/uniform' ).factory;
-var filledarrayBy = require( '@stdlib/array/filled-by' );
+var uniform = require( '@stdlib/random/array/uniform' );
 var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var pow = require( '@stdlib/math/base/special/pow' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -36,7 +35,9 @@ var sfill = tryRequire( resolve( __dirname, './../lib/sfill.native.js' ) );
 var opts = {
 	'skip': ( sfill instanceof Error )
 };
-var rand = uniform( -10.0, 10.0 );
+var options = {
+	'dtype': 'float32'
+};
 
 
 // FUNCTIONS //
@@ -49,7 +50,7 @@ var rand = uniform( -10.0, 10.0 );
 * @returns {Function} benchmark function
 */
 function createBenchmark( len ) {
-	var x = filledarrayBy( len, 'float32', rand );
+	var x = uniform( len, -10.0, 10.0, options );
 	return benchmark;
 
 	function benchmark( b ) {
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.js
index b72638c6a70..b4be6044299 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.js
@@ -21,8 +21,7 @@
 // MODULES //
 
 var bench = require( '@stdlib/bench' );
-var uniform = require( '@stdlib/random/base/uniform' ).factory;
-var filledarrayBy = require( '@stdlib/array/filled-by' );
+var uniform = require( '@stdlib/random/array/uniform' );
 var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var pow = require( '@stdlib/math/base/special/pow' );
 var pkg = require( './../package.json' ).name;
@@ -31,7 +30,9 @@ var sfill = require( './../lib/ndarray.js' );
 
 // VARIABLES //
 
-var rand = uniform( -10, 10 );
+var options = {
+	'dtype': 'float32'
+};
 
 
 // FUNCTIONS //
@@ -44,7 +45,7 @@ var rand = uniform( -10, 10 );
 * @returns {Function} benchmark function
 */
 function createBenchmark( len ) {
-	var x = filledarrayBy( len, 'float32', rand );
+	var x = uniform( len, -10.0, 10.0, options );
 	return benchmark;
 
 	function benchmark( b ) {
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.native.js
index 975eb4cf5be..55016e52841 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.native.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/benchmark.ndarray.native.js
@@ -22,8 +22,7 @@
 
 var resolve = require( 'path' ).resolve;
 var bench = require( '@stdlib/bench' );
-var uniform = require( '@stdlib/random/base/uniform' ).factory;
-var filledarrayBy = require( '@stdlib/array/filled-by' );
+var uniform = require( '@stdlib/random/array/uniform' );
 var isnanf = require( '@stdlib/math/base/assert/is-nanf' );
 var pow = require( '@stdlib/math/base/special/pow' );
 var tryRequire = require( '@stdlib/utils/try-require' );
@@ -36,7 +35,9 @@ var sfill = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) );
 var opts = {
 	'skip': ( sfill instanceof Error )
 };
-var rand = uniform( -10, 10 );
+var options = {
+	'dtype': 'float32'
+};
 
 
 // FUNCTIONS //
@@ -49,7 +50,7 @@ var rand = uniform( -10, 10 );
 * @returns {Function} benchmark function
 */
 function createBenchmark( len ) {
-	var x = filledarrayBy( len, 'float32', rand );
+	var x = uniform( len, -10.0, 10.0, options );
 	return benchmark;
 
 	function benchmark( b ) {
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/c/benchmark.length.c
index 68d68a5ca1e..9b02b5402ed 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/c/benchmark.length.c
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/benchmark/c/benchmark.length.c
@@ -94,7 +94,7 @@ static float rand_float( void ) {
 * @param len          array length
 * @return elapsed time in seconds
 */
-static double benchmark( int iterations, int len ) {
+static double benchmark1( int iterations, int len ) {
 	double elapsed;
 	float x[ len ];
 	double t;
@@ -118,6 +118,37 @@ static double benchmark( int iterations, int len ) {
 	return elapsed;
 }
 
+/**
+* Runs a benchmark.
+*
+* @param iterations   number of iterations
+* @param len          array length
+* @return elapsed time in seconds
+*/
+static double benchmark2( int iterations, int len ) {
+	double elapsed;
+	float x[ len ];
+	double t;
+	int i;
+
+	for ( i = 0; i < len; i++ ) {
+		x[ i ] = ( rand_float()*200.0f ) - 100.0f;
+	}
+	t = tic();
+	for ( i = 0; i < iterations; i++ ) {
+		c_sfill_ndarray( len, (float)i, x, 1, 0 );
+		if ( x[ 0 ] != x[ 0 ] ) {
+			printf( "should not return NaN\n" );
+			break;
+		}
+	}
+	elapsed = tic() - t;
+	if ( x[ 0 ] != x[ 0 ] ) {
+		printf( "should not return NaN\n" );
+	}
+	return elapsed;
+}
+
 /**
 * Main execution sequence.
 */
@@ -140,7 +171,14 @@ int main( void ) {
 		for ( j = 0; j < REPEATS; j++ ) {
 			count += 1;
 			printf( "# c::%s:len=%d\n", NAME, len );
-			elapsed = benchmark( iter, len );
+			elapsed = benchmark1( iter, len );
+			print_results( iter, elapsed );
+			printf( "ok %d benchmark finished\n", count );
+		}
+		for ( j = 0; j < REPEATS; j++ ) {
+			count += 1;
+			printf( "# c::%s:ndarray:len=%d\n", NAME, len );
+			elapsed = benchmark2( iter, len );
 			print_results( iter, elapsed );
 			printf( "ok %d benchmark finished\n", count );
 		}
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/repl.txt
index 1e9de741b8c..18044061a4b 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/repl.txt
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/repl.txt
@@ -1,10 +1,10 @@
 
-{{alias}}( N, alpha, x, stride )
+{{alias}}( N, alpha, x, strideX )
     Fills a single-precision floating-point strided array with a specified
     scalar value.
 
-    The `N` and `stride` parameters determine which elements in the strided 
-    array are accessed at runtime.
+    The `N` and stride parameters determine which elements in the strided array
+    are accessed at runtime.
 
     Indexing is relative to the first index. To introduce an offset, use typed
     array views.
@@ -22,7 +22,7 @@
     x: Float32Array
         Input array.
 
-    stride: integer
+    strideX: integer
         Index increment for `x`.
 
     Returns
@@ -37,7 +37,7 @@
     > {{alias}}( x.length, 5.0, x, 1 )
     <Float32Array>[ 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0 ]
 
-    // Using `N` and `stride` parameters:
+    // Using `N` and stride parameters:
     > x = new {{alias:@stdlib/array/float32}}( [ -2.0, 1.0, 3.0, -5.0, 4.0, -1.0, -3.0 ] );
     > {{alias}}( 3, 5.0, x, 2 )
     <Float32Array>[ 5.0, 1.0, 5.0, -5.0, 5.0, -1.0, -3.0 ]
@@ -51,13 +51,13 @@
     <Float32Array>[ 1.0, 5.0, 3.0, 5.0, 5.0, 5.0 ]
 
 
-{{alias}}.ndarray( N, alpha, x, stride, offset )
+{{alias}}.ndarray( N, alpha, x, strideX, offsetX )
     Fills a single-precision floating-point strided array with a specified
     scalar value using alternative indexing semantics.
 
     While typed array views mandate a view offset based on the underlying
-    buffer, the `offset` parameter supports indexing semantics based on a
-    starting index.
+    buffer, the offset parameter supports indexing semantics based on a starting
+    index.
 
     Parameters
     ----------
@@ -70,10 +70,10 @@
     x: Float32Array
         Input array.
 
-    stride: integer
+    strideX: integer
         Index increment for `x`.
 
-    offset: integer
+    offsetX: integer
         Starting index of `x`.
 
     Returns
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/types/index.d.ts
index 8b281176efb..33ff7babb6f 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/types/index.d.ts
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/docs/types/index.d.ts
@@ -28,7 +28,7 @@ interface Routine {
 	* @param N - number of indexed elements
 	* @param alpha - constant
 	* @param x - input array
-	* @param stride - stride length
+	* @param strideX - stride length
 	* @returns input array
 	*
 	* @example
@@ -39,7 +39,7 @@ interface Routine {
 	* sfill( x.length, 5.0, x, 1 );
 	* // x => <Float32Array>[ 5.0, 5.0, 5.0, 0.0, 5.0, 5.0, 5.0, 5.0 ]
 	*/
-	( N: number, alpha: number, x: Float32Array, stride: number ): Float32Array;
+	( N: number, alpha: number, x: Float32Array, strideX: number ): Float32Array;
 
 	/**
 	* Fills a single-precision floating-point strided array with a specified scalar value using alternative indexing semantics.
@@ -47,8 +47,8 @@ interface Routine {
 	* @param N - number of indexed elements
 	* @param alpha - constant
 	* @param x - input array
-	* @param stride - stride length
-	* @param offset - starting index
+	* @param strideX - stride length
+	* @param offsetX - starting index
 	* @returns input array
 	*
 	* @example
@@ -59,7 +59,7 @@ interface Routine {
 	* sfill.ndarray( x.length, 5.0, x, 1, 0 );
 	* // x => <Float32Array>[ 5.0, 5.0, 5.0, 0.0, 5.0, 5.0, 5.0, 5.0 ]
 	*/
-	ndarray( N: number, alpha: number, x: Float32Array, stride: number, offset: number ): Float32Array;
+	ndarray( N: number, alpha: number, x: Float32Array, strideX: number, offsetX: number ): Float32Array;
 }
 
 /**
@@ -68,7 +68,7 @@ interface Routine {
 * @param N - number of indexed elements
 * @param alpha - constant
 * @param x - input array
-* @param stride - stride length
+* @param strideX - stride length
 * @returns input array
 *
 * @example
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/c/example.c
index 7a1fee88d1a..1dad47ac586 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/c/example.c
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/c/example.c
@@ -21,13 +21,13 @@
 
 int main( void ) {
 	// Create a strided array:
-	float x[] = { 1.0, -2.0, 3.0, -4.0, 5.0, -6.0, 7.0, -8.0 };
+	float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f };
 
-	// Specify the number of elements:
-	int N = 8;
+	// Specify the number of indexed elements:
+	const int N = 8;
 
 	// Specify a stride:
-	int strideX = 1;
+	const int strideX = 1;
 
 	// Fill the array:
 	c_sfill( N, 5.0f, x, strideX );
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/index.js
index df91e67b06e..55b380ae292 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/index.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/examples/index.js
@@ -18,13 +18,12 @@
 
 'use strict';
 
-var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ).factory;
-var filledarrayBy = require( '@stdlib/array/filled-by' );
+var discreteUniform = require( '@stdlib/random/array/discrete-uniform' );
 var sfill = require( './../lib' );
 
-var rand = discreteUniform( -100, 100 );
-
-var x = filledarrayBy( 10, 'float32', rand );
+var x = discreteUniform( 10, -100, {
+	'dtype': 'float32'
+});
 console.log( x );
 
 sfill( x.length, 5.0, x, 1 );
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/include/stdlib/blas/ext/base/sfill.h b/lib/node_modules/@stdlib/blas/ext/base/sfill/include/stdlib/blas/ext/base/sfill.h
index c367ed21215..d084bc7dfab 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/include/stdlib/blas/ext/base/sfill.h
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/include/stdlib/blas/ext/base/sfill.h
@@ -19,7 +19,7 @@
 #ifndef STDLIB_BLAS_EXT_BASE_SFILL_H
 #define STDLIB_BLAS_EXT_BASE_SFILL_H
 
-#include <stdint.h>
+#include "stdlib/blas/base/shared.h"
 
 /*
 * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler.
@@ -31,7 +31,12 @@ extern "C" {
 /**
 * Fills a single-precision floating-point strided array with a specified scalar constant.
 */
-void c_sfill( const int64_t N, const float alpha, float *X, const int64_t stride );
+void API_SUFFIX(c_sfill)( const CBLAS_INT N, const float alpha, float *X, const CBLAS_INT strideX );
+
+/**
+* Fills a single-precision floating-point strided array with a specified scalar constant using alternative indexing semantics.
+*/
+void API_SUFFIX(c_sfill_ndarray)( const CBLAS_INT N, const float alpha, float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX );
 
 #ifdef __cplusplus
 }
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.js
index aa511a0b7a3..f309b9b5376 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.js
@@ -31,8 +31,8 @@ var M = 8;
 * @param {PositiveInteger} N - number of indexed elements
 * @param {number} alpha - scalar
 * @param {Float32Array} x - input array
-* @param {integer} stride - index increment
-* @param {NonNegativeInteger} offset - starting index
+* @param {integer} strideX - index increment
+* @param {NonNegativeInteger} offsetX - starting index
 * @returns {Float32Array} input array
 *
 * @example
@@ -43,7 +43,7 @@ var M = 8;
 * sfill( 3, 5.0, x, 1, x.length-3 );
 * // x => <Float32Array>[ 1.0, -2.0, 3.0, 5.0, 5.0, 5.0 ]
 */
-function sfill( N, alpha, x, stride, offset ) {
+function sfill( N, alpha, x, strideX, offsetX ) {
 	var ix;
 	var m;
 	var i;
@@ -51,17 +51,17 @@ function sfill( N, alpha, x, stride, offset ) {
 	if ( N <= 0 ) {
 		return x;
 	}
-	ix = offset;
+	ix = offsetX;
 
 	// Use loop unrolling if the stride is equal to `1`...
-	if ( stride === 1 ) {
+	if ( strideX === 1 ) {
 		m = N % M;
 
 		// If we have a remainder, run a clean-up loop...
 		if ( m > 0 ) {
 			for ( i = 0; i < m; i++ ) {
 				x[ ix ] = alpha;
-				ix += stride;
+				ix += strideX;
 			}
 		}
 		if ( N < M ) {
@@ -82,7 +82,7 @@ function sfill( N, alpha, x, stride, offset ) {
 	}
 	for ( i = 0; i < N; i++ ) {
 		x[ ix ] = alpha;
-		ix += stride;
+		ix += strideX;
 	}
 	return x;
 }
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.native.js
index 662a0ea7fb8..ffea0a872d4 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.native.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/ndarray.native.js
@@ -20,9 +20,7 @@
 
 // MODULES //
 
-var minViewBufferIndex = require( '@stdlib/strided/base/min-view-buffer-index' );
-var offsetView = require( '@stdlib/strided/base/offset-view' );
-var addon = require( './sfill.native.js' );
+var addon = require( './../src/addon.node' );
 
 
 // MAIN //
@@ -33,8 +31,8 @@ var addon = require( './sfill.native.js' );
 * @param {PositiveInteger} N - number of indexed elements
 * @param {number} alpha - scalar
 * @param {Float32Array} x - input array
-* @param {integer} stride - index increment
-* @param {NonNegativeInteger} offset - starting index
+* @param {integer} strideX - index increment
+* @param {NonNegativeInteger} offsetX - starting index
 * @returns {Float32Array} input array
 *
 * @example
@@ -45,11 +43,8 @@ var addon = require( './sfill.native.js' );
 * sfill( 3, 5.0, x, 1, x.length-3 );
 * // x => <Float32Array>[ 1.0, -2.0, 3.0, 5.0, 5.0, 5.0 ]
 */
-function sfill( N, alpha, x, stride, offset ) {
-	var view;
-	offset = minViewBufferIndex( N, stride, offset );
-	view = offsetView( x, offset );
-	addon( N, alpha, view, stride );
+function sfill( N, alpha, x, strideX, offsetX ) {
+	addon.ndarray( N, alpha, x, strideX, offsetX );
 	return x;
 }
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.js
index f70cf47c2a1..05616419522 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.js
@@ -32,7 +32,7 @@ var ndarray = require( './ndarray.js' );
 * @param {PositiveInteger} N - number of indexed elements
 * @param {number} alpha - scalar
 * @param {Float32Array} x - input array
-* @param {integer} stride - index increment
+* @param {integer} strideX - index increment
 * @returns {Float32Array} input array
 *
 * @example
@@ -43,8 +43,8 @@ var ndarray = require( './ndarray.js' );
 * sfill( x.length, 5.0, x, 1 );
 * // x => <Float32Array>[ 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0 ]
 */
-function sfill( N, alpha, x, stride ) {
-	return ndarray( N, alpha, x, stride, stride2offset( N, stride ) );
+function sfill( N, alpha, x, strideX ) {
+	return ndarray( N, alpha, x, strideX, stride2offset( N, strideX ) );
 }
 
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.native.js b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.native.js
index d52ad32151a..31e8917a398 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.native.js
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/lib/sfill.native.js
@@ -31,7 +31,7 @@ var addon = require( './../src/addon.node' );
 * @param {PositiveInteger} N - number of indexed elements
 * @param {number} alpha - scalar
 * @param {Float32Array} x - input array
-* @param {integer} stride - index increment
+* @param {integer} strideX - index increment
 * @returns {Float32Array} input array
 *
 * @example
@@ -42,8 +42,8 @@ var addon = require( './../src/addon.node' );
 * sfill( x.length, 5.0, x, 1 );
 * // x => <Float32Array>[ 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0 ]
 */
-function sfill( N, alpha, x, stride ) {
-	addon( N, alpha, x, stride );
+function sfill( N, alpha, x, strideX ) {
+	addon( N, alpha, x, strideX );
 	return x;
 }
 
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/sfill/manifest.json
index 983a02e89b2..9184cfd0377 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/manifest.json
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/manifest.json
@@ -28,16 +28,16 @@
     {
       "task": "build",
       "src": [
-        "./src/sfill.c"
+        "./src/main.c"
       ],
       "include": [
         "./include"
       ],
-      "libraries": [
-        "-lm"
-      ],
+      "libraries": [],
       "libpath": [],
       "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset",
         "@stdlib/napi/export",
         "@stdlib/napi/argv",
         "@stdlib/napi/argv-float",
@@ -48,30 +48,32 @@
     {
       "task": "benchmark",
       "src": [
-        "./src/sfill.c"
+        "./src/main.c"
       ],
       "include": [
         "./include"
       ],
-      "libraries": [
-        "-lm"
-      ],
+      "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     },
     {
       "task": "examples",
       "src": [
-        "./src/sfill.c"
+        "./src/main.c"
       ],
       "include": [
         "./include"
       ],
-      "libraries": [
-        "-lm"
-      ],
+      "libraries": [],
       "libpath": [],
-      "dependencies": []
+      "dependencies": [
+        "@stdlib/blas/base/shared",
+        "@stdlib/strided/base/stride2offset"
+      ]
     }
   ]
 }
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/sfill/src/addon.c
index 88f0f41d3f6..20994ed978c 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/src/addon.c
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/src/addon.c
@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/blas/ext/base/sfill.h"
+#include "stdlib/blas/base/shared.h"
 #include "stdlib/napi/export.h"
 #include "stdlib/napi/argv.h"
 #include "stdlib/napi/argv_float.h"
@@ -35,10 +36,28 @@ static napi_value addon( napi_env env, napi_callback_info info ) {
 	STDLIB_NAPI_ARGV( env, info, argv, argc, 4 );
 	STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 );
 	STDLIB_NAPI_ARGV_FLOAT( env, alpha, argv, 1 );
-	STDLIB_NAPI_ARGV_INT64( env, stride, argv, 3 );
-	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, stride, argv, 2 );
-	c_sfill( N, alpha, X, stride );
+	STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 );
+	API_SUFFIX(c_sfill)( N, alpha, X, strideX );
 	return NULL;
 }
 
-STDLIB_NAPI_MODULE_EXPORT_FCN( addon )
+/**
+* Receives JavaScript callback invocation data.
+*
+* @param env    environment under which the function is invoked
+* @param info   callback data
+* @return       Node-API value
+*/
+static napi_value addon_method( napi_env env, napi_callback_info info ) {
+	STDLIB_NAPI_ARGV( env, info, argv, argc, 5 );
+	STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 );
+	STDLIB_NAPI_ARGV_FLOAT( env, alpha, argv, 1 );
+	STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 );
+	STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 );
+	STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 );
+	API_SUFFIX(c_sfill_ndarray)( N, alpha, X, strideX, offsetX );
+	return NULL;
+}
+
+STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method )
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sfill/src/sfill.c b/lib/node_modules/@stdlib/blas/ext/base/sfill/src/main.c
similarity index 52%
rename from lib/node_modules/@stdlib/blas/ext/base/sfill/src/sfill.c
rename to lib/node_modules/@stdlib/blas/ext/base/sfill/src/main.c
index a1b68dcd95a..6c381394492 100644
--- a/lib/node_modules/@stdlib/blas/ext/base/sfill/src/sfill.c
+++ b/lib/node_modules/@stdlib/blas/ext/base/sfill/src/main.c
@@ -17,7 +17,8 @@
 */
 
 #include "stdlib/blas/ext/base/sfill.h"
-#include <stdint.h>
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/stride2offset.h"
 
 /**
 * Fills a single-precision floating-point strided array with a specified scalar constant.
@@ -25,49 +26,63 @@
 * @param N       number of indexed elements
 * @param alpha   scalar
 * @param X       input array
-* @param stride  index increment
+* @param strideX  index increment
 */
-void c_sfill( const int64_t N, const float alpha, float *X, const int64_t stride ) {
-	int64_t ix;
-	int64_t m;
-	int64_t i;
+
+void API_SUFFIX(c_sfill)( const CBLAS_INT N, const float alpha, float *X, const CBLAS_INT strideX ) {
+	CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX );
+	API_SUFFIX(c_sfill_ndarray)( N, alpha, X, strideX, ox );
+}
+
+/**
+* Fills a single-precision floating-point strided array with a specified scalar constant using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param alpha    scalar
+* @param X        input array
+* @param strideX  index increment
+* @param offsetX  starting index
+*/
+void API_SUFFIX(c_sfill_ndarray)( const CBLAS_INT N, const float alpha, float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) {
+	CBLAS_INT ix;
+	CBLAS_INT m;
+	CBLAS_INT i;
 
 	if ( N <= 0 ) {
 		return;
 	}
+	ix = offsetX;
+
 	// Use loop unrolling if the stride is equal to `1`...
-	if ( stride == 1 ) {
+	if ( strideX == 1 ) {
 		m = N % 8;
 
 		// If we have a remainder, run a clean-up loop...
 		if ( m > 0 ) {
 			for ( i = 0; i < m; i++ ) {
-				X[ i ] = alpha;
+				X[ ix ] = alpha;
+				ix += strideX;
 			}
 		}
 		if ( N < 8 ) {
 			return;
 		}
 		for ( i = m; i < N; i += 8 ) {
-			X[ i ] = alpha;
-			X[ i+1 ] = alpha;
-			X[ i+2 ] = alpha;
-			X[ i+3 ] = alpha;
-			X[ i+4 ] = alpha;
-			X[ i+5 ] = alpha;
-			X[ i+6 ] = alpha;
-			X[ i+7 ] = alpha;
+			X[ ix ] = alpha;
+			X[ ix+1 ] = alpha;
+			X[ ix+2 ] = alpha;
+			X[ ix+3 ] = alpha;
+			X[ ix+4 ] = alpha;
+			X[ ix+5 ] = alpha;
+			X[ ix+6 ] = alpha;
+			X[ ix+7 ] = alpha;
+			ix += 8;
 		}
 		return;
 	}
-	if ( stride < 0 ) {
-		ix = (1-N) * stride;
-	} else {
-		ix = 0;
-	}
 	for ( i = 0; i < N; i++ ) {
 		X[ ix ] = alpha;
-		ix += stride;
+		ix += strideX;
 	}
 	return;
 }

From d9655756945741bdcb75d1a83711c4ca1d639323 Mon Sep 17 00:00:00 2001
From: Gunj Joshi <gunjjoshi8372@gmail.com>
Date: Thu, 19 Sep 2024 03:52:20 +0530
Subject: [PATCH 5/5] build: add a script to find packages with structured meta
 data

PR-URL: https://github.com/stdlib-js/stdlib/pull/2922
Ref: https://github.com/stdlib-js/stdlib/issues/1147
Co-authored-by: Athan Reines <kgryte@gmail.com>
Reviewed-by: Athan Reines <kgryte@gmail.com>
---
 .../math-iter-unary/scripts/data.json         | 149 ++++++++++++++++++
 .../math-iter-unary/scripts/generate_data.js  |  79 ++++++++++
 2 files changed, 228 insertions(+)
 create mode 100644 lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/data.json
 create mode 100644 lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/generate_data.js

diff --git a/lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/data.json b/lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/data.json
new file mode 100644
index 00000000000..86110e6ba6d
--- /dev/null
+++ b/lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/data.json
@@ -0,0 +1,149 @@
+[
+  {
+    "$schema": "math/base@v1.0",
+    "base_alias": "cbrt",
+    "alias": "cbrt",
+    "pkg_desc": "compute the cube root",
+    "desc": "computes the cube root",
+    "short_desc": "cube root",
+    "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": [
+          64,
+          27,
+          0,
+          0,
+          -9,
+          8,
+          -1,
+          125,
+          -10.2,
+          11.3,
+          -12.4,
+          3.5,
+          -1.6,
+          15.7,
+          -16,
+          17.9,
+          -188,
+          19.11,
+          -200,
+          21.15
+        ]
+      }
+    ],
+    "returns": {
+      "desc": "cube root",
+      "type": {
+        "javascript": "number",
+        "jsdoc": "number",
+        "c": "double",
+        "dtype": "float64"
+      }
+    },
+    "keywords": [
+      "cube",
+      "root",
+      "cbrt",
+      "cubic",
+      "power"
+    ],
+    "extra_keywords": [
+      "math.cbrt"
+    ]
+  },
+  {
+    "$schema": "math/base@v1.0",
+    "base_alias": "exp",
+    "alias": "exp",
+    "pkg_desc": "evaluate the natural exponential function",
+    "desc": "evaluates the natural exponential function",
+    "short_desc": "natural exponential 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.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": [
+      "natural",
+      "exponential",
+      "exp",
+      "power"
+    ],
+    "extra_keywords": [
+      "math.exp"
+    ]
+  }
+]
diff --git a/lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/generate_data.js b/lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/generate_data.js
new file mode 100644
index 00000000000..1cada68fe58
--- /dev/null
+++ b/lib/node_modules/@stdlib/_tools/scaffold/math-iter-unary/scripts/generate_data.js
@@ -0,0 +1,79 @@
+#!/usr/bin/env node
+
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2024 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.
+* You may obtain a copy of the License at
+*
+*    http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+// MODULES //
+
+var resolve = require( 'path' ).resolve;
+var findPkgs = require( '@stdlib/_tools/pkgs/find' ).sync;
+var readJSON = require( '@stdlib/fs/read-json' ).sync;
+var writeFile = require( '@stdlib/fs/write-file' ).sync;
+var rootDir = require( '@stdlib/_tools/utils/root-dir' );
+
+
+// VARIABLES //
+
+var ROOT_DIR = rootDir();
+var opts = {
+	'encoding': 'utf8'
+};
+
+
+// MAIN //
+
+/**
+* Main execution sequence.
+*
+* @private
+*/
+function main() {
+	var path;
+	var json;
+	var pkgs;
+	var out;
+	var o;
+	var i;
+
+	// Resolve unary "special" math packages:
+	pkgs = findPkgs({
+		'dir': ROOT_DIR,
+		'pattern': '**/math/base/special/*/package.json'
+	});
+
+	// Filter for package's containing scaffold metadata and which have unary APIs...
+	out = [];
+	for ( i = 0; i < pkgs.length; i++ ) {
+		path = resolve( ROOT_DIR, pkgs[i] );
+		json = readJSON( resolve( path, 'package.json' ), opts );
+		if ( json instanceof Error ) {
+			console.error( 'Encountered an error when attempting to read package: %s. Error: %s.', pkgs[ i ], json.message );
+			continue;
+		}
+		o = json.__stdlib__; // eslint-disable-line no-underscore-dangle
+		if ( o && o.scaffold && o.scaffold.parameters.length === 1 ) {
+			out.push( o.scaffold );
+		}
+	}
+	// Write the metadata to a local file for subsequent consumption by other scripts:
+	writeFile( resolve( __dirname, 'data.json' ), JSON.stringify( out, null, '  ' )+'\n', opts );
+}
+
+main();