From 6c320f6e0882bdfbdff5ecb6721b937aeedc432d Mon Sep 17 00:00:00 2001 From: headlessNode Date: Sun, 8 Mar 2026 12:47:31 +0500 Subject: [PATCH 1/5] feat: add blas/ext/base/dcartesian-power --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: passed - task: lint_package_json status: passed - task: lint_repl_help status: passed - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: passed - task: lint_javascript_tests status: passed - task: lint_javascript_benchmarks status: passed - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: missing_dependencies - task: lint_c_benchmarks status: missing_dependencies - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: passed - task: lint_license_headers status: passed --- --- .../blas/ext/base/dcartesian-power/README.md | 332 ++++++++++++++++++ .../dcartesian-power/benchmark/benchmark.js | 104 ++++++ .../benchmark/benchmark.native.js | 109 ++++++ .../benchmark/benchmark.ndarray.js | 104 ++++++ .../benchmark/benchmark.ndarray.native.js | 109 ++++++ .../dcartesian-power/benchmark/c/Makefile | 146 ++++++++ .../benchmark/c/benchmark.length.c | 208 +++++++++++ .../ext/base/dcartesian-power/binding.gyp | 170 +++++++++ .../ext/base/dcartesian-power/docs/repl.txt | 111 ++++++ .../dcartesian-power/docs/types/index.d.ts | 109 ++++++ .../base/dcartesian-power/docs/types/test.ts | 310 ++++++++++++++++ .../base/dcartesian-power/examples/c/Makefile | 146 ++++++++ .../dcartesian-power/examples/c/example.c | 45 +++ .../base/dcartesian-power/examples/index.js | 30 ++ .../ext/base/dcartesian-power/include.gypi | 53 +++ .../stdlib/blas/ext/base/dcartesianpower.h | 45 +++ .../dcartesian-power/lib/dcartesianpower.js | 61 ++++ .../lib/dcartesianpower.native.js | 57 +++ .../ext/base/dcartesian-power/lib/index.js | 70 ++++ .../ext/base/dcartesian-power/lib/main.js | 35 ++ .../ext/base/dcartesian-power/lib/native.js | 35 ++ .../ext/base/dcartesian-power/lib/ndarray.js | 83 +++++ .../dcartesian-power/lib/ndarray.native.js | 61 ++++ .../ext/base/dcartesian-power/manifest.json | 82 +++++ .../ext/base/dcartesian-power/package.json | 74 ++++ .../ext/base/dcartesian-power/src/Makefile | 30 ++ .../ext/base/dcartesian-power/src/addon.c | 72 ++++ .../blas/ext/base/dcartesian-power/src/main.c | 81 +++++ .../test/test.dcartesianpower.js | 239 +++++++++++++ .../test/test.dcartesianpower.native.js | 248 +++++++++++++ .../ext/base/dcartesian-power/test/test.js | 82 +++++ .../dcartesian-power/test/test.ndarray.js | 254 ++++++++++++++ .../test/test.ndarray.native.js | 263 ++++++++++++++ 33 files changed, 3958 insertions(+) create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/binding.gyp create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include.gypi create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include/stdlib/blas/ext/base/dcartesianpower.h create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/main.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/package.json create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md new file mode 100644 index 000000000000..659ccbc5e23e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md @@ -0,0 +1,332 @@ + + +# dcartesianPower + +> Compute the Cartesian power for a double-precision floating-point strided array. + +
+ +
+ + + +
+ +## Usage + +```javascript +var dcartesianPower = require( '@stdlib/blas/ext/base/dcartesian-power' ); +``` + +#### dcartesianPower( N, k, x, strideX, out, strideOut1, strideOut2 ) + +Computes the Cartesian power for a double-precision floating-point strided array. + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, 2.0 ] ); +var out = new Float64Array( 8 ); + +dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); +// out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +``` + +The function has the following parameters: + +- **N**: number of indexed elements. +- **k**: power. +- **x**: input [`Float64Array`][@stdlib/array/float64]. +- **strideX**: stride length. +- **out**: output [`Float64Array`][@stdlib/array/float64]. +- **strideOut1**: stride length between consecutive output pairs. +- **strideOut2**: stride length between elements within each output pair. + +The `N`, `k`, and stride parameters determine which elements in the strided arrays are accessed at runtime. For example, to compute the Cartesian power of every other element: + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, 0.0, 2.0, 0.0 ] ); +var out = new Float64Array( 8 ); + +dcartesianPower( 2, 2, x, 2, out, 2, 1 ); +// out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +``` + +Note that indexing is relative to the first index. To introduce an offset, use [`typed array`][mdn-typed-array] views. + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +// Initial arrays... +var x0 = new Float64Array( [ 0.0, 1.0, 2.0 ] ); +var out0 = new Float64Array( 8 ); + +// Create offset views... +var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element + +dcartesianPower( 2, 2, x1, 1, out0, 2, 1 ); +// out0 => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +``` + + + +#### dcartesianPower.ndarray( N, k, x, strideX, offsetX, out, strideOut1, strideOut2, offsetOut ) + + + +Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, 2.0 ] ); +var out = new Float64Array( 8 ); + +dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ); +// out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +``` + +The function has the following additional parameters: + +- **offsetX**: starting index. +- **offsetOut**: starting index. + +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 compute the Cartesian power of every other value in the strided input array starting from the second value: + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 0.0, 1.0, 2.0 ] ); +var out = new Float64Array( 8 ); + +dcartesianPower.ndarray( 2, 2, x, 1, 1, out, 2, 1, 0 ); +// out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +``` + +
+ + + +
+ +## Notes + +- If `N <= 0` or `k <= 0`, both functions return `out` unchanged. +- The output array must contain at least `N^k` elements. + +
+ + + +
+ +## Examples + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var dcartesianPower = require( '@stdlib/blas/ext/base/dcartesian-power' ); + +var x = new Float64Array( [ 1.0, 2.0 ] ); +console.log( x ); + +var out = new Float64Array( 24 ); + +dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); +console.log( out ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/blas/ext/base/dcartesianpower.h" +``` + + + +#### stdlib_strided_dcartesian_power( N, k, \*X, strideX, \*Out, strideOut1, strideOut2 ) + + + +Computes the Cartesian power for a double-precision floating-point strided array. + +```c +const double X[] = { 1.0, 2.0 }; +double Out[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + +stdlib_strided_dcartesian_power( 2, 2, X, 1, Out, 2, 1 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **k**: `[in] CBLAS_INT` power. +- **X**: `[in] double*` input array. +- **strideX**: `[in] CBLAS_INT` stride length. +- **Out**: `[out] double*` output array. +- **strideOut1**: `[in] CBLAS_INT` stride length between consecutive output pairs. +- **strideOut2**: `[in] CBLAS_INT` stride length between elements within each output pair. + +```c +void stdlib_strided_dcartesian_power( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ); +``` + + + +#### stdlib_strided_dcartesian_power_ndarray( N, k, \*X, strideX, offsetX, \*Out, strideOut1, strideOut2, offsetOut ) + + + +Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. + +```c +const double X[] = { 1.0, 2.0 }; +double Out[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + +stdlib_strided_dcartesian_power_ndarray( 2, 2, X, 1, 0, Out, 2, 1, 0 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **k**: `[in] CBLAS_INT` power. +- **X**: `[in] double*` input array. +- **strideX**: `[in] CBLAS_INT` stride length. +- **offsetX**: `[in] CBLAS_INT` starting index. +- **Out**: `[out] double*` output array. +- **strideOut1**: `[in] CBLAS_INT` stride length between consecutive output pairs. +- **strideOut2**: `[in] CBLAS_INT` stride length between elements within each output pair. +- **offsetOut**: `[in] CBLAS_INT` starting index. + +```c +void stdlib_strided_dcartesian_power_ndarray( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2, const CBLAS_INT offsetOut ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/blas/ext/base/dcartesianpower.h" +#include + +int main( void ) { + // Create a strided input array: + const double X[] = { 1.0, 2.0 }; + + // Specify the number of indexed elements and power: + const int N = 2; + const int k = 2; + + // Create an output array: + double out[ 8 ] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + + // Specify strides: + const int strideX = 1; + const int strideOut1 = 2; + const int strideOut2 = 1; + + // Compute the Cartesian power: + stdlib_strided_dcartesian_power( N, k, X, strideX, out, strideOut1, strideOut2 ); + + // Print the result: + for ( int i = 0; i < N*N; i++ ) { + printf( "out[ %i ] = ( %lf, %lf )\n", i, out[ i*2 ], out[ i*2+1 ] ); + } +} +``` + +
+ + + +
+ + + +
+ +
+ + + + + +
+ +
+ + + + + +
+ +[@stdlib/array/float64]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/float64 + +[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray + +
+ + diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js new file mode 100644 index 000000000000..05d48cd0efca --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js @@ -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. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Float64Array = require( '@stdlib/array/float64' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var dcartesianPower = require( './../lib/dcartesianpower.js' ); + + +// VARIABLES // + +var K = 2; +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var out; + var x; + + x = uniform( len, -10.0, 10.0, options ); + out = new Float64Array( pow( len, K ) * K ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var z; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + z = dcartesianPower( x.length, K, x, 1, out, K, 1 ); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 3; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js new file mode 100644 index 000000000000..535efd55f221 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/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 resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var dcartesianPower = tryRequire( resolve( __dirname, './../lib/dcartesianpower.native.js' ) ); +var opts = { + 'skip': ( dcartesianPower instanceof Error ) +}; +var K = 2; +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var out; + var x; + + x = uniform( len, -10.0, 10.0, options ); + out = new Float64Array( pow( len, K ) * K ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var z; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + z = dcartesianPower( x.length, K, x, 1, out, K, 1 ); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 3; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s::native:len=%d', pkg, len ), opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.js new file mode 100644 index 000000000000..d07c9470b5ce --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.js @@ -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. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Float64Array = require( '@stdlib/array/float64' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var dcartesianPower = require( './../lib/ndarray.js' ); + + +// VARIABLES // + +var K = 2; +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var out; + var x; + + x = uniform( len, -10.0, 10.0, options ); + out = new Float64Array( pow( len, K ) * K ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var z; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + z = dcartesianPower( x.length, K, x, 1, 0, out, K, 1, 0 ); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 3; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s:ndarray:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.ndarray.native.js new file mode 100644 index 000000000000..cee20d158a90 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/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 resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var dcartesianPower = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( dcartesianPower instanceof Error ) +}; +var K = 2; +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var out; + var x; + + x = uniform( len, -10.0, 10.0, options ); + out = new Float64Array( pow( len, K ) * K ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var z; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + z = dcartesianPower( x.length, K, x, 1, 0, out, K, 1, 0 ); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( z[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 3; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s::native:ndarray:len=%d', pkg, len ), opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/Makefile new file mode 100644 index 000000000000..0756dc7da20a --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/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 := benchmark.length.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 benchmarks. +# +# @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/ext/base/dcartesian-power/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c new file mode 100644 index 000000000000..3a3901459268 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c @@ -0,0 +1,208 @@ +/** +* @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/ext/base/dcartesianpower.h" +#include +#include +#include +#include +#include + +#define NAME "dcartesian_power" +#define ITERATIONS 1000000 +#define REPEATS 3 +#define MIN 1 +#define MAX 3 +#define K 2 + +/** +* Prints the TAP version. +*/ +static void print_version( void ) { + printf( "TAP version 13\n" ); +} + +/** +* Prints the TAP summary. +* +* @param total total number of tests +* @param passing total number of passing tests +*/ +static void print_summary( int total, int passing ) { + printf( "#\n" ); + printf( "1..%d\n", total ); // TAP plan + printf( "# total %d\n", total ); + printf( "# pass %d\n", passing ); + printf( "#\n" ); + printf( "# ok\n" ); +} + +/** +* Prints benchmarks results. +* +* @param iterations number of iterations +* @param elapsed elapsed time in seconds +*/ +static void print_results( int iterations, double elapsed ) { + double rate = (double)iterations / elapsed; + printf( " ---\n" ); + printf( " iterations: %d\n", iterations ); + printf( " elapsed: %0.9f\n", elapsed ); + printf( " rate: %0.9f\n", rate ); + printf( " ...\n" ); +} + +/** +* Returns a clock time. +* +* @return clock time +*/ +static double tic( void ) { + struct timeval now; + gettimeofday( &now, NULL ); + return (double)now.tv_sec + (double)now.tv_usec/1.0e6; +} + +/** +* Generates a random number on the interval [0,1). +* +* @return random number +*/ +static double rand_double( void ) { + int r = rand(); + return (double)r / ( (double)RAND_MAX + 1.0 ); +} + +/** +* Runs a benchmark. +* +* @param iterations number of iterations +* @param len array length +* @return elapsed time in seconds +*/ +static double benchmark1( int iterations, int len ) { + double elapsed; + double *out; + double *x; + double t; + int outlen; + int i; + + x = (double *) malloc( len * sizeof( double ) ); + outlen = 1; + for ( i = 0; i < K; i++ ) { + outlen *= len; + } + out = (double *) malloc( outlen * K * sizeof( double ) ); + for ( i = 0; i < len; i++ ) { + x[ i ] = ( rand_double()*200.0 ) - 100.0; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + stdlib_strided_dcartesian_power( len, K, x, 1, out, K, 1 ); + if ( out[ 0 ] != out[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( out[ 0 ] != out[ 0 ] ) { + printf( "should not return NaN\n" ); + } + free( x ); + free( out ); + return elapsed; +} + +/** +* Runs a benchmark. +* +* @param iterations number of iterations +* @param len array length +* @return elapsed time in seconds +*/ +static double benchmark2( int iterations, int len ) { + double elapsed; + double *out; + double *x; + double t; + int outlen; + int i; + + x = (double *) malloc( len * sizeof( double ) ); + outlen = 1; + for ( i = 0; i < K; i++ ) { + outlen *= len; + } + out = (double *) malloc( outlen * K * sizeof( double ) ); + for ( i = 0; i < len; i++ ) { + x[ i ] = ( rand_double()*200.0 ) - 100.0; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + stdlib_strided_dcartesian_power_ndarray( len, K, x, 1, 0, out, K, 1, 0 ); + if ( out[ 0 ] != out[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( out[ 0 ] != out[ 0 ] ) { + printf( "should not return NaN\n" ); + } + free( x ); + free( out ); + return elapsed; +} + +/** +* Main execution sequence. +*/ +int main( void ) { + double elapsed; + int count; + int iter; + int len; + int i; + int j; + + // Use the current time to seed the random number generator: + srand( time( NULL ) ); + + print_version(); + count = 0; + for ( i = MIN; i <= MAX; i++ ) { + len = (int) pow( 10, i ); + iter = ITERATIONS / (int) pow( 10, i-1 ); + for ( j = 0; j < REPEATS; j++ ) { + count += 1; + printf( "# c::%s:len=%d\n", NAME, len ); + elapsed = benchmark1( iter, len ); + print_results( iter, elapsed ); + printf( "ok %d benchmark finished\n", count ); + } + for ( j = 0; j < REPEATS; j++ ) { + count += 1; + printf( "# c::%s:ndarray:len=%d\n", NAME, len ); + elapsed = benchmark2( iter, len ); + print_results( iter, elapsed ); + printf( "ok %d benchmark finished\n", count ); + } + } + print_summary( count, count ); +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/binding.gyp b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/binding.gyp new file mode 100644 index 000000000000..0d6508a12e99 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/binding.gyp @@ -0,0 +1,170 @@ +# @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', + + # 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 + }, # 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/ext/base/dcartesian-power/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt new file mode 100644 index 000000000000..8f036e35781f --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt @@ -0,0 +1,111 @@ +{{alias}}( N, k, x, strideX, out, strideOut1, strideOut2 ) + Computes the Cartesian power for a double-precision floating-point strided + array. + + The `N`, `k`, and stride parameters determine which elements in the strided + arrays are accessed at runtime. + + If `N <= 0` or `k <= 0`, the function returns `out` unchanged. + + Parameters + ---------- + N: integer + Number of indexed elements. + + k: integer + Power. + + x: Float64Array + Input array. + + strideX: integer + Stride length. + + out: Float64Array + Output array. + + strideOut1: integer + Stride length between consecutive output pairs. + + strideOut2: integer + Stride length between elements within each output pair. + + Returns + ------- + out: Float64Array + Output array. + + Examples + -------- + // Standard Usage: + > var x = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0 ] ); + > var out = new {{alias:@stdlib/array/float64}}( 8 ); + > {{alias}}( x.length, 2, x, 1, out, 2, 1 ) + [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] + + // Using `N` and stride parameters: + > x = new {{alias:@stdlib/array/float64}}( [ 1.0, 0.0, 2.0, 0.0 ] ); + > out = new {{alias:@stdlib/array/float64}}( 8 ); + > {{alias}}( 2, 2, x, 2, out, 2, 1 ) + [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] + + +{{alias}}.ndarray( N, k, x, strideX, offsetX, out, sO1, sO2, offsetOut ) + Computes the Cartesian power for a double-precision floating-point strided + array using alternative indexing semantics. + + While typed array views mandate a view offset based on the underlying + buffer, the offset parameters support indexing semantics based on starting + indices. + + Parameters + ---------- + N: integer + Number of indexed elements. + + k: integer + Power. + + x: Float64Array + Input array. + + strideX: integer + Stride length. + + offsetX: integer + Starting index. + + out: Float64Array + Output array. + + sO1: integer + Stride length between consecutive output pairs. + + sO2: integer + Stride length between elements within each output pair. + + offsetOut: integer + Starting index. + + Returns + ------- + out: Float64Array + Output array. + + Examples + -------- + // Standard Usage: + > var x = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0 ] ); + > var out = new {{alias:@stdlib/array/float64}}( 8 ); + > {{alias}}.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ) + [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] + + // Using an index offset: + > x = new {{alias:@stdlib/array/float64}}( [ 0.0, 1.0, 2.0 ] ); + > out = new {{alias:@stdlib/array/float64}}( 8 ); + > {{alias}}.ndarray( 2, 2, x, 1, 1, out, 2, 1, 0 ) + [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts new file mode 100644 index 000000000000..0240210bc806 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts @@ -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. +*/ + +// TypeScript Version: 4.1 + +/** +* Interface describing `dcartesianPower`. +*/ +interface Routine { + /** + * Computes the Cartesian power for a double-precision floating-point strided array. + * + * @param N - number of indexed elements + * @param k - power + * @param x - input array + * @param strideX - stride length + * @param out - output array + * @param strideOut1 - stride length between consecutive output pairs + * @param strideOut2 - stride length between elements within each output pair + * @returns output array + * + * @example + * var Float64Array = require( `@stdlib/array/float64` ); + * + * var x = new Float64Array( [ 1.0, 2.0 ] ); + * var out = new Float64Array( 8 ); + * + * dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] + */ + ( N: number, k: number, x: Float64Array, strideX: number, out: Float64Array, strideOut1: number, strideOut2: number ): Float64Array; + + /** + * Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. + * + * @param N - number of indexed elements + * @param k - power + * @param x - input array + * @param strideX - stride length + * @param offsetX - starting index + * @param out - output array + * @param strideOut1 - stride length between consecutive output pairs + * @param strideOut2 - stride length between elements within each output pair + * @param offsetOut - starting index + * @returns output array + * + * @example + * var Float64Array = require( `@stdlib/array/float64` ); + * + * var x = new Float64Array( [ 1.0, 2.0 ] ); + * var out = new Float64Array( 8 ); + * + * dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ); + * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] + */ + ndarray( N: number, k: number, x: Float64Array, strideX: number, offsetX: number, out: Float64Array, strideOut1: number, strideOut2: number, offsetOut: number ): Float64Array; +} + +/** +* Computes the Cartesian power for a double-precision floating-point strided array. +* +* @param N - number of indexed elements +* @param k - power +* @param x - input array +* @param strideX - stride length +* @param out - output array +* @param strideOut1 - stride length between consecutive output pairs +* @param strideOut2 - stride length between elements within each output pair +* @returns output array +* +* @example +* var Float64Array = require( `@stdlib/array/float64` ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +* +* @example +* var Float64Array = require( `@stdlib/array/float64` ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +*/ +declare var dcartesianPower: Routine; + + +// EXPORTS // + +export = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts new file mode 100644 index 000000000000..a3b37f63d693 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts @@ -0,0 +1,310 @@ +/* +* @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 dcartesianPower = require( './index' ); + + +// TESTS // + +// The function returns a Float64Array... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); // $ExpectType Float64Array +} + +// The compiler throws an error if the function is provided a first argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower( '10', 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( true, 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( false, 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( null, 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( undefined, 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( [], 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( {}, 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( ( x: number ): number => x, 2, x, 1, out, 2, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a second argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower( x.length, '10', x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, true, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, false, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, null, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, undefined, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, [], x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, {}, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, ( x: number ): number => x, x, 1, out, 2, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a third argument which is not a Float64Array... +{ + const out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, 10, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, '10', 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, true, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, false, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, null, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, undefined, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, [], 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, {}, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, ( x: number ): number => x, 1, out, 2, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fourth argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, '10', out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, true, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, false, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, null, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, undefined, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, [], out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, {}, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, ( x: number ): number => x, out, 2, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fifth argument which is not a Float64Array... +{ + const x = new Float64Array( 2 ); + + dcartesianPower( x.length, 2, x, 1, 10, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, '10', 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, true, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, false, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, null, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, undefined, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, [], 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, {}, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, ( x: number ): number => x, 2, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a sixth argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, out, '10', 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, true, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, false, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, null, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, undefined, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, [], 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, {}, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, ( x: number ): number => x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a seventh argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, out, 2, '10' ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, true ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, false ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, null ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, undefined ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, [] ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, {} ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower(); // $ExpectError + dcartesianPower( x.length ); // $ExpectError + dcartesianPower( x.length, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x ); // $ExpectError + dcartesianPower( x.length, 2, x, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, 1, 10 ); // $ExpectError +} + +// The `ndarray` method returns a Float64Array... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectType Float64Array +} + +// The compiler throws an error if the `ndarray` method is provided a first argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( '10', 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( true, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( false, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( null, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( undefined, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( [], 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( {}, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( ( x: number ): number => x, 2, x, 1, 0, out, 2, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a second argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, '10', x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, true, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, false, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, null, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, undefined, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, [], x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, {}, x, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, ( x: number ): number => x, x, 1, 0, out, 2, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a third argument which is not a Float64Array... +{ + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( 2, 2, 10, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, '10', 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, true, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, false, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, null, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, undefined, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, [], 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, {}, 1, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( 2, 2, ( x: number ): number => x, 1, 0, out, 2, 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 Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, 2, x, '10', 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, true, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, false, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, null, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, undefined, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, [], 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, {}, 0, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, ( x: number ): number => x, 0, out, 2, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, 2, x, 1, '10', out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, true, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, false, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, null, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, undefined, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, [], out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, {}, out, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, ( x: number ): number => x, out, 2, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a Float64Array... +{ + const x = new Float64Array( 2 ); + + dcartesianPower.ndarray( x.length, 2, x, 1, 0, 10, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, '10', 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, true, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, false, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, null, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, undefined, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, [], 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, {}, 2, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, ( x: number ): number => x, 2, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, '10', 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, true, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, false, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, null, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, undefined, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, [], 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, {}, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, ( x: number ): number => x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an eighth argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, '10', 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, true, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, false, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, null, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, undefined, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, [], 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, {}, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a ninth argument which is not a number... +{ + const x = new Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, '10' ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, true ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, false ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, null ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, undefined ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, [] ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, {} ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 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 Float64Array( 2 ); + const out = new Float64Array( 8 ); + + dcartesianPower.ndarray(); // $ExpectError + dcartesianPower.ndarray( x.length ); // $ExpectError + dcartesianPower.ndarray( x.length, 2 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1 ); // $ExpectError + dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0, 10 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/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/ext/base/dcartesian-power/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c new file mode 100644 index 000000000000..63c4351738b5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c @@ -0,0 +1,45 @@ +/** +* @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/ext/base/dcartesianpower.h" +#include + +int main( void ) { + // Create a strided input array: + const double X[] = { 1.0, 2.0 }; + + // Specify the number of indexed elements and power: + const int N = 2; + const int k = 2; + + // Create an output array: + double out[ 8 ] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + + // Specify strides: + const int strideX = 1; + const int strideOut1 = 2; + const int strideOut2 = 1; + + // Compute the Cartesian power: + stdlib_strided_dcartesian_power( N, k, X, strideX, out, strideOut1, strideOut2 ); + + // Print the result: + for ( int i = 0; i < N*N; i++ ) { + printf( "out[ %i ] = ( %lf, %lf )\n", i, out[ i*2 ], out[ i*2+1 ] ); + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js new file mode 100644 index 000000000000..6720a45a2655 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js @@ -0,0 +1,30 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var dcartesianPower = require( './../lib' ); + +var x = new Float64Array( [ 1.0, 2.0 ] ); +console.log( x ); + +var out = new Float64Array( 24 ); + +dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); +console.log( out ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include.gypi b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include.gypi new file mode 100644 index 000000000000..bee8d41a2caf --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include.gypi @@ -0,0 +1,53 @@ +# @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. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + ' [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +*/ +function dcartesianPower( N, k, x, strideX, out, strideOut1, strideOut2 ) { + var len = pow( N, k ); + var ox = stride2offset( N, strideX ); + var oo = stride2offset( len, strideOut1 ) + stride2offset( k, strideOut2 ); + return ndarray( N, k, x, strideX, ox, out, strideOut1, strideOut2, oo ); +} + + +// EXPORTS // + +module.exports = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js new file mode 100644 index 000000000000..836afe675913 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js @@ -0,0 +1,57 @@ +/** +* @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 addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Computes the Cartesian power for a double-precision floating-point strided array. +* +* @param {NonNegativeInteger} N - number of indexed elements +* @param {NonNegativeInteger} k - power +* @param {Float64Array} x - input array +* @param {integer} strideX - stride length +* @param {Float64Array} out - output array +* @param {integer} strideOut1 - stride length between consecutive output pairs +* @param {integer} strideOut2 - stride length between elements within each output pair +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower( 2, 2, x, 1, out, 2, 1 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +*/ +function dcartesianPower( N, k, x, strideX, out, strideOut1, strideOut2 ) { + addon( N, k, x, strideX, out, strideOut1, strideOut2 ); + return out; +} + + +// EXPORTS // + +module.exports = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js new file mode 100644 index 000000000000..302d47ad7a08 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js @@ -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. +*/ + +'use strict'; + +/** +* Compute the Cartesian power for a double-precision floating-point strided array. +* +* @module @stdlib/blas/ext/base/dcartesian-power +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var dcartesianPower = require( '@stdlib/blas/ext/base/dcartesian-power' ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var dcartesianPower = require( '@stdlib/blas/ext/base/dcartesian-power' ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.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 dcartesianPower; +var tmp = tryRequire( join( __dirname, './native.js' ) ); +if ( isError( tmp ) ) { + dcartesianPower = main; +} else { + dcartesianPower = tmp; +} + + +// EXPORTS // + +module.exports = dcartesianPower; + +// exports: { "ndarray": "dcartesianPower.ndarray" } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/main.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/main.js new file mode 100644 index 000000000000..8611af173e7f --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/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 dcartesianPower = require( './dcartesianpower.js' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +setReadOnly( dcartesianPower, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/native.js new file mode 100644 index 000000000000..b9025a1e6d77 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/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 dcartesianPower = require( './dcartesianpower.native.js' ); +var ndarray = require( './ndarray.native.js' ); + + +// MAIN // + +setReadOnly( dcartesianPower, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js new file mode 100644 index 000000000000..7e9634806077 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js @@ -0,0 +1,83 @@ +/** +* @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 pow = require( '@stdlib/math/base/special/pow' ); + + +// MAIN // + +/** +* Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. +* +* @param {NonNegativeInteger} N - number of indexed elements +* @param {NonNegativeInteger} k - power +* @param {Float64Array} x - input array +* @param {integer} strideX - stride length +* @param {integer} offsetX - starting index +* @param {Float64Array} out - output array +* @param {integer} strideOut1 - stride length between consecutive output pairs +* @param {integer} strideOut2 - stride length between elements within each output pair +* @param {integer} offsetOut - starting index +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower( 2, 2, x, 1, 0, out, 2, 1, 0 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +*/ +function dcartesianPower( N, k, x, strideX, offsetX, out, strideOut1, strideOut2, offsetOut ) { // eslint-disable-line max-len + var len; + var oj; + var io; + var ix; + var i; + var j; + var s; + var t; + + if ( N <= 0 || k <= 0 ) { + return out; + } + len = pow( N, k ); + io = offsetOut; + for ( i = 0; i < len; i++ ) { + t = i; + for ( j = k - 1; j >= 0; j-- ) { + s = t % N; + t = ( t - s ) / N; + ix = offsetX + ( s * strideX ); + oj = io + ( j * strideOut2 ); + out[ oj ] = x[ ix ]; + } + io += strideOut1; + } + return out; +} + + +// EXPORTS // + +module.exports = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js new file mode 100644 index 000000000000..b97f519aee49 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.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. +*/ + +/* eslint-disable max-len */ + +'use strict'; + +// MODULES // + +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. +* +* @param {NonNegativeInteger} N - number of indexed elements +* @param {NonNegativeInteger} k - power +* @param {Float64Array} x - input array +* @param {integer} strideX - stride length +* @param {integer} offsetX - starting index +* @param {Float64Array} out - output array +* @param {integer} strideOut1 - stride length between consecutive output pairs +* @param {integer} strideOut2 - stride length between elements within each output pair +* @param {integer} offsetOut - starting index +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, 2.0 ] ); +* var out = new Float64Array( 8 ); +* +* dcartesianPower( 2, 2, x, 1, 0, out, 2, 1, 0 ); +* // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] +*/ +function dcartesianPower( N, k, x, strideX, offsetX, out, strideOut1, strideOut2, offsetOut ) { + addon.ndarray( N, k, x, strideX, offsetX, out, strideOut1, strideOut2, offsetOut ); + return out; +} + + +// EXPORTS // + +module.exports = dcartesianPower; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json new file mode 100644 index 000000000000..f35e9860ebd5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json @@ -0,0 +1,82 @@ +{ + "options": { + "task": "build" + }, + "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", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/math/base/special/pow", + "@stdlib/strided/base/stride2offset", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-float64array", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-float64array" + ] + }, + { + "task": "benchmark", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/math/base/special/pow", + "@stdlib/strided/base/stride2offset" + ] + }, + { + "task": "examples", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/math/base/special/pow", + "@stdlib/strided/base/stride2offset" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/package.json b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/package.json new file mode 100644 index 000000000000..3c13a979f0e2 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/package.json @@ -0,0 +1,74 @@ +{ + "name": "@stdlib/blas/ext/base/dcartesian-power", + "version": "0.0.0", + "description": "Compute the Cartesian power for a double-precision floating-point strided array.", + "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", + "extended", + "cartesian", + "power", + "strided", + "array", + "ndarray", + "float64", + "double", + "float64array" + ], + "__stdlib__": { + "wasm": false + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/Makefile new file mode 100644 index 000000000000..57c89dbb96d4 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/Makefile @@ -0,0 +1,30 @@ +#/ +# @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. +#/ + +# RULES # + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c new file mode 100644 index 000000000000..715c84628687 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c @@ -0,0 +1,72 @@ +/** +* @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/ext/base/dcartesianpower.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/napi/argv_float64array.h" +#include "stdlib/napi/argv_strided_float64array.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_NAPI_ARGV( env, info, argv, argc, 7 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, k, argv, 1 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, strideOut1, argv, 5 ); + STDLIB_NAPI_ARGV_INT64( env, strideOut2, argv, 6 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_FLOAT64ARRAY( env, Out, out_len, argv, 4 ); + (void)out_len; + API_SUFFIX(stdlib_strided_dcartesian_power)( N, k, X, strideX, Out, strideOut1, strideOut2 ); + return NULL; +} + +/** +* Receives JavaScript callback invocation data. +* +* @param env environment under which the function is invoked +* @param info callback data +* @return Node-API value +*/ +static napi_value addon_method( napi_env env, napi_callback_info info ) { + STDLIB_NAPI_ARGV( env, info, argv, argc, 9 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, k, argv, 1 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 ); + STDLIB_NAPI_ARGV_INT64( env, strideOut1, argv, 6 ); + STDLIB_NAPI_ARGV_INT64( env, strideOut2, argv, 7 ); + STDLIB_NAPI_ARGV_INT64( env, offsetOut, argv, 8 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_FLOAT64ARRAY( env, Out, out_len, argv, 5 ); + (void)out_len; + API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, offsetX, Out, strideOut1, strideOut2, offsetOut ); + return NULL; +} + +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c new file mode 100644 index 000000000000..540f9e6f9990 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c @@ -0,0 +1,81 @@ +/** +* @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/ext/base/dcartesianpower.h" +#include "stdlib/math/base/special/pow.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" + +/** +* Computes the Cartesian power for a double-precision floating-point strided array. +* +* @param N number of indexed elements +* @param k power +* @param X input array +* @param strideX stride length +* @param Out output array +* @param strideOut1 stride length between consecutive output pairs +* @param strideOut2 stride length between elements within each output pair +*/ +void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ) { + CBLAS_INT len = stdlib_base_pow( N, k ); + CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); + CBLAS_INT oo = stdlib_strided_stride2offset( len, strideOut1 ) + stdlib_strided_stride2offset( k, strideOut2 ); + API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, ox, Out, strideOut1, strideOut2, oo ); +} + +/** +* Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. +* +* @param N number of indexed elements +* @param k power +* @param X input array +* @param strideX stride length +* @param offsetX starting index +* @param Out output array +* @param strideOut1 stride length between consecutive output pairs +* @param strideOut2 stride length between elements within each output pair +* @param offsetOut starting index +*/ +void API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2, const CBLAS_INT offsetOut ) { + CBLAS_INT len; + CBLAS_INT oj; + CBLAS_INT io; + CBLAS_INT ix; + CBLAS_INT i; + CBLAS_INT j; + CBLAS_INT s; + CBLAS_INT t; + + if ( N <= 0 || k <= 0 ) { + return; + } + len = stdlib_base_pow( N, k ); + io = offsetOut; + for ( i = 0; i < len; i++ ) { + t = i; + for ( j = k - 1; j >= 0; j-- ) { + s = t % N; + t = ( t - s ) / N; + ix = offsetX + ( s * strideX ); + oj = io + ( j * strideOut2 ); + Out[ oj ] = X[ ix ]; + } + io += strideOut1; + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js new file mode 100644 index 000000000000..76b74cb590a4 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js @@ -0,0 +1,239 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var dcartesianPower = require( './../lib/dcartesianpower.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dcartesianPower, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', function test( t ) { + t.strictEqual( dcartesianPower.length, 7, 'returns expected value' ); + t.end(); +}); + +tape( 'the function computes the Cartesian power', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 24 ); + + dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); + expected = new Float64Array([ + 1.0, + 1.0, + 1.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 1.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the output array', function test( t ) { + var out; + var x; + var v; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + v = dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + + t.strictEqual( v, out, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + expected = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dcartesianPower( 0, 2, x, 1, out, 2, 1 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + dcartesianPower( 2, 0, x, 1, out, 2, 1 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + t.end(); +}); + +tape( 'the function supports a strideX parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 1.0, // 0 + 0.0, + 2.0, // 1 + 0.0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 2, out, 2, 1 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut1 parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, out, 4, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 0.0, + 0.0, + 1.0, + 2.0, + 0.0, + 0.0, + 2.0, + 1.0, + 0.0, + 0.0, + 2.0, + 2.0, + 0.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut2 parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, out, 4, 2 ); + + expected = new Float64Array([ + 1.0, + 0.0, + 1.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports negative strides', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 2.0, // 1 + 1.0 // 0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, -1, out, -2, 1 ); + + expected = new Float64Array( [ 2.0, 2.0, 2.0, 1.0, 1.0, 2.0, 1.0, 1.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports view offsets', function test( t ) { + var expected; + var out0; + var x0; + var x1; + + x0 = new Float64Array( [ 0.0, 1.0, 2.0 ] ); + out0 = new Float64Array( 8 ); + + // Create an offset view... + x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element + + dcartesianPower( 2, 2, x1, 1, out0, 2, 1 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out0, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js new file mode 100644 index 000000000000..1348fb4b531a --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js @@ -0,0 +1,248 @@ +/** +* @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 resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var dcartesianPower = tryRequire( resolve( __dirname, './../lib/dcartesianpower.native.js' ) ); +var opts = { + 'skip': ( dcartesianPower instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dcartesianPower, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', opts, function test( t ) { + t.strictEqual( dcartesianPower.length, 7, 'returns expected value' ); + t.end(); +}); + +tape( 'the function computes the Cartesian power', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 24 ); + + dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); + expected = new Float64Array([ + 1.0, + 1.0, + 1.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 1.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the output array', opts, function test( t ) { + var out; + var x; + var v; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + v = dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + + t.strictEqual( v, out, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + expected = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dcartesianPower( 0, 2, x, 1, out, 2, 1 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + dcartesianPower( 2, 0, x, 1, out, 2, 1 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + t.end(); +}); + +tape( 'the function supports a strideX parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 1.0, // 0 + 0.0, + 2.0, // 1 + 0.0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 2, out, 2, 1 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut1 parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, out, 4, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 0.0, + 0.0, + 1.0, + 2.0, + 0.0, + 0.0, + 2.0, + 1.0, + 0.0, + 0.0, + 2.0, + 2.0, + 0.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut2 parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, out, 4, 2 ); + + expected = new Float64Array([ + 1.0, + 0.0, + 1.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports negative strides', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 2.0, // 1 + 1.0 // 0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, -1, out, -2, 1 ); + + expected = new Float64Array( [ 2.0, 2.0, 2.0, 1.0, 1.0, 2.0, 1.0, 1.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports view offsets', opts, function test( t ) { + var expected; + var out0; + var x0; + var x1; + + x0 = new Float64Array( [ 0.0, 1.0, 2.0 ] ); + out0 = new Float64Array( 8 ); + + // Create an offset view... + x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element + + dcartesianPower( 2, 2, x1, 1, out0, 2, 1 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out0, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.js new file mode 100644 index 000000000000..afe716d02b98 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/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 dcartesianPower = 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 dcartesianPower, '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 dcartesianPower.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 dcartesianPower = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( dcartesianPower, 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 dcartesianPower; + var main; + + main = require( './../lib/dcartesianpower.js' ); + + dcartesianPower = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( dcartesianPower, main, 'returns expected value' ); + t.end(); + + function tryRequire() { + return new Error( 'Cannot find module' ); + } +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js new file mode 100644 index 000000000000..6b012aaf55e5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js @@ -0,0 +1,254 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var dcartesianPower = require( './../lib/ndarray.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dcartesianPower, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 9', function test( t ) { + t.strictEqual( dcartesianPower.length, 9, 'returns expected value' ); + t.end(); +}); + +tape( 'the function computes the Cartesian power', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, 0, out, 2, 1, 0 ); + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 24 ); + + dcartesianPower( x.length, 3, x, 1, 0, out, 3, 1, 0 ); + expected = new Float64Array([ + 1.0, + 1.0, + 1.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 1.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the output array', function test( t ) { + var out; + var x; + var v; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + v = dcartesianPower( x.length, 2, x, 1, 0, out, 2, 1, 0 ); + + t.strictEqual( v, out, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + expected = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dcartesianPower( 0, 2, x, 1, 0, out, 2, 1, 0 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + dcartesianPower( 2, 0, x, 1, 0, out, 2, 1, 0 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + t.end(); +}); + +tape( 'the function supports a strideX parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 1.0, // 0 + 0.0, + 2.0, // 1 + 0.0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 2, 0, out, 2, 1, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut1 parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 4, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 0.0, + 0.0, + 1.0, + 2.0, + 0.0, + 0.0, + 2.0, + 1.0, + 0.0, + 0.0, + 2.0, + 2.0, + 0.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut2 parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 4, 2, 0 ); + + expected = new Float64Array([ + 1.0, + 0.0, + 1.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports negative strides', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 2.0, // 1 + 1.0 // 0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, -1, 1, out, -2, 1, 6 ); + + expected = new Float64Array( [ 2.0, 2.0, 2.0, 1.0, 1.0, 2.0, 1.0, 1.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports an offsetX parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 0.0, + 1.0, // 0 + 2.0 // 1 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 1, 1, out, 2, 1, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports an offsetOut parameter', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 12 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 2, 1, 4 ); + + expected = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js new file mode 100644 index 000000000000..7c1900d308bf --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js @@ -0,0 +1,263 @@ +/** +* @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 resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var dcartesianPower = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( dcartesianPower instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dcartesianPower, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 9', opts, function test( t ) { + t.strictEqual( dcartesianPower.length, 9, 'returns expected value' ); + t.end(); +}); + +tape( 'the function computes the Cartesian power', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + dcartesianPower( x.length, 2, x, 1, 0, out, 2, 1, 0 ); + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 24 ); + + dcartesianPower( x.length, 3, x, 1, 0, out, 3, 1, 0 ); + expected = new Float64Array([ + 1.0, + 1.0, + 1.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 1.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 1.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0, + 1.0, + 2.0, + 2.0, + 2.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the output array', opts, function test( t ) { + var out; + var x; + var v; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + v = dcartesianPower( x.length, 2, x, 1, 0, out, 2, 1, 0 ); + + t.strictEqual( v, out, 'same reference' ); + t.end(); +}); + +tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + expected = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dcartesianPower( 0, 2, x, 1, 0, out, 2, 1, 0 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + dcartesianPower( 2, 0, x, 1, 0, out, 2, 1, 0 ); + t.deepEqual( out, expected, 'returns `out` unchanged' ); + + t.end(); +}); + +tape( 'the function supports a strideX parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 1.0, // 0 + 0.0, + 2.0, // 1 + 0.0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 2, 0, out, 2, 1, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut1 parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 4, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 0.0, + 0.0, + 1.0, + 2.0, + 0.0, + 0.0, + 2.0, + 1.0, + 0.0, + 0.0, + 2.0, + 2.0, + 0.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports a strideOut2 parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 16 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 4, 2, 0 ); + + expected = new Float64Array([ + 1.0, + 0.0, + 1.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0, + 1.0, + 0.0, + 2.0, + 0.0, + 2.0, + 0.0 + ]); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports negative strides', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 2.0, // 1 + 1.0 // 0 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, -1, 1, out, -2, 1, 6 ); + + expected = new Float64Array( [ 2.0, 2.0, 2.0, 1.0, 1.0, 2.0, 1.0, 1.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports an offsetX parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array([ + 0.0, + 1.0, // 0 + 2.0 // 1 + ]); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 1, 1, out, 2, 1, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports an offsetOut parameter', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 12 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 2, 1, 4 ); + + expected = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); From 2e393b93231ebfecee48896c20bfe5c9f96ef95e Mon Sep 17 00:00:00 2001 From: headlessNode Date: Sun, 8 Mar 2026 23:45:41 +0500 Subject: [PATCH 2/5] fix: apply suggestions from code review --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: passed - task: lint_javascript_benchmarks status: na - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../blas/ext/base/dcartesian-power/lib/ndarray.js | 4 ++-- .../blas/ext/base/dcartesian-power/lib/ndarray.native.js | 4 ++-- .../@stdlib/blas/ext/base/dcartesian-power/manifest.json | 4 ++-- .../@stdlib/blas/ext/base/dcartesian-power/src/addon.c | 9 ++++----- .../base/dcartesian-power/test/test.dcartesianpower.js | 2 +- .../dcartesian-power/test/test.dcartesianpower.native.js | 2 +- .../blas/ext/base/dcartesian-power/test/test.ndarray.js | 2 +- .../base/dcartesian-power/test/test.ndarray.native.js | 2 +- 8 files changed, 14 insertions(+), 15 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js index 7e9634806077..9c7fa7b2fad8 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js @@ -32,11 +32,11 @@ var pow = require( '@stdlib/math/base/special/pow' ); * @param {NonNegativeInteger} k - power * @param {Float64Array} x - input array * @param {integer} strideX - stride length -* @param {integer} offsetX - starting index +* @param {NonNegativeInteger} offsetX - starting index * @param {Float64Array} out - output array * @param {integer} strideOut1 - stride length between consecutive output pairs * @param {integer} strideOut2 - stride length between elements within each output pair -* @param {integer} offsetOut - starting index +* @param {NonNegativeInteger} offsetOut - starting index * @returns {Float64Array} output array * * @example diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js index b97f519aee49..e437f3ea5e19 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.native.js @@ -34,11 +34,11 @@ var addon = require( './../src/addon.node' ); * @param {NonNegativeInteger} k - power * @param {Float64Array} x - input array * @param {integer} strideX - stride length -* @param {integer} offsetX - starting index +* @param {NonNegativeInteger} offsetX - starting index * @param {Float64Array} out - output array * @param {integer} strideOut1 - stride length between consecutive output pairs * @param {integer} strideOut2 - stride length between elements within each output pair -* @param {integer} offsetOut - starting index +* @param {NonNegativeInteger} offsetOut - starting index * @returns {Float64Array} output array * * @example diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json index f35e9860ebd5..a5be842aa2a3 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json @@ -41,9 +41,9 @@ "@stdlib/strided/base/stride2offset", "@stdlib/napi/export", "@stdlib/napi/argv", - "@stdlib/napi/argv-float64array", "@stdlib/napi/argv-int64", - "@stdlib/napi/argv-strided-float64array" + "@stdlib/napi/argv-strided-float64array", + "@stdlib/napi/argv-strided-float64array2d" ] }, { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c index 715c84628687..a9a1e7ac2bef 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c @@ -17,12 +17,13 @@ */ #include "stdlib/blas/ext/base/dcartesianpower.h" +#include "stdlib/math/base/special/pow.h" #include "stdlib/blas/base/shared.h" #include "stdlib/napi/export.h" #include "stdlib/napi/argv.h" #include "stdlib/napi/argv_int64.h" -#include "stdlib/napi/argv_float64array.h" #include "stdlib/napi/argv_strided_float64array.h" +#include "stdlib/napi/argv_strided_float64array2d.h" #include /** @@ -40,8 +41,7 @@ static napi_value addon( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_INT64( env, strideOut1, argv, 5 ); STDLIB_NAPI_ARGV_INT64( env, strideOut2, argv, 6 ); STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); - STDLIB_NAPI_ARGV_FLOAT64ARRAY( env, Out, out_len, argv, 4 ); - (void)out_len; + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY2D( env, Out, stdlib_base_pow( N, k ), k, strideOut1, strideOut2, argv, 4 ); API_SUFFIX(stdlib_strided_dcartesian_power)( N, k, X, strideX, Out, strideOut1, strideOut2 ); return NULL; } @@ -63,8 +63,7 @@ static napi_value addon_method( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_INT64( env, strideOut2, argv, 7 ); STDLIB_NAPI_ARGV_INT64( env, offsetOut, argv, 8 ); STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); - STDLIB_NAPI_ARGV_FLOAT64ARRAY( env, Out, out_len, argv, 5 ); - (void)out_len; + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY2D( env, Out, stdlib_base_pow( N, k ), k, strideOut1, strideOut2, argv, 5 ); API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, offsetX, Out, strideOut1, strideOut2, offsetOut ); return NULL; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js index 76b74cb590a4..7b168eb32cc2 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js @@ -99,7 +99,7 @@ tape( 'the function returns a reference to the output array', function test( t ) t.end(); }); -tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', function test( t ) { +tape( 'if provided an `N` or `k` parameter equal to `0`, the function returns `out` unchanged', function test( t ) { var expected; var out; var x; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js index 1348fb4b531a..5c9a01ebf48e 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js @@ -108,7 +108,7 @@ tape( 'the function returns a reference to the output array', opts, function tes t.end(); }); -tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', opts, function test( t ) { +tape( 'if provided an `N` or `k` parameter equal to `0`, the function returns `out` unchanged', opts, function test( t ) { var expected; var out; var x; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js index 6b012aaf55e5..8f4c354f4fbe 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js @@ -99,7 +99,7 @@ tape( 'the function returns a reference to the output array', function test( t ) t.end(); }); -tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', function test( t ) { +tape( 'if provided an `N` or `k` parameter equal to `0`, the function returns `out` unchanged', function test( t ) { var expected; var out; var x; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js index 7c1900d308bf..e7a490c7fab2 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js @@ -108,7 +108,7 @@ tape( 'the function returns a reference to the output array', opts, function tes t.end(); }); -tape( 'if provided an `N` or `k` parameter less than or equal to `0`, the function returns `out` unchanged', opts, function test( t ) { +tape( 'if provided an `N` or `k` parameter equal to `0`, the function returns `out` unchanged', opts, function test( t ) { var expected; var out; var x; From 85717ed6085822c72854d216218dd8828789fdb6 Mon Sep 17 00:00:00 2001 From: headlessNode Date: Mon, 9 Mar 2026 03:49:38 +0500 Subject: [PATCH 3/5] refactor: apply suggestions from code review --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: na - task: lint_javascript_benchmarks status: na - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../blas/ext/base/dcartesian-power/src/main.c | 15 ++++++++------- 1 file changed, 8 insertions(+), 7 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c index 540f9e6f9990..267e8a3536be 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c @@ -20,6 +20,7 @@ #include "stdlib/math/base/special/pow.h" #include "stdlib/blas/base/shared.h" #include "stdlib/strided/base/stride2offset.h" +#include /** * Computes the Cartesian power for a double-precision floating-point strided array. @@ -33,7 +34,7 @@ * @param strideOut2 stride length between elements within each output pair */ void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ) { - CBLAS_INT len = stdlib_base_pow( N, k ); + int64_t len = (int64_t)stdlib_base_pow( N, k ); CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); CBLAS_INT oo = stdlib_strided_stride2offset( len, strideOut1 ) + stdlib_strided_stride2offset( k, strideOut2 ); API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, ox, Out, strideOut1, strideOut2, oo ); @@ -53,19 +54,19 @@ void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS * @param offsetOut starting index */ void API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2, const CBLAS_INT offsetOut ) { - CBLAS_INT len; - CBLAS_INT oj; - CBLAS_INT io; CBLAS_INT ix; - CBLAS_INT i; CBLAS_INT j; CBLAS_INT s; - CBLAS_INT t; + int64_t len; + int64_t io; + int64_t oj; + int64_t i; + int64_t t; if ( N <= 0 || k <= 0 ) { return; } - len = stdlib_base_pow( N, k ); + len = (int64_t)stdlib_base_pow( N, k ); io = offsetOut; for ( i = 0; i < len; i++ ) { t = i; From d49a25e7447ceb672296b6429ae4ea85054029e9 Mon Sep 17 00:00:00 2001 From: headlessNode Date: Mon, 9 Mar 2026 03:53:01 +0500 Subject: [PATCH 4/5] fix: apply suggestions from code review --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: na - task: lint_package_json status: na - task: lint_repl_help status: na - task: lint_javascript_src status: na - task: lint_javascript_cli status: na - task: lint_javascript_examples status: na - task: lint_javascript_tests status: na - task: lint_javascript_benchmarks status: na - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: na - task: lint_c_benchmarks status: na - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: na - task: lint_license_headers status: passed --- --- .../@stdlib/blas/ext/base/dcartesian-power/src/main.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c index 267e8a3536be..974eff4b0a4f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c @@ -34,7 +34,7 @@ * @param strideOut2 stride length between elements within each output pair */ void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ) { - int64_t len = (int64_t)stdlib_base_pow( N, k ); + int64_t len = stdlib_base_pow( N, k ); CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); CBLAS_INT oo = stdlib_strided_stride2offset( len, strideOut1 ) + stdlib_strided_stride2offset( k, strideOut2 ); API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, ox, Out, strideOut1, strideOut2, oo ); @@ -66,7 +66,7 @@ void API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( const CBLAS_INT N, con if ( N <= 0 || k <= 0 ) { return; } - len = (int64_t)stdlib_base_pow( N, k ); + len = stdlib_base_pow( N, k ); io = offsetOut; for ( i = 0; i < len; i++ ) { t = i; From 77cbbf7690421e62305e56a188d4f71e3c2b5558 Mon Sep 17 00:00:00 2001 From: headlessNode Date: Tue, 10 Mar 2026 01:43:19 +0500 Subject: [PATCH 5/5] fix: apply suggestions from code review --- type: pre_commit_static_analysis_report description: Results of running static analysis checks when committing changes. report: - task: lint_filenames status: passed - task: lint_editorconfig status: passed - task: lint_markdown status: passed - task: lint_package_json status: na - task: lint_repl_help status: passed - task: lint_javascript_src status: passed - task: lint_javascript_cli status: na - task: lint_javascript_examples status: passed - task: lint_javascript_tests status: passed - task: lint_javascript_benchmarks status: passed - task: lint_python status: na - task: lint_r status: na - task: lint_c_src status: missing_dependencies - task: lint_c_examples status: missing_dependencies - task: lint_c_benchmarks status: missing_dependencies - task: lint_c_tests_fixtures status: na - task: lint_shell status: na - task: lint_typescript_declarations status: passed - task: lint_typescript_tests status: passed - task: lint_license_headers status: passed --- --- .../blas/ext/base/dcartesian-power/README.md | 49 ++++--- .../dcartesian-power/benchmark/benchmark.js | 2 +- .../benchmark/benchmark.native.js | 2 +- .../benchmark/c/benchmark.length.c | 6 +- .../ext/base/dcartesian-power/docs/repl.txt | 33 +++-- .../dcartesian-power/docs/types/index.d.ts | 26 ++-- .../base/dcartesian-power/docs/types/test.ts | 120 ++++++++---------- .../dcartesian-power/examples/c/example.c | 5 +- .../base/dcartesian-power/examples/index.js | 2 +- .../stdlib/blas/ext/base/dcartesianpower.h | 2 +- .../dcartesian-power/lib/dcartesianpower.js | 21 +-- .../lib/dcartesianpower.native.js | 9 +- .../ext/base/dcartesian-power/lib/index.js | 2 +- .../ext/base/dcartesian-power/lib/ndarray.js | 48 ++++--- .../ext/base/dcartesian-power/manifest.json | 3 + .../ext/base/dcartesian-power/src/addon.c | 14 +- .../blas/ext/base/dcartesian-power/src/main.c | 78 ++++++++---- .../test/test.dcartesianpower.js | 54 ++------ .../test/test.dcartesianpower.native.js | 54 ++------ .../dcartesian-power/test/test.ndarray.js | 15 +++ .../test/test.ndarray.native.js | 15 +++ 21 files changed, 276 insertions(+), 284 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md index 659ccbc5e23e..f59376eb861d 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/README.md @@ -36,7 +36,7 @@ limitations under the License. var dcartesianPower = require( '@stdlib/blas/ext/base/dcartesian-power' ); ``` -#### dcartesianPower( N, k, x, strideX, out, strideOut1, strideOut2 ) +#### dcartesianPower( N, k, x, strideX, out, LDO ) Computes the Cartesian power for a double-precision floating-point strided array. @@ -46,7 +46,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var x = new Float64Array( [ 1.0, 2.0 ] ); var out = new Float64Array( 8 ); -dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); +dcartesianPower( x.length, 2, x, 1, out, 2 ); // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ``` @@ -55,10 +55,9 @@ The function has the following parameters: - **N**: number of indexed elements. - **k**: power. - **x**: input [`Float64Array`][@stdlib/array/float64]. -- **strideX**: stride length. +- **strideX**: stride length for `x`. - **out**: output [`Float64Array`][@stdlib/array/float64]. -- **strideOut1**: stride length between consecutive output pairs. -- **strideOut2**: stride length between elements within each output pair. +- **LDO**: stride length for the leading dimension of `out`. The `N`, `k`, and stride parameters determine which elements in the strided arrays are accessed at runtime. For example, to compute the Cartesian power of every other element: @@ -68,7 +67,7 @@ var Float64Array = require( '@stdlib/array/float64' ); var x = new Float64Array( [ 1.0, 0.0, 2.0, 0.0 ] ); var out = new Float64Array( 8 ); -dcartesianPower( 2, 2, x, 2, out, 2, 1 ); +dcartesianPower( 2, 2, x, 2, out, 2 ); // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ``` @@ -86,7 +85,7 @@ var out0 = new Float64Array( 8 ); // Create offset views... var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element -dcartesianPower( 2, 2, x1, 1, out0, 2, 1 ); +dcartesianPower( 2, 2, x1, 1, out0, 2 ); // out0 => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ``` @@ -110,8 +109,10 @@ dcartesianPower.ndarray( x.length, 2, x, 1, 0, out, 2, 1, 0 ); The function has the following additional parameters: -- **offsetX**: starting index. -- **offsetOut**: starting index. +- **offsetX**: starting index for `x`. +- **strideOut1**: stride length for the first dimension of `out`. +- **strideOut2**: stride length for the second dimension of `out`. +- **offsetOut**: starting index for `out`. 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 compute the Cartesian power of every other value in the strided input array starting from the second value: @@ -134,7 +135,7 @@ dcartesianPower.ndarray( 2, 2, x, 1, 1, out, 2, 1, 0 ); ## Notes - If `N <= 0` or `k <= 0`, both functions return `out` unchanged. -- The output array must contain at least `N^k` elements. +- The output array must contain at least `N^k * k` elements. @@ -155,7 +156,7 @@ console.log( x ); var out = new Float64Array( 24 ); -dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); +dcartesianPower( x.length, 3, x, 1, out, 3 ); console.log( out ); ``` @@ -191,7 +192,7 @@ console.log( out ); -#### stdlib_strided_dcartesian_power( N, k, \*X, strideX, \*Out, strideOut1, strideOut2 ) +#### stdlib_strided_dcartesian_power( N, k, \*X, strideX, \*Out, LDO ) @@ -201,7 +202,7 @@ Computes the Cartesian power for a double-precision floating-point strided array const double X[] = { 1.0, 2.0 }; double Out[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; -stdlib_strided_dcartesian_power( 2, 2, X, 1, Out, 2, 1 ); +stdlib_strided_dcartesian_power( 2, 2, X, 1, Out, 2 ); ``` The function accepts the following arguments: @@ -209,13 +210,12 @@ The function accepts the following arguments: - **N**: `[in] CBLAS_INT` number of indexed elements. - **k**: `[in] CBLAS_INT` power. - **X**: `[in] double*` input array. -- **strideX**: `[in] CBLAS_INT` stride length. +- **strideX**: `[in] CBLAS_INT` stride length for `X`. - **Out**: `[out] double*` output array. -- **strideOut1**: `[in] CBLAS_INT` stride length between consecutive output pairs. -- **strideOut2**: `[in] CBLAS_INT` stride length between elements within each output pair. +- **LDO**: `[in] CBLAS_INT` stride length for the leading dimension of `Out`. ```c -void stdlib_strided_dcartesian_power( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ); +void stdlib_strided_dcartesian_power( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT LDO ); ``` @@ -238,12 +238,12 @@ The function accepts the following arguments: - **N**: `[in] CBLAS_INT` number of indexed elements. - **k**: `[in] CBLAS_INT` power. - **X**: `[in] double*` input array. -- **strideX**: `[in] CBLAS_INT` stride length. -- **offsetX**: `[in] CBLAS_INT` starting index. +- **strideX**: `[in] CBLAS_INT` stride length for `X`. +- **offsetX**: `[in] CBLAS_INT` starting index for `X`. - **Out**: `[out] double*` output array. -- **strideOut1**: `[in] CBLAS_INT` stride length between consecutive output pairs. -- **strideOut2**: `[in] CBLAS_INT` stride length between elements within each output pair. -- **offsetOut**: `[in] CBLAS_INT` starting index. +- **strideOut1**: `[in] CBLAS_INT` stride length for the first dimension of `Out`. +- **strideOut2**: `[in] CBLAS_INT` stride length for the second dimension of `Out`. +- **offsetOut**: `[in] CBLAS_INT` starting index for `Out`. ```c void stdlib_strided_dcartesian_power_ndarray( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2, const CBLAS_INT offsetOut ); @@ -284,11 +284,10 @@ int main( void ) { // Specify strides: const int strideX = 1; - const int strideOut1 = 2; - const int strideOut2 = 1; + const int LDO = 2; // Compute the Cartesian power: - stdlib_strided_dcartesian_power( N, k, X, strideX, out, strideOut1, strideOut2 ); + stdlib_strided_dcartesian_power( N, k, X, strideX, out, LDO ); // Print the result: for ( int i = 0; i < N*N; i++ ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js index 05d48cd0efca..aadb18d6488b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.js @@ -67,7 +67,7 @@ function createBenchmark( len ) { b.tic(); for ( i = 0; i < b.iterations; i++ ) { - z = dcartesianPower( x.length, K, x, 1, out, K, 1 ); + z = dcartesianPower( x.length, K, x, 1, out, K ); if ( isnan( z[ 0 ] ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js index 535efd55f221..93cb0878df8e 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/benchmark.native.js @@ -72,7 +72,7 @@ function createBenchmark( len ) { b.tic(); for ( i = 0; i < b.iterations; i++ ) { - z = dcartesianPower( x.length, K, x, 1, out, K, 1 ); + z = dcartesianPower( x.length, K, x, 1, out, K ); if ( isnan( z[ 0 ] ) ) { b.fail( 'should not return NaN' ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c index 3a3901459268..03fa06f99060 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/benchmark/c/benchmark.length.c @@ -98,9 +98,9 @@ static double rand_double( void ) { static double benchmark1( int iterations, int len ) { double elapsed; double *out; + int outlen; double *x; double t; - int outlen; int i; x = (double *) malloc( len * sizeof( double ) ); @@ -114,7 +114,7 @@ static double benchmark1( int iterations, int len ) { } t = tic(); for ( i = 0; i < iterations; i++ ) { - stdlib_strided_dcartesian_power( len, K, x, 1, out, K, 1 ); + stdlib_strided_dcartesian_power( len, K, x, 1, out, K ); if ( out[ 0 ] != out[ 0 ] ) { printf( "should not return NaN\n" ); break; @@ -139,9 +139,9 @@ static double benchmark1( int iterations, int len ) { static double benchmark2( int iterations, int len ) { double elapsed; double *out; + int outlen; double *x; double t; - int outlen; int i; x = (double *) malloc( len * sizeof( double ) ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt index 8f036e35781f..bc84cd5b2105 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/repl.txt @@ -1,4 +1,4 @@ -{{alias}}( N, k, x, strideX, out, strideOut1, strideOut2 ) +{{alias}}( N, k, x, strideX, out, LDO ) Computes the Cartesian power for a double-precision floating-point strided array. @@ -19,16 +19,13 @@ Input array. strideX: integer - Stride length. + Stride length for `x`. out: Float64Array Output array. - strideOut1: integer - Stride length between consecutive output pairs. - - strideOut2: integer - Stride length between elements within each output pair. + LDO: integer + Stride length for the leading dimension of `out`. Returns ------- @@ -40,17 +37,17 @@ // Standard Usage: > var x = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0 ] ); > var out = new {{alias:@stdlib/array/float64}}( 8 ); - > {{alias}}( x.length, 2, x, 1, out, 2, 1 ) + > {{alias}}( x.length, 2, x, 1, out, 2 ) [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] // Using `N` and stride parameters: > x = new {{alias:@stdlib/array/float64}}( [ 1.0, 0.0, 2.0, 0.0 ] ); > out = new {{alias:@stdlib/array/float64}}( 8 ); - > {{alias}}( 2, 2, x, 2, out, 2, 1 ) + > {{alias}}( 2, 2, x, 2, out, 2 ) [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] -{{alias}}.ndarray( N, k, x, strideX, offsetX, out, sO1, sO2, offsetOut ) +{{alias}}.ndarray( N, k, x, strideX, offsetX, out, so1, so2, oo ) Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. @@ -70,22 +67,22 @@ Input array. strideX: integer - Stride length. + Stride length for `x`. offsetX: integer - Starting index. + Starting index for `x`. out: Float64Array Output array. - sO1: integer - Stride length between consecutive output pairs. + so1: integer + Stride length for the first dimension of `out`. - sO2: integer - Stride length between elements within each output pair. + so2: integer + Stride length for the second dimension of `out`. - offsetOut: integer - Starting index. + oo: integer + Starting index for `out`. Returns ------- diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts index 0240210bc806..68f3a01537c5 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/index.d.ts @@ -28,10 +28,9 @@ interface Routine { * @param N - number of indexed elements * @param k - power * @param x - input array - * @param strideX - stride length + * @param strideX - stride length for `x` * @param out - output array - * @param strideOut1 - stride length between consecutive output pairs - * @param strideOut2 - stride length between elements within each output pair + * @param LDO - stride length for the leading dimension of `out` * @returns output array * * @example @@ -40,10 +39,10 @@ interface Routine { * var x = new Float64Array( [ 1.0, 2.0 ] ); * var out = new Float64Array( 8 ); * - * dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + * dcartesianPower( x.length, 2, x, 1, out, 2 ); * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] */ - ( N: number, k: number, x: Float64Array, strideX: number, out: Float64Array, strideOut1: number, strideOut2: number ): Float64Array; + ( N: number, k: number, x: Float64Array, strideX: number, out: Float64Array, LDO: number ): Float64Array; /** * Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. @@ -51,12 +50,12 @@ interface Routine { * @param N - number of indexed elements * @param k - power * @param x - input array - * @param strideX - stride length - * @param offsetX - starting index + * @param strideX - stride length for `x` + * @param offsetX - starting index for `x` * @param out - output array - * @param strideOut1 - stride length between consecutive output pairs - * @param strideOut2 - stride length between elements within each output pair - * @param offsetOut - starting index + * @param strideOut1 - stride length for the first dimension of `out` + * @param strideOut2 - stride length for the second dimension of `out` + * @param offsetOut - starting index for `out` * @returns output array * * @example @@ -77,10 +76,9 @@ interface Routine { * @param N - number of indexed elements * @param k - power * @param x - input array -* @param strideX - stride length +* @param strideX - stride length for `x` * @param out - output array -* @param strideOut1 - stride length between consecutive output pairs -* @param strideOut2 - stride length between elements within each output pair +* @param LDO - stride length for the leading dimension of `out` * @returns output array * * @example @@ -89,7 +87,7 @@ interface Routine { * var x = new Float64Array( [ 1.0, 2.0 ] ); * var out = new Float64Array( 8 ); * -* dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); +* dcartesianPower( x.length, 2, x, 1, out, 2 ); * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] * * @example diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts index a3b37f63d693..5495fede9427 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/docs/types/test.ts @@ -26,7 +26,7 @@ import dcartesianPower = require( './index' ); const x = new Float64Array( 2 ); const out = new Float64Array( 8 ); - dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); // $ExpectType Float64Array + dcartesianPower( x.length, 2, x, 1, out, 2 ); // $ExpectType Float64Array } // The compiler throws an error if the function is provided a first argument which is not a number... @@ -34,14 +34,14 @@ import dcartesianPower = require( './index' ); const x = new Float64Array( 2 ); const out = new Float64Array( 8 ); - dcartesianPower( '10', 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( true, 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( false, 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( null, 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( undefined, 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( [], 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( {}, 2, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( ( x: number ): number => x, 2, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( '10', 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( true, 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( false, 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( null, 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( undefined, 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( [], 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( {}, 2, x, 1, out, 2 ); // $ExpectError + dcartesianPower( ( x: number ): number => x, 2, x, 1, out, 2 ); // $ExpectError } // The compiler throws an error if the function is provided a second argument which is not a number... @@ -49,29 +49,29 @@ import dcartesianPower = require( './index' ); const x = new Float64Array( 2 ); const out = new Float64Array( 8 ); - dcartesianPower( x.length, '10', x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, true, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, false, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, null, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, undefined, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, [], x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, {}, x, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, ( x: number ): number => x, x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, '10', x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, true, x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, false, x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, null, x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, undefined, x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, [], x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, {}, x, 1, out, 2 ); // $ExpectError + dcartesianPower( x.length, ( x: number ): number => x, x, 1, out, 2 ); // $ExpectError } // The compiler throws an error if the function is provided a third argument which is not a Float64Array... { const out = new Float64Array( 8 ); - dcartesianPower( 2, 2, 10, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, '10', 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, true, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, false, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, null, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, undefined, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, [], 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, {}, 1, out, 2, 1 ); // $ExpectError - dcartesianPower( 2, 2, ( x: number ): number => x, 1, out, 2, 1 ); // $ExpectError + dcartesianPower( 2, 2, 10, 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, '10', 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, true, 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, false, 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, null, 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, undefined, 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, [], 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, {}, 1, out, 2 ); // $ExpectError + dcartesianPower( 2, 2, ( x: number ): number => x, 1, out, 2 ); // $ExpectError } // The compiler throws an error if the function is provided a fourth argument which is not a number... @@ -79,29 +79,29 @@ import dcartesianPower = require( './index' ); const x = new Float64Array( 2 ); const out = new Float64Array( 8 ); - dcartesianPower( x.length, 2, x, '10', out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, true, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, false, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, null, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, undefined, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, [], out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, {}, out, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, ( x: number ): number => x, out, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, '10', out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, true, out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, false, out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, null, out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, undefined, out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, [], out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, {}, out, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, ( x: number ): number => x, out, 2 ); // $ExpectError } // The compiler throws an error if the function is provided a fifth argument which is not a Float64Array... { const x = new Float64Array( 2 ); - dcartesianPower( x.length, 2, x, 1, 10, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, '10', 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, true, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, false, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, null, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, undefined, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, [], 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, {}, 2, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, ( x: number ): number => x, 2, 1 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, 10, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, '10', 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, true, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, false, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, null, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, undefined, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, [], 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, {}, 2 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, ( x: number ): number => x, 2 ); // $ExpectError } // The compiler throws an error if the function is provided a sixth argument which is not a number... @@ -109,29 +109,14 @@ import dcartesianPower = require( './index' ); const x = new Float64Array( 2 ); const out = new Float64Array( 8 ); - dcartesianPower( x.length, 2, x, 1, out, '10', 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, true, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, false, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, null, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, undefined, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, [], 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, {}, 1 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, ( x: number ): number => x, 1 ); // $ExpectError -} - -// The compiler throws an error if the function is provided a seventh argument which is not a number... -{ - const x = new Float64Array( 2 ); - const out = new Float64Array( 8 ); - - dcartesianPower( x.length, 2, x, 1, out, 2, '10' ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, true ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, false ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, null ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, undefined ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, [] ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, {} ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, ( x: number ): number => x ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, '10' ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, true ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, false ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, null ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, undefined ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, [] ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, {} ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, ( x: number ): number => x ); // $ExpectError } // The compiler throws an error if the function is provided an unsupported number of arguments... @@ -145,8 +130,7 @@ import dcartesianPower = require( './index' ); dcartesianPower( x.length, 2, x ); // $ExpectError dcartesianPower( x.length, 2, x, 1 ); // $ExpectError dcartesianPower( x.length, 2, x, 1, out ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2 ); // $ExpectError - dcartesianPower( x.length, 2, x, 1, out, 2, 1, 10 ); // $ExpectError + dcartesianPower( x.length, 2, x, 1, out, 2, 10 ); // $ExpectError } // The `ndarray` method returns a Float64Array... diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c index 63c4351738b5..fb32b7a21048 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/c/example.c @@ -32,11 +32,10 @@ int main( void ) { // Specify strides: const int strideX = 1; - const int strideOut1 = 2; - const int strideOut2 = 1; + const int LDO = 2; // Compute the Cartesian power: - stdlib_strided_dcartesian_power( N, k, X, strideX, out, strideOut1, strideOut2 ); + stdlib_strided_dcartesian_power( N, k, X, strideX, out, LDO ); // Print the result: for ( int i = 0; i < N*N; i++ ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js index 6720a45a2655..9d3580428548 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/examples/index.js @@ -26,5 +26,5 @@ console.log( x ); var out = new Float64Array( 24 ); -dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); +dcartesianPower( x.length, 3, x, 1, out, 3 ); console.log( out ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include/stdlib/blas/ext/base/dcartesianpower.h b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include/stdlib/blas/ext/base/dcartesianpower.h index 48c799b88db4..cc568df6794b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include/stdlib/blas/ext/base/dcartesianpower.h +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/include/stdlib/blas/ext/base/dcartesianpower.h @@ -31,7 +31,7 @@ extern "C" { /** * Computes the Cartesian power for a double-precision floating-point strided array. */ -void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ); +void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT LDO ); /** * Computes the Cartesian power for a double-precision floating-point strided array using alternative indexing semantics. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.js index 886b57e26b5d..256837234228 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.js @@ -33,10 +33,9 @@ var ndarray = require( './ndarray.js' ); * @param {NonNegativeInteger} N - number of indexed elements * @param {NonNegativeInteger} k - power * @param {Float64Array} x - input array -* @param {integer} strideX - stride length +* @param {integer} strideX - stride length for `x` * @param {Float64Array} out - output array -* @param {integer} strideOut1 - stride length between consecutive output pairs -* @param {integer} strideOut2 - stride length between elements within each output pair +* @param {integer} LDO - stride length for the leading dimension of `out` * @returns {Float64Array} output array * * @example @@ -45,14 +44,18 @@ var ndarray = require( './ndarray.js' ); * var x = new Float64Array( [ 1.0, 2.0 ] ); * var out = new Float64Array( 8 ); * -* dcartesianPower( 2, 2, x, 1, out, 2, 1 ); +* dcartesianPower( 2, 2, x, 1, out, 2 ); * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] */ -function dcartesianPower( N, k, x, strideX, out, strideOut1, strideOut2 ) { - var len = pow( N, k ); - var ox = stride2offset( N, strideX ); - var oo = stride2offset( len, strideOut1 ) + stride2offset( k, strideOut2 ); - return ndarray( N, k, x, strideX, ox, out, strideOut1, strideOut2, oo ); +function dcartesianPower( N, k, x, strideX, out, LDO ) { + var len; + var ox; + var oo; + + len = pow( N, k ); + ox = stride2offset( N, strideX ); + oo = stride2offset( len, LDO ); + return ndarray( N, k, x, strideX, ox, out, LDO, 1, oo ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js index 836afe675913..f9c2be9c4407 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/dcartesianpower.native.js @@ -33,8 +33,7 @@ var addon = require( './../src/addon.node' ); * @param {Float64Array} x - input array * @param {integer} strideX - stride length * @param {Float64Array} out - output array -* @param {integer} strideOut1 - stride length between consecutive output pairs -* @param {integer} strideOut2 - stride length between elements within each output pair +* @param {integer} LDO - leading dimension of the output array * @returns {Float64Array} output array * * @example @@ -43,11 +42,11 @@ var addon = require( './../src/addon.node' ); * var x = new Float64Array( [ 1.0, 2.0 ] ); * var out = new Float64Array( 8 ); * -* dcartesianPower( 2, 2, x, 1, out, 2, 1 ); +* dcartesianPower( 2, 2, x, 1, out, 2 ); * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] */ -function dcartesianPower( N, k, x, strideX, out, strideOut1, strideOut2 ) { - addon( N, k, x, strideX, out, strideOut1, strideOut2 ); +function dcartesianPower( N, k, x, strideX, out, LDO ) { + addon( N, k, x, strideX, out, LDO ); return out; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js index 302d47ad7a08..42f13a57d49f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/index.js @@ -30,7 +30,7 @@ * var x = new Float64Array( [ 1.0, 2.0 ] ); * var out = new Float64Array( 8 ); * -* dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); +* dcartesianPower( x.length, 2, x, 1, out, 2 ); * // out => [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] * * @example diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js index 9c7fa7b2fad8..9796f890554b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/lib/ndarray.js @@ -20,6 +20,7 @@ // MODULES // +var isRowMajor = require( '@stdlib/ndarray/base/assert/is-row-major' ); var pow = require( '@stdlib/math/base/special/pow' ); @@ -31,12 +32,12 @@ var pow = require( '@stdlib/math/base/special/pow' ); * @param {NonNegativeInteger} N - number of indexed elements * @param {NonNegativeInteger} k - power * @param {Float64Array} x - input array -* @param {integer} strideX - stride length -* @param {NonNegativeInteger} offsetX - starting index +* @param {integer} strideX - stride length for `x` +* @param {NonNegativeInteger} offsetX - starting index for `x` * @param {Float64Array} out - output array -* @param {integer} strideOut1 - stride length between consecutive output pairs -* @param {integer} strideOut2 - stride length between elements within each output pair -* @param {NonNegativeInteger} offsetOut - starting index +* @param {integer} strideOut1 - stride length for the first dimension of `out` +* @param {integer} strideOut2 - stride length for the second dimension of `out` +* @param {NonNegativeInteger} offsetOut - starting index for `out` * @returns {Float64Array} output array * * @example @@ -53,6 +54,7 @@ function dcartesianPower( N, k, x, strideX, offsetX, out, strideOut1, strideOut2 var oj; var io; var ix; + var sj; var i; var j; var s; @@ -62,17 +64,33 @@ function dcartesianPower( N, k, x, strideX, offsetX, out, strideOut1, strideOut2 return out; } len = pow( N, k ); - io = offsetOut; - for ( i = 0; i < len; i++ ) { - t = i; - for ( j = k - 1; j >= 0; j-- ) { - s = t % N; - t = ( t - s ) / N; - ix = offsetX + ( s * strideX ); - oj = io + ( j * strideOut2 ); - out[ oj ] = x[ ix ]; + if ( isRowMajor( [ strideOut1, strideOut2 ] ) ) { + // Row-major: outer loop over tuples, inner loop over tuple elements... + io = offsetOut; + for ( i = 0; i < len; i++ ) { + t = i; + for ( j = k - 1; j >= 0; j-- ) { + s = t % N; + t = ( t - s ) / N; + ix = offsetX + ( s * strideX ); + oj = io + ( j * strideOut2 ); + out[ oj ] = x[ ix ]; + } + io += strideOut1; + } + } else { + // Column-major: outer loop over tuple elements, inner loop over tuples... + sj = len; + for ( j = 0; j < k; j++ ) { + sj = sj / N; + io = offsetOut + ( j * strideOut2 ); + for ( i = 0; i < len; i++ ) { + s = ( ( i - ( i % sj ) ) / sj ) % N; + ix = offsetX + ( s * strideX ); + out[ io ] = x[ ix ]; + io += strideOut1; + } } - io += strideOut1; } return out; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json index a5be842aa2a3..7ddbc0d9fd9a 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/manifest.json @@ -38,6 +38,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/math/base/special/pow", + "@stdlib/ndarray/base/assert/is-row-major", "@stdlib/strided/base/stride2offset", "@stdlib/napi/export", "@stdlib/napi/argv", @@ -59,6 +60,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/math/base/special/pow", + "@stdlib/ndarray/base/assert/is-row-major", "@stdlib/strided/base/stride2offset" ] }, @@ -75,6 +77,7 @@ "dependencies": [ "@stdlib/blas/base/shared", "@stdlib/math/base/special/pow", + "@stdlib/ndarray/base/assert/is-row-major", "@stdlib/strided/base/stride2offset" ] } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c index a9a1e7ac2bef..8b6c5185c148 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/addon.c @@ -25,6 +25,7 @@ #include "stdlib/napi/argv_strided_float64array.h" #include "stdlib/napi/argv_strided_float64array2d.h" #include +#include /** * Receives JavaScript callback invocation data. @@ -34,15 +35,15 @@ * @return Node-API value */ static napi_value addon( napi_env env, napi_callback_info info ) { - STDLIB_NAPI_ARGV( env, info, argv, argc, 7 ); + STDLIB_NAPI_ARGV( env, info, argv, argc, 6 ); STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); STDLIB_NAPI_ARGV_INT64( env, k, argv, 1 ); STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); - STDLIB_NAPI_ARGV_INT64( env, strideOut1, argv, 5 ); - STDLIB_NAPI_ARGV_INT64( env, strideOut2, argv, 6 ); + STDLIB_NAPI_ARGV_INT64( env, LDO, argv, 5 ); STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); - STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY2D( env, Out, stdlib_base_pow( N, k ), k, strideOut1, strideOut2, argv, 4 ); - API_SUFFIX(stdlib_strided_dcartesian_power)( N, k, X, strideX, Out, strideOut1, strideOut2 ); + int64_t len = stdlib_base_pow( N, k ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY2D( env, Out, len, k, LDO, 1, argv, 4 ); + API_SUFFIX(stdlib_strided_dcartesian_power)( N, k, X, strideX, Out, LDO ); return NULL; } @@ -63,7 +64,8 @@ static napi_value addon_method( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_INT64( env, strideOut2, argv, 7 ); STDLIB_NAPI_ARGV_INT64( env, offsetOut, argv, 8 ); STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); - STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY2D( env, Out, stdlib_base_pow( N, k ), k, strideOut1, strideOut2, argv, 5 ); + int64_t len = stdlib_base_pow( N, k ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY2D( env, Out, len, k, strideOut1, strideOut2, argv, 5 ); API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, offsetX, Out, strideOut1, strideOut2, offsetOut ); return NULL; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c index 974eff4b0a4f..303414b9715a 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/src/main.c @@ -17,6 +17,7 @@ */ #include "stdlib/blas/ext/base/dcartesianpower.h" +#include "stdlib/ndarray/base/assert/is_row_major.h" #include "stdlib/math/base/special/pow.h" #include "stdlib/blas/base/shared.h" #include "stdlib/strided/base/stride2offset.h" @@ -25,19 +26,22 @@ /** * Computes the Cartesian power for a double-precision floating-point strided array. * -* @param N number of indexed elements -* @param k power -* @param X input array -* @param strideX stride length -* @param Out output array -* @param strideOut1 stride length between consecutive output pairs -* @param strideOut2 stride length between elements within each output pair +* @param N number of indexed elements +* @param k power +* @param X input array +* @param strideX stride length for X +* @param Out output array +* @param LDO stride length for the leading dimension of Out */ -void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2 ) { - int64_t len = stdlib_base_pow( N, k ); - CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); - CBLAS_INT oo = stdlib_strided_stride2offset( len, strideOut1 ) + stdlib_strided_stride2offset( k, strideOut2 ); - API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, ox, Out, strideOut1, strideOut2, oo ); +void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, double *Out, const CBLAS_INT LDO ) { + CBLAS_INT ox; + CBLAS_INT oo; + int64_t len; + + len = stdlib_base_pow( N, k ); + ox = stdlib_strided_stride2offset( N, strideX ); + oo = stdlib_strided_stride2offset( len, LDO ); + API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( N, k, X, strideX, ox, Out, LDO, 1, oo ); } /** @@ -46,20 +50,23 @@ void API_SUFFIX(stdlib_strided_dcartesian_power)( const CBLAS_INT N, const CBLAS * @param N number of indexed elements * @param k power * @param X input array -* @param strideX stride length -* @param offsetX starting index +* @param strideX stride length for X +* @param offsetX starting index for X * @param Out output array -* @param strideOut1 stride length between consecutive output pairs -* @param strideOut2 stride length between elements within each output pair -* @param offsetOut starting index +* @param strideOut1 stride length for the first dimension of Out +* @param strideOut2 stride length for the second dimension of Out +* @param offsetOut starting index for Out */ void API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( const CBLAS_INT N, const CBLAS_INT k, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Out, const CBLAS_INT strideOut1, const CBLAS_INT strideOut2, const CBLAS_INT offsetOut ) { + int64_t strides[ 2 ]; + CBLAS_INT isrm; CBLAS_INT ix; CBLAS_INT j; CBLAS_INT s; int64_t len; int64_t io; int64_t oj; + int64_t sj; int64_t i; int64_t t; @@ -67,16 +74,33 @@ void API_SUFFIX(stdlib_strided_dcartesian_power_ndarray)( const CBLAS_INT N, con return; } len = stdlib_base_pow( N, k ); - io = offsetOut; - for ( i = 0; i < len; i++ ) { - t = i; - for ( j = k - 1; j >= 0; j-- ) { - s = t % N; - t = ( t - s ) / N; - ix = offsetX + ( s * strideX ); - oj = io + ( j * strideOut2 ); - Out[ oj ] = X[ ix ]; + strides[ 0 ] = strideOut1; + strides[ 1 ] = strideOut2; + isrm = stdlib_ndarray_is_row_major( 2, strides ); + if ( isrm ) { + io = offsetOut; + for ( i = 0; i < len; i++ ) { + t = i; + for ( j = k - 1; j >= 0; j-- ) { + s = t % N; + t = ( t - s ) / N; + ix = offsetX + ( s * strideX ); + oj = io + ( j * strideOut2 ); + Out[ oj ] = X[ ix ]; + } + io += strideOut1; + } + } else { + sj = len; + for ( j = 0; j < k; j++ ) { + sj = sj / N; + io = offsetOut + ( j * strideOut2 ); + for ( i = 0; i < len; i++ ) { + s = ( ( i - ( i % sj ) ) / sj ) % N; + ix = offsetX + ( s * strideX ); + Out[ io ] = X[ ix ]; + io += strideOut1; + } } - io += strideOut1; } } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js index 7b168eb32cc2..7030b192c907 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.js @@ -33,8 +33,8 @@ tape( 'main export is a function', function test( t ) { t.end(); }); -tape( 'the function has an arity of 7', function test( t ) { - t.strictEqual( dcartesianPower.length, 7, 'returns expected value' ); +tape( 'the function has an arity of 6', function test( t ) { + t.strictEqual( dcartesianPower.length, 6, 'returns expected value' ); t.end(); }); @@ -46,14 +46,14 @@ tape( 'the function computes the Cartesian power', function test( t ) { x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 8 ); - dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + dcartesianPower( x.length, 2, x, 1, out, 2 ); expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); t.deepEqual( out, expected, 'returns expected value' ); x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 24 ); - dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); + dcartesianPower( x.length, 3, x, 1, out, 3 ); expected = new Float64Array([ 1.0, 1.0, @@ -93,7 +93,7 @@ tape( 'the function returns a reference to the output array', function test( t ) x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 8 ); - v = dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + v = dcartesianPower( x.length, 2, x, 1, out, 2 ); t.strictEqual( v, out, 'same reference' ); t.end(); @@ -108,10 +108,10 @@ tape( 'if provided an `N` or `k` parameter equal to `0`, the function returns `o out = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); expected = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); - dcartesianPower( 0, 2, x, 1, out, 2, 1 ); + dcartesianPower( 0, 2, x, 1, out, 2 ); t.deepEqual( out, expected, 'returns `out` unchanged' ); - dcartesianPower( 2, 0, x, 1, out, 2, 1 ); + dcartesianPower( 2, 0, x, 1, out, 2 ); t.deepEqual( out, expected, 'returns `out` unchanged' ); t.end(); @@ -130,7 +130,7 @@ tape( 'the function supports a strideX parameter', function test( t ) { ]); out = new Float64Array( 8 ); - dcartesianPower( 2, 2, x, 2, out, 2, 1 ); + dcartesianPower( 2, 2, x, 2, out, 2 ); expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); t.deepEqual( out, expected, 'returns expected value' ); @@ -145,7 +145,7 @@ tape( 'the function supports a strideOut1 parameter', function test( t ) { x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 16 ); - dcartesianPower( 2, 2, x, 1, out, 4, 1 ); + dcartesianPower( 2, 2, x, 1, out, 4 ); expected = new Float64Array([ 1.0, @@ -169,38 +169,6 @@ tape( 'the function supports a strideOut1 parameter', function test( t ) { t.end(); }); -tape( 'the function supports a strideOut2 parameter', function test( t ) { - var expected; - var out; - var x; - - x = new Float64Array( [ 1.0, 2.0 ] ); - out = new Float64Array( 16 ); - - dcartesianPower( 2, 2, x, 1, out, 4, 2 ); - - expected = new Float64Array([ - 1.0, - 0.0, - 1.0, - 0.0, - 1.0, - 0.0, - 2.0, - 0.0, - 2.0, - 0.0, - 1.0, - 0.0, - 2.0, - 0.0, - 2.0, - 0.0 - ]); - t.deepEqual( out, expected, 'returns expected value' ); - t.end(); -}); - tape( 'the function supports negative strides', function test( t ) { var expected; var out; @@ -212,7 +180,7 @@ tape( 'the function supports negative strides', function test( t ) { ]); out = new Float64Array( 8 ); - dcartesianPower( 2, 2, x, -1, out, -2, 1 ); + dcartesianPower( 2, 2, x, -1, out, -2 ); expected = new Float64Array( [ 2.0, 2.0, 2.0, 1.0, 1.0, 2.0, 1.0, 1.0 ] ); t.deepEqual( out, expected, 'returns expected value' ); @@ -231,7 +199,7 @@ tape( 'the function supports view offsets', function test( t ) { // Create an offset view... x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element - dcartesianPower( 2, 2, x1, 1, out0, 2, 1 ); + dcartesianPower( 2, 2, x1, 1, out0, 2 ); expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); t.deepEqual( out0, expected, 'returns expected value' ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js index 5c9a01ebf48e..8db8bfd9c125 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.dcartesianpower.native.js @@ -42,8 +42,8 @@ tape( 'main export is a function', opts, function test( t ) { t.end(); }); -tape( 'the function has an arity of 7', opts, function test( t ) { - t.strictEqual( dcartesianPower.length, 7, 'returns expected value' ); +tape( 'the function has an arity of 6', opts, function test( t ) { + t.strictEqual( dcartesianPower.length, 6, 'returns expected value' ); t.end(); }); @@ -55,14 +55,14 @@ tape( 'the function computes the Cartesian power', opts, function test( t ) { x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 8 ); - dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + dcartesianPower( x.length, 2, x, 1, out, 2 ); expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); t.deepEqual( out, expected, 'returns expected value' ); x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 24 ); - dcartesianPower( x.length, 3, x, 1, out, 3, 1 ); + dcartesianPower( x.length, 3, x, 1, out, 3 ); expected = new Float64Array([ 1.0, 1.0, @@ -102,7 +102,7 @@ tape( 'the function returns a reference to the output array', opts, function tes x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 8 ); - v = dcartesianPower( x.length, 2, x, 1, out, 2, 1 ); + v = dcartesianPower( x.length, 2, x, 1, out, 2 ); t.strictEqual( v, out, 'same reference' ); t.end(); @@ -117,10 +117,10 @@ tape( 'if provided an `N` or `k` parameter equal to `0`, the function returns `o out = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); expected = new Float64Array( [ 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 ] ); - dcartesianPower( 0, 2, x, 1, out, 2, 1 ); + dcartesianPower( 0, 2, x, 1, out, 2 ); t.deepEqual( out, expected, 'returns `out` unchanged' ); - dcartesianPower( 2, 0, x, 1, out, 2, 1 ); + dcartesianPower( 2, 0, x, 1, out, 2 ); t.deepEqual( out, expected, 'returns `out` unchanged' ); t.end(); @@ -139,7 +139,7 @@ tape( 'the function supports a strideX parameter', opts, function test( t ) { ]); out = new Float64Array( 8 ); - dcartesianPower( 2, 2, x, 2, out, 2, 1 ); + dcartesianPower( 2, 2, x, 2, out, 2 ); expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); t.deepEqual( out, expected, 'returns expected value' ); @@ -154,7 +154,7 @@ tape( 'the function supports a strideOut1 parameter', opts, function test( t ) { x = new Float64Array( [ 1.0, 2.0 ] ); out = new Float64Array( 16 ); - dcartesianPower( 2, 2, x, 1, out, 4, 1 ); + dcartesianPower( 2, 2, x, 1, out, 4 ); expected = new Float64Array([ 1.0, @@ -178,38 +178,6 @@ tape( 'the function supports a strideOut1 parameter', opts, function test( t ) { t.end(); }); -tape( 'the function supports a strideOut2 parameter', opts, function test( t ) { - var expected; - var out; - var x; - - x = new Float64Array( [ 1.0, 2.0 ] ); - out = new Float64Array( 16 ); - - dcartesianPower( 2, 2, x, 1, out, 4, 2 ); - - expected = new Float64Array([ - 1.0, - 0.0, - 1.0, - 0.0, - 1.0, - 0.0, - 2.0, - 0.0, - 2.0, - 0.0, - 1.0, - 0.0, - 2.0, - 0.0, - 2.0, - 0.0 - ]); - t.deepEqual( out, expected, 'returns expected value' ); - t.end(); -}); - tape( 'the function supports negative strides', opts, function test( t ) { var expected; var out; @@ -221,7 +189,7 @@ tape( 'the function supports negative strides', opts, function test( t ) { ]); out = new Float64Array( 8 ); - dcartesianPower( 2, 2, x, -1, out, -2, 1 ); + dcartesianPower( 2, 2, x, -1, out, -2 ); expected = new Float64Array( [ 2.0, 2.0, 2.0, 1.0, 1.0, 2.0, 1.0, 1.0 ] ); t.deepEqual( out, expected, 'returns expected value' ); @@ -240,7 +208,7 @@ tape( 'the function supports view offsets', opts, function test( t ) { // Create an offset view... x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element - dcartesianPower( 2, 2, x1, 1, out0, 2, 1 ); + dcartesianPower( 2, 2, x1, 1, out0, 2 ); expected = new Float64Array( [ 1.0, 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 2.0 ] ); t.deepEqual( out0, expected, 'returns expected value' ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js index 8f4c354f4fbe..3e526ce00afa 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.js @@ -201,6 +201,21 @@ tape( 'the function supports a strideOut2 parameter', function test( t ) { t.end(); }); +tape( 'the function supports a column-major memory layout', function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 1, 4, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 1.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + tape( 'the function supports negative strides', function test( t ) { var expected; var out; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js index e7a490c7fab2..70989cc0fa55 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dcartesian-power/test/test.ndarray.native.js @@ -210,6 +210,21 @@ tape( 'the function supports a strideOut2 parameter', opts, function test( t ) { t.end(); }); +tape( 'the function supports a column-major memory layout', opts, function test( t ) { + var expected; + var out; + var x; + + x = new Float64Array( [ 1.0, 2.0 ] ); + out = new Float64Array( 8 ); + + dcartesianPower( 2, 2, x, 1, 0, out, 1, 4, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 2.0, 2.0, 1.0, 2.0, 1.0, 2.0 ] ); + t.deepEqual( out, expected, 'returns expected value' ); + t.end(); +}); + tape( 'the function supports negative strides', opts, function test( t ) { var expected; var out;