@stdlib/blas-base-cswap 0.2.1 → 0.4.0

package/package.json

@@ -1,7 +1,7 @@
 {
   "name": "@stdlib/blas-base-cswap",
-  "version": "0.2.1",
-  "description": "Interchanges two complex single-precision floating-point vectors.",
+  "version": "0.4.0",
+  "description": "Interchange two complex single-precision floating-point vectors.",
   "license": "Apache-2.0",
   "author": {
     "name": "The Stdlib Authors",
@@ -34,15 +34,19 @@
     "url": "https://github.com/stdlib-js/stdlib/issues"
   },
   "dependencies": {
-    "@stdlib/assert-is-error": "^0.2.1",
-    "@stdlib/napi-argv": "^0.2.1",
-    "@stdlib/napi-argv-int64": "^0.2.1",
-    "@stdlib/napi-argv-strided-complex64array": "^0.2.1",
-    "@stdlib/napi-export": "^0.2.1",
-    "@stdlib/strided-base-reinterpret-complex64": "^0.2.0",
-    "@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.1",
-    "@stdlib/utils-library-manifest": "^0.2.1",
-    "@stdlib/utils-try-require": "^0.2.1"
+    "@stdlib/assert-is-error": "^0.2.2",
+    "@stdlib/blas-base-shared": "^0.2.0",
+    "@stdlib/complex-float32-ctor": "^0.0.2",
+    "@stdlib/napi-argv": "^0.2.2",
+    "@stdlib/napi-argv-int64": "^0.2.2",
+    "@stdlib/napi-argv-strided-complex64array": "^0.2.2",
+    "@stdlib/napi-export": "^0.3.0",
+    "@stdlib/strided-base-min-view-buffer-index": "^0.3.0",
+    "@stdlib/strided-base-reinterpret-complex64": "^0.2.1",
+    "@stdlib/strided-base-stride2offset": "^0.1.0",
+    "@stdlib/utils-define-nonenumerable-read-only-property": "^0.2.2",
+    "@stdlib/utils-library-manifest": "^0.2.3",
+    "@stdlib/utils-try-require": "^0.2.2"
   },
   "devDependencies": {},
   "engines": {
@@ -85,9 +89,6 @@
     "single",
     "float32array"
   ],
-  "__stdlib__": {
-    "wasm": false
-  },
   "funding": {
     "type": "opencollective",
     "url": "https://opencollective.com/stdlib"

package/src/addon.c

@@ -17,6 +17,7 @@
 */
 
 #include "stdlib/blas/base/cswap.h"
+#include "stdlib/blas/base/shared.h"
 #include "stdlib/napi/export.h"
 #include "stdlib/napi/argv.h"
 #include "stdlib/napi/argv_int64.h"
@@ -26,7 +27,6 @@
 /**
 * Receives JavaScript callback invocation data.
 *
-* @private
 * @param env    environment under which the function is invoked
 * @param info   callback data
 * @return       Node-API value
@@ -38,8 +38,28 @@ static napi_value addon( napi_env env, napi_callback_info info ) {
 	STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 4 );
 	STDLIB_NAPI_ARGV_STRIDED_COMPLEX64ARRAY( env, X, N, strideX, argv, 1 );
 	STDLIB_NAPI_ARGV_STRIDED_COMPLEX64ARRAY( env, Y, N, strideY, argv, 3 );
-	c_cswap( N, (void *)X, strideX, (void *)Y, strideY );
+	API_SUFFIX(c_cswap)( N, (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, 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_COMPLEX64ARRAY( env, X, N, strideX, argv, 1 );
+	STDLIB_NAPI_ARGV_STRIDED_COMPLEX64ARRAY( env, Y, N, strideY, argv, 4 );
+	API_SUFFIX(c_cswap_ndarray)( N, (void *)X, strideX, offsetX, (void *)Y, strideY, offsetY );
+	return NULL;
+}
+
+STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method )

package/src/cswap.c

@@ -17,6 +17,8 @@
 */
 
 #include "stdlib/blas/base/cswap.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/strided/base/stride2offset.h"
 
 /**
 * Interchanges two complex single-precision floating-point vectors.
@@ -27,52 +29,8 @@
 * @param Y        second input array
 * @param strideY  Y stride length
 */
-void c_cswap( const int N, void *X, const int strideX, void *Y, const int strideY ) {
-	float *x = (float *)X;
-	float *y = (float *)Y;
-	float tmp;
-	int ix;
-	int iy;
-	int i;
-	int j;
-
-	if ( N <= 0 ) {
-		return;
-	}
-	if ( strideX == 1 && strideY == 1 ) {
-		for ( i = 0; i < N*2; i += 2 ) {
-			tmp = x[ i ];
-			x[ i ] = y[ i ];
-			y[ i ] = tmp;
-
-			j = i + 1;
-			tmp = x[ j ];
-			x[ j ] = y[ j ];
-			y[ j ] = tmp;
-		}
-		return;
-	}
-	if ( strideX < 0 ) {
-		ix = 2 * (1-N) * strideX;
-	} else {
-		ix = 0;
-	}
-	if ( strideY < 0 ) {
-		iy = 2 * (1-N) * strideY;
-	} else {
-		iy = 0;
-	}
-	for ( i = 0; i < N; i++ ) {
-		tmp = x[ ix ];
-		x[ ix ] = y[ iy ];
-		y[ iy ] = tmp;
-
-		tmp = x[ ix+1 ];
-		x[ ix+1 ] = y[ iy+1 ];
-		y[ iy+1 ] = tmp;
-
-		ix += strideX * 2;
-		iy += strideY * 2;
-	}
-	return;
+void API_SUFFIX(c_cswap)( const CBLAS_INT N, 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_cswap_ndarray)( N, X, strideX, ox, Y, strideY, oy );
 }

package/src/cswap.f

@@ -96,4 +96,4 @@ subroutine cswap( N, cx, strideX, cy, strideY )
     end do
   end if
   return
-end subroutine cswap
+end subroutine cswap

package/src/cswap_cblas.c

@@ -18,6 +18,9 @@
 
 #include "stdlib/blas/base/cswap.h"
 #include "stdlib/blas/base/cswap_cblas.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/complex/float32/ctor.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Interchanges two complex single-precision floating-point vectors.
@@ -28,6 +31,26 @@
 * @param Y        second input array
 * @param strideY  Y stride length
 */
-void c_cswap( const int N, void *X, const int strideX, void *Y, const int strideY ) {
-	cblas_cswap( N, X, strideX, Y, strideY );
+void API_SUFFIX(c_cswap)( const CBLAS_INT N, void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
+	API_SUFFIX(cblas_cswap)( N, X, strideX, Y, strideY );
+}
+
+/**
+* Interchanges two complex single-precision floating-point vectors using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param X        first input array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second input array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+*/
+void API_SUFFIX(c_cswap_ndarray)( const CBLAS_INT N, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	stdlib_complex64_t *x = (stdlib_complex64_t *)X;
+	stdlib_complex64_t *y = (stdlib_complex64_t *)Y;
+
+	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_cswap)( N, (void *)x, strideX, (void *)y, strideY );
 }

package/src/cswap_f.c

@@ -18,6 +18,9 @@
 
 #include "stdlib/blas/base/cswap.h"
 #include "stdlib/blas/base/cswap_fortran.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/complex/float32/ctor.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Interchanges two complex single-precision floating-point vectors.
@@ -28,6 +31,26 @@
 * @param Y        second input array
 * @param strideY  Y stride length
 */
-void c_cswap( const int N, void *X, const int strideX, void *Y, const int strideY ) {
+void API_SUFFIX(c_cswap)( const CBLAS_INT N, void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
 	cswap( &N, X, &strideX, Y, &strideY );
 }
+
+/**
+* Interchanges two complex single-precision floating-point vectors using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @param X        first input array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second input array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+*/
+void API_SUFFIX(c_cswap_ndarray)( const CBLAS_INT N, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	stdlib_complex64_t *x = (stdlib_complex64_t *)X;
+	stdlib_complex64_t *y = (stdlib_complex64_t *)Y;
+
+	x += stdlib_strided_min_view_buffer_index( N, strideX, offsetX );
+	y += stdlib_strided_min_view_buffer_index( N, strideY, offsetY );
+	cswap( &N, (void *)x, &strideX, (void *)y, &strideY );
+}

package/src/cswap_ndarray.c

@@ -0,0 +1,63 @@
+/**
+* @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/cswap.h"
+#include "stdlib/blas/base/shared.h"
+
+/**
+* Interchanges two complex single-precision floating-point vectors.
+*
+* @param N        number of indexed elements
+* @param X        first input array
+* @param strideX  X stride length
+* @param offsetX  starting index for X
+* @param Y        second input array
+* @param strideY  Y stride length
+* @param offsetY  starting index for Y
+*/
+void API_SUFFIX(c_cswap_ndarray)( const CBLAS_INT N, void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	float *x = (float *)X;
+	float *y = (float *)Y;
+	float tmp;
+	CBLAS_INT ix;
+	CBLAS_INT iy;
+	CBLAS_INT sx;
+	CBLAS_INT sy;
+	CBLAS_INT i;
+
+	if ( N <= 0 ) {
+		return;
+	}
+	ix = offsetX * 2;
+	iy = offsetY * 2;
+	sx = strideX * 2;
+	sy = strideY * 2;
+	for ( i = 0; i < N; i++ ) {
+		tmp = x[ ix ];
+		x[ ix ] = y[ iy ];
+		y[ iy ] = tmp;
+
+		tmp = x[ ix+1 ];
+		x[ ix+1 ] = y[ iy+1 ];
+		y[ iy+1 ] = tmp;
+
+		ix += sx;
+		iy += sy;
+	}
+	return;
+}

package/include.gypi

@@ -1,70 +0,0 @@
-# @license Apache-2.0
-#
-# Copyright (c) 2020 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.
-
-# A GYP include file for building a Node.js native add-on.
-#
-# Note that nesting variables is required due to how GYP processes a configuration. Any variables defined within a nested 'variables' section is defined in the outer scope. Thus, conditions in the outer variable scope are free to use these variables without running into "variable undefined" errors.
-#
-# Main documentation:
-#
-# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md
-# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md
-#
-# Variable nesting hacks:
-#
-# [3]: https://chromium.googlesource.com/external/skia/gyp/+/master/common_variables.gypi
-# [4]: https://src.chromium.org/viewvc/chrome/trunk/src/build/common.gypi?revision=127004
-{
-  # Define variables to be used throughout the configuration for all targets:
-  'variables': {
-    'variables': {
-      # Host BLAS library (to override -Dblas=<blas>):
-      'blas%': '',
-
-      # Path to BLAS library (to override -Dblas_dir=<path>):
-      'blas_dir%': '',
-    }, # end variables
-
-    # Source directory:
-    'src_dir': './src',
-
-    # Include directories:
-    'include_dirs': [
-      '<@(blas_dir)',
-      '<!@(node -e "var arr = require(\'@stdlib/utils-library-manifest\')(\'./manifest.json\',{\'os\':\'<(OS)\',\'blas\':\'<(blas)\'},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).include; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
-    ],
-
-    # Add-on destination directory:
-    'addon_output_dir': './src',
-
-    # Source files:
-    'src_files': [
-      '<(src_dir)/addon.c',
-      '<!@(node -e "var arr = require(\'@stdlib/utils-library-manifest\')(\'./manifest.json\',{\'os\':\'<(OS)\',\'blas\':\'<(blas)\'},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).src; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
-    ],
-
-    # Library dependencies:
-    'libraries': [
-      '<!@(node -e "var arr = require(\'@stdlib/utils-library-manifest\')(\'./manifest.json\',{\'os\':\'<(OS)\',\'blas\':\'<(blas)\'},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).libraries; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
-    ],
-
-    # Library directories:
-    'library_dirs': [
-      '<@(blas_dir)',
-      '<!@(node -e "var arr = require(\'@stdlib/utils-library-manifest\')(\'./manifest.json\',{\'os\':\'<(OS)\',\'blas\':\'<(blas)\'},{\'basedir\':process.cwd(),\'paths\':\'posix\'}).libpath; for ( var i = 0; i < arr.length; i++ ) { console.log( arr[ i ] ); }")',
-    ],
-  }, # end variables
-}