From 0806fda5a25ee0044b19d47c4539c1ec4f7adfb9 Mon Sep 17 00:00:00 2001 From: vivek maurya Date: Thu, 6 Mar 2025 14:07:45 +0530 Subject: [PATCH 1/4] feat: add blas/ext/base/dnancusumkbn --- .../blas/ext/base/dnancusumkbn/README.md | 361 ++++++++++++++++++ .../base/dnancusumkbn/benchmark/benchmark.js | 117 ++++++ .../benchmark/benchmark.native.js | 122 ++++++ .../benchmark/benchmark.ndarray.js | 117 ++++++ .../benchmark/benchmark.ndarray.native.js | 122 ++++++ .../base/dnancusumkbn/benchmark/c/Makefile | 146 +++++++ .../benchmark/c/benchmark.length.c | 201 ++++++++++ .../blas/ext/base/dnancusumkbn/binding.gyp | 170 +++++++++ .../blas/ext/base/dnancusumkbn/docs/repl.txt | 120 ++++++ .../base/dnancusumkbn/docs/types/index.d.ts | 106 +++++ .../ext/base/dnancusumkbn/docs/types/test.ts | 280 ++++++++++++++ .../ext/base/dnancusumkbn/examples/c/Makefile | 146 +++++++ .../base/dnancusumkbn/examples/c/example.c | 41 ++ .../ext/base/dnancusumkbn/examples/index.js | 44 +++ .../blas/ext/base/dnancusumkbn/include.gypi | 53 +++ .../stdlib/blas/ext/base/dnancusumkbn.h | 46 +++ .../ext/base/dnancusumkbn/lib/dnancusumkbn.js | 65 ++++ .../dnancusumkbn/lib/dnancusumkbn.native.js | 56 +++ .../blas/ext/base/dnancusumkbn/lib/index.js | 70 ++++ .../blas/ext/base/dnancusumkbn/lib/main.js | 35 ++ .../blas/ext/base/dnancusumkbn/lib/native.js | 35 ++ .../blas/ext/base/dnancusumkbn/lib/ndarray.js | 92 +++++ .../base/dnancusumkbn/lib/ndarray.native.js | 58 +++ .../blas/ext/base/dnancusumkbn/manifest.json | 85 +++++ .../blas/ext/base/dnancusumkbn/package.json | 83 ++++ .../blas/ext/base/dnancusumkbn/src/Makefile | 70 ++++ .../blas/ext/base/dnancusumkbn/src/addon.c | 68 ++++ .../blas/ext/base/dnancusumkbn/src/main.c | 105 +++++ .../dnancusumkbn/test/test.dnancusumkbn.js | 304 +++++++++++++++ .../test/test.dnancusumkbn.native.js | 313 +++++++++++++++ .../blas/ext/base/dnancusumkbn/test/test.js | 82 ++++ .../base/dnancusumkbn/test/test.ndarray.js | 336 ++++++++++++++++ .../dnancusumkbn/test/test.ndarray.native.js | 345 +++++++++++++++++ 33 files changed, 4394 insertions(+) create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/manifest.json create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js create mode 100644 lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md new file mode 100644 index 000000000000..9c51ab628e47 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md @@ -0,0 +1,361 @@ + + +# dnancusumkbn + +> Calculate the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. + +
+ +
+ + + +
+ +## Usage + +```javascript +var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); +``` + +#### dnancusumkbn( N, sum, x, strideX, y, strideY ) + +Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, -2.0, NaN ] ); +var y = new Float64Array( x.length ); + +dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); +// y => [ 1.0, -1.0, -1.0 ] + +x = new Float64Array( [ 1.0, -2.0, NaN ] ); +y = new Float64Array( x.length ); + +dnancusumkbn( x.length, 10.0, x, 1, y, 1 ); +// y => [ 11.0, 9.0, 9.0 ] +``` + +The function has the following parameters: + +- **N**: number of indexed elements. +- **sum**: initial sum. +- **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 sum of every other element: + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, 2.0, 2.0, -7.0, -2.0, 3.0, NaN, NaN ] ); +var y = new Float64Array( x.length ); + +var v = dnancusumkbn( 4, 0.0, x, 2, y, 1 ); +// y => [ 1.0, 3.0, 1.0, 1.0, 0.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( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, NaN, NaN ] ); +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 + +dnancusumkbn( 4, 0.0, x1, -2, y1, 1 ); +// y0 => [ 0.0, 0.0, 0.0, 0.0, 2.0, 0.0, 1.0, 0.0 ] +``` + +#### dnancusumkbn.ndarray( N, sum, x, strideX, offsetX, y, strideY, offsetY ) + +Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 1.0, -2.0, NaN ] ); +var y = new Float64Array( x.length ); + +dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0 ); +// y => [ 1.0, -1.0, -1.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 sum of every other value in the strided input array starting from the second value and to store in the last `N` elements of the strided output array starting from the last element: + +```javascript +var Float64Array = require( '@stdlib/array/float64' ); + +var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, NaN, NaN ] ); +var y = new Float64Array( x.length ); + +dnancusumkbn.ndarray( 4, 0.0, x, 2, 1, y, -1, y.length-1 ); +// y => [ 0.0, 0.0, 0.0, 0.0, 1.0, 1.0, -1.0, 1.0 ] +``` + +
+ + + +
+ +## Notes + +- If `N <= 0`, both functions return `y` unchanged. + +
+ + + +
+ +## Examples + + + +```javascript +var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); +var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); + +function rand() { + if ( bernoulli( 0.7 ) > 0 ) { + return discreteUniform( 0, 100 ); + } + return NaN; +} + +function constant() { + return 0.0; +} + +var x = filledarrayBy( 10, 'float64', rand ); +console.log( x ); + +var y = filledarrayBy( 10, 'float64', constant ); +console.log( y ); + +dnancusumkbn( x.length, 0.0, x, 1, y, -1 ); +console.log( y ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/blas/ext/base/dnancusumkbn.h" +``` + +#### stdlib_strided_dnancusumkbn( N, sum, \*X, strideX, \*Y, strideY ) + +Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. + +```c +const double x[] = { 1.0, 2.0, 3.0, 4.0 } +double y[] = { 0.0, 0.0, 0.0, 0.0 } + +stdlib_strided_dnancusumkbn( 4, 0.0, x, 1, y, 1 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **sum**: `[in] double` initial sum. +- **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_dnancusumkbn( const CBLAS_INT N, const double sum, const double *X, const CBLAS_INT strideX, double *Y, const CBLAS_INT strideY ); +``` + + + +#### stdlib_strided_dnancusumkbn_ndarray( N, sum, \*X, strideX, offsetX, \*Y, strideY, offsetY ) + + + +Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. + +```c +const double x[] = { 1.0, 2.0, 3.0, 4.0 } +double y[] = { 0.0, 0.0, 0.0, 0.0 } + +stdlib_strided_dnancusumkbn_ndarray( 4, 0.0, x, 1, 0, y, 1, 0 ); +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **sum**: `[in] double` initial sum. +- **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_dnancusumkbn_ndarray( const CBLAS_INT N, const double sum, 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/blas/ext/base/dnancusumkbn.h" +#include + +int main( void ) { + // Create strided arrays: + const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 0.0/0.0, 0.0/0.0 }; + double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + + // Specify the number of elements: + const int N = 4; + + // Specify stride lengths: + const int strideX = 2; + const int strideY = -2; + + // Compute the cumulative sum: + stdlib_strided_dnancusumkbn( N, 0.0, x, strideX, y, strideY ); + + // Print the result: + for ( int i = 0; i < 8; i++ ) { + printf( "y[ %d ] = %lf\n", i, y[ i ] ); + } +} +``` + +
+ + + +
+ + + +
+ +## References + +- Neumaier, Arnold. 1974. "Rounding Error Analysis of Some Methods for Summing Finite Sums." _Zeitschrift Für Angewandte Mathematik Und Mechanik_ 54 (1): 39–51. doi:[10.1002/zamm.19740540106][@neumaier:1974a]. + +
+ + + + + +
+ +* * * + +## See Also + +- [`@stdlib/blas/ext/base/dcusum`][@stdlib/blas/ext/base/dcusum]: calculate the cumulative sum of double-precision floating-point strided array elements. +- [`@stdlib/blas/ext/base/gcusumkbn`][@stdlib/blas/ext/base/gcusumkbn]: calculate the cumulative sum of strided array elements using an improved Kahan–Babuška algorithm. +- [`@stdlib/blas/ext/base/scusumkbn`][@stdlib/blas/ext/base/scusumkbn]: calculate the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm. + +
+ + + + + +
+ +[@stdlib/array/float64]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/float64 + +[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray + +[@neumaier:1974a]: https://doi.org/10.1002/zamm.19740540106 + + + +[@stdlib/blas/ext/base/dcusum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dcusum + +[@stdlib/blas/ext/base/gcusumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/gcusumkbn + +[@stdlib/blas/ext/base/scusumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/scusumkbn + + + +
+ + diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js new file mode 100644 index 000000000000..3c0ab69e2d54 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js @@ -0,0 +1,117 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var pkg = require( './../package.json' ).name; +var dnancusumkbn = require( './../lib/dnancusumkbn.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number +*/ +function rand() { + if ( bernoulli( 0.8 ) > 0 ) { + return uniform( -100.0, 100.0 ); + } + return NaN; +} + +function constant() { + return 0.0; +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, options.dtype, rand ); + var y = filledarrayBy( len, options.dtype, constant ); + return benchmark; + + function benchmark( b ) { + var v; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + x[ 0 ] += 1.0; + v = dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+':len='+len, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js new file mode 100644 index 000000000000..d37500e9e1d5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js @@ -0,0 +1,122 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var options = { + 'dtype': 'float64' +}; +var dnancusumkbn = tryRequire( resolve( __dirname, './../lib/dnancusumkbn.native.js' ) ); +var opts = { + 'skip': ( dnancusumkbn instanceof Error ) +}; + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number +*/ +function rand() { + if ( bernoulli( 0.8 ) > 0 ) { + return uniform( -100.0, 100.0 ); + } + return NaN; +} + +function constant() { + return 0.0; +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, options.dtype, rand ); + var y = filledarrayBy( len, options.dtype, constant ); + return benchmark; + + function benchmark( b ) { + var v; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + x[ 0 ] += 1.0; + v = dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+'::native:len='+len, opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js new file mode 100644 index 000000000000..f0c4d98c2f9d --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js @@ -0,0 +1,117 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var pkg = require( './../package.json' ).name; +var dnancusumkbn = require( './../lib/ndarray.js' ); + + +// VARIABLES // + +var options = { + 'dtype': 'float64' +}; + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number +*/ +function rand() { + if ( bernoulli( 0.8 ) > 0 ) { + return uniform( -100.0, 100.0 ); + } + return NaN; +} + +function constant() { + return 0.0; +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, options.dtype, rand ); + var y = filledarrayBy( len, options.dtype, constant ); + return benchmark; + + function benchmark( b ) { + var v; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + x[ 0 ] += 1.0; + v = dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+':ndarray:len='+len, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js new file mode 100644 index 000000000000..9801d4136007 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js @@ -0,0 +1,122 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var bench = require( '@stdlib/bench' ); +var uniform = require( '@stdlib/random/base/uniform' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var tryRequire = require( '@stdlib/utils/try-require' ); +var pkg = require( './../package.json' ).name; + + +// VARIABLES // + +var options = { + 'dtype': 'float64' +}; +var dnancusumkbn = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( dnancusumkbn instanceof Error ) +}; + + +// FUNCTIONS // + +/** +* Returns a random number. +* +* @private +* @returns {number} random number +*/ +function rand() { + if ( bernoulli( 0.8 ) > 0 ) { + return uniform( -100.0, 100.0 ); + } + return NaN; +} + +function constant() { + return 0.0; +} + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var x = filledarrayBy( len, options.dtype, rand ); + var y = filledarrayBy( len, options.dtype, constant ); + return benchmark; + + function benchmark( b ) { + var v; + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + x[ 0 ] += 1.0; + v = dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( v[ i % len ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var min; + var max; + var f; + var i; + + min = 1; // 10^min + max = 6; // 10^max + + for ( i = min; i <= max; i++ ) { + len = pow( 10, i ); + f = createBenchmark( len ); + bench( pkg+'::native:ndarray:len='+len, opts, f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile new file mode 100644 index 000000000000..7280962b4c4d --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2020 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := benchmark.length.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled benchmarks. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c new file mode 100644 index 000000000000..a07bf06d6c4c --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c @@ -0,0 +1,201 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/dnancusumkbn.h" +#include +#include +#include +#include +#include + +#define NAME "dnancusumkbn" +#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[ len ]; + double y[ len ]; + double t; + int i; + + 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++ ) { + x[ 0 ] += 1.0; + stdlib_strided_dnancusumkbn( len, 0.0, x, 1, y, 1 ); + if ( y[ 0 ] != y[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( y[ len-1 ] != y[ len-1 ] ) { + printf( "should not return NaN\n" ); + } + 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[ len ]; + double y[ len ]; + double t; + int i; + + for ( i = 0; i < len; i++ ) { + x[ i ] = ( rand_double() * 20000.0 ) - 10000.0; + y[ i ] = 0.0; + } + t = tic(); + for ( i = 0; i < iterations; i++ ) { + x[ 0 ] += 1.0; + stdlib_strided_dnancusumkbn_ndarray( len, 0.0, x, 1, 0, y, 1, 0 ); + if ( y[ 0 ] != y[ 0 ] ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( y[ len-1 ] != y[ len-1 ] ) { + printf( "should not return NaN\n" ); + } + 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/blas/ext/base/dnancusumkbn/binding.gyp b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp new file mode 100644 index 000000000000..7d0005b2e390 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp @@ -0,0 +1,170 @@ +# @license Apache-2.0 +# +# Copyright (c) 2020 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. + +# A `.gyp` file for building a Node.js native add-on. +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # List of files to include in this file: + 'includes': [ + './include.gypi', + ], + + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Target name should match the add-on export name: + 'addon_target_name%': 'addon', + + # Set variables based on the host OS: + 'conditions': [ + [ + 'OS=="win"', + { + # Define the object file suffix: + 'obj': 'obj', + }, + { + # Define the object file suffix: + 'obj': 'o', + } + ], # end condition (OS=="win") + ], # end conditions + }, # end variables + + # Define compile targets: + 'targets': [ + + # Target to generate an add-on: + { + # The target name should match the add-on export name: + 'target_name': '<(addon_target_name)', + + # Define dependencies: + 'dependencies': [], + + # Define directories which contain relevant include headers: + 'include_dirs': [ + # Local include directory: + '<@(include_dirs)', + ], + + # List of source files: + 'sources': [ + '<@(src_files)', + ], + + # Settings which should be applied when a target's object files are used as linker input: + 'link_settings': { + # Define libraries: + 'libraries': [ + '<@(libraries)', + ], + + # Define library directories: + 'library_dirs': [ + '<@(library_dirs)', + ], + }, + + # C/C++ compiler flags: + 'cflags': [ + # Enable commonly used warning options: + '-Wall', + + # Aggressive optimization: + '-O3', + ], + + # C specific compiler flags: + 'cflags_c': [ + # Specify the C standard to which a program is expected to conform: + '-std=c99', + ], + + # C++ specific compiler flags: + 'cflags_cpp': [ + # Specify the C++ standard to which a program is expected to conform: + '-std=c++11', + ], + + # Linker flags: + 'ldflags': [], + + # Apply conditions based on the host OS: + 'conditions': [ + [ + 'OS=="mac"', + { + # Linker flags: + 'ldflags': [ + '-undefined dynamic_lookup', + '-Wl,-no-pie', + '-Wl,-search_paths_first', + ], + }, + ], # end condition (OS=="mac") + [ + 'OS!="win"', + { + # C/C++ flags: + 'cflags': [ + # Generate platform-independent code: + '-fPIC', + ], + }, + ], # end condition (OS!="win") + ], # end conditions + }, # end target <(addon_target_name) + + # Target to copy a generated add-on to a standard location: + { + 'target_name': 'copy_addon', + + # Declare that the output of this target is not linked: + 'type': 'none', + + # Define dependencies: + 'dependencies': [ + # Require that the add-on be generated before building this target: + '<(addon_target_name)', + ], + + # Define a list of actions: + 'actions': [ + { + 'action_name': 'copy_addon', + 'message': 'Copying addon...', + + # Explicitly list the inputs in the command-line invocation below: + 'inputs': [], + + # Declare the expected outputs: + 'outputs': [ + '<(addon_output_dir)/<(addon_target_name).node', + ], + + # Define the command-line invocation: + 'action': [ + 'cp', + '<(PRODUCT_DIR)/<(addon_target_name).node', + '<(addon_output_dir)/<(addon_target_name).node', + ], + }, + ], # end actions + }, # end target copy_addon + ], # end targets +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt new file mode 100644 index 000000000000..bd2b8d818ec3 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt @@ -0,0 +1,120 @@ + +{{alias}}( N, sum, x, strideX, y, strideY ) + Computes the cumulative sum of double-precision floating-point strided array + elements using an improved Kahan–Babuška algorithm, 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. + + sum: number + Initial sum. + + 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 ] ); + > var y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}( x.length, 0.0, x, 1, y, 1 ) + [ 1.0, 1.0, 3.0 ] + + // Using `N` and stride parameters: + > x = new {{alias:@stdlib/array/float64}}( [ -2.0, 1.0, 1.0, NaN, 2.0, -1.0 ] ); + > y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}( 3, 0.0, x, 2, y, 2 ) + [ -2.0, 0.0, -1.0, 0.0, 1.0, 0.0 ] + + // Using view offsets: + > var x0 = new {{alias:@stdlib/array/float64}}( [ 1.0, -2.0, 3.0, NaN, 5.0, -1.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, 0.0, x1, 2, y1, 1 ) + [ -2.0, -2.0, -3.0 ] + > y0 + [ 0.0, 0.0, 0.0, -2.0, -2.0, -3.0 ] + + +{{alias}}.ndarray( N, sum, x, strideX, offsetX, y, strideY, offsetY ) + Computes the cumulative sum of double-precision floating-point strided array + elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values + and alternative indexing semantics. + + While typed array views mandate a view offset based on the underlying + buffer, the offset parameters support indexing semantics based on starting + indices. + + Parameters + ---------- + N: integer + Number of indexed elements. + + sum: number + Initial sum. + + 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, -2.0, NaN ] ); + > var y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0 ) + [ 1.0, -1.0, -1.0 ] + + // Advanced indexing: + > x = new {{alias:@stdlib/array/float64}}( [ NaN, -2.0, 3.0, 2.0, 5.0, -1.0 ] ); + > y = new {{alias:@stdlib/array/float64}}( x.length ); + > {{alias}}.ndarray( 3, 0.0, x, 2, 1, y, -1, y.length-1 ) + [ 0.0, 0.0, 0.0, -1.0, 0.0, -2.0 ] + + See Also + -------- + diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts new file mode 100644 index 000000000000..f2dced279706 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts @@ -0,0 +1,106 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +// TypeScript Version: 4.1 + +/** +* Interface describing `dnancusumkbn`. +*/ +interface Routine { + /** + * Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. + * + * @param N - number of indexed elements + * @param sum - initial sum + * @param x - input array + * @param strideX - stride length for `x` + * @param y - output array + * @param strideY - stride length for `y` + * @returns output array + * + * @example + * var Float64Array = require( '@stdlib/array/float64' ); + * + * var x = new Float64Array( [ 1.0, NaN, 2.0 ] ); + * var y = new Float64Array( x.length ); + * + * dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + * // y => [ 1.0, 1.0, 3.0 ] + */ + ( N: number, sum: number, x: Float64Array, strideX: number, y: Float64Array, strideY: number ): Float64Array; + + /** + * Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. + * + * @param N - number of indexed elements + * @param sum - initial sum + * @param x - input array + * @param strideX - stride length for `x` + * @param offsetX - starting index for `x` + * @param y - output array + * @param strideY - stride length for `y` + * @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 ] ); + * var y = new Float64Array( x.length ); + * + * dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0 ); + * // y => [ 1.0, -1.0, -1.0 ] + */ + ndarray( N: number, sum: number, x: Float64Array, strideX: number, offsetX: number, y: Float64Array, strideY: number, offsetY: number ): Float64Array; +} + +/** +* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* +* @param N - number of indexed elements +* @param sum - initial sum +* @param x - input array +* @param strideX - stride length for `x` +* @param y - output array +* @param strideY - stride length for `y` +* @returns output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN ] ); +* var y = new Float64Array( x.length ); +* +* dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); +* // y => [ 1.0, -1.0, -1.0 ] +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN ] ); +* var y = new Float64Array( x.length ); +* +* dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0 ); +* // y => [ 1.0, -1.0, -1.0 ] +*/ +declare var dnancusumkbn: Routine; + + +// EXPORTS // + +export = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts new file mode 100644 index 000000000000..d97f43361438 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts @@ -0,0 +1,280 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +import dnancusumkbn = require( './index' ); + + +// TESTS // + +// The function returns a Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn( x.length, 0.0, 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 ); + + dnancusumkbn( '10', 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( true, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( false, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( null, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( undefined, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( [], 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( {}, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( ( x: number ): number => x, 0.0, x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a second argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn( x.length, '10', 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, true, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, false, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, null, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, undefined, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, [], 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, {}, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, ( x: number ): number => x, 0.0, x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a third argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn( x.length, 0.0, 10, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, '10', 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, true, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, false, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, null, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, undefined, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, [ '1' ], 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, {}, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, ( x: number ): number => x, 1, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fourth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn( x.length, 0.0, x, '10', y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, true, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, false, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, null, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, undefined, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, [], y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, {}, y, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, ( x: number ): number => x, y, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fifth argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + + dnancusumkbn( x.length, 0.0, x, 1, 10, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, '10', 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, true, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, false, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, null, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, undefined, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, [ '1' ], 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, {}, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, ( x: number ): number => x, 1 ); // $ExpectError +} + +// The compiler throws an error if the function is provided a sixth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn( x.length, 0.0, x, 1, y, '10' ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y, true ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y, false ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y, null ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y, undefined ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y, [] ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y, {} ); // $ExpectError + dnancusumkbn( x.length, 0.0, 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 ); + + dnancusumkbn(); // $ExpectError + dnancusumkbn( x.length ); // $ExpectError + dnancusumkbn( x.length, 0.0 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1 ); // $ExpectError + dnancusumkbn( x.length, 0.0, x, 1, y ); // $ExpectError + dnancusumkbn( x.length, 0.0, 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 ); + + dnancusumkbn.ndarray( x.length, 0.0, 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 ); + + dnancusumkbn.ndarray( '10', 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( true, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( false, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( null, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( undefined, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( [], 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( {}, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( ( x: number ): number => x, 0.0, 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 number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn.ndarray( x.length, '10', 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, true, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, false, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, null, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, undefined, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, [], 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, {}, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, ( x: number ): number => x, 0.0, 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 Float64Array... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn.ndarray( x.length, 0.0, 10, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, '10', 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, true, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, false, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, null, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, undefined, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, [ '1' ], 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, {}, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, ( x: number ): number => x, 1, 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 ); + + dnancusumkbn.ndarray( x.length, 0.0, x, '10', 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, true, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, false, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, null, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, undefined, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, [], 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, {}, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, ( x: number ): number => x, 0, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn.ndarray( x.length, 0.0, x, 1, '10', y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, true, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, false, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, null, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, undefined, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, [], y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, {}, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, ( x: number ): number => x, y, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a Float64Array... +{ + const x = new Float64Array( 10 ); + + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, 10, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, '10', 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, true, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, false, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, null, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, undefined, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, [ '1' ], 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, {}, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, ( x: number ): number => x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, '10', 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, true, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, false, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, null, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, undefined, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, [], 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, {}, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an eighth argument which is not a number... +{ + const x = new Float64Array( 10 ); + const y = new Float64Array( 10 ); + + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, '10' ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, true ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, false ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, null ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, undefined ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, [] ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, {} ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, 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 ); + + dnancusumkbn.ndarray(); // $ExpectError + dnancusumkbn.ndarray( x.length ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1 ); // $ExpectError + dnancusumkbn.ndarray( x.length, 0.0, x, 1, 0, y, 1, 0, 10 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile new file mode 100644 index 000000000000..ff5293d3059f --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile @@ -0,0 +1,146 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2020 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + +# Define the program used for compiling C source files: +ifdef C_COMPILER + CC := $(C_COMPILER) +else + CC := gcc +endif + +# Define the command-line options when compiling C files: +CFLAGS ?= \ + -std=c99 \ + -O3 \ + -Wall \ + -pedantic + +# Determine whether to generate position independent code ([1][1], [2][2]). +# +# [1]: https://gcc.gnu.org/onlinedocs/gcc/Code-Gen-Options.html#Code-Gen-Options +# [2]: http://stackoverflow.com/questions/5311515/gcc-fpic-option +ifeq ($(OS), WINNT) + fPIC ?= +else + fPIC ?= -fPIC +endif + +# List of includes (e.g., `-I /foo/bar -I /beep/boop/include`): +INCLUDE ?= + +# List of source files: +SOURCE_FILES ?= + +# List of libraries (e.g., `-lopenblas -lpthread`): +LIBRARIES ?= + +# List of library paths (e.g., `-L /foo/bar -L /beep/boop`): +LIBPATH ?= + +# List of C targets: +c_targets := example.out + + +# RULES # + +#/ +# Compiles source files. +# +# @param {string} [C_COMPILER] - C compiler (e.g., `gcc`) +# @param {string} [CFLAGS] - C compiler options +# @param {(string|void)} [fPIC] - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} [INCLUDE] - list of includes (e.g., `-I /foo/bar -I /beep/boop/include`) +# @param {string} [SOURCE_FILES] - list of source files +# @param {string} [LIBPATH] - list of library paths (e.g., `-L /foo/bar -L /beep/boop`) +# @param {string} [LIBRARIES] - list of libraries (e.g., `-lopenblas -lpthread`) +# +# @example +# make +# +# @example +# make all +#/ +all: $(c_targets) + +.PHONY: all + +#/ +# Compiles C source files. +# +# @private +# @param {string} CC - C compiler (e.g., `gcc`) +# @param {string} CFLAGS - C compiler options +# @param {(string|void)} fPIC - compiler flag determining whether to generate position independent code (e.g., `-fPIC`) +# @param {string} INCLUDE - list of includes (e.g., `-I /foo/bar`) +# @param {string} SOURCE_FILES - list of source files +# @param {string} LIBPATH - list of library paths (e.g., `-L /foo/bar`) +# @param {string} LIBRARIES - list of libraries (e.g., `-lopenblas`) +#/ +$(c_targets): %.out: %.c + $(QUIET) $(CC) $(CFLAGS) $(fPIC) $(INCLUDE) -o $@ $(SOURCE_FILES) $< $(LIBPATH) -lm $(LIBRARIES) + +#/ +# Runs compiled examples. +# +# @example +# make run +#/ +run: $(c_targets) + $(QUIET) ./$< + +.PHONY: run + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: + $(QUIET) -rm -f *.o *.out + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c new file mode 100644 index 000000000000..5268dbd01293 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c @@ -0,0 +1,41 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/dnancusumkbn.h" +#include + +int main( void ) { + // Create strided arrays: + const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 0.0/0.0, 0.0/0.0, 8.0 }; + double y[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 }; + + // Specify the number of elements: + const int N = 4; + + // Specify stride lengths: + const int strideX = 2; + const int strideY = -2; + + // Compute the cumulative sum: + stdlib_strided_dnancusumkbn( N, 0.0, x, strideX, y, strideY ); + + // Print the result: + for ( int i = 0; i < 8; i++ ) { + printf( "y[ %d ] = %lf\n", i, y[ i ] ); + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js new file mode 100644 index 000000000000..82e02721694b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js @@ -0,0 +1,44 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ); +var bernoulli = require( '@stdlib/random/base/bernoulli' ); +var filledarrayBy = require( '@stdlib/array/filled-by' ); +var dnancusumkbn = require( './../lib' ); + +function rand() { + if ( bernoulli( 0.7 ) > 0 ) { + return discreteUniform( 0, 100 ); + } + return NaN; +} + +function constant() { + return 0.0; +} + +var x = filledarrayBy( 10, 'float64', rand ); +console.log( x ); + +var y = filledarrayBy( 10, 'float64', constant ); +console.log( y ); + +dnancusumkbn( x.length, 0.0, x, 1, y, -1 ); +console.log( y ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi new file mode 100644 index 000000000000..26476a8c2655 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi @@ -0,0 +1,53 @@ +# @license Apache-2.0 +# +# Copyright (c) 2020 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. + +# A GYP include file for building a Node.js native add-on. +# +# Main documentation: +# +# [1]: https://gyp.gsrc.io/docs/InputFormatReference.md +# [2]: https://gyp.gsrc.io/docs/UserDocumentation.md +{ + # Define variables to be used throughout the configuration for all targets: + 'variables': { + # Source directory: + 'src_dir': './src', + + # Include directories: + 'include_dirs': [ + '[ 1.0, -1.0, -1.0 ] +*/ +function dnancusumkbn( N, sum, x, strideX, y, strideY ) { + ndarray( N, sum, x, strideX, stride2offset( N, strideX ), y, strideY, stride2offset( N, strideY ) ); // eslint-disable-line max-len + return y; +} + + +// EXPORTS // + +module.exports = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js new file mode 100644 index 000000000000..1e68a2cf287d --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js @@ -0,0 +1,56 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} sum - initial sum +* @param {Float64Array} x - input array +* @param {integer} strideX - stride length for `x` +* @param {Float64Array} y - output array +* @param {integer} strideY - stride length for `y` +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN ] ); +* var y = new Float64Array( x.length ); +* +* var v = dnancusumkbn( 3, 0.0, x, 1, y, 1 ); +* // returns [ 1.0, -1.0, -1.0 ] +*/ +function dnancusumkbn( N, sum, x, strideX, y, strideY ) { + addon( N, sum, x, strideX, y, strideY ); + return y; +} + + +// EXPORTS // + +module.exports = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js new file mode 100644 index 000000000000..43943848a3cf --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js @@ -0,0 +1,70 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +/** +* Compute the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* +* @module @stdlib/blas/ext/base/dnancusumkbn +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); +* +* var x = new Float64Array( [ 1.0, -2.0, NaN ] ); +* var y = new Float64Array( x.length ); +* +* dnancusumkbn( 3, 0.0, x, 1, y, 1 ); +* // y => [ 1.0, -1.0, -1.0 ] +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); +* +* var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, NaN, 3.0, NaN ] ); +* var y = new Float64Array( x.length ); +* +* dnancusumkbn.ndarray( 4, 0.0, x, 2, 1, y, 1, 0 ); +* // y => [ 1.0, -1.0, -1.0, -1.0, 0.0, 0.0, 0.0, 0.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 dnancusumkbn; +var tmp = tryRequire( join( __dirname, './native.js' ) ); +if ( isError( tmp ) ) { + dnancusumkbn = main; +} else { + dnancusumkbn = tmp; +} + + +// EXPORTS // + +module.exports = dnancusumkbn; + +// exports: { "ndarray": "dnancusumkbn.ndarray" } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js new file mode 100644 index 000000000000..2bd38e738fe7 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var dnancusumkbn = require( './dnancusumkbn.js' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +setReadOnly( dnancusumkbn, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js new file mode 100644 index 000000000000..cf37a0d329c6 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var dnancusumkbn = require( './dnancusumkbn.native.js' ); +var ndarray = require( './ndarray.native.js' ); + + +// MAIN // + +setReadOnly( dnancusumkbn, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js new file mode 100644 index 000000000000..626a3e6915c0 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js @@ -0,0 +1,92 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var abs = require( '@stdlib/math/base/special/abs' ); + + +// MAIN // + +/** +* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} sum - initial sum +* @param {Float64Array} x - input array +* @param {integer} strideX - stride length for `x` +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @param {Float64Array} y - output array +* @param {integer} strideY - stride length `y` +* @param {NonNegativeInteger} offsetY - starting index for `y` +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, NaN, NaN ] ); +* var y = new Float64Array( x.length ); +* +* var v = dnancusumkbn( 4, 0.0, x, 2, 1, y, 1, 0 ); +* // returns [ 1.0, -1.0, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0 ] +*/ +function dnancusumkbn( N, sum, x, strideX, offsetX, y, strideY, offsetY ) { + var ix; + var iy; + var s; + var v; + var t; + var c; + var i; + + if ( N <= 0 ) { + return y; + } + ix = offsetX; + iy = offsetY; + s = sum; + c = 0.0; + for ( i = 0; i < N; i++ ) { + v = x[ ix ]; + if ( isnan( v ) ) { + y[ iy ] = s + c; + ix += strideX; + iy += strideY; + continue; + } + t = s + v; + if ( abs( s ) >= abs( v ) ) { + c += (s-t) + v; + } else { + c += (v-t) + s; + } + s = t; + y[ iy ] = s + c; + ix += strideX; + iy += strideY; + } + return y; +} + + +// EXPORTS // + +module.exports = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js new file mode 100644 index 000000000000..a49bf1163d8e --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js @@ -0,0 +1,58 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var addon = require( './../src/addon.node' ); + + +// MAIN // + +/** +* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* +* @param {PositiveInteger} N - number of indexed elements +* @param {number} sum - initial sum +* @param {Float64Array} x - input array +* @param {integer} strideX - stride length for `x` +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @param {Float64Array} y - output array +* @param {integer} strideY - stride length for `y` +* @param {NonNegativeInteger} offsetY - starting index for `y` +* @returns {Float64Array} output array +* +* @example +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, NaN, NaN ] ); +* var y = new Float64Array( x.length ); +* +* var v = dnancusumkbn( 4, 0.0, x, 2, 1, y, 1, 0 ); +* // returns [ 1.0, -1.0, 1.0, 1.0, 0.0, 0.0, 0.0, 0.0 ] +*/ +function dnancusumkbn( N, sum, x, strideX, offsetX, y, strideY, offsetY ) { + addon.ndarray( N, sum, x, strideX, offsetX, y, strideY, offsetY ); + return y; +} + + +// EXPORTS // + +module.exports = dnancusumkbn; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/manifest.json new file mode 100644 index 000000000000..a9daff554a13 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/manifest.json @@ -0,0 +1,85 @@ +{ + "options": { + "task": "build" + }, + "fields": [ + { + "field": "src", + "resolve": true, + "relative": true + }, + { + "field": "include", + "resolve": true, + "relative": true + }, + { + "field": "libraries", + "resolve": false, + "relative": false + }, + { + "field": "libpath", + "resolve": true, + "relative": false + } + ], + "confs": [ + { + "task": "build", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/special/abs", + "@stdlib/math/base/assert/is-nan", + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-double", + "@stdlib/napi/argv-strided-float64array", + "@stdlib/napi/export" + ] + }, + { + "task": "benchmark", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/special/abs", + "@stdlib/math/base/assert/is-nan", + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset" + ] + }, + { + "task": "examples", + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/math/base/special/abs", + "@stdlib/math/base/assert/is-nan", + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset" + ] + } + ] +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json new file mode 100644 index 000000000000..6187ec3dd3bd --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json @@ -0,0 +1,83 @@ +{ + "name": "@stdlib/blas/ext/base/dnancusumkbn", + "version": "0.0.0", + "description": "Calculate the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, 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", + "blas", + "extended", + "sum", + "total", + "summation", + "cumulative", + "accumulate", + "kbn", + "kahan", + "compensated", + "strided", + "strided array", + "typed", + "array", + "float64", + "double", + "float64array" + ], + "__stdlib__": { + "wasm": false + } +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile new file mode 100644 index 000000000000..dd720a3de8f2 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile @@ -0,0 +1,70 @@ +#/ +# @license Apache-2.0 +# +# Copyright (c) 2020 The Stdlib Authors. +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. +#/ + +# VARIABLES # + +ifndef VERBOSE + QUIET := @ +else + QUIET := +endif + +# Determine the OS ([1][1], [2][2]). +# +# [1]: https://en.wikipedia.org/wiki/Uname#Examples +# [2]: http://stackoverflow.com/a/27776822/2225624 +OS ?= $(shell uname) +ifneq (, $(findstring MINGW,$(OS))) + OS := WINNT +else +ifneq (, $(findstring MSYS,$(OS))) + OS := WINNT +else +ifneq (, $(findstring CYGWIN,$(OS))) + OS := WINNT +else +ifneq (, $(findstring Windows_NT,$(OS))) + OS := WINNT +endif +endif +endif +endif + + +# RULES # + +#/ +# Removes generated files for building an add-on. +# +# @example +# make clean-addon +#/ +clean-addon: + $(QUIET) -rm -f *.o *.node + +.PHONY: clean-addon + +#/ +# Removes generated files. +# +# @example +# make clean +#/ +clean: clean-addon + +.PHONY: clean diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c new file mode 100644 index 000000000000..129e04213423 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c @@ -0,0 +1,68 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/dnancusumkbn.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/napi/argv_double.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, 6 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_DOUBLE( env, sum, argv, 1 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 5 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, Y, N, strideY, argv, 4 ); + API_SUFFIX(stdlib_strided_dnancusumkbn)( N, sum, 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, 8 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_DOUBLE( env, sum, argv, 1 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 ); + STDLIB_NAPI_ARGV_INT64( env, strideY, argv, 6 ); + STDLIB_NAPI_ARGV_INT64( env, offsetY, argv, 7 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, Y, N, strideY, argv, 5 ); + API_SUFFIX(stdlib_strided_dnancusumkbn_ndarray)( N, sum, 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/blas/ext/base/dnancusumkbn/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c new file mode 100644 index 000000000000..25f51baecba8 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c @@ -0,0 +1,105 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/blas/ext/base/dnancusumkbn.h" +#include "stdlib/strided/base/stride2offset.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/math/base/special/abs.h" +#include "stdlib/math/base/assert/is_nan.h" + +/** +* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* +* ## Method +* +* - This implementation uses an "improved Kahan–Babuška algorithm", as described by Neumaier (1974). +* +* ## References +* +* - Neumaier, Arnold. 1974. "Rounding Error Analysis of Some Methods for Summing Finite Sums." _Zeitschrift Für Angewandte Mathematik Und Mechanik_ 54 (1): 39–51. doi:[10.1002/zamm.19740540106](https://doi.org/10.1002/zamm.19740540106). +* +* @param N number of indexed elements +* @param sum initial sum +* @param X input array +* @param strideX stride length for X +* @param Y output array +* @param strideY stride length for Y +*/ +void API_SUFFIX(stdlib_strided_dnancusumkbn)( const CBLAS_INT N, const double sum, 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_dnancusumkbn_ndarray)( N, sum, X, strideX, ox, Y, strideY, oy ); +} + +/** +* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. +* +* ## Method +* +* - This implementation uses an "improved Kahan–Babuška algorithm", as described by Neumaier (1974). +* +* ## References +* +* - Neumaier, Arnold. 1974. "Rounding Error Analysis of Some Methods for Summing Finite Sums." _Zeitschrift Für Angewandte Mathematik Und Mechanik_ 54 (1): 39–51. doi:[10.1002/zamm.19740540106](https://doi.org/10.1002/zamm.19740540106). +* +* @param N number of indexed elements +* @param sum initial sum +* @param X input array +* @param strideX stride length for X +* @param offsetX starting index for X +* @param Y output array +* @param strideY stride length for Y +* @param offsetY starting index for Y +*/ +void API_SUFFIX(stdlib_strided_dnancusumkbn_ndarray)( const CBLAS_INT N, const double sum, 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 s; + double v; + double t; + double c; + + if ( N <= 0 ) { + return; + } + ix = offsetX; + iy = offsetY; + s = sum; + c = 0.0; + for ( i = 0; i < N; i++ ) { + v = X[ ix ]; + if ( stdlib_base_is_nan( v ) ) { + Y[ iy ] = s + c; + ix += strideX; + iy += strideY; + continue; + } + t = s + v; + if ( stdlib_base_abs( s ) >= stdlib_base_abs( v ) ) { + c += (s-t) + v; + } else { + c += (v-t) + s; + } + s = t; + Y[ iy ] = s + c; + ix += strideX; + iy += strideY; + } + return; +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js new file mode 100644 index 000000000000..6e56cc715959 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js @@ -0,0 +1,304 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var dnancusumkbn = require( './../lib/dnancusumkbn.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancusumkbn, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 6', function test( t ) { + t.strictEqual( dnancusumkbn.length, 6, 'returns expected value' ); + t.end(); +}); + +tape( 'the function computes the cumulative sum (ignoring NaN values)', function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, 2.0, NaN, 4.0, NaN ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + expected = new Float64Array([ + 1.0, + 3.0, + 3.0, + 7.0, + 7.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 10.0, x, 1, y, 1 ); + expected = new Float64Array([ + 11.0, + 13.0, + 13.0, + 13.0, + 18.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + expected = new Float64Array([ + 0.0, + 0.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, NaN, 3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 4.0, + 4.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 1.0e100, NaN, -1.0e100 ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0e100, + 1.0e100, + 1.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( 1e3 ); + y = new Float64Array( x.length ); + expected = new Float64Array( x.length ); + for ( i = 0; i < x.length; i++ ) { + x[ i ] = i + 1; + if ( i === 0 ) { + expected[ i ] += x[ i ]; + } else { + expected[ i ] += expected[ i-1 ] + x[ i ]; + } + } + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns expected value' ); + + 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 = dnancusumkbn( x.length, 0.0, 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 ] ); + + dnancusumkbn( -1, 0.0, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancusumkbn( 0, 0.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, + NaN, // 1 + 4.0, + 5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, y, 1 ); + + expected = new Float64Array( [ 1.0, 1.0, 6.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 + NaN, // 2 + NaN, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancusumkbn( 3, 0.0, x, 1, y, 2 ); + + expected = new Float64Array( [ 1.0, 0.0, 3.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([ + NaN, // 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 + ]); + + dnancusumkbn( 3, 0.0, x, -2, y, -1 ); + + expected = new Float64Array( [ 8.0, 8.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, + NaN, // 1 + NaN, + 5.0, // 2 + 6.0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, y, -1 ); + + expected = new Float64Array( [ 6.0, 1.0, 1.0, 0.0, 0.0, 0.0 ] ); + + t.deepEqual( y, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports view offsets', function test( t ) { + var expected; + var x0; + var y0; + var x1; + var y1; + + // Initial arrays... + x0 = new Float64Array([ + 1.0, + NaN, // 2 + NaN, + 4.0, // 1 + 5.0, + 6.0 // 0 + ]); + y0 = new Float64Array([ + 0.0, + 0.0, + 0.0, + 0.0, // 0 + 0.0, // 1 + 0.0 // 2 + ]); + + // Create offset views... + x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // begin at 2nd element + y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // begin at the 4th element + + dnancusumkbn( 3, 0.0, x1, -2, y1, 1 ); + expected = new Float64Array( [ 0.0, 0.0, 0.0, 6.0, 10.0, 10.0 ] ); + + t.deepEqual( y0, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js new file mode 100644 index 000000000000..e2795ec6ce28 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js @@ -0,0 +1,313 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var dnancusumkbn = tryRequire( resolve( __dirname, './../lib/dnancusumkbn.native.js' ) ); +var opts = { + 'skip': ( dnancusumkbn instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancusumkbn, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 6', opts, function test( t ) { + t.strictEqual( dnancusumkbn.length, 6, 'returns expected value' ); + t.end(); +}); + +tape( 'the function computes the cumulative sum (ignoring NaN values)', opts, function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + expected = new Float64Array([ + 1.0, + 3.0, + 3.0, + 3.0, + 8.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 10.0, x, 1, y, 1 ); + expected = new Float64Array([ + 11.0, + 13.0, + 13.0, + 13.0, + 18.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + expected = new Float64Array([ + 0.0, + 0.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, NaN, 3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 4.0, + 4.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 1.0e100, NaN, -1.0e100 ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + + expected = new Float64Array([ + 1.0, + 1.0e100, + 1.0e100, + 1.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( 1e3 ); + y = new Float64Array( x.length ); + expected = new Float64Array( x.length ); + for ( i = 0; i < x.length; i++ ) { + x[ i ] = i + 1; + if ( i === 0 ) { + expected[ i ] += x[ i ]; + } else { + expected[ i ] += expected[ i-1 ] + x[ i ]; + } + } + dnancusumkbn( x.length, 0.0, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the output array', opts, function test( t ) { + var out; + var x; + var 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 = dnancusumkbn( x.length, 0.0, 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 ] ); + + dnancusumkbn( -1, 0.0, x, 1, y, 1 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancusumkbn( 0, 0.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([ + NaN, // 0 + 2.0, + NaN, // 1 + 4.0, + 5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, y, 1 ); + + expected = new Float64Array( [ 0.0, 0.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 + NaN, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancusumkbn( 3, 0.0, x, 1, y, 2 ); + + expected = new Float64Array( [ 1.0, 0.0, 3.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([ + NaN, // 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 + ]); + + dnancusumkbn( 3, 0.0, x, -2, y, -1 ); + + expected = new Float64Array( [ 8.0, 8.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, + NaN, // 1 + NaN, + 5.0, // 2 + 6.0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, y, -1 ); + + expected = new Float64Array( [ 6.0, 1.0, 1.0, 0.0, 0.0, 0.0 ] ); + + t.deepEqual( y, expected, 'returns expected value' ); + t.end(); +}); + +tape( 'the function supports view offsets', opts, function test( t ) { + var expected; + var x0; + var y0; + var x1; + var y1; + + // Initial arrays... + x0 = new Float64Array([ + 1.0, + NaN, // 2 + 3.0, + 4.0, // 1 + 5.0, + 6.0 // 0 + ]); + y0 = new Float64Array([ + 0.0, + 0.0, + 0.0, + 0.0, // 0 + 0.0, // 1 + 0.0 // 2 + ]); + + // Create offset views... + x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // begin at 2nd element + y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // begin at the 4th element + + dnancusumkbn( 3, 0.0, x1, -2, y1, 1 ); + expected = new Float64Array( [ 0.0, 0.0, 0.0, 6.0, 10.0, 10.0 ] ); + + t.deepEqual( y0, expected, 'returns expected value' ); + t.end(); +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js new file mode 100644 index 000000000000..f5e7629c4e1b --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js @@ -0,0 +1,82 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var proxyquire = require( 'proxyquire' ); +var IS_BROWSER = require( '@stdlib/assert/is-browser' ); +var dnancusumkbn = require( './../lib' ); + + +// VARIABLES // + +var opts = { + 'skip': IS_BROWSER +}; + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancusumkbn, '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 dnancusumkbn.ndarray, 'function', 'method is a function' ); + t.end(); +}); + +tape( 'if a native implementation is available, the main export is the native implementation', opts, function test( t ) { + var dnancusumkbn = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( dnancusumkbn, mock, 'returns expected value' ); + t.end(); + + function tryRequire() { + return mock; + } + + function mock() { + // Mock... + } +}); + +tape( 'if a native implementation is not available, the main export is a JavaScript implementation', opts, function test( t ) { + var dnancusumkbn; + var main; + + main = require( './../lib/dnancusumkbn.js' ); + + dnancusumkbn = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( dnancusumkbn, main, 'returns expected value' ); + t.end(); + + function tryRequire() { + return new Error( 'Cannot find module' ); + } +}); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js new file mode 100644 index 000000000000..9b8ff13a29e1 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js @@ -0,0 +1,336 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var dnancusumkbn = require( './../lib/ndarray.js' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancusumkbn, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 8', function test( t ) { + t.strictEqual( dnancusumkbn.length, 8, 'returns expected value' ); + t.end(); +}); + +tape( 'the function calculates the cumulative sum', function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + expected = new Float64Array([ + 1.0, + 3.0, + 3.0, + 3.0, + 8.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 10.0, x, 1, 0, y, 1, 0 ); + expected = new Float64Array([ + 11.0, + 13.0, + 13.0, + 13.0, + 18.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + expected = new Float64Array([ + 0.0, + 0.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, NaN, 3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 4.0, + 4.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 1.0e100, NaN, -1.0e100 ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0e100, + 1.0e100, + 1.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( 1e3 ); + y = new Float64Array( x.length ); + expected = new Float64Array( x.length ); + for ( i = 0; i < x.length; i++ ) { + x[ i ] = i + 1; + if ( i === 0 ) { + expected[ i ] += x[ i ]; + } else { + expected[ i ] += expected[ i-1 ] + x[ i ]; + } + } + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns expected value' ); + + 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 = dnancusumkbn( x.length, 0.0, 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 ] ); + + dnancusumkbn( -1, 0.0, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancusumkbn( 0, 0.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, + NaN, // 1 + NaN, + 5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, 0, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 6.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 + NaN, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancusumkbn( 3, 0.0, x, 1, 0, y, 2, 0 ); + + expected = new Float64Array( [ 1.0, 0.0, 3.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([ + NaN, // 2 + NaN, + 3.0, // 1 + 4.0, + 5.0 // 0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, -2, x.length-1, y, -1, 2 ); + + expected = new Float64Array( [ 8.0, 8.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, + NaN, // 0 + NaN, + -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 + ]); + + dnancusumkbn( 4, 0.0, x, 2, 1, y, 1, 0 ); + + expected = new Float64Array( [ 0.0, -2.0, 0.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 + NaN, // 3 + NaN, + 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 + ]); + + dnancusumkbn( 4, 0.0, x, 1, 0, y, 2, 1 ); + + expected = new Float64Array( [ 0.0, 2.0, 0.0, 3.0, 0.0, 5.0, 0.0, 5.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, + NaN, // 2 + NaN + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, 0, y, -1, 2 ); + + expected = new Float64Array( [ 4.0, 4.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/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js new file mode 100644 index 000000000000..2b8915038bba --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js @@ -0,0 +1,345 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2020 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var resolve = require( 'path' ).resolve; +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var tryRequire = require( '@stdlib/utils/try-require' ); + + +// VARIABLES // + +var dnancusumkbn = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); +var opts = { + 'skip': ( dnancusumkbn instanceof Error ) +}; + + +// TESTS // + +tape( 'main export is a function', opts, function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dnancusumkbn, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 8', opts, function test( t ) { + t.strictEqual( dnancusumkbn.length, 8, 'returns expected value' ); + t.end(); +}); + +tape( 'the function calculates the cumulative sum', opts, function test( t ) { + var expected; + var x; + var y; + var i; + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + expected = new Float64Array([ + 1.0, + 3.0, + 3.0, + 3.0, + 8.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 2.0, NaN, NaN, 5.0 ] ); + y = new Float64Array( x.length ); + + dnancusumkbn( x.length, 10.0, x, 1, 0, y, 1, 0 ); + expected = new Float64Array([ + 11.0, + 13.0, + 13.0, + 13.0, + 18.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ NaN, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + expected = new Float64Array([ + 0.0, + 0.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, NaN, 3.0, NaN ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0, + 4.0, + 4.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( [ 1.0, 1.0e100, NaN, -1.0e100 ] ); + y = new Float64Array( x.length ); + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + + expected = new Float64Array([ + 1.0, + 1.0e100, + 1.0e100, + 1.0 + ]); + t.deepEqual( y, expected, 'returns expected value' ); + + x = new Float64Array( 1e3 ); + y = new Float64Array( x.length ); + expected = new Float64Array( x.length ); + for ( i = 0; i < x.length; i++ ) { + x[ i ] = i + 1; + if ( i === 0 ) { + expected[ i ] += x[ i ]; + } else { + expected[ i ] += expected[ i-1 ] + x[ i ]; + } + } + dnancusumkbn( x.length, 0.0, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a reference to the output array', opts, function test( t ) { + var out; + var x; + var 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 = dnancusumkbn( x.length, 0.0, 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 ] ); + + dnancusumkbn( -1, 0.0, x, 1, 0, y, 1, 0 ); + t.deepEqual( y, expected, 'returns `y` unchanged' ); + + dnancusumkbn( 0, 0.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, + NaN, // 1 + NaN, + 5.0 // 2 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, 0, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, 1.0, 6.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 + NaN, // 2 + 4.0, + 5.0 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, + 0.0, // 1 + 0.0, + 0.0 // 2 + ]); + + dnancusumkbn( 3, 0.0, x, 1, 0, y, 2, 0 ); + + expected = new Float64Array( [ 1.0, 0.0, 3.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([ + NaN, // 2 + NaN, + 3.0, // 1 + 4.0, + 5.0 // 0 + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, -2, x.length-1, y, -1, 2 ); + + expected = new Float64Array( [ 8.0, 8.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, + NaN // 3 + ]); + y = new Float64Array([ + 0.0, // 0 + 0.0, // 1 + 0.0, // 2 + 0.0, // 3 + 0.0, + 0.0, + 0.0, + 0.0 + ]); + + dnancusumkbn( 4, 0.0, x, 2, 1, y, 1, 0 ); + + expected = new Float64Array( [ 1.0, -1.0, 1.0, 1.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 + NaN, // 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 + ]); + + dnancusumkbn( 4, 0.0, x, 1, 0, y, 2, 1 ); + + expected = new Float64Array( [ 0.0, 2.0, 0.0, 3.0, 0.0, 5.0, 0.0, 5.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 + NaN + ]); + y = new Float64Array([ + 0.0, // 2 + 0.0, // 1 + 0.0, // 0 + 0.0, + 0.0, + 0.0 + ]); + + dnancusumkbn( 3, 0.0, x, 2, 0, y, -1, 2 ); + + expected = new Float64Array( [ 9.0, 4.0, 1.0, 0.0, 0.0, 0.0 ] ); + + t.deepEqual( y, expected, 'returns expected value' ); + t.end(); +}); From 2dd5bff80a965ede26dee68b42090da6fb85b084 Mon Sep 17 00:00:00 2001 From: stdlib-bot <82920195+stdlib-bot@users.noreply.github.com> Date: Thu, 6 Mar 2025 08:41:59 +0000 Subject: [PATCH 2/4] chore: update copyright years --- lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js | 2 +- .../blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js | 2 +- .../blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js | 2 +- .../ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile | 2 +- .../blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c | 2 +- lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/examples/index.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/include.gypi | 2 +- .../dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js | 2 +- .../blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/lib/index.js | 2 +- lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/lib/native.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/src/Makefile | 2 +- lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c | 2 +- lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c | 2 +- .../blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js | 2 +- .../blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/test/test.js | 2 +- .../@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js | 2 +- .../blas/ext/base/dnancusumkbn/test/test.ndarray.native.js | 2 +- 30 files changed, 30 insertions(+), 30 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md index 9c51ab628e47..b3ac46e9491b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md @@ -2,7 +2,7 @@ @license Apache-2.0 -Copyright (c) 2020 The Stdlib Authors. +Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js index 3c0ab69e2d54..76d73cffb5ef 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js index d37500e9e1d5..8bcb8f2df231 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js index f0c4d98c2f9d..ab27306bbd05 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js index 9801d4136007..6f3b88009a1f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile index 7280962b4c4d..cce2c865d7ad 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2020 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c index a07bf06d6c4c..08e6ff29cc8c 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp index 7d0005b2e390..68a1ca11d160 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/binding.gyp @@ -1,6 +1,6 @@ # @license Apache-2.0 # -# Copyright (c) 2020 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts index f2dced279706..3c0037bcc20f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts @@ -1,7 +1,7 @@ /* * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts index d97f43361438..393dcb62e470 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts @@ -1,7 +1,7 @@ /* * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile index ff5293d3059f..25ced822f96a 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2020 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c index 5268dbd01293..019e16e6d8f7 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/c/example.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js index 82e02721694b..affbf15c747c 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi index 26476a8c2655..ecfaf82a3279 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include.gypi @@ -1,6 +1,6 @@ # @license Apache-2.0 # -# Copyright (c) 2020 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h index a990e5560b5c..80f5133bcd3e 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js index 8083f5532a1a..61bc53625786 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js index 1e68a2cf287d..284ee54da979 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js index 43943848a3cf..fa186f3a0684 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js index 2bd38e738fe7..a74c95825542 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/main.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js index cf37a0d329c6..d2e94799e6eb 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js index 626a3e6915c0..eaa06ded1094 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js index a49bf1163d8e..4609e6af5f10 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile index dd720a3de8f2..7733b6180cb4 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/Makefile @@ -1,7 +1,7 @@ #/ # @license Apache-2.0 # -# Copyright (c) 2020 The Stdlib Authors. +# Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c index 129e04213423..4ab8661d6501 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/addon.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c index 25f51baecba8..e3ad44ff099b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js index 6e56cc715959..9f767da6a2e4 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js index e2795ec6ce28..162adb4a9c6e 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js index f5e7629c4e1b..7956fcf860ca 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js index 9b8ff13a29e1..5aec3d82abd2 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js index 2b8915038bba..ea409883b845 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2025 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. From 115b6655ceaea78d36ad5dea3ff9d1ca62cd53f0 Mon Sep 17 00:00:00 2001 From: Athan Date: Wed, 14 May 2025 02:27:22 -0700 Subject: [PATCH 3/4] docs: fix description Signed-off-by: Athan --- lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md index b3ac46e9491b..fec3839980af 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md @@ -98,7 +98,7 @@ dnancusumkbn( 4, 0.0, x1, -2, y1, 1 ); #### dnancusumkbn.ndarray( N, sum, x, strideX, offsetX, y, strideY, offsetY ) -Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. +Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm and alternative indexing semantics. ```javascript var Float64Array = require( '@stdlib/array/float64' ); From 80a583b2f48e55930f7a02c0d5f92d13d579f1d9 Mon Sep 17 00:00:00 2001 From: Athan Date: Wed, 14 May 2025 02:51:26 -0700 Subject: [PATCH 4/4] chore: clean-up --- 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: na - 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 --- --- type: pre_push_report description: Results of running various checks prior to pushing changes. report: - task: run_javascript_examples status: na - task: run_c_examples status: na - task: run_cpp_examples status: na - task: run_javascript_readme_examples status: na - task: run_c_benchmarks status: na - task: run_cpp_benchmarks status: na - task: run_fortran_benchmarks status: na - task: run_javascript_benchmarks status: na - task: run_julia_benchmarks status: na - task: run_python_benchmarks status: na - task: run_r_benchmarks status: na - task: run_javascript_tests status: na --- --- .../blas/ext/base/dnancusumkbn/README.md | 37 ++++--------------- .../base/dnancusumkbn/benchmark/benchmark.js | 7 +--- .../benchmark/benchmark.native.js | 7 +--- .../benchmark/benchmark.ndarray.js | 7 +--- .../benchmark/benchmark.ndarray.native.js | 7 +--- .../benchmark/c/benchmark.length.c | 6 ++- .../blas/ext/base/dnancusumkbn/docs/repl.txt | 5 ++- .../base/dnancusumkbn/docs/types/index.d.ts | 6 +-- .../ext/base/dnancusumkbn/docs/types/test.ts | 32 ++++++++-------- .../ext/base/dnancusumkbn/examples/index.js | 7 +--- .../stdlib/blas/ext/base/dnancusumkbn.h | 4 +- .../ext/base/dnancusumkbn/lib/dnancusumkbn.js | 2 +- .../dnancusumkbn/lib/dnancusumkbn.native.js | 2 +- .../blas/ext/base/dnancusumkbn/lib/index.js | 2 +- .../blas/ext/base/dnancusumkbn/lib/ndarray.js | 2 +- .../base/dnancusumkbn/lib/ndarray.native.js | 2 +- .../blas/ext/base/dnancusumkbn/package.json | 2 +- .../blas/ext/base/dnancusumkbn/src/main.c | 4 +- .../dnancusumkbn/test/test.dnancusumkbn.js | 2 +- .../test/test.dnancusumkbn.native.js | 2 +- .../base/dnancusumkbn/test/test.ndarray.js | 2 +- .../dnancusumkbn/test/test.ndarray.native.js | 2 +- 22 files changed, 59 insertions(+), 90 deletions(-) diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md index fec3839980af..e3b4adb02c28 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/README.md @@ -20,7 +20,7 @@ limitations under the License. # dnancusumkbn -> Calculate the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +> Calculate the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm.
@@ -38,7 +38,7 @@ var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); #### dnancusumkbn( N, sum, x, strideX, y, strideY ) -Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. ```javascript var Float64Array = require( '@stdlib/array/float64' ); @@ -115,7 +115,7 @@ 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 sum of every other value in the strided input array starting from the second value and to store in the last `N` elements of the strided output array starting from the last element: +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 sum of every other element in the strided input array starting from the second element and to store in the last `N` elements of the strided output array starting from the last element: ```javascript var Float64Array = require( '@stdlib/array/float64' ); @@ -148,10 +148,11 @@ dnancusumkbn.ndarray( 4, 0.0, x, 2, 1, y, -1, y.length-1 ); ```javascript -var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ); var bernoulli = require( '@stdlib/random/base/bernoulli' ); var filledarrayBy = require( '@stdlib/array/filled-by' ); +var zeros = require( '@stdlib/array/zeros' ); +var dnancusumkbn = require( '@stdlib/blas/ext/base/dnancusumkbn' ); function rand() { if ( bernoulli( 0.7 ) > 0 ) { @@ -160,14 +161,10 @@ function rand() { return NaN; } -function constant() { - return 0.0; -} - var x = filledarrayBy( 10, 'float64', rand ); console.log( x ); -var y = filledarrayBy( 10, 'float64', constant ); +var y = zeros( 10, 'float64' ); console.log( y ); dnancusumkbn( x.length, 0.0, x, 1, y, -1 ); @@ -206,7 +203,7 @@ console.log( y ); #### stdlib_strided_dnancusumkbn( N, sum, \*X, strideX, \*Y, strideY ) -Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. ```c const double x[] = { 1.0, 2.0, 3.0, 4.0 } @@ -234,7 +231,7 @@ void stdlib_strided_dnancusumkbn( const CBLAS_INT N, const double sum, const dou -Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. +Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm and alternative indexing semantics. ```c const double x[] = { 1.0, 2.0, 3.0, 4.0 } @@ -324,14 +321,6 @@ int main( void ) {
-* * * - -## See Also - -- [`@stdlib/blas/ext/base/dcusum`][@stdlib/blas/ext/base/dcusum]: calculate the cumulative sum of double-precision floating-point strided array elements. -- [`@stdlib/blas/ext/base/gcusumkbn`][@stdlib/blas/ext/base/gcusumkbn]: calculate the cumulative sum of strided array elements using an improved Kahan–Babuška algorithm. -- [`@stdlib/blas/ext/base/scusumkbn`][@stdlib/blas/ext/base/scusumkbn]: calculate the cumulative sum of single-precision floating-point strided array elements using an improved Kahan–Babuška algorithm. -
@@ -346,16 +335,6 @@ int main( void ) { [@neumaier:1974a]: https://doi.org/10.1002/zamm.19740540106 - - -[@stdlib/blas/ext/base/dcusum]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/dcusum - -[@stdlib/blas/ext/base/gcusumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/gcusumkbn - -[@stdlib/blas/ext/base/scusumkbn]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/blas/ext/base/scusumkbn - - -
diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js index 76d73cffb5ef..d7a655ccbb7b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.js @@ -24,6 +24,7 @@ var bench = require( '@stdlib/bench' ); var uniform = require( '@stdlib/random/base/uniform' ); var bernoulli = require( '@stdlib/random/base/bernoulli' ); var filledarrayBy = require( '@stdlib/array/filled-by' ); +var zeros = require( '@stdlib/array/zeros' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var pkg = require( './../package.json' ).name; @@ -52,10 +53,6 @@ function rand() { return NaN; } -function constant() { - return 0.0; -} - /** * Creates a benchmark function. * @@ -65,7 +62,7 @@ function constant() { */ function createBenchmark( len ) { var x = filledarrayBy( len, options.dtype, rand ); - var y = filledarrayBy( len, options.dtype, constant ); + var y = zeros( len, options.dtype ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js index 8bcb8f2df231..30d8fc4b22e2 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.native.js @@ -25,6 +25,7 @@ var bench = require( '@stdlib/bench' ); var uniform = require( '@stdlib/random/base/uniform' ); var bernoulli = require( '@stdlib/random/base/bernoulli' ); var filledarrayBy = require( '@stdlib/array/filled-by' ); +var zeros = require( '@stdlib/array/zeros' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var tryRequire = require( '@stdlib/utils/try-require' ); @@ -57,10 +58,6 @@ function rand() { return NaN; } -function constant() { - return 0.0; -} - /** * Creates a benchmark function. * @@ -70,7 +67,7 @@ function constant() { */ function createBenchmark( len ) { var x = filledarrayBy( len, options.dtype, rand ); - var y = filledarrayBy( len, options.dtype, constant ); + var y = zeros( len, options.dtype ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js index ab27306bbd05..44ab9c4b9fda 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.js @@ -24,6 +24,7 @@ var bench = require( '@stdlib/bench' ); var uniform = require( '@stdlib/random/base/uniform' ); var bernoulli = require( '@stdlib/random/base/bernoulli' ); var filledarrayBy = require( '@stdlib/array/filled-by' ); +var zeros = require( '@stdlib/array/zeros' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var pkg = require( './../package.json' ).name; @@ -52,10 +53,6 @@ function rand() { return NaN; } -function constant() { - return 0.0; -} - /** * Creates a benchmark function. * @@ -65,7 +62,7 @@ function constant() { */ function createBenchmark( len ) { var x = filledarrayBy( len, options.dtype, rand ); - var y = filledarrayBy( len, options.dtype, constant ); + var y = zeros( len, options.dtype ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js index 6f3b88009a1f..85f21847dca6 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/benchmark.ndarray.native.js @@ -25,6 +25,7 @@ var bench = require( '@stdlib/bench' ); var uniform = require( '@stdlib/random/base/uniform' ); var bernoulli = require( '@stdlib/random/base/bernoulli' ); var filledarrayBy = require( '@stdlib/array/filled-by' ); +var zeros = require( '@stdlib/array/zeros' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var tryRequire = require( '@stdlib/utils/try-require' ); @@ -57,10 +58,6 @@ function rand() { return NaN; } -function constant() { - return 0.0; -} - /** * Creates a benchmark function. * @@ -70,7 +67,7 @@ function constant() { */ function createBenchmark( len ) { var x = filledarrayBy( len, options.dtype, rand ); - var y = filledarrayBy( len, options.dtype, constant ); + var y = zeros( len, options.dtype ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c index 08e6ff29cc8c..b9b1f65f7e4f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/benchmark/c/benchmark.length.c @@ -140,7 +140,11 @@ static double benchmark2( int iterations, int len ) { int i; for ( i = 0; i < len; i++ ) { - x[ i ] = ( rand_double() * 20000.0 ) - 10000.0; + 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(); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt index bd2b8d818ec3..e9404a1c6b9b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/repl.txt @@ -1,7 +1,8 @@ {{alias}}( N, sum, x, strideX, y, strideY ) Computes the cumulative sum of double-precision floating-point strided array - elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. + elements ignoring `NaN` values and using an improved Kahan–Babuška + algorithm. The `N` and stride parameters determine which elements in the strided arrays are accessed at runtime. @@ -63,7 +64,7 @@ {{alias}}.ndarray( N, sum, x, strideX, offsetX, y, strideY, offsetY ) Computes the cumulative sum of double-precision floating-point strided array - elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values + elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm and alternative indexing semantics. While typed array views mandate a view offset based on the underlying diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts index 3c0037bcc20f..bb6b235b389e 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/index.d.ts @@ -23,7 +23,7 @@ */ interface Routine { /** - * Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. + * Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * @param N - number of indexed elements * @param sum - initial sum @@ -45,7 +45,7 @@ interface Routine { ( N: number, sum: number, x: Float64Array, strideX: number, y: Float64Array, strideY: number ): Float64Array; /** - * Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. + * Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm and alternative indexing semantics. * * @param N - number of indexed elements * @param sum - initial sum @@ -70,7 +70,7 @@ interface Routine { } /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * @param N - number of indexed elements * @param sum - initial sum diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts index 393dcb62e470..eacfac0a61d8 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/docs/types/test.ts @@ -49,14 +49,14 @@ import dnancusumkbn = require( './index' ); const x = new Float64Array( 10 ); const y = new Float64Array( 10 ); - dnancusumkbn( x.length, '10', 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, true, 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, false, 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, null, 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, undefined, 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, [], 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, {}, 0.0, x, 1, y, 1 ); // $ExpectError - dnancusumkbn( x.length, ( x: number ): number => x, 0.0, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, '10', x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, true, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, false, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, null, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, undefined, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, [], x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, {}, x, 1, y, 1 ); // $ExpectError + dnancusumkbn( x.length, ( x: number ): number => x, x, 1, y, 1 ); // $ExpectError } // The compiler throws an error if the function is provided a third argument which is not a Float64Array... @@ -162,14 +162,14 @@ import dnancusumkbn = require( './index' ); const x = new Float64Array( 10 ); const y = new Float64Array( 10 ); - dnancusumkbn.ndarray( x.length, '10', 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, true, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, false, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, null, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, undefined, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, [], 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, {}, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError - dnancusumkbn.ndarray( x.length, ( x: number ): number => x, 0.0, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, '10', x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, true, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, false, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, null, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, undefined, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, [], x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, {}, x, 1, 0, y, 1, 0 ); // $ExpectError + dnancusumkbn.ndarray( x.length, ( x: number ): number => x, 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 Float64Array... diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js index affbf15c747c..8a9939435f46 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/examples/index.js @@ -21,6 +21,7 @@ var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ); var bernoulli = require( '@stdlib/random/base/bernoulli' ); var filledarrayBy = require( '@stdlib/array/filled-by' ); +var zeros = require( '@stdlib/array/zeros' ); var dnancusumkbn = require( './../lib' ); function rand() { @@ -30,14 +31,10 @@ function rand() { return NaN; } -function constant() { - return 0.0; -} - var x = filledarrayBy( 10, 'float64', rand ); console.log( x ); -var y = filledarrayBy( 10, 'float64', constant ); +var y = zeros( x.length, 'float64' ); console.log( y ); dnancusumkbn( x.length, 0.0, x, 1, y, -1 ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h index 80f5133bcd3e..9a51b7848dae 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/include/stdlib/blas/ext/base/dnancusumkbn.h @@ -29,12 +29,12 @@ extern "C" { #endif /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. */ void API_SUFFIX(stdlib_strided_dnancusumkbn)( const CBLAS_INT N, const double sum, const double *X, const CBLAS_INT strideX, double *Y, const CBLAS_INT strideY ); /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm and alternative indexing semantics. */ void API_SUFFIX(stdlib_strided_dnancusumkbn_ndarray)( const CBLAS_INT N, const double sum, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, double *Y, const CBLAS_INT strideY, const CBLAS_INT offsetY ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js index 61bc53625786..185a1ee6a2f6 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.js @@ -27,7 +27,7 @@ var ndarray = require( './ndarray.js' ); // MAIN // /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * ## Method * diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js index 284ee54da979..009a026f1fa3 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/dnancusumkbn.native.js @@ -26,7 +26,7 @@ var addon = require( './../src/addon.node' ); // MAIN // /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * @param {PositiveInteger} N - number of indexed elements * @param {number} sum - initial sum diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js index fa186f3a0684..451a8b9f325b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/index.js @@ -19,7 +19,7 @@ 'use strict'; /** -* Compute the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Compute the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * @module @stdlib/blas/ext/base/dnancusumkbn * diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js index eaa06ded1094..42dfbfc1aeaa 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.js @@ -27,7 +27,7 @@ var abs = require( '@stdlib/math/base/special/abs' ); // MAIN // /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * @param {PositiveInteger} N - number of indexed elements * @param {number} sum - initial sum diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js index 4609e6af5f10..de708dbb6048 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/lib/ndarray.native.js @@ -26,7 +26,7 @@ var addon = require( './../src/addon.node' ); // MAIN // /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * @param {PositiveInteger} N - number of indexed elements * @param {number} sum - initial sum diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json index 6187ec3dd3bd..29b16ba02003 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/package.json @@ -1,7 +1,7 @@ { "name": "@stdlib/blas/ext/base/dnancusumkbn", "version": "0.0.0", - "description": "Calculate the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values.", + "description": "Calculate the cumulative sum of double-precision floating-point strided array elements ignoring NaN values and using an improved Kahan–Babuška algorithm.", "license": "Apache-2.0", "author": { "name": "The Stdlib Authors", diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c index e3ad44ff099b..530d61d5fe21 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/src/main.c @@ -23,7 +23,7 @@ #include "stdlib/math/base/assert/is_nan.h" /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm. * * ## Method * @@ -47,7 +47,7 @@ void API_SUFFIX(stdlib_strided_dnancusumkbn)( const CBLAS_INT N, const double su } /** -* Computes the cumulative sum of double-precision floating-point strided array elements using an improved Kahan–Babuška algorithm, ignoring `NaN` values and alternative indexing semantics. +* Computes the cumulative sum of double-precision floating-point strided array elements ignoring `NaN` values and using an improved Kahan–Babuška algorithm and alternative indexing semantics. * * ## Method * diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js index 9f767da6a2e4..8984a9d086d9 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.js @@ -294,7 +294,7 @@ tape( 'the function supports view offsets', function test( t ) { // Create offset views... x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // begin at 2nd element - y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // begin at the 4th element + y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // begin at 4th element dnancusumkbn( 3, 0.0, x1, -2, y1, 1 ); expected = new Float64Array( [ 0.0, 0.0, 0.0, 6.0, 10.0, 10.0 ] ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js index 162adb4a9c6e..b62573a51ac3 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.dnancusumkbn.native.js @@ -303,7 +303,7 @@ tape( 'the function supports view offsets', opts, function test( t ) { // Create offset views... x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // begin at 2nd element - y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // begin at the 4th element + y1 = new Float64Array( y0.buffer, y0.BYTES_PER_ELEMENT*3 ); // begin at 4th element dnancusumkbn( 3, 0.0, x1, -2, y1, 1 ); expected = new Float64Array( [ 0.0, 0.0, 0.0, 6.0, 10.0, 10.0 ] ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js index 5aec3d82abd2..0f6aa70868a6 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.js @@ -38,7 +38,7 @@ tape( 'the function has an arity of 8', function test( t ) { t.end(); }); -tape( 'the function calculates the cumulative sum', function test( t ) { +tape( 'the function calculates the cumulative sum (ignoring NaN values)', function test( t ) { var expected; var x; var y; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js index ea409883b845..4393c38a4eae 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dnancusumkbn/test/test.ndarray.native.js @@ -47,7 +47,7 @@ tape( 'the function has an arity of 8', opts, function test( t ) { t.end(); }); -tape( 'the function calculates the cumulative sum', opts, function test( t ) { +tape( 'the function calculates the cumulative sum (ignoring NaN values)', opts, function test( t ) { var expected; var x; var y;