diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/README.md b/lib/node_modules/@stdlib/stats/base/dmskmin/README.md index 7f15599fe02c..7886fb7cb0e3 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/README.md +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/README.md @@ -38,7 +38,7 @@ var dmskmin = require( '@stdlib/stats/base/dmskmin' ); #### dmskmin( N, x, strideX, mask, strideMask ) -Computes the minimum value of a double-precision floating-point strided array `x` according to a `mask`. +Computes the minimum value of a double-precision floating-point strided array according to a mask. ```javascript var Float64Array = require( '@stdlib/array/float64' ); @@ -55,22 +55,20 @@ The function has the following parameters: - **N**: number of indexed elements. - **x**: input [`Float64Array`][@stdlib/array/float64]. -- **strideX**: index increment for `x`. +- **strideX**: stride length for `x`. - **mask**: mask [`Uint8Array`][@stdlib/array/uint8]. If a `mask` array element is `0`, the corresponding element in `x` is considered valid and **included** in computation. If a `mask` array element is `1`, the corresponding element in `x` is considered invalid/missing and **excluded** from computation. -- **strideMask**: index increment for `mask`. +- **strideMask**: stride length for `mask`. -The `N` and `stride` parameters determine which elements are accessed at runtime. For example, to compute the minimum value of every other element in `x`, +The `N` and stride parameters determine which elements in the strided arrays are accessed at runtime. For example, to compute the minimum value of every other element in `x`, ```javascript var Float64Array = require( '@stdlib/array/float64' ); var Uint8Array = require( '@stdlib/array/uint8' ); -var floor = require( '@stdlib/math/base/special/floor' ); var x = new Float64Array( [ 1.0, 2.0, 7.0, -2.0, -4.0, 3.0, -5.0, -6.0 ] ); var mask = new Uint8Array( [ 0, 0, 0, 0, 0, 0, 1, 1 ] ); -var N = floor( x.length / 2 ); -var v = dmskmin( N, x, 2, mask, 2 ); +var v = dmskmin( 4, x, 2, mask, 2 ); // returns -4.0 ``` @@ -81,7 +79,6 @@ Note that indexing is relative to the first index. To introduce offsets, use [`t ```javascript var Float64Array = require( '@stdlib/array/float64' ); var Uint8Array = require( '@stdlib/array/uint8' ); -var floor = require( '@stdlib/math/base/special/floor' ); var x0 = new Float64Array( [ 2.0, 1.0, -2.0, -2.0, 3.0, 4.0, -5.0, -6.0 ] ); var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element @@ -89,15 +86,13 @@ var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd var mask0 = new Uint8Array( [ 0, 0, 0, 0, 0, 0, 1, 1 ] ); var mask1 = new Uint8Array( mask0.buffer, mask0.BYTES_PER_ELEMENT*1 ); // start at 2nd element -var N = floor( x0.length / 2 ); - -var v = dmskmin( N, x1, 2, mask1, 2 ); +var v = dmskmin( 4, x1, 2, mask1, 2 ); // returns -2.0 ``` #### dmskmin.ndarray( N, x, strideX, offsetX, mask, strideMask, offsetMask ) -Computes the minimum value of a double-precision floating-point strided array according to a `mask` and using alternative indexing semantics. +Computes the minimum value of a double-precision floating-point strided array according to a mask and using alternative indexing semantics. ```javascript var Float64Array = require( '@stdlib/array/float64' ); @@ -115,18 +110,16 @@ The function has the following additional parameters: - **offsetX**: starting index for `x`. - **offsetMask**: starting index for `mask`. -While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying `buffer`, the `offset` parameter supports indexing semantics based on a starting index. For example, to calculate the minimum value for every other value in `x` starting from the second value +While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on a starting indices. For example, to calculate the minimum value for every other element in `x` starting from the second element ```javascript var Float64Array = require( '@stdlib/array/float64' ); var Uint8Array = require( '@stdlib/array/uint8' ); -var floor = require( '@stdlib/math/base/special/floor' ); var x = new Float64Array( [ 2.0, 1.0, -2.0, -2.0, 3.0, 4.0, -5.0, -6.0 ] ); var mask = new Uint8Array( [ 0, 0, 0, 0, 0, 0, 1, 1 ] ); -var N = floor( x.length / 2 ); -var v = dmskmin.ndarray( N, x, 2, 1, mask, 2, 1 ); +var v = dmskmin.ndarray( 4, x, 2, 1, mask, 2, 1 ); // returns -2.0 ``` @@ -151,26 +144,19 @@ var v = dmskmin.ndarray( N, x, 2, 1, mask, 2, 1 ); ```javascript -var randu = require( '@stdlib/random/base/randu' ); -var round = require( '@stdlib/math/base/special/round' ); -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var bernoulli = require( '@stdlib/random/array/bernoulli' ); var dmskmin = require( '@stdlib/stats/base/dmskmin' ); -var mask; -var x; -var i; - -x = new Float64Array( 10 ); -mask = new Uint8Array( x.length ); -for ( i = 0; i < x.length; i++ ) { - if ( randu() < 0.2 ) { - mask[ i ] = 1; - } else { - mask[ i ] = 0; - } - x[ i ] = round( (randu()*100.0) - 50.0 ); -} +var uniformOptions = { + 'dtype': 'float64' +}; +var bernoulliOptions = { + 'dtype': 'uint8' +}; + +var x = uniform( 10, -50.0, 50.0, uniformOptions ); +var mask = bernoulli( x.length, 0.2, bernoulliOptions ); console.log( x ); console.log( mask ); @@ -182,6 +168,125 @@ console.log( v ); + + +
+ +### Usage + +```c +#include "stdlib/stats/base/dmskmin.h" +``` + +#### stdlib_strided_dmskmin( N, \*X, strideX, \*Mask, strideMask ) + +Computes the minimum value of a double-precision floating-point strided array according to a mask. + +```c +#include + +const double x[] = { 1.0, -2.0, 2.0 }; +const uint8_t mask[] = { 0, 1, 0 }; + +double v = stdlib_strided_dmskmin( 3, x, 1, mask, 1 ); +// returns 1.0 +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **X**: `[in] double*` input array. +- **strideX**: `[in] CBLAS_INT` stride length for `X`. +- **Mask**: `[in] uint8_t*` mask array. If a `Mask` array element is `0`, the corresponding element in `X` is considered valid and included in computation. If a `Mask` array element is `1`, the corresponding element in `X` is considered invalid/missing and excluded from computation. +- **strideMask**: `[in] CBLAS_INT` stride length for `Mask`. + +```c +double stdlib_strided_dmskmin( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const uint8_t *Mask, const CBLAS_INT strideMask ); +``` + + + +#### stdlib_strided_dmskmin_ndarray( N, \*X, strideX, offsetX, \*Mask, strideMask, offsetMask ) + +Computes the minimum value of a double-precision floating-point strided array according to a mask and using alternative indexing semantics. + +```c +#include + +const double x[] = { 1.0, -2.0, 2.0 }; +const uint8_t mask[] = { 0, 1, 0 }; + +double v = stdlib_strided_dmskmin( 3, x, 1, 0, mask, 1, 0 ); +// returns 1.0 +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **X**: `[in] double*` input array. +- **strideX**: `[in] CBLAS_INT` stride length for `X`. +- **offsetX**: `[in] CBLAS_INT` starting index for `X`. +- **Mask**: `[in] uint8_t*` mask array. If a `Mask` array element is `0`, the corresponding element in `X` is considered valid and included in computation. If a `Mask` array element is `1`, the corresponding element in `X` is considered invalid/missing and excluded from computation. +- **strideMask**: `[in] CBLAS_INT` stride length for `Mask`. +- **offsetMask**: `[in] CBLAS_INT` starting index for `Mask`. + +```c +double stdlib_strided_dmskmin_ndarray( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const uint8_t *Mask, const CBLAS_INT strideMask, const CBLAS_INT offsetMask ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/stats/base/dmskmin.h" +#include +#include + +int main( void ) { + // Create a strided array: + const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 }; + + // Create a mask array: + const uint8_t mask[] = { 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 }; + + // Specify the number of elements: + const int N = 5; + + // Specify the stride lengths: + const int strideX = 2; + const int strideMask = 2; + + // Compute the minimum value: + double v = stdlib_strided_dmskmin( N, x, strideX, mask, strideMask ); + + // Print the result: + printf( "min: %lf\n", v ); +} +``` + +
+ + + + + + +
diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.js b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.js index 50b44e6d5266..3dd78bbe6c98 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.js @@ -21,15 +21,24 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var bernoulli = require( '@stdlib/random/array/bernoulli' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); var pkg = require( './../package.json' ).name; var dmskmin = require( './../lib/dmskmin.js' ); +// VARIABLES // + +var uniformOptions = { + 'dtype': 'float64' +}; +var bernoulliOptions = { + 'dtype': 'uint8' +}; + + // FUNCTIONS // /** @@ -40,20 +49,8 @@ var dmskmin = require( './../lib/dmskmin.js' ); * @returns {Function} benchmark function */ function createBenchmark( len ) { - var mask; - var x; - var i; - - x = new Float64Array( len ); - mask = new Uint8Array( len ); - for ( i = 0; i < x.length; i++ ) { - if ( randu() < 0.2 ) { - mask[ i ] = 1; - } else { - mask[ i ] = 0; - } - x[ i ] = ( randu()*20.0 ) - 10.0; - } + var mask = bernoulli( len, 0.2, bernoulliOptions ); + var x = uniform( len, -10.0, 10.0, uniformOptions ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.native.js index b757c6df95d4..b5774740b75e 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.native.js @@ -22,11 +22,10 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var bernoulli = require( '@stdlib/random/array/bernoulli' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); var tryRequire = require( '@stdlib/utils/try-require' ); var pkg = require( './../package.json' ).name; @@ -37,6 +36,12 @@ var dmskmin = tryRequire( resolve( __dirname, './../lib/dmskmin.native.js' ) ); var opts = { 'skip': ( dmskmin instanceof Error ) }; +var uniformOptions = { + 'dtype': 'float64' +}; +var bernoulliOptions = { + 'dtype': 'uint8' +}; // FUNCTIONS // @@ -49,20 +54,8 @@ var opts = { * @returns {Function} benchmark function */ function createBenchmark( len ) { - var mask; - var x; - var i; - - x = new Float64Array( len ); - mask = new Uint8Array( len ); - for ( i = 0; i < x.length; i++ ) { - if ( randu() < 0.2 ) { - mask[ i ] = 1; - } else { - mask[ i ] = 0; - } - x[ i ] = ( randu()*20.0 ) - 10.0; - } + var mask = bernoulli( len, 0.2, bernoulliOptions ); + var x = uniform( len, -10.0, 10.0, uniformOptions ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.js index 902a67af63eb..da95dc91adf4 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.js @@ -21,15 +21,24 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var bernoulli = require( '@stdlib/random/array/bernoulli' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); var pkg = require( './../package.json' ).name; var dmskmin = require( './../lib/ndarray.js' ); +// VARIABLES // + +var uniformOptions = { + 'dtype': 'float64' +}; +var bernoulliOptions = { + 'dtype': 'uint8' +}; + + // FUNCTIONS // /** @@ -40,20 +49,8 @@ var dmskmin = require( './../lib/ndarray.js' ); * @returns {Function} benchmark function */ function createBenchmark( len ) { - var mask; - var x; - var i; - - x = new Float64Array( len ); - mask = new Uint8Array( len ); - for ( i = 0; i < x.length; i++ ) { - if ( randu() < 0.2 ) { - mask[ i ] = 1; - } else { - mask[ i ] = 0; - } - x[ i ] = ( randu()*20.0 ) - 10.0; - } + var mask = bernoulli( len, 0.2, bernoulliOptions ); + var x = uniform( len, -10.0, 10.0, uniformOptions ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.native.js index 6d8897c6dd31..6e7b504fd865 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.native.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/benchmark.ndarray.native.js @@ -22,11 +22,10 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); -var randu = require( '@stdlib/random/base/randu' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var bernoulli = require( '@stdlib/random/array/bernoulli' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); var tryRequire = require( '@stdlib/utils/try-require' ); var pkg = require( './../package.json' ).name; @@ -37,6 +36,12 @@ var dmskmin = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) ); var opts = { 'skip': ( dmskmin instanceof Error ) }; +var uniformOptions = { + 'dtype': 'float64' +}; +var bernoulliOptions = { + 'dtype': 'uint8' +}; // FUNCTIONS // @@ -49,20 +54,8 @@ var opts = { * @returns {Function} benchmark function */ function createBenchmark( len ) { - var mask; - var x; - var i; - - x = new Float64Array( len ); - mask = new Uint8Array( len ); - for ( i = 0; i < x.length; i++ ) { - if ( randu() < 0.2 ) { - mask[ i ] = 1; - } else { - mask[ i ] = 0; - } - x[ i ] = ( randu()*20.0 ) - 10.0; - } + var mask = bernoulli( len, 0.2, bernoulliOptions ); + var x = uniform( len, -10.0, 10.0, uniformOptions ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/c/benchmark.length.c index 81bbf03e14a5..eb44fa95f589 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/benchmark/c/benchmark.length.c @@ -94,7 +94,7 @@ static double rand_double( void ) { * @param len array length * @return elapsed time in seconds */ -static double benchmark( int iterations, int len ) { +static double benchmark1( int iterations, int len ) { unsigned char mask[ len ]; double elapsed; double x[ len ]; @@ -113,6 +113,7 @@ static double benchmark( int iterations, int len ) { v = 0.0; t = tic(); for ( i = 0; i < iterations; i++ ) { + // cppcheck-suppress uninitvar v = stdlib_strided_dmskmin( len, x, 1, mask, 1 ); if ( v != v ) { printf( "should not return NaN\n" ); @@ -126,6 +127,46 @@ static double benchmark( int iterations, int len ) { 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 ) { + unsigned char mask[ len ]; + double elapsed; + double x[ len ]; + double v; + double t; + int i; + + for ( i = 0; i < len; i++ ) { + if ( rand_double() < 0.2 ) { + mask[ i ] = 1; // missing + } else { + mask[ i ] = 0; + } + x[ i ] = ( rand_double() * 20000.0 ) - 10000.0; + } + v = 0.0; + t = tic(); + for ( i = 0; i < iterations; i++ ) { + // cppcheck-suppress uninitvar + v = stdlib_strided_dmskmin_ndarray( len, x, 1, 0, mask, 1, 0 ); + if ( v != v ) { + printf( "should not return NaN\n" ); + break; + } + } + elapsed = tic() - t; + if ( v != v ) { + printf( "should not return NaN\n" ); + } + return elapsed; +} + /** * Main execution sequence. */ @@ -148,7 +189,18 @@ int main( void ) { for ( j = 0; j < REPEATS; j++ ) { count += 1; printf( "# c::%s:len=%d\n", NAME, len ); - elapsed = benchmark( iter, 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 ); } diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/docs/repl.txt b/lib/node_modules/@stdlib/stats/base/dmskmin/docs/repl.txt index 8205792a6759..573087eb0ff6 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/docs/repl.txt +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/docs/repl.txt @@ -3,8 +3,8 @@ Computes the minimum value of a double-precision floating-point strided array according to a mask. - The `N` and `stride` parameters determine which elements are accessed at - runtime. + 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 offsets, use a typed array views. @@ -26,13 +26,13 @@ Input array. strideX: integer - Index increment for `x`. + Stride Length for `x`. mask: Uint8Array Mask array. strideMask: integer - Index increment for `mask`. + Stride Length for `mask`. Returns ------- @@ -47,11 +47,10 @@ > {{alias}}( x.length, x, 1, mask, 1 ) -2.0 - // Using `N` and `stride` parameters: + // Using `N` and stride parameters: > x = new {{alias:@stdlib/array/float64}}( [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0, -4.0 ] ); > mask = new {{alias:@stdlib/array/uint8}}( [ 0, 0, 0, 0, 0, 0, 1 ] ); - > var N = {{alias:@stdlib/math/base/special/floor}}( x.length / 2 ); - > {{alias}}( N, x, 2, mask, 2 ) + > {{alias}}( 3, x, 2, mask, 2 ) -2.0 // Using view offsets: @@ -59,10 +58,10 @@ > var x1 = new {{alias:@stdlib/array/float64}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); > var mask0 = new {{alias:@stdlib/array/uint8}}( [ 0, 0, 0, 0, 0, 0, 1 ] ); > var mask1 = new {{alias:@stdlib/array/uint8}}( mask0.buffer, mask0.BYTES_PER_ELEMENT*1 ); - > N = {{alias:@stdlib/math/base/special/floor}}( x0.length / 2 ); - > {{alias}}( N, x1, 2, mask1, 2 ) + > {{alias}}( 3, x1, 2, mask1, 2 ) -2.0 + {{alias}}.ndarray( N, x, strideX, offsetX, mask, strideMask, offsetMask ) Computes the minimum value of a double-precision floating-point strided array according to a mask and using alternative indexing semantics. @@ -80,7 +79,7 @@ Input array. strideX: integer - Index increment for `x`. + Stride Length for `x`. offsetX: integer Starting index for `x`. @@ -89,7 +88,7 @@ Mask array. strideMask: integer - Index increment for `mask`. + Stride Length for `mask`. offsetMask: integer Starting index for `mask`. @@ -110,8 +109,7 @@ // Using offset parameter: > x = new {{alias:@stdlib/array/float64}}( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0, -4.0 ] ); > mask = new {{alias:@stdlib/array/uint8}}( [ 0, 0, 0, 0, 0, 0, 1 ] ); - > var N = {{alias:@stdlib/math/base/special/floor}}( x.length / 2 ); - > {{alias}}.ndarray( N, x, 2, 1, mask, 2, 1 ) + > {{alias}}.ndarray( 3, x, 2, 1, mask, 2, 1 ) -2.0 See Also diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/examples/c/example.c b/lib/node_modules/@stdlib/stats/base/dmskmin/examples/c/example.c index 5d83cdd2c079..9fa5cafa30f6 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/examples/c/example.c +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/examples/c/example.c @@ -22,17 +22,17 @@ int main( void ) { // Create a strided array: - double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 }; + const double x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0 }; // Create a mask array: - uint8_t mask[] = { 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 }; + const uint8_t mask[] = { 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 }; // Specify the number of elements: - int64_t N = 5; + const int N = 5; // Specify the stride lengths: - int64_t strideX = 2; - int64_t strideMask = 2; + const int strideX = 2; + const int strideMask = 2; // Compute the minimum value: double v = stdlib_strided_dmskmin( N, x, strideX, mask, strideMask ); diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/examples/index.js b/lib/node_modules/@stdlib/stats/base/dmskmin/examples/index.js index ccdc244a1c2b..cd17ef2d8fd1 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/examples/index.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/examples/index.js @@ -18,26 +18,19 @@ 'use strict'; -var randu = require( '@stdlib/random/base/randu' ); -var round = require( '@stdlib/math/base/special/round' ); -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); +var uniform = require( '@stdlib/random/array/uniform' ); +var bernoulli = require( '@stdlib/random/array/bernoulli' ); var dmskmin = require( './../lib' ); -var mask; -var x; -var i; +var uniformOptions = { + 'dtype': 'float64' +}; +var bernoulliOptions = { + 'dtype': 'uint8' +}; -x = new Float64Array( 10 ); -mask = new Uint8Array( x.length ); -for ( i = 0; i < x.length; i++ ) { - if ( randu() < 0.2 ) { - mask[ i ] = 1; - } else { - mask[ i ] = 0; - } - x[ i ] = round( (randu()*100.0) - 50.0 ); -} +var x = uniform( 10, -50.0, 50.0, uniformOptions ); +var mask = bernoulli( x.length, 0.2, bernoulliOptions ); console.log( x ); console.log( mask ); diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/include.gypi b/lib/node_modules/@stdlib/stats/base/dmskmin/include.gypi index 868c5c12e852..26476a8c2655 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/include.gypi +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/include.gypi @@ -36,7 +36,7 @@ # Source files: 'src_files': [ - '<(src_dir)/addon.cpp', + '<(src_dir)/addon.c', ' /* @@ -31,7 +32,12 @@ extern "C" { /** * Computes the minimum value of a double-precision floating-point strided array according to a mask. */ -double stdlib_strided_dmskmin( const int64_t N, const double *X, const int64_t strideX, const uint8_t *Mask, const int64_t strideMask ); +double API_SUFFIX(stdlib_strided_dmskmin)( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const uint8_t *Mask, const CBLAS_INT strideMask ); + +/** +* Computes the minimum value of a double-precision floating-point strided array according to a mask and using alternative indexing semantics. +*/ +double API_SUFFIX(stdlib_strided_dmskmin_ndarray)( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const uint8_t *Mask, const CBLAS_INT strideMask, const CBLAS_INT offsetMask ); #ifdef __cplusplus } diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/dmskmin.js b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/dmskmin.js index 86d2e78c4626..c3d59d082840 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/dmskmin.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/dmskmin.js @@ -20,8 +20,8 @@ // MODULES // -var isnan = require( '@stdlib/math/base/assert/is-nan' ); -var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); +var stride2offset = require( '@stdlib/strided/base/stride2offset' ); +var ndarray = require( './ndarray.js' ); // MAIN // @@ -47,55 +47,9 @@ var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); * // returns -2.0 */ function dmskmin( N, x, strideX, mask, strideMask ) { - var min; - var ix; - var im; - var v; - var i; - - if ( N <= 0 ) { - return NaN; - } - if ( strideX < 0 ) { - ix = (1-N) * strideX; - } else { - ix = 0; - } - if ( strideMask < 0 ) { - im = (1-N) * strideMask; - } else { - im = 0; - } - for ( i = 0; i < N; i++ ) { - if ( mask[ im ] === 0 ) { - break; - } - ix += strideX; - im += strideMask; - } - if ( i === N ) { - return NaN; - } - min = x[ ix ]; - if ( isnan( min ) ) { - return min; - } - i += 1; - for ( i; i < N; i++ ) { - ix += strideX; - im += strideMask; - if ( mask[ im ] ) { - continue; - } - v = x[ ix ]; - if ( isnan( v ) ) { - return v; - } - if ( v < min || ( v === min && isNegativeZero( v ) ) ) { - min = v; - } - } - return min; + var ox = stride2offset( N, strideX ); + var om = stride2offset( N, strideMask ); + return ndarray( N, x, strideX, ox, mask, strideMask, om ); } diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/index.js b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/index.js index 12a4a30d4a40..8bfd1833449b 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/index.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/index.js @@ -37,14 +37,12 @@ * @example * var Float64Array = require( '@stdlib/array/float64' ); * var Uint8Array = require( '@stdlib/array/uint8' ); -* var floor = require( '@stdlib/math/base/special/floor' ); * var dmskmin = require( '@stdlib/stats/base/dmskmin' ); * * var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0, -5.0, -6.0 ] ); * var mask = new Uint8Array( [ 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 ] ); -* var N = floor( x.length / 2 ); * -* var v = dmskmin.ndarray( N, x, 2, 1, mask, 2, 1 ); +* var v = dmskmin.ndarray( 5, x, 2, 1, mask, 2, 1 ); * // returns -2.0 */ diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.js b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.js index 6bcbe3bbe996..20258d249ac7 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.js @@ -45,9 +45,8 @@ var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); * * var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0, -5.0, -6.0 ] ); * var mask = new Uint8Array( [ 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 ] ); -* var N = floor( x.length / 2 ); * -* var v = dmskmin( N, x, 2, 1, mask, 2, 1 ); +* var v = dmskmin( 5, x, 2, 1, mask, 2, 1 ); * // returns -2.0 */ function dmskmin( N, x, strideX, offsetX, mask, strideMask, offsetMask ) { diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.native.js b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.native.js index 3a19aef5489e..46405c90df59 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/lib/ndarray.native.js @@ -20,9 +20,7 @@ // MODULES // -var Float64Array = require( '@stdlib/array/float64' ); -var Uint8Array = require( '@stdlib/array/uint8' ); -var addon = require( './dmskmin.native.js' ); +var addon = require( './../src/addon.node' ); // MAIN // @@ -42,27 +40,15 @@ var addon = require( './dmskmin.native.js' ); * @example * var Float64Array = require( '@stdlib/array/float64' ); * var Uint8Array = require( '@stdlib/array/uint8' ); -* var floor = require( '@stdlib/math/base/special/floor' ); * * var x = new Float64Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0, -5.0, -6.0 ] ); * var mask = new Uint8Array( [ 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 ] ); -* var N = floor( x.length / 2 ); * -* var v = dmskmin( N, x, 2, 1, mask, 2, 1 ); +* var v = dmskmin( 5, x, 2, 1, mask, 2, 1 ); * // returns -2.0 */ function dmskmin( N, x, strideX, offsetX, mask, strideMask, offsetMask ) { - var viewX; - var viewM; - if ( strideX < 0 ) { - offsetX += (N-1) * strideX; - } - if ( strideMask < 0 ) { - offsetMask += (N-1) * strideMask; - } - viewX = new Float64Array( x.buffer, x.byteOffset+(x.BYTES_PER_ELEMENT*offsetX), x.length-offsetX ); // eslint-disable-line max-len - viewM = new Uint8Array( mask.buffer, mask.byteOffset+(mask.BYTES_PER_ELEMENT*offsetMask), mask.length-offsetMask ); // eslint-disable-line max-len - return addon( N, viewX, strideX, viewM, strideMask ); + return addon.ndarray( N, x, strideX, offsetX, mask, strideMask, offsetMask ); // eslint-disable-line max-len } diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/manifest.json b/lib/node_modules/@stdlib/stats/base/dmskmin/manifest.json index 6481051a91a9..3f0b4cbbc446 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/manifest.json +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/manifest.json @@ -1,5 +1,8 @@ { - "options": {}, + "options": { + "task": "build", + "wasm": false + }, "fields": [ { "field": "src", @@ -24,17 +27,79 @@ ], "confs": [ { + "task": "build", + "wasm": false, "src": [ - "./src/dmskmin.c" + "./src/main.c" ], "include": [ "./include" ], - "libraries": [ - "-lm" + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-negative-zero", + "@stdlib/math/base/assert/is-nan", + "@stdlib/napi/export", + "@stdlib/napi/argv", + "@stdlib/napi/argv-int64", + "@stdlib/napi/argv-strided-float64array", + "@stdlib/napi/argv-strided-uint8array", + "@stdlib/napi/create-double" + ] + }, + { + "task": "benchmark", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/assert/is-negative-zero" + ] + }, + { + "task": "examples", + "wasm": false, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" + ], + "libraries": [], + "libpath": [], + "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/assert/is-nan", + "@stdlib/math/base/assert/is-negative-zero" + ] + }, + { + "task": "", + "wasm": true, + "src": [ + "./src/main.c" + ], + "include": [ + "./include" ], + "libraries": [], "libpath": [], "dependencies": [ + "@stdlib/blas/base/shared", + "@stdlib/strided/base/stride2offset", "@stdlib/math/base/assert/is-nan", "@stdlib/math/base/assert/is-negative-zero" ] diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/src/addon.c b/lib/node_modules/@stdlib/stats/base/dmskmin/src/addon.c new file mode 100644 index 000000000000..d56a9c13acd5 --- /dev/null +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/src/addon.c @@ -0,0 +1,66 @@ +/** +* @license Apache-2.0 +* +* 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. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +#include "stdlib/stats/base/dmskmin.h" +#include "stdlib/napi/export.h" +#include "stdlib/napi/argv.h" +#include "stdlib/napi/argv_int64.h" +#include "stdlib/napi/argv_strided_float64array.h" +#include "stdlib/napi/argv_strided_uint8array.h" +#include "stdlib/napi/create_double.h" +#include + +/** +* Receives JavaScript callback invocation data. +* +* @param env environment under which the function is invoked +* @param info callback data +* @return Node-API value +*/ +static napi_value addon( napi_env env, napi_callback_info info ) { + STDLIB_NAPI_ARGV( env, info, argv, argc, 5 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_INT64( env, strideMask, argv, 4 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 1 ); + STDLIB_NAPI_ARGV_STRIDED_UINT8ARRAY( env, Mask, N, strideMask, argv, 3 ); + STDLIB_NAPI_CREATE_DOUBLE( env, stdlib_strided_dmskmin( N, X, strideX, Mask, strideMask ), v ); + return v; +} + +/** +* Receives JavaScript callback invocation data. +* +* @param env environment under which the function is invoked +* @param info callback data +* @return Node-API value +*/ +static napi_value addon_method( napi_env env, napi_callback_info info ) { + STDLIB_NAPI_ARGV( env, info, argv, argc, 7 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 2 ); + STDLIB_NAPI_ARGV_INT64( env, strideMask, argv, 5 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, offsetMask, argv, 6 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 1 ); + STDLIB_NAPI_ARGV_STRIDED_UINT8ARRAY( env, Mask, N, strideMask, argv, 4 ); + STDLIB_NAPI_CREATE_DOUBLE( env, stdlib_strided_dmskmin_ndarray( N, X, strideX, offsetX, Mask, strideMask, offsetMask ), v ); + return v; +} + +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ) diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/src/addon.cpp b/lib/node_modules/@stdlib/stats/base/dmskmin/src/addon.cpp deleted file mode 100644 index 86acb4265ca2..000000000000 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/src/addon.cpp +++ /dev/null @@ -1,153 +0,0 @@ -/** -* @license Apache-2.0 -* -* Copyright (c) 2020 The Stdlib Authors. -* -* Licensed under the Apache License, Version 2.0 (the "License"); -* you may not use this file except in compliance with the License. -* You may obtain a copy of the License at -* -* http://www.apache.org/licenses/LICENSE-2.0 -* -* Unless required by applicable law or agreed to in writing, software -* distributed under the License is distributed on an "AS IS" BASIS, -* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. -* See the License for the specific language governing permissions and -* limitations under the License. -*/ - -#include "stdlib/stats/base/dmskmin.h" -#include -#include -#include -#include -#include - -/** -* Add-on namespace. -*/ -namespace stdlib_stats_base_dmskmin { - - /** - * Computes the minimum value of a double-precision floating-point strided array according to a mask. - * - * ## Notes - * - * - When called from JavaScript, the function expects five arguments: - * - * - `N`: number of indexed elements - * - `X`: input array - * - `strideX`: `X` stride length - * - `Mask`: mask array - * - `strideMask`: `Mask` stride length - */ - napi_value node_dmskmin( napi_env env, napi_callback_info info ) { - napi_status status; - - size_t argc = 5; - napi_value argv[ 5 ]; - status = napi_get_cb_info( env, info, &argc, argv, nullptr, nullptr ); - assert( status == napi_ok ); - - if ( argc < 5 ) { - napi_throw_error( env, nullptr, "invalid invocation. Must provide 5 arguments." ); - return nullptr; - } - - napi_valuetype vtype0; - status = napi_typeof( env, argv[ 0 ], &vtype0 ); - assert( status == napi_ok ); - if ( vtype0 != napi_number ) { - napi_throw_type_error( env, nullptr, "invalid argument. First argument must be a number." ); - return nullptr; - } - - bool res1; - status = napi_is_typedarray( env, argv[ 1 ], &res1 ); - assert( status == napi_ok ); - if ( res1 == false ) { - napi_throw_type_error( env, nullptr, "invalid argument. Second argument must be a Float64Array." ); - return nullptr; - } - - napi_valuetype vtype2; - status = napi_typeof( env, argv[ 2 ], &vtype2 ); - assert( status == napi_ok ); - if ( vtype2 != napi_number ) { - napi_throw_type_error( env, nullptr, "invalid argument. Third argument must be a number." ); - return nullptr; - } - - bool res3; - status = napi_is_typedarray( env, argv[ 3 ], &res3 ); - assert( status == napi_ok ); - if ( res3 == false ) { - napi_throw_type_error( env, nullptr, "invalid argument. Fourth argument must be a Uint8Array." ); - return nullptr; - } - - napi_valuetype vtype4; - status = napi_typeof( env, argv[ 4 ], &vtype4 ); - assert( status == napi_ok ); - if ( vtype4 != napi_number ) { - napi_throw_type_error( env, nullptr, "invalid argument. Fifth argument must be a number." ); - return nullptr; - } - - int64_t N; - status = napi_get_value_int64( env, argv[ 0 ], &N ); - assert( status == napi_ok ); - - int64_t strideX; - status = napi_get_value_int64( env, argv[ 2 ], &strideX ); - assert( status == napi_ok ); - - int64_t strideMask; - status = napi_get_value_int64( env, argv[ 4 ], &strideMask ); - assert( status == napi_ok ); - - napi_typedarray_type vtype1; - size_t xlen; - void *X; - status = napi_get_typedarray_info( env, argv[ 1 ], &vtype1, &xlen, &X, nullptr, nullptr ); - assert( status == napi_ok ); - if ( vtype1 != napi_float64_array ) { - napi_throw_type_error( env, nullptr, "invalid argument. Second argument must be a Float64Array." ); - return nullptr; - } - if ( (N-1)*llabs(strideX) >= (int64_t)xlen ) { - napi_throw_range_error( env, nullptr, "invalid argument. Second argument has insufficient elements based on the associated stride and the number of indexed elements." ); - return nullptr; - } - - napi_typedarray_type vtype3; - size_t mlen; - void *Mask; - status = napi_get_typedarray_info( env, argv[ 3 ], &vtype3, &mlen, &Mask, nullptr, nullptr ); - assert( status == napi_ok ); - if ( vtype3 != napi_uint8_array ) { - napi_throw_type_error( env, nullptr, "invalid argument. Fourth argument must be a Uint8Array." ); - return nullptr; - } - if ( (N-1)*llabs(strideMask) >= (int64_t)mlen ) { - napi_throw_range_error( env, nullptr, "invalid argument. Fourth argument has insufficient elements based on the associated stride and the number of indexed elements." ); - return nullptr; - } - - napi_value v; - status = napi_create_double( env, stdlib_strided_dmskmin( N, (double *)X, strideX, (uint8_t *)Mask, strideMask ), &v ); - assert( status == napi_ok ); - - return v; - } - - napi_value Init( napi_env env, napi_value exports ) { - napi_status status; - napi_value fcn; - status = napi_create_function( env, "exports", NAPI_AUTO_LENGTH, node_dmskmin, NULL, &fcn ); - assert( status == napi_ok ); - return fcn; - } - - NAPI_MODULE( NODE_GYP_MODULE_NAME, Init ) -} // end namespace stdlib_stats_base_dmskmin diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/src/dmskmin.c b/lib/node_modules/@stdlib/stats/base/dmskmin/src/main.c similarity index 56% rename from lib/node_modules/@stdlib/stats/base/dmskmin/src/dmskmin.c rename to lib/node_modules/@stdlib/stats/base/dmskmin/src/main.c index dab8ccf3c9eb..5d1dd57d6b8d 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/src/dmskmin.c +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/src/main.c @@ -1,7 +1,7 @@ /** * @license Apache-2.0 * -* Copyright (c) 2020 The Stdlib Authors. +* Copyright (c) 2024 The Stdlib Authors. * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. @@ -19,6 +19,8 @@ #include "stdlib/stats/base/dmskmin.h" #include "stdlib/math/base/assert/is_nan.h" #include "stdlib/math/base/assert/is_negative_zero.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" #include /** @@ -31,26 +33,36 @@ * @param strideMask Mask stride length * @return output value */ -double stdlib_strided_dmskmin( const int64_t N, const double *X, const int64_t strideX, const uint8_t *Mask, const int64_t strideMask ) { +double API_SUFFIX(stdlib_strided_dmskmin)( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const uint8_t *Mask, const CBLAS_INT strideMask ) { + const CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); + const CBLAS_INT om = stdlib_strided_stride2offset( N, strideMask ); + return API_SUFFIX(stdlib_strided_dmskmin_ndarray)( N, X, strideX, ox, Mask, strideMask, om ); +} + +/** +* Computes the minimum value of a double-precision floating-point strided array according to a mask and using alternative indexing semantics. +* +* @param N number of indexed elements +* @param X input array +* @param strideX X stride length +* @param offsetX starting index for X +* @param Mask mask array +* @param strideMask Mask stride length +* @param offsetMask starting index for Mask +* @return output value +*/ +double API_SUFFIX(stdlib_strided_dmskmin_ndarray)( const CBLAS_INT N, const double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX, const uint8_t *Mask, const CBLAS_INT strideMask, const CBLAS_INT offsetMask ) { + CBLAS_INT ix; + CBLAS_INT im; + CBLAS_INT i; double min; - int64_t ix; - int64_t im; - int64_t i; double v; if ( N <= 0 ) { return 0.0 / 0.0; // NaN } - if ( strideX < 0 ) { - ix = (1-N) * strideX; - } else { - ix = 0; - } - if ( strideMask < 0 ) { - im = (1-N) * strideMask; - } else { - im = 0; - } + ix = offsetX; + im = offsetMask; for ( i = 0; i < N; i++ ) { if ( Mask[ im ] == 0 ) { break; diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.js b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.js index 0ef22637413f..5fbaf3ed47b2 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.js @@ -21,7 +21,6 @@ // MODULES // var tape = require( 'tape' ); -var floor = require( '@stdlib/math/base/special/floor' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); var Float64Array = require( '@stdlib/array/float64' ); @@ -133,7 +132,6 @@ tape( 'if provided an `N` parameter equal to `1`, the function returns the first tape( 'the function supports `stride` parameters', function test( t ) { var mask; - var N; var x; var v; @@ -162,8 +160,7 @@ tape( 'the function supports `stride` parameters', function test( t ) { 1 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, 2, mask, 2 ); + v = dmskmin( 5, x, 2, mask, 2 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -171,7 +168,6 @@ tape( 'the function supports `stride` parameters', function test( t ) { tape( 'the function supports negative `stride` parameters', function test( t ) { var mask; - var N; var x; var v; @@ -200,8 +196,7 @@ tape( 'the function supports negative `stride` parameters', function test( t ) { 0 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, -2, mask, -2 ); + v = dmskmin( 5, x, -2, mask, -2 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -212,7 +207,6 @@ tape( 'the function supports view offsets', function test( t ) { var mask1; var x0; var x1; - var N; var v; x0 = new Float64Array([ @@ -244,9 +238,8 @@ tape( 'the function supports view offsets', function test( t ) { x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element mask1 = new Uint8Array( mask0.buffer, mask0.BYTES_PER_ELEMENT*1 ); // start at 2nd element - N = floor( x1.length/2 ); - v = dmskmin( N, x1, 2, mask1, 2 ); + v = dmskmin( 5, x1, 2, mask1, 2 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.native.js b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.native.js index dfb405f17387..a31ba1a8d3b8 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.native.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.dmskmin.native.js @@ -22,7 +22,6 @@ var resolve = require( 'path' ).resolve; var tape = require( 'tape' ); -var floor = require( '@stdlib/math/base/special/floor' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); var Float64Array = require( '@stdlib/array/float64' ); @@ -142,7 +141,6 @@ tape( 'if provided an `N` parameter equal to `1`, the function returns the first tape( 'the function supports `stride` parameters', opts, function test( t ) { var mask; - var N; var x; var v; @@ -171,8 +169,7 @@ tape( 'the function supports `stride` parameters', opts, function test( t ) { 1 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, 2, mask, 2 ); + v = dmskmin( 5, x, 2, mask, 2 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -180,7 +177,6 @@ tape( 'the function supports `stride` parameters', opts, function test( t ) { tape( 'the function supports negative `stride` parameters', opts, function test( t ) { var mask; - var N; var x; var v; @@ -209,8 +205,7 @@ tape( 'the function supports negative `stride` parameters', opts, function test( 0 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, -2, mask, -2 ); + v = dmskmin( 5, x, -2, mask, -2 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -221,7 +216,6 @@ tape( 'the function supports view offsets', opts, function test( t ) { var mask1; var x0; var x1; - var N; var v; x0 = new Float64Array([ @@ -253,9 +247,8 @@ tape( 'the function supports view offsets', opts, function test( t ) { x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element mask1 = new Uint8Array( mask0.buffer, mask0.BYTES_PER_ELEMENT*1 ); // start at 2nd element - N = floor( x1.length/2 ); - v = dmskmin( N, x1, 2, mask1, 2 ); + v = dmskmin( 5, x1, 2, mask1, 2 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.js b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.js index 2756d56e6da2..cb0c228316bf 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.js @@ -21,7 +21,6 @@ // MODULES // var tape = require( 'tape' ); -var floor = require( '@stdlib/math/base/special/floor' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); var Float64Array = require( '@stdlib/array/float64' ); @@ -133,7 +132,6 @@ tape( 'if provided an `N` parameter equal to `1`, the function returns the first tape( 'the function supports `stride` parameters', function test( t ) { var mask; - var N; var x; var v; @@ -162,8 +160,7 @@ tape( 'the function supports `stride` parameters', function test( t ) { 1 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, 2, 0, mask, 2, 0 ); + v = dmskmin( 5, x, 2, 0, mask, 2, 0 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -171,7 +168,6 @@ tape( 'the function supports `stride` parameters', function test( t ) { tape( 'the function supports negative `stride` parameters', function test( t ) { var mask; - var N; var x; var v; @@ -200,8 +196,7 @@ tape( 'the function supports negative `stride` parameters', function test( t ) { 0 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, -2, 8, mask, -2, 8 ); + v = dmskmin( 5, x, -2, 8, mask, -2, 8 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -209,7 +204,6 @@ tape( 'the function supports negative `stride` parameters', function test( t ) { tape( 'the function supports `offset` parameters', function test( t ) { var mask; - var N; var x; var v; @@ -237,9 +231,8 @@ tape( 'the function supports `offset` parameters', function test( t ) { 1, 1 // 4 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, 2, 1, mask, 2, 1 ); + v = dmskmin( 5, x, 2, 1, mask, 2, 1 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); diff --git a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.native.js b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.native.js index 2c8edd097d19..ab53ffece127 100644 --- a/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/stats/base/dmskmin/test/test.ndarray.native.js @@ -22,7 +22,6 @@ var resolve = require( 'path' ).resolve; var tape = require( 'tape' ); -var floor = require( '@stdlib/math/base/special/floor' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var isNegativeZero = require( '@stdlib/math/base/assert/is-negative-zero' ); var Float64Array = require( '@stdlib/array/float64' ); @@ -120,7 +119,7 @@ tape( 'if provided an `N` parameter less than or equal to `0`, the function retu v = dmskmin( 0, x, 1, 0, mask, 1, 0 ); t.strictEqual( isnan( v ), true, 'returns expected value' ); - v = dmskmin( -1, x, 1, 0, mask, 1 ); + v = dmskmin( -1, x, 1, 0, mask, 1, 0 ); t.strictEqual( isnan( v ), true, 'returns expected value' ); t.end(); @@ -142,7 +141,6 @@ tape( 'if provided an `N` parameter equal to `1`, the function returns the first tape( 'the function supports `stride` parameters', opts, function test( t ) { var mask; - var N; var x; var v; @@ -171,8 +169,7 @@ tape( 'the function supports `stride` parameters', opts, function test( t ) { 1 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, 2, 0, mask, 2, 0 ); + v = dmskmin( 5, x, 2, 0, mask, 2, 0 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -180,7 +177,6 @@ tape( 'the function supports `stride` parameters', opts, function test( t ) { tape( 'the function supports negative `stride` parameters', opts, function test( t ) { var mask; - var N; var x; var v; @@ -209,8 +205,7 @@ tape( 'the function supports negative `stride` parameters', opts, function test( 0 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, -2, 8, mask, -2, 8 ); + v = dmskmin( 5, x, -2, 8, mask, -2, 8 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end(); @@ -218,7 +213,6 @@ tape( 'the function supports negative `stride` parameters', opts, function test( tape( 'the function supports `offset` parameters', opts, function test( t ) { var mask; - var N; var x; var v; @@ -246,9 +240,8 @@ tape( 'the function supports `offset` parameters', opts, function test( t ) { 1, 1 // 4 ]); - N = floor( x.length / 2 ); - v = dmskmin( N, x, 2, 1, mask, 2, 1 ); + v = dmskmin( 5, x, 2, 1, mask, 2, 1 ); t.strictEqual( v, -2.0, 'returns expected value' ); t.end();