refactor: add ndarray implementation

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

Author: ShabiShett07

lib/node_modules/@stdlib/blas/base/zaxpy/src/addon.c

@@ -37,11 +37,32 @@ static napi_value addon( napi_env env, napi_callback_info info ) {
 	STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 );
 	STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 );
 	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 5 );
-	STDLIB_NAPI_ARGV_COMPLEX128( env, za, argv, 1 );
-	STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, ZX, N, strideX, argv, 2 );
-	STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, ZY, N, strideY, argv, 4 );
-	API_SUFFIX(c_zaxpy)( N, za, (void *)ZX, strideX, (void *)ZY, strideY );
+	STDLIB_NAPI_ARGV_COMPLEX128( env, alpha, argv, 1 );
+	STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, X, N, strideX, argv, 2 );
+	STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, Y, N, strideY, argv, 4 );
+	API_SUFFIX(c_zaxpy)( N, alpha, (void *)X, strideX, (void *)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_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_COMPLEX128( env, alpha, argv, 1 );
+	STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, X, N, strideX, argv, 2 );
+	STDLIB_NAPI_ARGV_STRIDED_COMPLEX128ARRAY( env, Y, N, strideY, argv, 5 );
+	API_SUFFIX(c_zaxpy_ndarray)( N, alpha, (void *)X, strideX, offsetX, (void *)Y, strideY, offsetY );
+	return NULL;
+}
+
+STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method )

lib/node_modules/@stdlib/blas/base/zaxpy/src/zaxpy.c

@@ -19,35 +19,20 @@
 #include "stdlib/blas/base/zaxpy.h"
 #include "stdlib/blas/base/shared.h"
 #include "stdlib/complex/float64/ctor.h"
-#include "stdlib/complex/float64/base/mul.h"
-#include "stdlib/complex/float64/base/add.h"
+#include "stdlib/strided/base/stride2offset.h"
 
 /**
 * Scales a double-precision complex floating-point vector by a double-precision complex floating-point constant and adds the result to a double-precision complex floating-point vector.
 *
 * @param N        number of indexed elements
-* @param za       scalar constant
-* @param ZX       first input array
-* @param strideX  ZX stride length
-* @param ZY       second input array
-* @param strideY  ZY stride length
+* @param alpha    scalar constant
+* @param X        input array
+* @param strideX  X stride length
+* @param Y        output array
+* @param strideY  Y stride length
 */
-void API_SUFFIX(c_zaxpy)( const CBLAS_INT N, const stdlib_complex128_t za, void *ZX, const CBLAS_INT strideX, void *ZY, const CBLAS_INT strideY ) {
-	stdlib_complex128_t zx;
-	stdlib_complex128_t zy;
-	CBLAS_INT i;
-
-	uint8_t *ip1 = (uint8_t *)ZX;
-	uint8_t *ip2 = (uint8_t *)ZY;
-	int64_t is1 = 16 * (int64_t)strideX;
-	int64_t is2 = 16 * (int64_t)strideY;
-	if ( N <= 0 || strideX <= 0 || strideY <= 0 ) {
-		return;
-	}
-	for ( i = 0; i < N; i++, ip1 += is1, ip2 += is2 ) {
-		zx = *(stdlib_complex128_t *)ip1;
-		zy = *(stdlib_complex128_t *)ip2;
-		*(stdlib_complex128_t *)ip2 = stdlib_base_complex128_add( stdlib_base_complex128_mul( za, zx ), zy );
-	}
-	return;
+void API_SUFFIX(c_zaxpy)( const CBLAS_INT N, const stdlib_complex128_t alpha, const void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
+	CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX );
+	CBLAS_INT oy = stdlib_strided_stride2offset( N, strideY );
+	API_SUFFIX(c_zaxpy_ndarray)( N, alpha, X, strideX, ox, Y, strideY, oy );
 }

lib/node_modules/@stdlib/blas/base/zaxpy/src/zaxpy_cblas.c

