@stdlib/blas-base-ccopy 0.2.1 → 0.4.0

package/src/ccopy.f

@@ -48,19 +48,19 @@
 ! > * We will gladly answer any questions regarding the software. If a modification is done, however, it is the responsibility of the person who modified the routine to provide support.
 !
 ! @param {integer} N - number of indexed elements
-! @param {Array<complex>} cx - input array
-! @param {integer} strideX - `cx` stride length
-! @param {Array<complex>} cy - output array
-! @param {integer} strideY - `cy` stride length
+! @param {Array<complex>} x - input array
+! @param {integer} strideX - `x` stride length
+! @param {Array<complex>} y - output array
+! @param {integer} strideY - `y` stride length
 !<
-subroutine ccopy( N, cx, strideX, cy, strideY )
+subroutine ccopy( N, x, strideX, y, strideY )
   implicit none
   ! ..
   ! Scalar arguments:
   integer :: strideX, strideY, N
   ! ..
   ! Array arguments:
-  complex :: cx(*), cy(*)
+  complex :: x(*), y(*)
   ! ..
   ! Local scalars:
   integer :: ix, iy, i
@@ -71,7 +71,7 @@ subroutine ccopy( N, cx, strideX, cy, strideY )
   ! ..
   if ( strideX == 1 .AND. strideY == 1 ) then
     do i = 1, N
-      cy( i ) = cx( i )
+      y( i ) = x( i )
     end do
   else
     if ( strideX < 0 ) then
@@ -85,10 +85,10 @@ subroutine ccopy( N, cx, strideX, cy, strideY )
       iy = 1
     end if
     do i = 1, N
-      cy( iy ) = cx( ix )
+      y( iy ) = x( ix )
       ix = ix + strideX
       iy = iy + strideY
     end do
   end if
   return
-end subroutine ccopy
+end subroutine ccopy

package/src/ccopy_cblas.c

@@ -18,6 +18,9 @@
 
 #include "stdlib/blas/base/ccopy.h"
 #include "stdlib/blas/base/ccopy_cblas.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/complex/float32/ctor.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Copies values from one complex single-precision floating-point vector to another complex single-precision floating-point vector.
@@ -28,6 +31,27 @@
 * @param Y        output array
 * @param strideY  Y stride length
 */
-void c_ccopy( const int N, const void *X, const int strideX, void *Y, const int strideY ) {
-	cblas_ccopy( N, X, strideX, Y, strideY );
+void API_SUFFIX(c_ccopy)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
+	API_SUFFIX(cblas_ccopy)( N, X, strideX, Y, strideY );
+}
+
+/**
+* Copies values from one complex single-precision floating-point vector to another complex single-precision floating-point vector using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @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_ccopy_ndarray)( const CBLAS_INT N, const 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_ccopy)( N, (void *)x, strideX, (void *)y, strideY );
 }

package/src/ccopy_f.c

@@ -18,6 +18,9 @@
 
 #include "stdlib/blas/base/ccopy.h"
 #include "stdlib/blas/base/ccopy_fortran.h"
+#include "stdlib/blas/base/shared.h"
+#include "stdlib/complex/float32/ctor.h"
+#include "stdlib/strided/base/min_view_buffer_index.h"
 
 /**
 * Copies values from one complex single-precision floating-point vector to another complex single-precision floating-point vector.
@@ -28,6 +31,27 @@
 * @param Y        output array
 * @param strideY  Y stride length
 */
-void c_ccopy( const int N, const void *X, const int strideX, void *Y, const int strideY ) {
+void API_SUFFIX(c_ccopy)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, void *Y, const CBLAS_INT strideY ) {
 	ccopy( &N, X, &strideX, Y, &strideY );
 }
+
+/**
+* Copies values from one complex single-precision floating-point vector to another complex single-precision floating-point vector using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @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_ccopy_ndarray)( const CBLAS_INT N, const 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 );
+
+	ccopy( &N, (void *)x, &strideX, (void *)y, &strideY );
+}

package/src/ccopy_ndarray.c

@@ -0,0 +1,56 @@
+/**
+* @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/ccopy.h"
+#include "stdlib/blas/base/shared.h"
+
+/**
+* Copies values from one complex single-precision floating-point vector to another complex single-precision floating-point vector using alternative indexing semantics.
+*
+* @param N        number of indexed elements
+* @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_ccopy_ndarray)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) {
+	const float *x = (float *)X;
+	float *y = (float *)Y;
+	CBLAS_INT ix;
+	CBLAS_INT iy;
+	CBLAS_INT sx;
+	CBLAS_INT sy;
+	CBLAS_INT i;
+
+	if ( N <= 0 ) {
+		return;
+	}
+	sx = strideX * 2;
+	sy = strideY * 2;
+	ix = offsetX * 2;
+	iy = offsetY * 2;
+	for ( i = 0; i < N; i++ ) {
+		y[ iy ] = x[ ix ];
+		y[ iy+1 ] = x[ ix+1 ];
+		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
-}