From da4326ba5c066145e0672b07e01a253f871c911c Mon Sep 17 00:00:00 2001 From: Om-A-osc Date: Mon, 2 Mar 2026 01:02:56 +0530 Subject: [PATCH 1/2] feat: add `stats/strided/dnancumaxabs` --- 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: passed - task: lint_c_examples status: passed - task: lint_c_benchmarks status: passed - 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 --- --- .../stats/strided/dnancumaxabs/README.md | 345 +++++++++++++++++ .../dnancumaxabs/benchmark/benchmark.js | 112 ++++++ .../benchmark/benchmark.native.js | 121 ++++++ .../benchmark/benchmark.ndarray.js | 112 ++++++ .../benchmark/benchmark.ndarray.native.js | 121 ++++++ .../strided/dnancumaxabs/benchmark/c/Makefile | 146 ++++++++ .../benchmark/c/benchmark.length.c | 213 +++++++++++ .../stats/strided/dnancumaxabs/binding.gyp | 170 +++++++++ .../stats/strided/dnancumaxabs/docs/repl.txt | 113 ++++++ .../dnancumaxabs/docs/types/index.d.ts | 94 +++++ .../strided/dnancumaxabs/docs/types/test.ts | 248 +++++++++++++ .../strided/dnancumaxabs/examples/c/Makefile | 146 ++++++++ .../strided/dnancumaxabs/examples/c/example.c | 41 ++ .../strided/dnancumaxabs/examples/index.js | 45 +++ .../stats/strided/dnancumaxabs/include.gypi | 53 +++ .../stdlib/stats/strided/dnancumaxabs.h | 45 +++ .../strided/dnancumaxabs/lib/dnancumaxabs.js | 58 +++ .../dnancumaxabs/lib/dnancumaxabs.native.js | 55 +++ .../stats/strided/dnancumaxabs/lib/index.js | 70 ++++ .../stats/strided/dnancumaxabs/lib/main.js | 35 ++ .../stats/strided/dnancumaxabs/lib/native.js | 35 ++ .../stats/strided/dnancumaxabs/lib/ndarray.js | 100 +++++ .../dnancumaxabs/lib/ndarray.native.js | 57 +++ .../stats/strided/dnancumaxabs/manifest.json | 137 +++++++ .../stats/strided/dnancumaxabs/package.json | 84 +++++ .../stats/strided/dnancumaxabs/src/Makefile | 70 ++++ .../stats/strided/dnancumaxabs/src/addon.c | 65 ++++ .../stats/strided/dnancumaxabs/src/main.c | 97 +++++ .../dnancumaxabs/test/test.dnancumaxabs.js | 273 ++++++++++++++ .../test/test.dnancumaxabs.native.js | 282 ++++++++++++++ .../stats/strided/dnancumaxabs/test/test.js | 38 ++ .../strided/dnancumaxabs/test/test.ndarray.js | 341 +++++++++++++++++ .../dnancumaxabs/test/test.ndarray.native.js | 350 ++++++++++++++++++ 33 files changed, 4272 insertions(+) create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.native.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/c/benchmark.length.c create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/binding.gyp create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/index.d.ts create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/test.ts create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/index.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/include.gypi create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/include/stdlib/stats/strided/dnancumaxabs.h create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/dnancumaxabs.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/dnancumaxabs.native.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/index.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/main.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/native.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.native.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/manifest.json create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/package.json create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/Makefile create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/addon.c create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/main.c create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.native.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.js create mode 100644 lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.native.js diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md new file mode 100644 index 000000000000..fd2598315e31 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md @@ -0,0 +1,345 @@ + + + + +# dnancumaxabs + +> Calculate the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. + +
+ +
+ + + +
+ +## Usage + +```javascript +var dnancumaxabs = require( '@stdlib/stats/strided/dnancumaxabs' ); +``` + +#### dnancumaxabs( N, x, strideX, y, strideY ) + +Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, NaN, -2.0, 3.0, NaN, 2.0 ] ); +var y = new Float64Array( x.length ); + +dnancumaxabs( x.length, x, 1, y, 1 ); +// y => [ 1.0, 1.0, 2.0, 3.0, 3.0, 3.0 ] +``` + +The function has the following parameters: + +- **N**: number of indexed elements. +- **x**: input [`Float64Array`][@stdlib/array/float64]. +- **strideX**: stride length for `x`. +- **y**: output [`Float64Array`][@stdlib/array/float64]. +- **strideY**: stride length for `y`. + +The `N` and stride parameters determine which elements in the strided arrays are accessed at runtime. For example, to compute the cumulative maximum absolute value of every other element in `x`, + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, 2.0, -3.0, 4.0, 5.0, 6.0 ] ); +var y = new Float64Array( x.length ); + +dnancumaxabs( 3, x, 2, y, 1 ); +// y => [ 1.0, 3.0, 5.0, 0.0, 0.0, 0.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( [ 1.0, NaN, -2.0, 3.0, 2.0, 5.0 ] ); +var y0 = new Float64Array( x0.length ); + +// Create offset views... +var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element +var y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // start at 4th element + +dnancumaxabs( 3, x1, 2, y1, 1 ); +// y0 => [ 0.0, 0.0, 0.0, NaN, 3.0, 5.0 ] +``` + +#### dnancumaxabs.ndarray( N, x, strideX, offsetX, y, strideY, offsetY ) + +Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values, using alternative indexing semantics. + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, NaN, -2.0, 3.0, NaN, 2.0 ] ); +var y = new Float64Array( x.length ); + +dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, 0 ); +// y => [ 1.0, 1.0, 2.0, 3.0, 3.0, 3.0 ] +``` + +The function has the following additional parameters: + +- **offsetX**: starting index for `x`. +- **offsetY**: starting index for `y`. + +While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example, to calculate the cumulative maximum absolute value of every other element in `x` starting from the second element and to store in the last `N` elements of `y` starting from the last element + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, NaN, -2.0, 3.0, 2.0, 5.0, -1.0 ] ); +var y = new Float64Array( x.length ); + +dnancumaxabs.ndarray( 3, x, 2, 1, y, -1, y.length-1 ); +// y => [ 0.0, 0.0, 0.0, 0.0, 5.0, 3.0, NaN ] +``` + +
+ + + +
+ +## Notes + +- If `N <= 0`, both functions return `y` unchanged. + +
+ + + +
+ +## Examples + + + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var dnancumaxabs = require( '@stdlib/stats/strided/dnancumaxabs' ); + +function rand() { + if ( bernoulli( 0.8 ) < 1 ) { + return NaN; + } + return uniform( -50.0, 50.0 ); +} + +var x = filledarrayBy( 10, 'float64', rand ); +console.log( x ); + +var y = new Float64Array( x.length ); +dnancumaxabs( x.length, x, 1, y, 1 ); +console.log( y ); + +dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, 0 ); +console.log( y ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/stats/strided/dnancumaxabs.h" +``` + +#### stdlib_strided_dnancumaxabs( N, \*X, strideX, \*Y, strideY ) + +Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. + +```c +const double x[] = { 1.0, 0.0/0.0, -2.0, 3.0, 0.0/0.0, 2.0 }; +double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + +stdlib_strided_dnancumaxabs( 6, x, 1, y, 1 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **X**: `[in] double*` input array. +- **strideX**: `[in] CBLAS_INT` stride length for `X`. +- **Y**: `[out] double*` output array. +- **strideY**: `[in] CBLAS_INT` stride length for `Y`. + +```c +void stdlib_strided_dnancumaxabs( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, double *Y, const CBLAS_INT strideY ); +``` + +#### stdlib_strided_dnancumaxabs_ndarray( N, \*X, strideX, offsetX, \*Y, strideY, offsetY ) + +Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values, using alternative indexing semantics. + +```c +const double x[] = { 1.0, 0.0/0.0, -2.0, 3.0, 0.0/0.0, 2.0 }; +double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + +stdlib_strided_dnancumaxabs_ndarray( 6, x, 1, 0, y, 1, 0 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **X**: `[in] double*` input array. +- **strideX**: `[in] CBLAS_INT` stride length for `X`. +- **offsetX**: `[in] CBLAS_INT` starting index for `X`. +- **Y**: `[out] double*` output array. +- **strideY**: `[in] CBLAS_INT` stride length for `Y`. +- **offsetY**: `[in] CBLAS_INT` starting index for `Y`. + +```c +void stdlib_strided_dnancumaxabs_ndarray( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/stats/strided/dnancumaxabs.h" +#include + +int main( void ) { + // Create strided arrays: + const double x[] = { 1.0, 0.0, 2.0, 0.0, 0.0/0.0, 0.0, -3.0, 0.0, 4.0, 0.0, -5.0, 0.0, 6.0, 0.0, 7.0, 0.0 }; + double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + + // Specify the number of elements: + const int N = 8; + + // Specify stride lengths: + const int strideX = 2; + const int strideY = 2; + + // Compute the cumulative maximum absolute value: + stdlib_strided_dnancumaxabs( N, x, strideX, y, strideY ); + + // Print the result: + for ( int i = 0; i < 16; i++ ) { + printf( "y[ %d ] = %lf\n", i, y[ i ] ); + } +} +``` + +
+ + + +
+ + + +
+ +
+ + + + + + + + + + + + + + diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.js new file mode 100644 index 000000000000..f816b6b3b116 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.js @@ -0,0 +1,112 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var format = require( '@stdlib/string/format' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var pkg = require( './../package.json' ).name; +var dnancumaxabs = require( './../lib' ); + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number or `NaN` +*/ +function rand() { + if ( bernoulli( 0.8 ) < 1 ) { + return NaN; + } + return uniform( -10.0, 10.0 ); +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, 'float64', rand ); + var y = new Float64Array( len ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + dnancumaxabs( x.length, x, 1, y, 1 ); + if ( isnan( y[ i % x.length ] ) && !isnan( x[ i % x.length ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y[ 0 ] ) && !isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.native.js new file mode 100644 index 000000000000..c5987302f2ea --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.native.js @@ -0,0 +1,121 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var format = require( '@stdlib/string/format' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var dnancumaxabs = tryRequire( resolve( __dirname, './../lib/dnancumaxabs.native.js' ) ); +var opts = { + 'skip': ( dnancumaxabs instanceof Error ) +}; + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number or `NaN` +*/ +function rand() { + if ( bernoulli( 0.8 ) < 1 ) { + return NaN; + } + return uniform( -10.0, 10.0 ); +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, 'float64', rand ); + var y = new Float64Array( len ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + dnancumaxabs( x.length, x, 1, y, 1 ); + if ( isnan( y[ i % x.length ] ) && !isnan( x[ i % x.length ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y[ 0 ] ) && !isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s::native:len=%d', pkg, len ), opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.js new file mode 100644 index 000000000000..399d324055a6 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.js @@ -0,0 +1,112 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var format = require( '@stdlib/string/format' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var pkg = require( './../package.json' ).name; +var dnancumaxabs = require( './../lib/ndarray.js' ); + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number or `NaN` +*/ +function rand() { + if ( bernoulli( 0.8 ) < 1 ) { + return NaN; + } + return uniform( -10.0, 10.0 ); +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, 'float64', rand ); + var y = new Float64Array( len ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + if ( isnan( y[ i % x.length ] ) && !isnan( x[ i % x.length ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y[ 0 ] ) && !isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s:ndarray:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.native.js new file mode 100644 index 000000000000..506e8d5e87dc --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/benchmark.ndarray.native.js @@ -0,0 +1,121 @@ +/** +* @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 Float64Array = require( '@stdlib/array/float64' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var format = require( '@stdlib/string/format' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var dnancumaxabs = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( dnancumaxabs instanceof Error ) +}; + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number or `NaN` +*/ +function rand() { + if ( bernoulli( 0.8 ) < 1 ) { + return NaN; + } + return uniform( -10.0, 10.0 ); +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, 'float64', rand ); + var y = new Float64Array( len ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + if ( isnan( y[ i % x.length ] ) && !isnan( x[ i % x.length ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( y[ 0 ] ) && !isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( format( '%s::native:ndarray:len=%d', pkg, len ), opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/c/Makefile b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/c/Makefile new file mode 100644 index 000000000000..0756dc7da20a --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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/stats/strided/dnancumaxabs/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/c/benchmark.length.c new file mode 100644 index 000000000000..06ea07c3d673 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/benchmark/c/benchmark.length.c @@ -0,0 +1,213 @@ +/** +* @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/stats/strided/dnancumaxabs.h" +#include +#include +#include +#include +#include + +#define NAME "dnancumaxabs" +#define ITERATIONS 1000000 +#define REPEATS 3 +#define MIN 1 +#define MAX 6 + +/** +* 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 *x; + double *y; + double t; + int i; + + x = (double *)malloc( len * sizeof(double) ); + y = (double *)malloc( len * sizeof(double) ); + + for ( i = 0; i < len; i++ ) { + if ( rand_double() < 0.2 ) { + x[ i ] = 0.0 / 0.0; // NaN + } else { + x[ i ] = ( rand_double() * 20000.0 ) - 10000.0; + } + y[ i ] = 0.0; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + stdlib_strided_dnancumaxabs( len, x, 1, y, 1 ); + if ( y[ 0 ] != y[ 0 ] && x[ 0 ] == x[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( y[ len-1 ] != y[ len-1 ] && x[ len-1 ] == x[ len-1 ] ) { + printf( "should not return NaN\n" ); + } + free( x ); + free( y ); + 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 *x; + double *y; + double t; + int i; + + x = (double *)malloc( len * sizeof(double) ); + y = (double *)malloc( len * sizeof(double) ); + + for ( i = 0; i < len; i++ ) { + if ( rand_double() < 0.2 ) { + x[ i ] = 0.0 / 0.0; // NaN + } else { + x[ i ] = ( rand_double() * 20000.0 ) - 10000.0; + } + y[ i ] = 0.0; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + stdlib_strided_dnancumaxabs_ndarray( len, x, 1, 0, y, 1, 0 ); + if ( y[ 0 ] != y[ 0 ] && x[ 0 ] == x[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( y[ len-1 ] != y[ len-1 ] && x[ len-1 ] == x[ len-1 ] ) { + printf( "should not return NaN\n" ); + } + free( x ); + free( y ); + 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 = pow( 10, i ); + iter = ITERATIONS / 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 ( i = MIN; i <= MAX; i++ ) { + len = pow( 10, i ); + iter = ITERATIONS / pow( 10, i-1 ); + 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/stats/strided/dnancumaxabs/binding.gyp b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/binding.gyp new file mode 100644 index 000000000000..0d6508a12e99 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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/stats/strided/dnancumaxabs/docs/repl.txt b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt new file mode 100644 index 000000000000..16ee77049280 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt @@ -0,0 +1,113 @@ +{{alias}}( N, x, strideX, y, strideY ) + Computes the cumulative maximum absolute value of double-precision floating- + point strided array elements, ignoring `NaN` values. + + The `N` and stride parameters determine which elements in the strided arrays + are accessed at runtime. + + Indexing is relative to the first index. To introduce an offset, use a typed + array view. + + If `N <= 0`, the function returns `y` unchanged. + + Parameters + ---------- + N: integer + Number of indexed elements. + + x: Float64Array + Input array. + + strideX: integer + Stride length for `x`. + + y: Float64Array + Output array. + + strideY: integer + Stride length for `y`. + + Returns + ------- + out: Float64Array + Output array. + + Examples + -------- + // Standard Usage: + > var x = new {{alias:@stdlib/array/float64}}( [ 1.0, NaN, -2.0, 3.0, NaN, 2.0 ] ); + > var y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}( x.length, x, 1, y, 1 ) + [ 1.0, 1.0, 2.0, 3.0, 3.0, 3.0 ] + + // Using `N` and stride parameters: + > x = new {{alias:@stdlib/array/float64}}( [ -2.0, 1.0, NaN, -5.0, 2.0, -1.0 ] ); + > y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}( 3, x, 2, y, 2 ) + [ 2.0, 0.0, 2.0, 0.0, 2.0, 0.0 ] + + // Using view offsets: + > var x0 = new {{alias:@stdlib/array/float64}}( [ 1.0, NaN, -2.0, 3.0, 2.0, 5.0 ] ); + > var y0 = new {{alias:@stdlib/array/float64}}( x0.length ); + > var x1 = new {{alias:@stdlib/array/float64}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); + > var y1 = new {{alias:@stdlib/array/float64}}( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); + > {{alias}}( 3, x1, 2, y1, 1 ) + [ NaN, 3.0, 5.0 ] + > y0 + [ 0.0, 0.0, 0.0, NaN, 3.0, 5.0 ] + + +{{alias}}.ndarray( N, x, strideX, offsetX, y, strideY, offsetY ) + Computes the cumulative maximum absolute value of double-precision floating- + point strided array elements, ignoring `NaN` values, using alternative + indexing semantics. + + While typed array views mandate a view offset based on the underlying + buffer, the `offset` parameter supports indexing semantics based on a + starting index. + + Parameters + ---------- + N: integer + Number of indexed elements. + + x: Float64Array + Input array. + + strideX: integer + Stride length for `x`. + + offsetX: integer + Starting index for `x`. + + y: Float64Array + Output array. + + strideY: integer + Stride length for `y`. + + offsetY: integer + Starting index for `y`. + + Returns + ------- + out: Float64Array + Output array. + + Examples + -------- + // Standard Usage: + > var x = new {{alias:@stdlib/array/float64}}( [ 1.0, NaN, -2.0, 2.0 ] ); + > var y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}.ndarray( x.length, x, 1, 0, y, 1, 0 ) + [ 1.0, 1.0, 2.0, 2.0 ] + + // Advanced indexing: + > x = new {{alias:@stdlib/array/float64}}( [ 1.0, NaN, -2.0, 3.0, 2.0, 5.0, -1.0 ] ); + > y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}.ndarray( 3, x, 2, 1, y, -1, y.length-1 ) + [ 0.0, 0.0, 0.0, 0.0, 5.0, 3.0, NaN ] + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/index.d.ts b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/index.d.ts new file mode 100644 index 000000000000..9ede1925f87a --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/index.d.ts @@ -0,0 +1,94 @@ +/** +* @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 `dnancumaxabs`. +*/ +interface Routine { + /** + * Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. + * + * @param N - number of indexed elements + * @param x - input array + * @param strideX - `x` stride length + * @param y - output array + * @param strideY - `y` stride length + * @returns output array + * + * @example + * var Float64Array = require( '@stdlib/array/float64' ); + * + * var x = new Float64Array( [ 1.0, -2.0, NaN, 3.0, -4.0 ] ); + * var y = new Float64Array( x.length ); + * + * dnancumaxabs( x.length, x, 1, y, 1 ); + * // y => [ 1.0, 2.0, 2.0, 3.0, 4.0 ] + */ + ( N: number, x: Float64Array, strideX: number, y: Float64Array, strideY: number ): Float64Array; + + /** + * Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values, using alternative indexing semantics. + * + * @param N - number of indexed elements + * @param x - input array + * @param strideX - `x` stride length + * @param offsetX - starting index for `x` + * @param y - output array + * @param strideY - `y` stride length + * @param offsetY - starting index for `y` + * @returns output array + * + * @example + * var Float64Array = require( '@stdlib/array/float64' ); + * + * var x = new Float64Array( [ 1.0, -2.0, NaN, 3.0, -4.0 ] ); + * var y = new Float64Array( x.length ); + * + * dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, 0 ); + * // y => [ 1.0, 2.0, 2.0, 3.0, 4.0 ] + */ + ndarray( N: number, x: Float64Array, strideX: number, offsetX: number, y: Float64Array, strideY: number, offsetY: number ): Float64Array; +} + +/** +* Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. +* +* @param N - number of indexed elements +* @param x - input array +* @param strideX - `x` stride length +* @param y - output array +* @param strideY - `y` stride length +* @returns output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN, 3.0, -4.0 ] ); +* var y = new Float64Array( x.length ); +* +* dnancumaxabs( x.length, x, 1, y, 1 ); +* // y => [ 1.0, 2.0, 2.0, 3.0, 4.0 ] +*/ +declare var dnancumaxabs: Routine; + + +// EXPORTS // + +export = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/test.ts b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/test.ts new file mode 100644 index 000000000000..672d7dbb7057 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/types/test.ts @@ -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. +*/ + +import dnancumaxabs = require( './index' ); + + +// TESTS // + +// The function returns a Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs( x.length, x, 1, y, 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( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs( '10', x, 1, y, 1 ); // $ExpectError + dnancumaxabs( true, x, 1, y, 1 ); // $ExpectError + dnancumaxabs( false, x, 1, y, 1 ); // $ExpectError + dnancumaxabs( null, x, 1, y, 1 ); // $ExpectError + dnancumaxabs( undefined, x, 1, y, 1 ); // $ExpectError + dnancumaxabs( [], x, 1, y, 1 ); // $ExpectError + dnancumaxabs( {}, x, 1, y, 1 ); // $ExpectError + dnancumaxabs( ( x: number ): number => x, x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a second argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs( x.length, 10, 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, '10', 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, true, 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, false, 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, null, 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, undefined, 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, [ '1' ], 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, {}, 1, y, 1 ); // $ExpectError + dnancumaxabs( x.length, ( x: number ): number => x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a third argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs( x.length, x, '10', y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, true, y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, false, y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, null, y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, undefined, y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, [], y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, {}, y, 1 ); // $ExpectError + dnancumaxabs( x.length, x, ( x: number ): number => x, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fourth argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + + dnancumaxabs( x.length, x, 1, 10, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, '10', 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, true, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, false, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, null, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, undefined, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, [ '1' ], 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, {}, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, ( x: number ): number => x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fifth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs( x.length, x, 1, y, '10' ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, true ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, false ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, null ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, undefined ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, [] ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, {} ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs(); // $ExpectError + dnancumaxabs( x.length ); // $ExpectError + dnancumaxabs( x.length, x ); // $ExpectError + dnancumaxabs( x.length, x, 1 ); // $ExpectError + dnancumaxabs( x.length, x, 1, y ); // $ExpectError + dnancumaxabs( x.length, x, 1, y, 1, 10 ); // $ExpectError +} + +// Attached to main export is an `ndarray` method which returns a Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 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( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( '10', x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( true, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( false, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( null, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( undefined, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( [], x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( {}, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( ( x: number ): number => x, x, 1, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a second argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, 10, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, '10', 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, true, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, false, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, null, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, undefined, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, [ '1' ], 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, {}, 1, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, ( x: number ): number => x, 1, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a third argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, x, '10', 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, true, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, false, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, null, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, undefined, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, [], 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, {}, 0, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, ( x: number ): number => x, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fourth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, x, 1, '10', y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, true, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, false, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, null, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, undefined, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, [], y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, {}, y, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, ( x: number ): number => x, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, x, 1, 0, 10, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, '10', 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, true, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, false, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, null, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, undefined, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, [ '1' ], 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, {}, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, ( x: number ): number => x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, x, 1, 0, y, '10', 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, true, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, false, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, null, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, undefined, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, [], 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, {}, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, '10' ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, true ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, false ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, null ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, undefined ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, [] ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, {} ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancumaxabs.ndarray(); // $ExpectError + dnancumaxabs.ndarray( x.length ); // $ExpectError + dnancumaxabs.ndarray( x.length, x ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1 ); // $ExpectError + dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, 0, 10 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/c/Makefile b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/c/Makefile new file mode 100644 index 000000000000..c8f8e9a1517b --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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/stats/strided/dnancumaxabs/examples/c/example.c b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/c/example.c new file mode 100644 index 000000000000..22efaf08a2a1 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/c/example.c @@ -0,0 +1,41 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/stats/strided/dnancumaxabs.h" +#include + +int main( void ) { + // Create strided arrays: + const double x[] = { 1.0, -2.0, -3.0, 4.0, 0.0/0.0, -6.0, 7.0, 8.0, 0.0/0.0, 10.0 }; + double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + + // Specify the number of elements: + const int N = 5; + + // Specify stride lengths: + const int strideX = 2; + const int strideY = -2; + + // Compute the cumulative maximum absolute value: + stdlib_strided_dnancumaxabs( N, x, strideX, y, strideY ); + + // Print the result: + for ( int i = 0; i < 10; i++ ) { + printf( "y[ %d ] = %lf\n", i, y[ i ] ); + } +} diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/index.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/index.js new file mode 100644 index 000000000000..53226cc2ce88 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/examples/index.js @@ -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. +*/ + +'use strict'; + +var Float64Array = require( '@stdlib/array/float64' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var dnancumaxabs = require( './../lib' ); + +/** +* Returns a random number. +* +* @private +* @returns {number} random number or `NaN` +*/ +function rand() { + if ( bernoulli( 0.8 ) < 1 ) { + return NaN; + } + return uniform( -50.0, 50.0 ); +} + +var x = filledarrayBy( 10, 'float64', rand ); +console.log( x ); + +var y = new Float64Array( x.length ); +dnancumaxabs( x.length, x, 1, y, 1 ); +console.log( y ); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/include.gypi b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/include.gypi new file mode 100644 index 000000000000..bee8d41a2caf --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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, 2.0, 3.0, 3.0, 3.0 ] +*/ +function dnancumaxabs( N, x, strideX, y, strideY ) { + var ox = stride2offset( N, strideX ); + var oy = stride2offset( N, strideY ); + return ndarray( N, x, strideX, ox, y, strideY, oy ); +} + + +// EXPORTS // + +module.exports = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/dnancumaxabs.native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/dnancumaxabs.native.js new file mode 100644 index 000000000000..3cdc3f6e5f0d --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/dnancumaxabs.native.js @@ -0,0 +1,55 @@ +/** +* @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 cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {Float64Array} x - input array +* @param {integer} strideX - `x` stride length +* @param {Float64Array} y - output array +* @param {integer} strideY - `y` stride length +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN, 3.0, -4.0 ] ); +* var y = new Float64Array( x.length ); +* +* dnancumaxabs( x.length, x, 1, y, 1 ); +* // y => [ 1.0, 2.0, 2.0, 3.0, 4.0 ] +*/ +function dnancumaxabs( N, x, strideX, y, strideY ) { + addon( N, x, strideX, y, strideY ); + return y; +} + + +// EXPORTS // + +module.exports = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/index.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/index.js new file mode 100644 index 000000000000..e084471f3885 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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 cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. +* +* @module @stdlib/stats/strided/dnancumaxabs +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var dnancumaxabs = require( '@stdlib/stats/strided/dnancumaxabs' ); +* +* var x = new Float64Array( [ 1.0, NaN, -2.0, 3.0, NaN, 2.0 ] ); +* var y = new Float64Array( x.length ); +* +* dnancumaxabs( x.length, x, 1, y, 1 ); +* // y => [ 1.0, 1.0, 2.0, 3.0, 3.0, 3.0 ] +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var dnancumaxabs = require( '@stdlib/stats/strided/dnancumaxabs' ); +* +* var x = new Float64Array( [ 1.0, NaN, -2.0, 3.0, NaN, 2.0 ] ); +* var y = new Float64Array( x.length ); +* +* dnancumaxabs.ndarray( 6, x, 1, 0, y, 1, 0 ); +* // y => [ 1.0, 1.0, 2.0, 3.0, 3.0, 3.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 dnancumaxabs; +var tmp = tryRequire( join( __dirname, './native.js' ) ); +if ( isError( tmp ) ) { + dnancumaxabs = main; +} else { + dnancumaxabs = tmp; +} + + +// EXPORTS // + +module.exports = dnancumaxabs; + +// exports: { "ndarray": "dnancumaxabs.ndarray" } diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/main.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/main.js new file mode 100644 index 000000000000..4cbe7d905be2 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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 dnancumaxabs = require( './dnancumaxabs.js' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +setReadOnly( dnancumaxabs, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/native.js new file mode 100644 index 000000000000..d75b17319089 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/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 dnancumaxabs = require( './dnancumaxabs.native.js' ); +var ndarray = require( './ndarray.native.js' ); + + +// MAIN // + +setReadOnly( dnancumaxabs, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.js new file mode 100644 index 000000000000..04cc07a02113 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.js @@ -0,0 +1,100 @@ +/** +* @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 abs = require( '@stdlib/math/base/special/abs' ); + + +// MAIN // + +/** +* Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {Float64Array} x - input array +* @param {integer} strideX - `x` stride length +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @param {Float64Array} y - output array +* @param {integer} strideY - `y` stride length +* @param {NonNegativeInteger} offsetY - starting index for `y` +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, NaN, -2.0, 3.0, NaN, 2.0 ] ); +* var y = new Float64Array( x.length ); +* +* dnancumaxabs( 6, x, 1, 0, y, 1, 0 ); +* // y => [ 1.0, 1.0, 2.0, 3.0, 3.0, 3.0 ] +*/ +function dnancumaxabs( N, x, strideX, offsetX, y, strideY, offsetY ) { + var max; + var ix; + var iy; + var v; + var i; + + if ( N <= 0 ) { + return y; + } + ix = offsetX; + iy = offsetY; + + // Find the first non-NaN element: + for ( i = 0; i < N; i++ ) { + v = x[ ix ]; + if ( v === v ) { + break; + } + y[ iy ] = NaN; + ix += strideX; + iy += strideY; + } + if ( i === N ) { + return y; + } + max = abs( v ); + y[ iy ] = max; + + i += 1; + ix += strideX; + iy += strideY; + + for ( i; i < N; i++ ) { + v = x[ ix ]; + if ( v === v ) { + v = abs( v ); + if ( v > max ) { + max = v; + } + } + y[ iy ] = max; + ix += strideX; + iy += strideY; + } + return y; +} + + +// EXPORTS // + +module.exports = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.native.js new file mode 100644 index 000000000000..c0d5608681e1 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/lib/ndarray.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 cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring `NaN` values, using alternative indexing semantics. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {Float64Array} x - input array +* @param {integer} strideX - `x` stride length +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @param {Float64Array} y - output array +* @param {integer} strideY - `y` stride length +* @param {NonNegativeInteger} offsetY - starting index for `y` +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN, 3.0, -4.0 ] ); +* var y = new Float64Array( x.length ); +* +* dnancumaxabs( 5, x, 1, 0, y, 1, 0 ); +* // y => [ 1.0, 2.0, 2.0, 3.0, 4.0 ] +*/ +function dnancumaxabs( N, x, strideX, offsetX, y, strideY, offsetY ) { + addon.ndarray( N, x, strideX, offsetX, y, strideY, offsetY ); + return y; +} + + +// EXPORTS // + +module.exports = dnancumaxabs; diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/manifest.json b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/manifest.json new file mode 100644 index 000000000000..3a6cdf8a726d --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/manifest.json @@ -0,0 +1,137 @@ +{ + "name": "dnancumaxabs", + "alias": "dnancumaxabs", + "description": "Calculate the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring NaN values.", + "keywords": [ + "statistics", + "stats", + "mathematics", + "math", + "maximum", + "max", + "abs", + "absolute", + "absolute value", + "accumulate", + "cumulative", + "strided", + "strided array", + "float64", + "double", + "float64array", + "nan", + "ignorenan" + ], + "dependencies": {}, + "extensions": { + "dag": { + "type": "function", + "arity": 5, + "returns": "Float64Array" + } + }, + "options": { + "task": "build", + "wasm": false + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true + }, + { + "field": "include", + "resolve": true, + "relative": true + }, + { + "field": "libraries", + "resolve": false, + "relative": false + }, + { + "field": "libpath", + "resolve": true, + "relative": false + } + ], + "confs": [ + { + "task": "build", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/assert/is-nan", + "@stdlib/blas/base/shared", + "@stdlib/math/base/special/abs", + "@stdlib/strided/base/stride2offset", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-float64array" + ] + }, + { + "task": "benchmark", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/abs" + ] + }, + { + "task": "examples", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/abs" + ] + }, + { + "task": "", + "wasm": true, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/special/abs" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/package.json b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/package.json new file mode 100644 index 000000000000..dd185e48644a --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/package.json @@ -0,0 +1,84 @@ +{ + "name": "@stdlib/stats/strided/dnancumaxabs", + "version": "0.0.0", + "description": "Calculate the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring NaN values.", + "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", + "statistics", + "stats", + "mathematics", + "math", + "maximum", + "max", + "abs", + "absolute", + "absolute value", + "range", + "extremes", + "domain", + "extent", + "accumulate", + "cumulative", + "strided", + "strided array", + "typed", + "array", + "float64", + "double", + "float64array", + "nan", + "ignorenan" + ], + "__stdlib__": {} +} diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/Makefile b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/Makefile new file mode 100644 index 000000000000..2caf905cedbe --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2026 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + + +# RULES # + +#/ +# Removes generated files for building an add-on. +# +# @example +# make clean-addon +#/ +clean-addon: + $(QUIET) -rm -f *.o *.node + +.PHONY: clean-addon + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: clean-addon + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/addon.c b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/addon.c new file mode 100644 index 000000000000..13da33def999 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/addon.c @@ -0,0 +1,65 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/stats/strided/dnancumaxabs.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/blas/base/shared.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, 5 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 4 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 1 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, Y, N, strideY, argv, 3 ); + API_SUFFIX(stdlib_strided_dnancumaxabs)( N, X, strideX, Y, strideY ); + 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, 7 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 5 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, offsetY, argv, 6 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 1 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, Y, N, strideY, argv, 4 ); + API_SUFFIX(stdlib_strided_dnancumaxabs_ndarray)( N, X, strideX, offsetX, Y, strideY, offsetY ); + return NULL; +} + +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ) diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/main.c b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/main.c new file mode 100644 index 000000000000..7c5bba5688fb --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/src/main.c @@ -0,0 +1,97 @@ +/** +* @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/stats/strided/dnancumaxabs.h" +#include "stdlib/math/base/assert/is_nan.h" +#include "stdlib/math/base/special/abs.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" + +/** +* Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring NaN values. +* +* @param N number of indexed elements +* @param X input array +* @param strideX X stride length +* @param Y output array +* @param strideY Y stride length +*/ +void API_SUFFIX(stdlib_strided_dnancumaxabs)( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, double *Y, const CBLAS_INT strideY ) { + const CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); + const CBLAS_INT oy = stdlib_strided_stride2offset( N, strideY ); + API_SUFFIX(stdlib_strided_dnancumaxabs_ndarray)( N, X, strideX, ox, Y, strideY, oy ); + return; +} + +/** +* Computes the cumulative maximum absolute value of double-precision floating-point strided array elements, ignoring NaN values, using alternative indexing semantics. +* +* @param N number of indexed elements +* @param X input array +* @param strideX X stride length +* @param offsetX starting index for X +* @param Y output array +* @param strideY Y stride length +* @param offsetY starting index for Y +*/ +void API_SUFFIX(stdlib_strided_dnancumaxabs_ndarray)( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ) { + CBLAS_INT ix; + CBLAS_INT iy; + CBLAS_INT i; + double max; + double v; + + if ( N <= 0 ) { + return; + } + ix = offsetX; + iy = offsetY; + + for ( i = 0; i < N; i++ ) { + v = X[ ix ]; + if ( !stdlib_base_is_nan( v ) ) { + break; + } + Y[ iy ] = v; // NaN + ix += strideX; + iy += strideY; + } + if ( i == N ) { + return; + } + max = stdlib_base_abs( v ); + Y[ iy ] = max; + + i += 1; + ix += strideX; + iy += strideY; + + for (; i < N; i++ ) { + v = X[ ix ]; + if ( !stdlib_base_is_nan( v ) ) { + v = stdlib_base_abs( v ); + if ( v > max ) { + max = v; + } + } + Y[ iy ] = max; + ix += strideX; + iy += strideY; + } + return; +} diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.js new file mode 100644 index 000000000000..3a75bcbdeada --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.js @@ -0,0 +1,273 @@ +/** +* @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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var isPositiveZero = require( '@stdlib/math/base/assert/is-positive-zero' ); +var Float64Array = require( '@stdlib/array/float64' ); +var dnancumaxabs = require( './../lib/dnancumaxabs.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancumaxabs, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 5', function test( t ) { + t.strictEqual( dnancumaxabs.length, 5, 'has expected arity' ); + t.end(); +}); + +tape( 'the function computes the cumulative maximum absolute value, ignoring NaNs', function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, -2.0, 3.0, -4.0, 5.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 2.0, + 3.0, + 4.0, + 5.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ -0.0, 0.0, -0.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + 0.0, + 0.0, + 0.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isPositiveZero( expected[ i ] ) ) { + t.strictEqual( isPositiveZero( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], true, 'returns expected value. i: ' + i ); + } + } + + x = new Float64Array( [ NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 3.0, + 3.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + NaN, + 1.0, + 1.0, + 3.0, + 3.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isnan( expected[ i ] ) ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], 'returns expected value. i: ' + i ); + } + } + t.end(); +}); + +tape( 'the function returns a reference to the output array', function test( t ) { + var out; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + + out = dnancumaxabs( x.length, x, 1, y, 1 ); + t.strictEqual( out, y, 'same reference' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns `y` unchanged', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + expected = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dnancumaxabs( -1, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancumaxabs( 0, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + t.end(); +}); + +tape( 'the function supports an `x` stride', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, y, 1 ); + + expected = new Float64Array( [ 1.0, 3.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` stride', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, // 1 + 3.0, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancumaxabs( 3, x, 1, y, 2 ); + + expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports negative strides', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 2 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, -2, y, -1 ); + + expected = new Float64Array( [ 5.0, 5.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports complex access patterns', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + 2.0, + -3.0, // 1 + 4.0, + 5.0, // 2 + 6.0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, y, -1 ); + + expected = new Float64Array( [ 5.0, 3.0, 1.0, 0.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.native.js new file mode 100644 index 000000000000..3c9ec4d8435e --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.dnancumaxabs.native.js @@ -0,0 +1,282 @@ +/** +* @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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var isPositiveZero = require( '@stdlib/math/base/assert/is-positive-zero' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var dnancumaxabs = tryRequire( resolve( __dirname, './../lib/dnancumaxabs.native.js' ) ); +var opts = { + 'skip': ( dnancumaxabs instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancumaxabs, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 5', opts, function test( t ) { + t.strictEqual( dnancumaxabs.length, 5, 'has expected arity' ); + t.end(); +}); + +tape( 'the function computes the cumulative maximum absolute value, ignoring NaNs', opts, function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, -2.0, 3.0, -4.0, 5.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 2.0, + 3.0, + 4.0, + 5.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ -0.0, 0.0, -0.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + 0.0, + 0.0, + 0.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isPositiveZero( expected[ i ] ) ) { + t.strictEqual( isPositiveZero( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], true, 'returns expected value. i: ' + i ); + } + } + + x = new Float64Array( [ NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 3.0, + 3.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, y, 1 ); + + expected = new Float64Array([ + NaN, + 1.0, + 1.0, + 3.0, + 3.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isnan( expected[ i ] ) ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], 'returns expected value. i: ' + i ); + } + } + t.end(); +}); + +tape( 'the function returns a reference to the output array', opts, function test( t ) { + var out; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + + out = dnancumaxabs( x.length, x, 1, y, 1 ); + t.strictEqual( out, y, 'same reference' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns `y` unchanged', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + expected = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dnancumaxabs( -1, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancumaxabs( 0, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + t.end(); +}); + +tape( 'the function supports an `x` stride', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, y, 1 ); + + expected = new Float64Array( [ 1.0, 3.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` stride', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, // 1 + 3.0, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancumaxabs( 3, x, 1, y, 2 ); + + expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports negative strides', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 2 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, -2, y, -1 ); + + expected = new Float64Array( [ 5.0, 5.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports complex access patterns', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + 2.0, + -3.0, // 1 + 4.0, + 5.0, // 2 + 6.0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, y, -1 ); + + expected = new Float64Array( [ 5.0, 3.0, 1.0, 0.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.js new file mode 100644 index 000000000000..77c83e1868a9 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.js @@ -0,0 +1,38 @@ +/** +* @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 dnancumaxabs = require( './../lib' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancumaxabs, '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 dnancumaxabs.ndarray, 'function', 'method is a function' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.js new file mode 100644 index 000000000000..0367aaed7ca8 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.js @@ -0,0 +1,341 @@ +/** +* @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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var isPositiveZero = require( '@stdlib/math/base/assert/is-positive-zero' ); +var Float64Array = require( '@stdlib/array/float64' ); +var dnancumaxabs = require( './../lib/ndarray.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancumaxabs, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', function test( t ) { + t.strictEqual( dnancumaxabs.length, 7, 'has expected arity' ); + t.end(); +}); + +tape( 'the function calculates the cumulative maximum absolute value, ignoring NaNs', function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, -2.0, 3.0, -4.0, 5.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 2.0, + 3.0, + 4.0, + 5.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ -0.0, 0.0, -0.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 0.0, + 0.0, + 0.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isPositiveZero( expected[ i ] ) ) { + t.strictEqual( isPositiveZero( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], true, 'returns expected value. i: ' + i ); + } + } + + x = new Float64Array( [ NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 3.0, + 3.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + NaN, + 1.0, + 1.0, + 3.0, + 3.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isnan( expected[ i ] ) ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], 'returns expected value. i: ' + i ); + } + } + t.end(); +}); + +tape( 'the function returns a reference to the output array', function test( t ) { + var out; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + + out = dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + t.strictEqual( out, y, 'same reference' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns `y` unchanged', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + expected = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dnancumaxabs( -1, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancumaxabs( 0, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + t.end(); +}); + +tape( 'the function supports an `x` stride', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, 0, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, 3.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` stride', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, // 1 + 3.0, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancumaxabs( 3, x, 1, 0, y, 2, 0 ); + + expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports negative strides', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 2 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, -2, x.length-1, y, -1, 2 ); + + expected = new Float64Array( [ 5.0, 5.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports an `x` offset', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 2.0, + 1.0, // 0 + 2.0, + -2.0, // 1 + -2.0, + 2.0, // 2 + 3.0, + 4.0 // 3 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, // 3 + 0.0, + 0.0, + 0.0, + 0.0 + ]); + + dnancumaxabs( 4, x, 2, 1, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, 2.0, 2.0, 4.0, 0.0, 0.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` offset', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 2.0, // 0 + 1.0, // 1 + 2.0, // 2 + -2.0, // 3 + -2.0, + 2.0, + 3.0, + 4.0 + ]); + y = new Float64Array([ + 0.0, + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0, // 2 + 0.0, + 0.0 // 3 + ]); + + dnancumaxabs( 4, x, 1, 0, y, 2, 1 ); + + expected = new Float64Array( [ 0.0, 2.0, 0.0, 2.0, 0.0, 2.0, 0.0, 2.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports complex access patterns', function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + 2.0, + -3.0, // 1 + 4.0, + 5.0, // 2 + 6.0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, 0, y, -1, 2 ); + + expected = new Float64Array( [ 5.0, 3.0, 1.0, 0.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.native.js b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.native.js new file mode 100644 index 000000000000..e203731845b7 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/test/test.ndarray.native.js @@ -0,0 +1,350 @@ +/** +* @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 isnan = require( '@stdlib/math/base/assert/is-nan' ); +var isPositiveZero = require( '@stdlib/math/base/assert/is-positive-zero' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var dnancumaxabs = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( dnancumaxabs instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancumaxabs, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 7', opts, function test( t ) { + t.strictEqual( dnancumaxabs.length, 7, 'has expected arity' ); + t.end(); +}); + +tape( 'the function calculates the cumulative maximum absolute value, ignoring NaNs', opts, function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, -2.0, 3.0, -4.0, 5.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 2.0, + 3.0, + 4.0, + 5.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ -0.0, 0.0, -0.0 ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 0.0, + 0.0, + 0.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isPositiveZero( expected[ i ] ) ) { + t.strictEqual( isPositiveZero( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], true, 'returns expected value. i: ' + i ); + } + } + + x = new Float64Array( [ NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + for ( i = 0; i < y.length; i++ ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } + + x = new Float64Array( [ 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 3.0, + 3.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, 1.0, NaN, -3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + NaN, + 1.0, + 1.0, + 3.0, + 3.0 + ]); + for ( i = 0; i < y.length; i++ ) { + if ( isnan( expected[ i ] ) ) { + t.strictEqual( isnan( y[ i ] ), true, 'returns expected value. i: ' + i ); + } else { + t.strictEqual( y[ i ], expected[ i ], 'returns expected value. i: ' + i ); + } + } + t.end(); +}); + +tape( 'the function returns a reference to the output array', opts, function test( t ) { + var out; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + + out = dnancumaxabs( x.length, x, 1, 0, y, 1, 0 ); + t.strictEqual( out, y, 'same reference' ); + + t.end(); +}); + +tape( 'if provided an `N` parameter less than or equal to `0`, the function returns `y` unchanged', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + y = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + expected = new Float64Array( [ 6.0, 7.0, 8.0, 9.0, 10.0 ] ); + + dnancumaxabs( -1, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancumaxabs( 0, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + t.end(); +}); + +tape( 'the function supports an `x` stride', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, 0, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, 3.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` stride', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + -2.0, // 1 + 3.0, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancumaxabs( 3, x, 1, 0, y, 2, 0 ); + + expected = new Float64Array( [ 1.0, 0.0, 2.0, 0.0, 3.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports negative strides', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 2 + -2.0, + 3.0, // 1 + 4.0, + -5.0 // 0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, -2, x.length-1, y, -1, 2 ); + + expected = new Float64Array( [ 5.0, 5.0, 5.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports an `x` offset', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 2.0, + 1.0, // 0 + 2.0, + -2.0, // 1 + -2.0, + 2.0, // 2 + 3.0, + 4.0 // 3 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, // 3 + 0.0, + 0.0, + 0.0, + 0.0 + ]); + + dnancumaxabs( 4, x, 2, 1, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, 2.0, 2.0, 4.0, 0.0, 0.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports a `y` offset', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 2.0, // 0 + 1.0, // 1 + 2.0, // 2 + -2.0, // 3 + -2.0, + 2.0, + 3.0, + 4.0 + ]); + y = new Float64Array([ + 0.0, + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0, // 2 + 0.0, + 0.0 // 3 + ]); + + dnancumaxabs( 4, x, 1, 0, y, 2, 1 ); + + expected = new Float64Array( [ 0.0, 2.0, 0.0, 2.0, 0.0, 2.0, 0.0, 2.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports complex access patterns', opts, function test( t ) { + var expected; + var x; + var y; + + x = new Float64Array([ + 1.0, // 0 + 2.0, + -3.0, // 1 + 4.0, + 5.0, // 2 + 6.0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancumaxabs( 3, x, 2, 0, y, -1, 2 ); + + expected = new Float64Array( [ 5.0, 3.0, 1.0, 0.0, 0.0, 0.0 ] ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); From 7c88f8048316f0832b57200e82acbf51e8e94652 Mon Sep 17 00:00:00 2001 From: Om-A-osc Date: Wed, 4 Mar 2026 23:48:09 +0530 Subject: [PATCH 2/2] fix: fixing markdown files --- 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: 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: na - 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 --- --- .../stats/strided/dnancumaxabs/README.md | 30 ++++--------------- .../stats/strided/dnancumaxabs/docs/repl.txt | 1 + 2 files changed, 6 insertions(+), 25 deletions(-) diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md index fd2598315e31..d25763034d64 100644 --- a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/README.md @@ -162,9 +162,6 @@ console.log( x ); var y = new Float64Array( x.length ); dnancumaxabs( x.length, x, 1, y, 1 ); console.log( y ); - -dnancumaxabs.ndarray( x.length, x, 1, 0, y, 1, 0 ); -console.log( y ); ``` @@ -269,21 +266,21 @@ void stdlib_strided_dnancumaxabs_ndarray( const CBLAS_INT N, const double *X, co int main( void ) { // Create strided arrays: - const double x[] = { 1.0, 0.0, 2.0, 0.0, 0.0/0.0, 0.0, -3.0, 0.0, 4.0, 0.0, -5.0, 0.0, 6.0, 0.0, 7.0, 0.0 }; - double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + const double x[] = { 1.0, -2.0, -3.0, 4.0, 0.0/0.0, -6.0, 7.0, 8.0, 0.0/0.0, 10.0 }; + double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; // Specify the number of elements: - const int N = 8; + const int N = 5; // Specify stride lengths: const int strideX = 2; - const int strideY = 2; + const int strideY = -2; // Compute the cumulative maximum absolute value: stdlib_strided_dnancumaxabs( N, x, strideX, y, strideY ); // Print the result: - for ( int i = 0; i < 16; i++ ) { + for ( int i = 0; i < 10; i++ ) { printf( "y[ %d ] = %lf\n", i, y[ i ] ); } } @@ -307,15 +304,6 @@ int main( void ) { @@ -330,14 +318,6 @@ int main( void ) { -[@stdlib/stats/base/dnancumaxabs]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/stats/base/dnancumaxabs - -[@stdlib/stats/strided/dcumaxabs]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/stats/strided/dcumaxabs - -[@stdlib/stats/strided/dnancumax]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/stats/strided/dnancumax - -[@stdlib/stats/strided/dnancuminabs]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/stats/strided/dnancuminabs - diff --git a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt index 16ee77049280..9373bfd72fa0 100644 --- a/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt +++ b/lib/node_modules/@stdlib/stats/strided/dnancumaxabs/docs/repl.txt @@ -1,3 +1,4 @@ + {{alias}}( N, x, strideX, y, strideY ) Computes the cumulative maximum absolute value of double-precision floating- point strided array elements, ignoring `NaN` values.