diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/README.md b/lib/node_modules/@stdlib/blas/base/cdotc/README.md new file mode 100644 index 000000000000..2f7d6947b6ae --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/README.md @@ -0,0 +1,335 @@ + + +# cdotc + +> Calculate the dot product `x^H * y` of `x` and `y`. + +
+ +The [dot product][dot-product] (or scalar product) is defined as + +
+ + + +
+ +## Usage + +```javascript +var cdotc = require( '@stdlib/blas/base/cdotc' ); +``` + +#### cdotc( N, x, strideX, y, strideY ) + +Calculates the dot product `x^H * y` of `x` and `y`. + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); + +var out = cdotc( x.length, x, 1, y, 1 ); +// returns [ 54.0, -80.0 ] +``` + +The function has the following parameters: + +- **N**: number of indexed elements. +- **x**: input [`Complex64Array`][@stdlib/array/complex64]. +- **strideX**: index increment for `x`. +- **y**: input [`Complex64Array`][@stdlib/array/complex64]. +- **strideY**: index increment for `y`. + +The `N` and stride parameters determine which elements in the strided arrays are accessed at runtime. For example, to calculate the dot product of every other value in `x` and the first `N` elements of `y` in reverse order, + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ -1.0, -9.0, 2.0, -8.0 ] ); +var y = new Complex64Array( [ -5.0, 1.0, -6.0, 7.0 ] ); + +var out = cdotc( x.length, x, 1, y, -1 ); +// returns [ -75.0, -99.0 ] +``` + +#### cdotc.ndarray( N, x, strideX, offsetX, y, strideY, offsetY ) + +Calculates the dot product `x^H * y` of `x` and `y` using alternative indexing semantics. + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); + +var out = cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); +// returns [ 54.0, -80.0 ] +``` + +The function has the following additional parameters: + +- **offsetX**: starting index for `x`. +- **offsetY**: starting index for `y`. + +While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example, to calculate the dot product of every other value in `x` starting from the second value with the last 2 elements in `y` in reverse order + +```javascript +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + +var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); + +var out = cdotc.ndarray( x.length, x, 1, 0, y, -1, y.length-1 ); +// returns [ -55.0, 23.0 ] +``` + +
+ + + +
+ +## Notes + +- If `N <= 0`, both functions return complex64 with real and imaginary as `0.0`. +- `cdotc()` corresponds to the [BLAS][blas] level 1 function [`cdotc`][cdotc]. + +
+ + + +
+ +## Examples + + + +```javascript +var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var cdotc = require( '@stdlib/blas/base/cdotc' ); + +function rand() { + return new Complex64( discreteUniform( 0, 10 ), discreteUniform( 1, 5 ) ); +} + +var x = filledarrayBy( 10, 'complex64', rand ); +console.log( x.toString() ); + +var y = filledarrayBy( 10, 'complex64', rand ); +console.log( y.toString() ); + +// Perform dot product of x and y +var out = cdotc.ndarray( x.length, x, 1, 0, y, -1, y.length-1 ); +console.log( out ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/blas/base/cdotc.h" +``` + +#### c_cdotc( N, \*X, strideX, \*Y, strideY ) + +Calculates the dot product `x^H * y` of `x` and `y`. + +```c +#include "stdlib/complex/float32/ctor.h" + +float x[] = { 7.0f, -8.0f, -1.0f, -9.0f }; +float y[] = { 6.0f, -6.0f, -9.0f, 5.0f }; + +stdlib_complex64_t out = c_cdotc( 2, (void *)x, 1, (void *)y, 1 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **X**: `[in] void*` first input array. +- **strideX**: `[in] CBLAS_INT` index increment for `X`. +- **Y**: `[in] void*` second input array. +- **strideY**: `[in] CBLAS_INT` index increment for `Y`. + +```c +stdlib_complex64_t c_cdotc( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const void *Y, const CBLAS_INT strideY ); +``` + +#### c_cdotc_ndarray( N, \*X, strideX, offsetX, \*Y, strideY, offsetY ) + +Calculates the dot product `x^H * y` of `x` and `y` using alternative indexing semantics. + +```c +#include "stdlib/complex/float32/ctor.h" + +float x[] = { 7.0f, -8.0f, -1.0f, -9.0f }; +float y[] = { 6.0f, -6.0f, -9.0f, 5.0f }; + +stdlib_complex64_t out = c_cdotc_ndarray( 2, (void *)x, 1, 0, (void *)y, 1, 0 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **X**: `[in] void*` first input array. +- **strideX**: `[in] CBLAS_INT` index increment for `X`. +- **offsetX**: `[in] CBLAS_INT` starting index for `X`. +- **Y**: `[in] void*` second input array. +- **strideY**: `[in] CBLAS_INT` index increment for `Y`. +- **offsetY**: `[in] CBLAS_INT` starting index for `Y`. + +```c +stdlib_complex64_t c_cdotc_ndarray( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/blas/base/cdotc.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/reim.h" +#include + +int main( void ) { + // Create strided arrays of interleaved real and imaginary components: + float x[] = { 7.0f, -8.0f, -1.0f, -9.0f }; + float y[] = { 6.0f, -6.0f, -9.0f, 5.0f }; + + // Specify the number of elements: + const int N = 2; + + // Specify stride lengths: + const int strideX = 1; + const int strideY = 1; + + // Compute the dot product: + stdlib_complex64_t dot = c_cdotc( N, (void *)x, strideX, (void *)y, strideY ); + + // Print the result: + float re; + float im; + stdlib_complex64_reim( dot, &re, &im ); + printf( "cdotc( x, y ) = %f + %fi\n", re, im ); + + // Compute the dot product using alternative indexing semantics: + dot = c_cdotc_ndarray( N, (void *)x, -strideX, 1, (void *)y, strideY, 0 ); + + // Print the result: + stdlib_complex64_reim( dot, &re, &im ); + printf( "cdotc( x, y ) = %f + %fi\n", re, im ); +} +``` + +
+ + + +
+ + + + + +
+ +* * * + +## See Also + +- [`@stdlib/blas/base/cdotu`][@stdlib/blas/base/cdotu]: calculate the dot product of two single-precision complex floating-point vectors. +- [`@stdlib/blas/base/zdotc`][@stdlib/blas/base/zdotc]: calculate the dot product of two double-precision complex floating-point vectors, conjugating the first vector. + +
+ + + + + +
+ +[dot-product]: https://en.wikipedia.org/wiki/Dot_product + +[blas]: http://www.netlib.org/blas + +[cdotc]: http://www.netlib.org/lapack/explore-html/df/d28/group__single__blas__level1.html + +[@stdlib/array/complex64]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/complex64 + +[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray + + + +[@stdlib/blas/base/cdotu]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/base/cdotu + +[@stdlib/blas/base/zdotc]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/base/zdotc + + + +
+ + diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.js new file mode 100644 index 000000000000..0d269158e792 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.js @@ -0,0 +1,109 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var realf = require( '@stdlib/complex/float32/real' ); +var pkg = require( './../package.json' ).name; +var cdotc = require( './../lib/cdotc.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x; + var y; + + x = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + y = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var out; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + out = cdotc( x.length, x, 1, y, 1 ); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+':len='+len, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.native.js new file mode 100644 index 000000000000..83dbd451257e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.native.js @@ -0,0 +1,109 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var realf = require( '@stdlib/complex/float32/real' ); +var pkg = require( './../package.json' ).name; +var cdotc = require( './../lib/cdotc.native.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x; + var y; + + x = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + y = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var out; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + out = cdotc( x.length, x, 1, y, 1 ); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+':native:len='+len, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.ndarray.js new file mode 100644 index 000000000000..31c9151eaf0f --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.ndarray.js @@ -0,0 +1,109 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var realf = require( '@stdlib/complex/float32/real' ); +var pkg = require( './../package.json' ).name; +var cdotc = require( './../lib/ndarray.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x; + var y; + + x = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + y = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var out; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + out = cdotc( x.length, x, 1, 0, y, 1, 0 ); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+':len='+len, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.ndarray.native.js new file mode 100644 index 000000000000..9282b76b99bd --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/benchmark/benchmark.ndarray.native.js @@ -0,0 +1,109 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var isnanf = require( '@stdlib/math/base/assert/is-nanf' ); +var realf = require( '@stdlib/complex/float32/real' ); +var pkg = require( './../package.json' ).name; +var cdotc = require( './../lib/ndarray.native.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float32' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x; + var y; + + x = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + y = new Complex64Array( uniform( len*2, -100.0, 100.0, options ) ); + + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var out; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + out = cdotc( x.length, x, 1, 0, y, 1, 0 ); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnanf( realf( out ) ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+':ndarray:native:len='+len, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/binding.gyp b/lib/node_modules/@stdlib/blas/base/cdotc/binding.gyp new file mode 100644 index 000000000000..60dce9d0b31a --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/binding.gyp @@ -0,0 +1,265 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 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` file for building a Node.js native add-on. +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # List of files to include in this file: + 'includes': [ + './include.gypi', + ], + + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Target name should match the add-on export name: + 'addon_target_name%': 'addon', + + # Fortran compiler (to override -Dfortran_compiler=): + 'fortran_compiler%': 'gfortran', + + # Fortran compiler flags: + 'fflags': [ + # Specify the Fortran standard to which a program is expected to conform: + '-std=f95', + + # Indicate that the layout is free-form source code: + '-ffree-form', + + # Aggressive optimization: + '-O3', + + # Enable commonly used warning options: + '-Wall', + + # Warn if source code contains problematic language features: + '-Wextra', + + # Warn if a procedure is called without an explicit interface: + '-Wimplicit-interface', + + # Do not transform names of entities specified in Fortran source files by appending underscores (i.e., don't mangle names, thus allowing easier usage in C wrappers): + '-fno-underscoring', + + # Warn if source code contains Fortran 95 extensions and C-language constructs: + '-pedantic', + + # Compile but do not link (output is an object file): + '-c', + ], + + # Set variables based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + { + # Define the object file suffix: + 'obj': 'obj', + }, + { + # Define the object file suffix: + 'obj': 'o', + } + ], # end condition (OS=="win") + ], # end conditions + }, # end variables + + # Define compile targets: + 'targets': [ + + # Target to generate an add-on: + { + # The target name should match the add-on export name: + 'target_name': '<(addon_target_name)', + + # Define dependencies: + 'dependencies': [], + + # Define directories which contain relevant include headers: + 'include_dirs': [ + # Local include directory: + '<@(include_dirs)', + ], + + # List of source files: + 'sources': [ + '<@(src_files)', + ], + + # Settings which should be applied when a target's object files are used as linker input: + 'link_settings': { + # Define libraries: + 'libraries': [ + '<@(libraries)', + ], + + # Define library directories: + 'library_dirs': [ + '<@(library_dirs)', + ], + }, + + # C/C++ compiler flags: + 'cflags': [ + # Enable commonly used warning options: + '-Wall', + + # Aggressive optimization: + '-O3', + ], + + # C specific compiler flags: + 'cflags_c': [ + # Specify the C standard to which a program is expected to conform: + '-std=c99', + ], + + # C++ specific compiler flags: + 'cflags_cpp': [ + # Specify the C++ standard to which a program is expected to conform: + '-std=c++11', + ], + + # Linker flags: + 'ldflags': [], + + # Apply conditions based on the host OS: + 'conditions': [ + [ + 'OS=="mac"', + { + # Linker flags: + 'ldflags': [ + '-undefined dynamic_lookup', + '-Wl,-no-pie', + '-Wl,-search_paths_first', + ], + }, + ], # end condition (OS=="mac") + [ + 'OS!="win"', + { + # C/C++ flags: + 'cflags': [ + # Generate platform-independent code: + '-fPIC', + ], + }, + ], # end condition (OS!="win") + ], # end conditions + + # Define custom build actions for particular inputs: + 'rules': [ + { + # Define a rule for processing Fortran files: + 'extension': 'f', + + # Define the pathnames to be used as inputs when performing processing: + 'inputs': [ + # Full path of the current input: + '<(RULE_INPUT_PATH)' + ], + + # Define the outputs produced during processing: + 'outputs': [ + # Store an output object file in a directory for placing intermediate results (only accessible within a single target): + '<(INTERMEDIATE_DIR)/<(RULE_INPUT_ROOT).<(obj)' + ], + + # Define the rule for compiling Fortran based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + + # Rule to compile Fortran on Windows: + { + 'rule_name': 'compile_fortran_windows', + 'message': 'Compiling Fortran file <(RULE_INPUT_PATH) on Windows...', + + 'process_outputs_as_sources': 0, + + # Define the command-line invocation: + 'action': [ + '<(fortran_compiler)', + '<@(fflags)', + '<@(_inputs)', + '-o', + '<@(_outputs)', + ], + }, + + # Rule to compile Fortran on non-Windows: + { + 'rule_name': 'compile_fortran_linux', + 'message': 'Compiling Fortran file <(RULE_INPUT_PATH) on Linux...', + + 'process_outputs_as_sources': 1, + + # Define the command-line invocation: + 'action': [ + '<(fortran_compiler)', + '<@(fflags)', + '-fPIC', # generate platform-independent code + '<@(_inputs)', + '-o', + '<@(_outputs)', + ], + } + ], # end condition (OS=="win") + ], # end conditions + }, # end rule (extension=="f") + ], # end rules + }, # end target <(addon_target_name) + + # Target to copy a generated add-on to a standard location: + { + 'target_name': 'copy_addon', + + # Declare that the output of this target is not linked: + 'type': 'none', + + # Define dependencies: + 'dependencies': [ + # Require that the add-on be generated before building this target: + '<(addon_target_name)', + ], + + # Define a list of actions: + 'actions': [ + { + 'action_name': 'copy_addon', + 'message': 'Copying addon...', + + # Explicitly list the inputs in the command-line invocation below: + 'inputs': [], + + # Declare the expected outputs: + 'outputs': [ + '<(addon_output_dir)/<(addon_target_name).node', + ], + + # Define the command-line invocation: + 'action': [ + 'cp', + '<(PRODUCT_DIR)/<(addon_target_name).node', + '<(addon_output_dir)/<(addon_target_name).node', + ], + }, + ], # end actions + }, # end target copy_addon + ], # end targets +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/docs/repl.txt b/lib/node_modules/@stdlib/blas/base/cdotc/docs/repl.txt new file mode 100644 index 000000000000..cade653ad7c5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/docs/repl.txt @@ -0,0 +1,107 @@ + +{{alias}}( N, x, strideX, y, strideY ) + Computes the dot product of two single-precision complex vectors. + + The `N` and stride parameters determine which elements in the strided arrays + are accessed at runtime. + + Indexing is relative to the first index. To introduce an offset, use a typed + array view. + + If `N <= 0`, the function returns `0`. + + Parameters + ---------- + N: integer + Number of indexed elements. + + x: Complex64Array + First input array. + + strideX: integer + Index increment for `x`. + + y: Complex64Array + Second input array. + + strideY: integer + Index increment for `y`. + + Returns + ------- + z: Complex64 + The dot product. + + Examples + -------- + // Standard usage: + > var x = new {{alias:@stdlib/array/complex64}}( [ 7.0, -8.0, -1.0, -9.0 ] ); + > var y = new {{alias:@stdlib/array/complex64}}( [ 6.0, -6.0, -9.0, 5.0 ] ); + > var z = {{alias}}( x.length, x, 1, y, 1 ) + [ 54.0, -80.0 ] + + // Strides: + > x = new {{alias:@stdlib/array/complex64}}( [ 7.0, -8.0, -4.0, -7.0, -1.0, -9.0 ] ); + > y = new {{alias:@stdlib/array/complex64}}( [ 6.0, -6.0, -9.0, 5.0, 7.0, -6.0 ] ); + > z = {{alias}}( 2, x, 2, y, 1 ) + [ 54.0, -80.0 ] + + +{{alias}}.ndarray( N, x, strideX, offsetX, y, strideY, offsetY ) + Computes the dot product of two single-precision complex vectors + using alternative indexing semantics. + + While typed array views mandate a view offset based on the underlying + buffer, the offset parameters support indexing based on a starting index. + + Parameters + ---------- + N: integer + Number of indexed elements. + + x: Complex64Array + First input array. + + strideX: integer + Index increment for `x`. + + offsetX: integer + Starting index for `x`. + + y: Complex64Array + Second input array. + + strideY: integer + Index increment for `y`. + + offsetY: integer + Starting index for `y`. + + Returns + ------- + z: Complex64 + The dot product. + + Examples + -------- + // Standard usage: + > var x = new {{alias:@stdlib/array/complex64}}( [ 7.0, -8.0, -1.0, -9.0 ] ); + > var y = new {{alias:@stdlib/array/complex64}}( [ 6.0, -6.0, -9.0, 5.0 ] ); + > var z = {{alias}}.ndarray( x.length, x, 1, 0, y, 1, 0 ) + [ 54.0, -80.0 ] + + // Strides: + > x = new {{alias:@stdlib/array/complex64}}( [ 7.0, -8.0, -4.0, -7.0, -1.0, -9.0 ] ); + > y = new {{alias:@stdlib/array/complex64}}( [ 6.0, -6.0, -9.0, 5.0, 7.0, -6.0 ] ); + > z = {{alias}}.ndarray( 2, x, 2, 0, y, 1, 0 ) + [ 54.0, -80.0 ] + + // Using offset indices: + > x = new {{alias:@stdlib/array/complex64}}( [ 7.0, -8.0, -4.0, -7.0, -1.0, -9.0 ] ); + > y = new {{alias:@stdlib/array/complex64}}( [ 6.0, -6.0, -9.0, 5.0, 7.0, -6.0 ] ); + > z = {{alias}}.ndarray( 2, x, -2, x.length-1, y, 1, 1 ) + [ 61.0, -72.0 ] + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/base/cdotc/docs/types/index.d.ts new file mode 100644 index 000000000000..77cbf7d4c7ea --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/docs/types/index.d.ts @@ -0,0 +1,107 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 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. +*/ + +import { Complex64 } from '@stdlib/types/complex'; +import { Complex64Array } from '@stdlib/types/array'; + + +// TypeScript Version: 4.1 + +/** +* Interface describing `cdotc`. +*/ +interface Rzine { + /** + * Computes the dot product `x^H * y` of `x` and `y`. + * + * @param N - number of indexed elements + * @param x - first input complex array + * @param strideX - `x` stride length + * @param y - second input complex array + * @param strideY - `y` stride length + * @returns dot product + * + * @example + * var Complex64Array = require( '@stdlib/array/complex64' ); + * + * var x = new Complex64Array( [ 7, -8, -1, -9 ] ); + * var y = new Complex64Array( [ 6, -6, -9, 5 ] ); + * + * var z = cdotc( x.length, x, 1, y, 1 ); + * // returns [ 54.0, -80.0 ] + */ + ( N: number, x: Complex64Array, strideX: number, y: Complex64Array, strideY: number ): Complex64; + + /** + * Computes the dot product `x^H * y` of `x` and `y` 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` + * @returns dot product + * + * @example + * var Complex64Array = require( '@stdlib/array/complex64' ); + * + * var x = new Complex64Array( [ 7, -8, -1, -9 ] ); + * var y = new Complex64Array( [ 6, -6, -9, 5 ] ); + * + * var z = cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); + * // returns [ 54.0, -80.0 ] + */ + ndarray( N: number, x: Complex64Array, strideX: number, offsetX: number, y: Complex64Array, strideY: number, offsetY: number ): Complex64; +} + +/** +* Computes the dot product `x^H * y` of `x` and `y`. +* +* @param N - number of indexed elements +* @param x - first input array +* @param strideX - `x` stride length +* @param y - second input array +* @param strideY - `y` stride length +* @returns dot product +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* +* var x = new Complex64Array( [ 7, -8, -1, -9 ] ); +* var y = new Complex64Array( [ 6, -6, -9, 5 ] ); +* +* var z = cdotc( x.length, x, 1, y, 1 ); +* // returns [ 54.0, -80.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* +* var x = new Complex64Array( [ 7, -8, -1, -9 ] ); +* var y = new Complex64Array( [ 6, -6, -9, 5 ] ); +* +* var z = cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); +* // returns [ 54.0, -80.0 ] +*/ +declare var cdotc: Rzine; + + +// EXPORTS // + +export = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/docs/types/test.ts b/lib/node_modules/@stdlib/blas/base/cdotc/docs/types/test.ts new file mode 100644 index 000000000000..6aa5fb9b0cbf --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/docs/types/test.ts @@ -0,0 +1,249 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 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. +*/ + +import Complex64Array = require( '@stdlib/array/complex64' ); +import cdotc = require( './index' ); + + +// TESTS // + +// The function returns a Complex64... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc( x.length, x, 1, y, 1 ); // $ExpectType Complex64 +} + +// The compiler throws an error if the function is provided a first argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc( '10', x, 1, y, 1 ); // $ExpectError + cdotc( true, x, 1, y, 1 ); // $ExpectError + cdotc( false, x, 1, y, 1 ); // $ExpectError + cdotc( null, x, 1, y, 1 ); // $ExpectError + cdotc( undefined, x, 1, y, 1 ); // $ExpectError + cdotc( [], x, 1, y, 1 ); // $ExpectError + cdotc( {}, x, 1, y, 1 ); // $ExpectError + cdotc( ( x: number ): number => x, x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a second argument which is not a Complex64Array.. +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc( x.length, 10, 1, y, 1 ); // $ExpectError + cdotc( x.length, '10', 1, y, 1 ); // $ExpectError + cdotc( x.length, true, 1, y, 1 ); // $ExpectError + cdotc( x.length, false, 1, y, 1 ); // $ExpectError + cdotc( x.length, null, 1, y, 1 ); // $ExpectError + cdotc( x.length, undefined, 1, y, 1 ); // $ExpectError + cdotc( x.length, [], 1, y, 1 ); // $ExpectError + cdotc( x.length, {}, 1, y, 1 ); // $ExpectError + cdotc( x.length, ( x: number ): number => x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a third argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc( x.length, x, '10', y, 1 ); // $ExpectError + cdotc( x.length, x, true, y, 1 ); // $ExpectError + cdotc( x.length, x, false, y, 1 ); // $ExpectError + cdotc( x.length, x, null, y, 1 ); // $ExpectError + cdotc( x.length, x, undefined, y, 1 ); // $ExpectError + cdotc( x.length, x, [], y, 1 ); // $ExpectError + cdotc( x.length, x, {}, y, 1 ); // $ExpectError + cdotc( x.length, x, ( x: number ): number => x, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fourth argument which is not a Complex64Array... +{ + const x = new Complex64Array( 10 ); + + cdotc( x.length, x, 1, 10, 1 ); // $ExpectError + cdotc( x.length, x, 1, '10', 1 ); // $ExpectError + cdotc( x.length, x, 1, true, 1 ); // $ExpectError + cdotc( x.length, x, 1, false, 1 ); // $ExpectError + cdotc( x.length, x, 1, null, 1 ); // $ExpectError + cdotc( x.length, x, 1, undefined, 1 ); // $ExpectError + cdotc( x.length, x, 1, [], 1 ); // $ExpectError + cdotc( x.length, x, 1, {}, 1 ); // $ExpectError + cdotc( x.length, x, 1, ( x: number ): number => x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fifth argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc( x.length, x, 1, y, '10' ); // $ExpectError + cdotc( x.length, x, 1, y, true ); // $ExpectError + cdotc( x.length, x, 1, y, false ); // $ExpectError + cdotc( x.length, x, 1, y, null ); // $ExpectError + cdotc( x.length, x, 1, y, undefined ); // $ExpectError + cdotc( x.length, x, 1, y, [] ); // $ExpectError + cdotc( x.length, x, 1, y, {} ); // $ExpectError + cdotc( x.length, x, 1, y, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc(); // $ExpectError + cdotc( x.length ); // $ExpectError + cdotc( x.length, x ); // $ExpectError + cdotc( x.length, x, 1 ); // $ExpectError + cdotc( x.length, x, 1, y ); // $ExpectError + cdotc( x.length, x, 1, y, 1, 10 ); // $ExpectError +} + +// Attached to main export is an `ndarray` method which returns a Complex64... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); // $ExpectType Complex64 +} + +// The compiler throws an error if the `ndarray` method is provided a first argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( '10', x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( true, x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( false, x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( null, x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( undefined, x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( [], x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( {}, x, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( ( x: number ): number => x, x, 1, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a second argument which is not a Complex64Array... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, 10, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, '10', 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, true, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, false, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, null, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, undefined, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, [], 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, {}, 1, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, ( x: number ): number => x, 1, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a third argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, x, '10', 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, true, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, false, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, null, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, undefined, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, [], 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, {}, 0, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, ( x: number ): number => x, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fourth argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, x, 1, '10', y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, true, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, false, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, null, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, undefined, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, [], y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, {}, y, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, ( x: number ): number => x, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a Complex64Array... +{ + const x = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, x, 1, 0, 10, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, '10', 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, true, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, false, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, null, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, undefined, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, [], 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, {}, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, ( x: number ): number => x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, x, 1, 0, y, '10', 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, true, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, false, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, null, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, undefined, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, [], 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, {}, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray( x.length, x, 1, 0, y, 1, '10' ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, true ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, false ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, null ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, undefined ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, [] ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, {} ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments... +{ + const x = new Complex64Array( 10 ); + const y = new Complex64Array( 10 ); + + cdotc.ndarray(); // $ExpectError + cdotc.ndarray( x.length ); // $ExpectError + cdotc.ndarray( x.length, x ); // $ExpectError + cdotc.ndarray( x.length, x, 1 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1 ); // $ExpectError + cdotc.ndarray( x.length, x, 1, 0, y, 1, 0, 10 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/examples/c/Makefile b/lib/node_modules/@stdlib/blas/base/cdotc/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := example.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled examples. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/examples/c/example.c b/lib/node_modules/@stdlib/blas/base/cdotc/examples/c/example.c new file mode 100644 index 000000000000..929a3ce20d3e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/examples/c/example.c @@ -0,0 +1,44 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/cdotc.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/real.h" +#include "stdlib/complex/float32/imag.h" +#include + +int main( void ) { + stdlib_complex64_t x[ 2 ]; + stdlib_complex64_t y[ 2 ]; + stdlib_complex64_t dot; + + // Create complex arrays: + x[ 0 ] = stdlib_complex64( 7.0f, -8.0f ); + x[ 1 ] = stdlib_complex64( -1.0f, -9.0f ); + + y[ 0 ] = stdlib_complex64( 6.0f, -6.0f ); + y[ 1 ] = stdlib_complex64( -9.0f, 5.0f ); + + // Compute the dot product: + dot = c_cdotc( 2, (void *)x, 1, (void *)y, 1 ); + + // Print the result: + printf( "cdotc( x, y ) = %f + %fi\n", stdlib_complex64_real( dot ), stdlib_complex64_imag( dot ) ); + + return 0; +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/examples/index.js b/lib/node_modules/@stdlib/blas/base/cdotc/examples/index.js new file mode 100644 index 000000000000..9851062d3758 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/examples/index.js @@ -0,0 +1,41 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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'; + +var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var cdotc = require( './../lib' ); + +function rand() { + return new Complex64( discreteUniform( 0, 10 ), discreteUniform( 1, 5 ) ); +} + +var x = filledarrayBy( 10, 'complex64', rand ); +console.log( x.toString() ); + +var y = filledarrayBy( 10, 'complex64', rand ); +console.log( y.toString() ); + +// Perform dot product of x and y +var z = cdotc( x.length, x, 1, y, 1 ); +console.log( z ); + +z = cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); +console.log( z ); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/include.gypi b/lib/node_modules/@stdlib/blas/base/cdotc/include.gypi new file mode 100644 index 000000000000..dcb556d250e8 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/include.gypi @@ -0,0 +1,70 @@ +# @license Apache-2.0 +# +# Copyright (c) 2026 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%': '', + + # Path to BLAS library (to override -Dblas_dir=): + 'blas_dir%': '', + }, # end variables + + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + '<@(blas_dir)', + '[ 54.0, -80.0 ] +*/ +function cdotc( N, x, strideX, y, strideY ) { + var ix; + var iy; + + ix = stride2offset( N, strideX ); + iy = stride2offset( N, strideY ); + return ndarray( N, x, strideX, ix, y, strideY, iy ); +} + + +// EXPORTS // + +module.exports = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/lib/cdotc.native.js b/lib/node_modules/@stdlib/blas/base/cdotc/lib/cdotc.native.js new file mode 100644 index 000000000000..4bf6c3e44e7a --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/lib/cdotc.native.js @@ -0,0 +1,61 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 reinterpret = require( '@stdlib/strided/base/reinterpret-complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Computes the dot product `x^H * y` of `x` and `y`. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {Complex64Array} x - first input array +* @param {integer} strideX - `x` stride length +* @param {Complex64Array} y - second input array +* @param {integer} strideY - `y` stride length +* @returns {Complex64} dot product +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var cdotc = require( '@stdlib/blas/base/cdotc' ); +* +* var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +* var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); +* +* var z = cdotc( x.length, x, 1, y, 1 ); +* // returns [ 54.0, -80.0 ] +*/ +function cdotc( N, x, strideX, y, strideY ) { + var viewX = reinterpret( x, 0 ); + var viewY = reinterpret( y, 0 ); + var v = addon( N, viewX, strideX, viewY, strideY ); + return new Complex64( v.re, v.im ); +} + + +// EXPORTS // + +module.exports = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/lib/index.js b/lib/node_modules/@stdlib/blas/base/cdotc/lib/index.js new file mode 100644 index 000000000000..0043d2da0022 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/lib/index.js @@ -0,0 +1,71 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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'; + +/** +* BLAS level 1 routine to compute the dot product `x^H * y` of `x` and `y`. +* +* @module @stdlib/blas/base/cdotc +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var cdotc = require( '@stdlib/blas/base/cdotc' ); +* +* var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +* var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); +* var z = cdotc( x.length, x, 1, y, 1 ); +* // returns [ 54.0, -80.0 ] +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var cdotc = require( '@stdlib/blas/base/cdotc' ); +* +* var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +* var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); +* +* var z = cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); +* // returns [ 54.0, -80.0 ] +*/ + +// MODULES // + +var join = require( 'path' ).join; +var tryRequire = require( '@stdlib/utils/try-require' ); +var isError = require( '@stdlib/assert/is-error' ); +var main = require( './main.js' ); + + +// MAIN // + +var cdotc; +var tmp = tryRequire( join( __dirname, './native.js' ) ); +if ( isError( tmp ) ) { + cdotc = main; +} else { + cdotc = tmp; +} + + +// EXPORTS // + +module.exports = cdotc; + +// exports: { "ndarray": "cdotc.ndarray" } diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/lib/main.js b/lib/node_modules/@stdlib/blas/base/cdotc/lib/main.js new file mode 100644 index 000000000000..57740ab97494 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/lib/main.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var cdotc = require( './cdotc.js' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +setReadOnly( cdotc, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/lib/native.js b/lib/node_modules/@stdlib/blas/base/cdotc/lib/native.js new file mode 100644 index 000000000000..0702e321c187 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/lib/native.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var cdotc = require( './cdotc.native.js' ); +var ndarray = require( './ndarray.native.js' ); + + +// MAIN // + +setReadOnly( cdotc, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/lib/ndarray.js b/lib/node_modules/@stdlib/blas/base/cdotc/lib/ndarray.js new file mode 100644 index 000000000000..5d1eee66169e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/lib/ndarray.js @@ -0,0 +1,89 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 reinterpret = require( '@stdlib/strided/base/reinterpret-complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); + + +// MAIN // + +/** +* Computes the dot product `x^H * y` of `x` and `y`. +* +* @param {integer} N - number of indexed elements +* @param {Complex64Array} x - first input array +* @param {integer} strideX - `x` stride length +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @param {Complex64Array} y - second input array +* @param {integer} strideY - `y` stride length +* @param {NonNegativeInteger} offsetY - starting index for `y` +* @returns {Complex64} dot product +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* +* var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +* var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); +* +* var z = cdotc( x.length, x, 1, 0, y, 1, 0 ); +* // returns [ 54.0, -80.0 ] +*/ +function cdotc( N, x, strideX, offsetX, y, strideY, offsetY ) { + var viewX; + var viewY; + var tr; + var ti; + var ix; + var iy; + var sx; + var sy; + var i; + + if ( N <= 0 ) { + return new Complex64( 0.0, 0.0 ); + } + viewX = reinterpret( x, 0 ); + viewY = reinterpret( y, 0 ); + + ix = offsetX * 2; + iy = offsetY * 2; + sx = strideX * 2; + sy = strideY * 2; + + tr = 0.0; + ti = 0.0; + for ( i = 0; i < N; i++ ) { + tr += ( viewX[ ix ] * viewY[ iy ] ) + + ( viewX[ ix + 1 ] * viewY[ iy + 1 ] ); + ti += ( viewX[ ix ] * viewY[ iy + 1 ] ) - + ( viewX[ ix + 1 ] * viewY[ iy ] ); + ix += sx; + iy += sy; + } + return new Complex64( tr, ti ); +} + + +// EXPORTS // + +module.exports = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/base/cdotc/lib/ndarray.native.js new file mode 100644 index 000000000000..a6023526c218 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/lib/ndarray.native.js @@ -0,0 +1,63 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 reinterpret = require( '@stdlib/strided/base/reinterpret-complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Computes the dot product `x^H * y` of `x` and `y` using alternative indexing semantics. +* +* @param {integer} N - number of indexed elements +* @param {Complex64Array} x - first input array +* @param {integer} strideX - `x` stride length +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @param {Complex64Array} y - second input array +* @param {integer} strideY - `y` stride length +* @param {NonNegativeInteger} offsetY - starting index for `y` +* @returns {Complex64} dot product +* +* @example +* var Complex64Array = require( '@stdlib/array/complex64' ); +* var Complex64 = require( '@stdlib/complex/float32/ctor' ); +* var cdotc = require( '@stdlib/blas/base/cdotc' ); +* +* var x = new Complex64Array( [ 7.0, -8.0, -1.0, -9.0 ] ); +* var y = new Complex64Array( [ 6.0, -6.0, -9.0, 5.0 ] ); +* +* var z = cdotc.ndarray( x.length, x, 1, 0, y, 1, 0 ); +* // returns [ 54.0, -80.0 ] +*/ +function cdotc( N, x, strideX, offsetX, y, strideY, offsetY ) { + var viewX = reinterpret( x, 0 ); + var viewY = reinterpret( y, 0 ); + var v = addon.ndarray( N, viewX, strideX, offsetX, viewY, strideY, offsetY ); + return new Complex64( v.re, v.im ); +} + + +// EXPORTS // + +module.exports = cdotc; diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/manifest.json b/lib/node_modules/@stdlib/blas/base/cdotc/manifest.json new file mode 100644 index 000000000000..fdedb20eca6f --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/manifest.json @@ -0,0 +1,523 @@ +{ + "options": { + "task": "build", + "os": "linux", + "blas": "", + "wasm": false + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true + }, + { + "field": "include", + "resolve": true, + "relative": true + }, + { + "field": "libraries", + "resolve": false, + "relative": false + }, + { + "field": "libpath", + "resolve": true, + "relative": false + } + ], + "confs": [ + { + "task": "build", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/napi/create-complex-like", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + }, + { + "task": "benchmark", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + }, + { + "task": "examples", + "os": "linux", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + + { + "task": "build", + "os": "linux", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/napi/create-complex-like", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + { + "task": "benchmark", + "os": "linux", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + { + "task": "examples", + "os": "linux", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + + { + "task": "build", + "os": "mac", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/napi/create-complex-like", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + }, + { + "task": "benchmark", + "os": "mac", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + }, + { + "task": "examples", + "os": "mac", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + + { + "task": "build", + "os": "mac", + "blas": "apple_accelerate_framework", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lblas" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/napi/create-complex-like", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + { + "task": "benchmark", + "os": "mac", + "blas": "apple_accelerate_framework", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lblas" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + { + "task": "examples", + "os": "mac", + "blas": "apple_accelerate_framework", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lblas" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + + { + "task": "build", + "os": "mac", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/napi/create-complex-like", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + { + "task": "benchmark", + "os": "mac", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + { + "task": "examples", + "os": "mac", + "blas": "openblas", + "wasm": false, + "src": [ + "./src/cdotc_cblas.c" + ], + "include": [ + "./include" + ], + "libraries": [ + "-lopenblas", + "-lpthread" + ], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/min-view-buffer-index", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + + { + "task": "build", + "os": "win", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-complex64array", + "@stdlib/napi/create-complex-like", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + }, + { + "task": "benchmark", + "os": "win", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + }, + { + "task": "examples", + "os": "win", + "blas": "", + "wasm": false, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj", + "@stdlib/complex/float32/real", + "@stdlib/complex/float32/imag" + ] + }, + + { + "task": "build", + "os": "", + "blas": "", + "wasm": true, + "src": [ + "./src/cdotc.c", + "./src/cdotc_ndarray.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/complex/float32/ctor", + "@stdlib/complex/float32/base/mul", + "@stdlib/complex/float32/base/add", + "@stdlib/complex/float32/conj" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/package.json b/lib/node_modules/@stdlib/blas/base/cdotc/package.json new file mode 100644 index 000000000000..dacf8f6025fc --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/package.json @@ -0,0 +1,86 @@ +{ + "name": "@stdlib/blas/base/cdotc", + "version": "0.0.0", + "description": "Calculate the dot product of two single-precision complex vectors, conjugating the first vector.", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "browser": "./lib/main.js", + "gypfile": true, + "directories": { + "benchmark": "./benchmark", + "doc": "./docs", + "example": "./examples", + "include": "./include", + "lib": "./lib", + "src": "./src", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "mathematics", + "math", + "blas", + "level 1", + "cdotc", + "dot", + "product", + "dot product", + "scalar", + "scalar product", + "inner", + "inner product", + "linear", + "algebra", + "subroutines", + "vector", + "array", + "ndarray", + "complex", + "complex64", + "float32", + "float", + "single", + "float32array" + ], + "__stdlib__": { + "wasm": false + } +} \ No newline at end of file diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/Makefile b/lib/node_modules/@stdlib/blas/base/cdotc/src/Makefile new file mode 100644 index 000000000000..2caf905cedbe --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 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. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + + +# RULES # + +#/ +# Removes generated files for building an add-on. +# +# @example +# make clean-addon +#/ +clean-addon: + $(QUIET) -rm -f *.o *.node + +.PHONY: clean-addon + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: clean-addon + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/addon.c b/lib/node_modules/@stdlib/blas/base/cdotc/src/addon.c new file mode 100644 index 000000000000..2187a0e53483 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/addon.c @@ -0,0 +1,84 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/cdotc.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/napi/argv_strided_complex64array.h" +#include "stdlib/napi/create_complex_like.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/real.h" +#include "stdlib/complex/float32/imag.h" +#include + +/** +* 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( napi_env env, napi_callback_info info ) { + stdlib_complex64_t v; + CBLAS_INT ix; + CBLAS_INT iy; + + STDLIB_NAPI_ARGV( env, info, argv, argc, 5 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 2 ); + 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 ); + + ix = stdlib_strided_stride2offset( N, strideX ); + iy = stdlib_strided_stride2offset( N, strideY ); + v = API_SUFFIX(c_cdotc_ndarray)( N, (void *)X, strideX, ix, (void *)Y, strideY, iy ); + STDLIB_NAPI_CREATE_COMPLEX_LIKE( env, stdlib_complex64_real( v ), stdlib_complex64_imag( v ), out ); + + return out; +} + +/** +* 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_complex64_t v; + + 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 ); + + v = API_SUFFIX(c_cdotc_ndarray)( N, (void *)X, strideX, offsetX, (void *)Y, strideY, offsetY ); + STDLIB_NAPI_CREATE_COMPLEX_LIKE( env, stdlib_complex64_real( v ), stdlib_complex64_imag( v ), out ); + + return out; +} + +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ) diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc.c b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc.c new file mode 100644 index 000000000000..1014dd49993c --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc.c @@ -0,0 +1,37 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/cdotc.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" + +/** +* Computes the dot product `x^H * y` of two single-precision complex floating-point vectors. +* +* @param N number of indexed elements +* @param X first input array +* @param strideX X stride length +* @param Y second input array +* @param strideY Y stride length +* @return dot product +*/ +stdlib_complex64_t API_SUFFIX(c_cdotc)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const void *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_cdotc_ndarray)( N, X, strideX, ox, Y, strideY, oy ); +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc.f b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc.f new file mode 100644 index 000000000000..408d9b94a9e5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc.f @@ -0,0 +1,104 @@ +!> +! @license Apache-2.0 +! +! Copyright (c) 2026 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. +!< + +!> Computes the dot product `x^H * y` of two single-precision complex vectors. +! +! ## Notes +! +! * Modified version of reference BLAS level1 routine (version 3.12.0). Updated to "free form" Fortran 95. +! +! ## Authors +! +! * Univ. of Tennessee +! * Univ. of California Berkeley +! * Univ. of Colorado Denver +! * NAG Ltd. +! +! ## History +! +! * Jack Dongarra, linpack, 3/11/78. +! +! - modified 12/3/93, array(1) declarations changed to array(*) +! +! ## License +! +! From : +! +! > The reference BLAS is a freely-available software package. It is available from netlib via anonymous ftp and the World Wide Web. Thus, it can be included in commercial software packages (and has been). We only ask that proper credit be given to the authors. +! > +! > Like all software, it is copyrighted. It is not trademarked, but we do ask the following: +! > +! > * If you modify the source for these routines we ask that you change the name of the routine and comment the changes made to the original. +! > +! > * 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} cx - first array +! @param {integer} strideX - `cx` stride length +! @param {Array} cy - second array +! @param {integer} strideY - `cy` stride length +! @returns {complex} the dot product +!< +complex function cdotc( N, cx, strideX, cy, strideY ) + implicit none + ! .. + ! Scalar arguments: + integer :: strideX, strideY, N + ! .. + ! Array arguments: + complex, intent(in) :: cx(*), cy(*) + ! .. + ! Local scalars: + complex :: ctemp + integer :: ix, iy, i + ! .. + ! Intrinsic functions: + intrinsic conjg + ! .. + ctemp = (0.0, 0.0) + cdotc = (0.0, 0.0) + ! .. + if ( N <= 0 ) then + return + end if + ! .. + ! If both strides are equal to `1`, use unrolled loops... + if ( strideX == 1 .AND. strideY == 1 ) then + do i = 1, N + ctemp = ctemp + conjg( cx( i ) ) * cy( i ) + end do + else + if ( strideX < 0 ) then + ix = ((1-N)*strideX) + 1 + else + ix = 1 + endif + if ( strideY < 0 ) then + iy = ((1-N)*strideY) + 1 + else + iy = 1 + endif + do i = 1, N + ctemp = ctemp + conjg( cx( ix ) ) * cy( iy ) + ix = ix + strideX + iy = iy + strideY + end do + endif + cdotc = ctemp + return +end function cdotc diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_cblas.c b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_cblas.c new file mode 100644 index 000000000000..9c395b1a5a43 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_cblas.c @@ -0,0 +1,62 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/cdotc.h" +#include "stdlib/blas/base/cdotc_cblas.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/min_view_buffer_index.h" +#include "stdlib/complex/float32/ctor.h" + +/** +* Computes the dot product `x^H * y` of two single-precision complex floating-point vectors. +* +* @param N number of indexed elements +* @param X first input array +* @param strideX X stride length +* @param Y second input array +* @param strideY Y stride length +* @return dot product +*/ +stdlib_complex64_t API_SUFFIX(c_cdotc)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const void *Y, const CBLAS_INT strideY ) { + stdlib_complex64_t dot; + cblas_cdotc_sub( N, X, strideX, Y, strideY, &dot ); + return dot; +} + +/** +* Computes the dot product `x^H * y` of two single-precision complex 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 +*/ +stdlib_complex64_t API_SUFFIX(c_cdotc_ndarray)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) { + stdlib_complex64_t dot; + stdlib_complex64_t *ip1 = (stdlib_complex64_t *)X; + stdlib_complex64_t *ip2 = (stdlib_complex64_t *)Y; + + ip1 += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); + ip2 += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); + cblas_cdotc_sub( N, ip1, strideX, ip2, strideY, &dot ); + return dot; +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_f.c b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_f.c new file mode 100644 index 000000000000..69d3647099d0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_f.c @@ -0,0 +1,62 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/cdotc.h" +#include "stdlib/blas/base/cdotc_fortran.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/min_view_buffer_index.h" +#include "stdlib/complex/float32/ctor.h" + +/** +* Computes the dot product `x^H * y` of two single-precision complex floating-point vectors. +* +* @param N number of indexed elements +* @param X first input array +* @param strideX X stride length +* @param Y second input array +* @param strideY Y stride length +* @return dot product +*/ +stdlib_complex64_t API_SUFFIX(c_cdotc)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const void *Y, const CBLAS_INT strideY ) { + stdlib_complex64_t dot; + cdotcsub( &N, X, &strideX, Y, &strideY, &dot ); + return dot; +} + +/** +* Computes the dot product `x^H * y` of two single-precision complex 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 +*/ +stdlib_complex64_t API_SUFFIX(c_cdotc_ndarray)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) { + stdlib_complex64_t dot; + stdlib_complex64_t *ip1 = (stdlib_complex64_t *)X; + stdlib_complex64_t *ip2 = (stdlib_complex64_t *)Y; + + ip1 += stdlib_strided_min_view_buffer_index( N, strideX, offsetX ); + ip2 += stdlib_strided_min_view_buffer_index( N, strideY, offsetY ); + cdotcsub( &N, ip1, &strideX, ip2, &strideY, &dot ); + return dot; +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_ndarray.c b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_ndarray.c new file mode 100644 index 000000000000..f9253e4db8a6 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotc_ndarray.c @@ -0,0 +1,65 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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/cdotc.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/complex/float32/ctor.h" +#include "stdlib/complex/float32/base/mul.h" +#include "stdlib/complex/float32/base/add.h" +#include "stdlib/complex/float32/conj.h" +#include + +/** +* Computes the dot product `x^H * y` of two single-precision complex 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 +* @return dot product +*/ +stdlib_complex64_t API_SUFFIX(c_cdotc_ndarray)( const CBLAS_INT N, const void *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const void *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) { + stdlib_complex64_t dot; + stdlib_complex64_t conjX; + stdlib_complex64_t xy; + CBLAS_INT ix; + CBLAS_INT iy; + CBLAS_INT i; + + dot = stdlib_complex64_from_float32( 0.0f ); + if ( N <= 0 ) { + return dot; + } + ix = offsetX; + iy = offsetY; + + stdlib_complex64_t *ip1 = (stdlib_complex64_t *)X; + stdlib_complex64_t *ip2 = (stdlib_complex64_t *)Y; + + for ( i = 0; i < N; i++ ) { + conjX = stdlib_complex64_conj( ip1[ ix ] ); + xy = stdlib_base_complex64_mul( conjX, ip2[ iy ] ); + dot = stdlib_base_complex64_add( dot, xy ); + ix += strideX; + iy += strideY; + } + return dot; +} diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotcsub.f b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotcsub.f new file mode 100644 index 000000000000..9439db50e2c7 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/src/cdotcsub.f @@ -0,0 +1,49 @@ +!> +! @license Apache-2.0 +! +! Copyright (c) 2026 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. +!< + +!> Wraps `cdotc` as a subroutine. +! +! @param {integer} N - number of indexed elements +! @param {Array} cx - first array +! @param {integer} strideX - `cx` stride length +! @param {Array} cy - second array +! @param {integer} strideY - `cy` stride length +! @param {complex} dot - output variable reference +!< +subroutine cdotcsub( N, cx, strideX, cy, strideY, dot ) + implicit none + ! .. + ! External functions: + interface + complex function cdotc( N, cx, strideX, cy, strideY ) + complex, intent(in) :: cx(*), cy(*) + integer :: strideX, strideY, N + end function cdotc + end interface + ! .. + ! Scalar arguments: + integer :: strideX, strideY, N + complex :: dot + ! .. + ! Array arguments: + complex, intent(in) :: cx(*), cy(*) + ! .. + ! Compute the dot product: + dot = cdotc( N, cx, strideX, cy, strideY ) + return +end subroutine cdotcsub diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/test/test.cdotc.js b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.cdotc.js new file mode 100644 index 000000000000..ca53b9381440 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.cdotc.js @@ -0,0 +1,280 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var cdotc = require( './../lib/cdotc.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.deepEqual( typeof cdotc, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 5', function test( t ) { + t.deepEqual( cdotc.length, 5, 'returns expected value' ); + t.end(); +}); + +tape( 'the function calculates the dot product of vectors `x` and `y`', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 0.7, // 0 + -0.8, // 0 + -0.4, // 1 + -0.7 // 1 + ]); + y = new Complex64Array([ + 0.6, // 0 + -0.6, // 0 + -0.9, // 1 + 0.5 // 1 + ]); + expected = new Complex64( 0.91, -0.77 ); + + dot = cdotc( x.length, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns Complex( 0.0, 0.0 )', function test( t ) { + var expected; + var dot; + var x; + var y; + + expected = new Complex64( 0.0, 0.0 ); + x = new Complex64Array( [-0.1, -0.9, 0.2, -0.8 ] ); + y = new Complex64Array( [ 0.7, -0.6, 0.1, -0.5 ] ); + + dot = cdotc( 0, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + dot = cdotc( -4, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, 2, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `x` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + -1, // 1 + -9, // 1 + 7, // 0 + -8, // 0 + -4, + -7, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, -1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `y` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, + 5, + 7, // 1 + -6, // 1 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, y, 2 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `y` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 7, // 1 + -6, // 1 + 6, // 0 + -6, // 0 + -9, + 5, + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, y, -1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports negative strides', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 1 + -8, // 1 + -4, // 0 + -7, // 0 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, + 5, + 7, // 0 + -6, // 0 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, -1, y, -2 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports complex access patterns', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, // 0 + 5, // 0 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( -55, 23 ); + + dot = cdotc( 2, x, 2, y, -1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/test/test.cdotc.native.js b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.cdotc.native.js new file mode 100644 index 000000000000..0cb6f31ae216 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.cdotc.native.js @@ -0,0 +1,288 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var IS_BROWSER = require( '@stdlib/assert/is-browser' ); +var cdotc = require( './../lib/cdotc.native.js' ); + + +// VARIABLES // + +var opts = { + 'skip': IS_BROWSER +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.deepEqual( typeof cdotc, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 5', opts, function test( t ) { + t.deepEqual( cdotc.length, 5, 'returns expected value' ); + t.end(); +}); + +tape( 'the function calculates the dot product of vectors `x` and `y`', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 0.7, // 0 + -0.8, // 0 + -0.4, // 1 + -0.7 // 1 + ]); + y = new Complex64Array([ + 0.6, // 0 + -0.6, // 0 + -0.9, // 1 + 0.5 // 1 + ]); + expected = new Complex64( 0.91, -0.77 ); + + dot = cdotc( x.length, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns Complex( 0.0, 0.0 )', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + expected = new Complex64( 0.0, 0.0 ); + x = new Complex64Array( [ -0.1, -0.9, 0.2, -0.8 ] ); + y = new Complex64Array( [ 0.7, -0.6, 0.1, -0.5 ] ); + + dot = cdotc( 0, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + dot = cdotc( -4, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, 2, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `x` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + -1, // 1 + -9, // 1 + 7, // 0 + -8, // 0 + -4, + -7, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, -1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `y` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, + 5, + 7, // 1 + -6, // 1 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, y, 2 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `y` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 7, // 1 + -6, // 1 + 6, // 0 + -6, // 0 + -9, + 5, + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, y, -1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports negative strides', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 1 + -8, // 1 + -4, // 0 + -7, // 0 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, + 5, + 7, // 0 + -6, // 0 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, -1, y, -2 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports complex access patterns', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, // 0 + 5, // 0 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( -55, 23 ); + + dot = cdotc( 2, x, 2, y, -1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/test/test.js b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.js new file mode 100644 index 000000000000..d5cbd0c19aeb --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.js @@ -0,0 +1,82 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var proxyquire = require( 'proxyquire' ); +var IS_BROWSER = require( '@stdlib/assert/is-browser' ); +var cdotc = require( './../lib' ); + + +// VARIABLES // + +var opts = { + 'skip': IS_BROWSER +}; + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof cdotc, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'attached to the main export is a method providing an ndarray interface', function test( t ) { + t.strictEqual( typeof cdotc.ndarray, 'function', 'method is a function' ); + t.end(); +}); + +tape( 'if a native implementation is available, the main export is the native implementation', opts, function test( t ) { + var cdotc = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( cdotc, mock, 'returns expected value' ); + t.end(); + + function tryRequire() { + return mock; + } + + function mock() { + // Mock... + } +}); + +tape( 'if a native implementation is not available, the main export is a JavaScript implementation', opts, function test( t ) { + var cdotc; + var main; + + main = require( './../lib/cdotc.js' ); + + cdotc = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( cdotc, main, 'returns expected value' ); + t.end(); + + function tryRequire() { + return new Error( 'Cannot find module' ); + } +}); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.ndarray.js new file mode 100644 index 000000000000..df5427d3d651 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.ndarray.js @@ -0,0 +1,336 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var cdotc = require( './../lib/ndarray.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.deepEqual( typeof cdotc, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', function test( t ) { + t.deepEqual( cdotc.length, 7, 'returns expected value' ); + t.end(); +}); + +tape( 'the function calculates the dot product of complex vectors `x` and `y`', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 0.7, // 0 + -0.8, // 0 + -0.4, // 1 + -0.7 // 1 + ]); + y = new Complex64Array([ + 0.6, // 0 + -0.6, // 0 + -0.9, // 1 + 0.5 // 1 + ]); + expected = new Complex64( 0.91, -0.77 ); + + dot = cdotc( x.length, x, 1, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns Complex( 0.0, 0.0 )', function test( t ) { + var expected; + var dot; + var x; + var y; + + expected = new Complex64( 0.0, 0.0 ); + x = new Complex64Array( [ -0.1, -0.9, 0.2, -0.8 ] ); + y = new Complex64Array( [ 0.7, -0.6, 0.1, -0.5 ] ); + + dot = cdotc( 0, x, 1, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + dot = cdotc( -4, x, 1, y, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, 2, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `x` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + -1, // 1 + -9, // 1 + 7, // 0 + -8, // 0 + -4, + -7, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, -1, 1, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` offset', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 0.0, + 0.0, + 1.0, // 0 + 2.0, // 0 + 3.0, // 1 + 4.0 // 1 + ]); + y = new Complex64Array([ + -5.0, // 0 + 1.0, // 0 + -6.0, // 1 + 7.0 // 1 + ]); + expected = new Complex64( 7.0, 56.0 ); + + dot = cdotc( 2, x, 1, 1, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, + 5, + 7, // 1 + -6, // 1 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, 0, y, 2, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `y` stride', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 7, // 1 + -6, // 1 + 6, // 0 + -6, // 0 + -9, + 5, + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, 0, y, -1, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `y` offset', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 1.0, // 0 + 2.0, // 0 + 3.0, // 1 + 4.0 // 1 + ]); + y = new Complex64Array([ + 0.0, + 0.0, + -5.0, // 0 + 1.0, // 0 + -6.0, // 1 + 7.0 // 1 + ]); + expected = new Complex64( 7.0, 56.0 ); + + dot = cdotc( 2, x, 1, 0, y, 1, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports negative strides', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 1 + -8, // 1 + -4, // 0 + -7, // 0 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, + 5, + 7, // 0 + -6, // 0 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, -1, 1, y, -2, 2 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports complex access patterns', function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, // 0 + 5, // 0 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( -55, 23 ); + + dot = cdotc( 2, x, 2, 0, y, -1, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/base/cdotc/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.ndarray.native.js new file mode 100644 index 000000000000..4ca8e5a6be1b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/base/cdotc/test/test.ndarray.native.js @@ -0,0 +1,344 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 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 tape = require( 'tape' ); +var Complex64Array = require( '@stdlib/array/complex64' ); +var Complex64 = require( '@stdlib/complex/float32/ctor' ); +var IS_BROWSER = require( '@stdlib/assert/is-browser' ); +var cdotc = require( './../lib/ndarray.native.js' ); + + +// VARIABLES // + +var opts = { + 'skip': IS_BROWSER +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.deepEqual( typeof cdotc, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', opts, function test( t ) { + t.deepEqual( cdotc.length, 7, 'returns expected value' ); + t.end(); +}); + +tape( 'the function calculates the dot product of complex vectors `x` and `y`', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 0.7, // 0 + -0.8, // 0 + -0.4, // 1 + -0.7 // 1 + ]); + y = new Complex64Array([ + 0.6, // 0 + -0.6, // 0 + -0.9, // 1 + 0.5 // 1 + ]); + expected = new Complex64( 0.91, -0.77 ); + + dot = cdotc( x.length, x, 1, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns Complex( 0.0, 0.0 )', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + expected = new Complex64( 0.0, 0.0 ); + x = new Complex64Array( [ -0.1, -0.9, 0.2, -0.8 ] ); + y = new Complex64Array( [ 0.7, -0.6, 0.1, -0.5 ] ); + + dot = cdotc( 0, x, 1, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + dot = cdotc( -4, x, 1, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, 2, 0, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `x` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + -1, // 1 + -9, // 1 + 7, // 0 + -8, // 0 + -4, + -7, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, // 1 + 5, // 1 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( 54, -80 ); + dot = cdotc( 2, x, -1, 1, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `x` offset', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 0.0, + 0.0, + 1.0, // 0 + 2.0, // 0 + 3.0, // 1 + 4.0 // 1 + ]); + y = new Complex64Array([ + -5.0, // 0 + 1.0, // 0 + -6.0, // 1 + 7.0 // 1 + ]); + expected = new Complex64( 7.0, 56.0 ); + + dot = cdotc( 2, x, 1, 1, y, 1, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 0 + -6, // 0 + -9, + 5, + 7, // 1 + -6, // 1 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, 0, y, 2, 0 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a negative `y` stride', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, // 1 + -7, // 1 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 7, // 1 + -6, // 1 + 6, // 0 + -6, // 0 + -9, + 5, + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, 1, 0, y, -1, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a `y` offset', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 1.0, // 0 + 2.0, // 0 + 3.0, // 1 + 4.0 // 1 + ]); + y = new Complex64Array([ + 0.0, + 0.0, + -5.0, // 0 + 1.0, // 0 + -6.0, // 1 + 7.0 // 1 + ]); + expected = new Complex64( 7.0, 56.0 ); + + dot = cdotc( 2, x, 1, 0, y, 1, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports negative strides', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 1 + -8, // 1 + -4, // 0 + -7, // 0 + -1, + -9, + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, + 5, + 7, // 0 + -6, // 0 + 1, + -5 + ]); + expected = new Complex64( 104, 79 ); + + dot = cdotc( 2, x, -1, 1, y, -2, 2 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports complex access patterns', opts, function test( t ) { + var expected; + var dot; + var x; + var y; + + x = new Complex64Array([ + 7, // 0 + -8, // 0 + -4, + -7, + -1, // 1 + -9, // 1 + 2, + -8 + ]); + y = new Complex64Array([ + 6, // 1 + -6, // 1 + -9, // 0 + 5, // 0 + 7, + -6, + 1, + -5 + ]); + expected = new Complex64( -55, 23 ); + + dot = cdotc( 2, x, 2, 0, y, -1, 1 ); + t.deepEqual( dot, expected, 'returns expected value' ); + t.end(); +});