@@ -25,12 +25,12 @@
 * Scales a double-precision complex floating-point vector by a double-precision complex floating-point constant and adds the result to a double-precision complex floating-point vector.
 *
 * @param N        number of indexed elements
-* @param za       scalar constant
-* @param ZX       first input array
-* @param strideX  ZX stride length
-* @param ZY       second input array
-* @param strideY  ZY stride length
+* @param alpha       scalar constant
+* @param X       first input array
+* @param strideX  X stride length
+* @param Y       second input array
+* @param strideY  Y stride length
 */
-void API_SUFFIX(c_zaxpy)( const CBLAS_INT N, const stdlib_complex128_t za, void *ZX, const CBLAS_INT strideX, void *ZY, const CBLAS_INT strideY ) {
-	API_SUFFIX(cblas_zaxpy)( N, za, ZX, strideX, ZY, strideY );
+void API_SUFFIX(c_zaxpy)( const CBLAS_INT N, const stdlib_complex128_t alpha, void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
+	API_SUFFIX(cblas_zaxpy)( N, alpha, X, strideX, Y, strideY );
 }

lib/node_modules/@stdlib/blas/base/zaxpy/src/zaxpy_f.c

@@ -25,12 +25,12 @@
 * Scales a double-precision complex floating-point vector by a double-precision complex floating-point constant and adds the result to a double-precision complex floating-point vector.
 *
 * @param N        number of indexed elements
-* @param za       scalar constant
-* @param ZX       first input array
-* @param strideX  ZX stride length
-* @param ZY       second input array
-* @param strideY  ZY stride length
+* @param alpha    scalar constant
+* @param X        first input array
+* @param strideX  X stride length
+* @param Y        second input array
+* @param strideY  Y stride length
 */
-void API_SUFFIX(c_zaxpy)( const CBLAS_INT N, const stdlib_complex128_t za, void *ZX, const CBLAS_INT strideX, void *ZY, const CBLAS_INT strideY ) {
-	zaxpy( &N, &za, ZX, &strideX, ZY, &strideY );
+void API_SUFFIX(c_zaxpy)( const CBLAS_INT N, const stdlib_complex128_t alpha, void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
+	zaxpy( &N, &alpha, X, &strideX, Y, &strideY );
 }

lib/node_modules/@stdlib/blas/base/zaxpy/src/zaxpy_ndarray.c

@@ -0,0 +1,64 @@
+/**
+* @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.
+*/
+
+#include "stdlib/blas/base/zaxpy.h"
+#include "stdlib/blas/base/scabs1.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/complex/float64/ctor.h"
+#include "stdlib/complex/float64/base/add.h"
+#include "stdlib/complex/float64/base/mul.h"
+#include <stdint.h>
+
+/**
+* Scales a double-precision complex floating-point vector by a double-precision complex floating-point constant and adds the result to a double-precision complex floating-point vector using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param alpha    scalar constant
+* @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_zaxpy_ndarray)( const CBLAS_INT N, const stdlib_complex128_t alpha, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	stdlib_complex128_t x;
+	stdlib_complex128_t y;
+	int64_t is1;
+	int64_t is2;
+	int64_t i;
+
+	if ( N <= 0 ) {
+		return;
+	}
+	if ( c_scabs1( alpha ) == 0.0 ) {
+		return;
+	}
+	stdlib_complex128_t *ip1 = (stdlib_complex128_t *)X;
+	stdlib_complex128_t *ip2 = (stdlib_complex128_t *)Y;
+	is1 = (int64_t)strideX;
+	is2 = (int64_t)strideY;
+	ip1 += offsetX;
+	ip2 += offsetY;
+	for ( i = 0; i < N; i++, ip1 += is1, ip2 += is2 ) {
+		x = *ip1;
+		y = *ip2;
+		*ip2 = stdlib_base_complex128_add( stdlib_base_complex128_mul( alpha, x ), y );
+	}
+	return;
+}