diff --git a/lib/node_modules/@stdlib/array/base/docs/types/index.d.ts b/lib/node_modules/@stdlib/array/base/docs/types/index.d.ts index 524e624f03b9..512946c21756 100644 --- a/lib/node_modules/@stdlib/array/base/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/array/base/docs/types/index.d.ts @@ -3365,10 +3365,10 @@ interface Namespace { * * @example * var x = [ 1, 2, 3, 4 ]; - * var idx = [ 5, 6, 7, 8 ]; + * var y = [ 5, 6, 7, 8 ]; * var mask = [ 0, 1, 0, 1 ]; * - * var out = ns.mskfilter2( x, idx, mask ); + * var out = ns.mskfilter2( x, y, mask ); * // returns [ [ 2, 4 ], [ 6, 8 ] ] */ mskfilter2: typeof mskfilter2; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/README.md b/lib/node_modules/@stdlib/blas/ext/base/dapx/README.md index f0a97cb8c937..c6a3afbd39f4 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/README.md @@ -48,7 +48,7 @@ The function has the following parameters: - **N**: number of indexed elements. - **alpha**: scalar constant. - **x**: input [`Float64Array`][@stdlib/array/float64]. -- **strideX**: stride length for `x`. +- **strideX**: stride length. The `N` and stride parameters determine which elements in the strided array are accessed at runtime. For example, to add a constant to every other element: @@ -168,14 +168,14 @@ console.log( x ); #include "stdlib/blas/ext/base/dapx.h" ``` -#### c_dapx( N, alpha, \*X, strideX ) +#### stdlib_strided_dapx( N, alpha, \*X, strideX ) Adds a scalar constant to each element in a double-precision floating-point strided array. ```c double x[] = { 1.0, 2.0, 3.0, 4.0 }; -c_dapx( 4, 5.0, x, 1 ); +stdlib_strided_dapx( 4, 5.0, x, 1 ); ``` @@ -187,17 +187,17 @@ The function accepts the following arguments: - **strideX**: `[in] CBLAS_INT` stride length for `X`. ```c -void c_dapx( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX ); +void stdlib_strided_dapx( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX ); ``` -#### c_dapx_ndarray( N, alpha, \*X, strideX, offsetX ) +#### stdlib_strided_dapx_ndarray( N, alpha, \*X, strideX, offsetX ) Adds a scalar constant to each element in a double-precision floating-point strided array using alternative indexing semantics. ```c double x[] = { 1.0, 2.0, 3.0, 4.0 }; -c_dapx_ndarray( 4, 5.0, x, 1, 0 ); +stdlib_strided_dapx_ndarray( 4, 5.0, x, 1, 0 ); ``` The function accepts the following arguments: @@ -205,11 +205,11 @@ The function accepts the following arguments: - **N**: `[in] CBLAS_INT` number of indexed elements. - **alpha**: `[in] double` scalar constant. - **X**: `[inout] double*` input array. -- **strideX**: `[in] CBLAS_INT` stride length for `X`. -- **offsetX**: `[in] CBLAS_INT` starting index for `X`. +- **strideX**: `[in] CBLAS_INT` stride length. +- **offsetX**: `[in] CBLAS_INT` starting index. ```c -void c_dapx_ndarray( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); +void stdlib_strided_dapx_ndarray( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); ``` @@ -245,7 +245,7 @@ int main( void ) { const int strideX = 1; // Fill the array: - c_dapx( N, 5.0, x, strideX ); + stdlib_strided_dapx( N, 5.0, x, strideX ); // Print the result: for ( int i = 0; i < 8; i++ ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/dapx/benchmark/c/benchmark.length.c index a4e9dbd372b0..678f2c5df6c1 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/benchmark/c/benchmark.length.c @@ -106,7 +106,7 @@ static double benchmark1( int iterations, int len ) { t = tic(); for ( i = 0; i < iterations; i++ ) { // cppcheck-suppress uninitvar - c_dapx( len, 5.0, x, 1 ); + stdlib_strided_dapx( len, 5.0, x, 1 ); if ( x[ 0 ] != x[ 0 ] ) { printf( "should not return NaN\n" ); break; @@ -131,7 +131,7 @@ static double benchmark2( int iterations, int len ) { t = tic(); for ( i = 0; i < iterations; i++ ) { // cppcheck-suppress uninitvar - c_dapx_ndarray( len, 5.0, x, 1, 0 ); + stdlib_strided_dapx_ndarray( len, 5.0, x, 1, 0 ); if ( x[ 0 ] != x[ 0 ] ) { printf( "should not return NaN\n" ); break; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/dapx/examples/c/example.c index ec1ac38026ef..e77cbe0d09ba 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/examples/c/example.c @@ -30,7 +30,7 @@ int main( void ) { const int strideX = 1; // Add a constant to each element: - c_dapx( N, 5.0, x, strideX ); + stdlib_strided_dapx( N, 5.0, x, strideX ); // Print the result: for ( int i = 0; i < 8; i++ ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/include/stdlib/blas/ext/base/dapx.h b/lib/node_modules/@stdlib/blas/ext/base/dapx/include/stdlib/blas/ext/base/dapx.h index d292078b1ca3..7298dbee5a8f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/include/stdlib/blas/ext/base/dapx.h +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/include/stdlib/blas/ext/base/dapx.h @@ -31,12 +31,12 @@ extern "C" { /** * Adds a scalar constant to each element in a double-precision floating-point strided array. */ -void API_SUFFIX(c_dapx)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX ); +void API_SUFFIX(stdlib_strided_dapx)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX ); /** * Adds a scalar constant to each element in a double-precision floating-point strided array using alternative indexing semantics. */ -void API_SUFFIX(c_dapx_ndarray)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); +void API_SUFFIX(stdlib_strided_dapx_ndarray)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); #ifdef __cplusplus } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.js index efe6621ddb44..d13e30fd1128 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.js @@ -39,9 +39,8 @@ var M = 5; * var Float64Array = require( '@stdlib/array/float64' ); * * var x = new Float64Array( [ 1.0, -2.0, 3.0, -4.0, 5.0, -6.0 ] ); -* var alpha = 5.0; * -* dapx( 3, alpha, x, 1, x.length-3 ); +* dapx( 3, 5.0, x, 1, x.length-3 ); * // x => [ 1.0, -2.0, 3.0, 1.0, 10.0, -1.0 ] */ function dapx( N, alpha, x, strideX, offsetX ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.native.js index f743b6fac10a..c9291a1a3c20 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/lib/ndarray.native.js @@ -39,9 +39,8 @@ var addon = require( './../src/addon.node' ); * var Float64Array = require( '@stdlib/array/float64' ); * * var x = new Float64Array( [ 1.0, -2.0, 3.0, -4.0, 5.0, -6.0 ] ); -* var alpha = 5.0; * -* dapx( 3, alpha, x, 1, x.length-3 ); +* dapx( 3, 5.0, x, 1, x.length-3 ); * // x => [ 1.0, -2.0, 3.0, 1.0, 10.0, -1.0 ] */ function dapx( N, alpha, x, strideX, offsetX ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/dapx/src/addon.c index 455f74e295db..88e2c6bb251c 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/src/addon.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/src/addon.c @@ -38,7 +38,7 @@ static napi_value addon( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_DOUBLE( env, alpha, argv, 1 ); STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); - API_SUFFIX(c_dapx)( N, alpha, X, strideX ); + API_SUFFIX(stdlib_strided_dapx)( N, alpha, X, strideX ); return NULL; } @@ -56,7 +56,7 @@ static napi_value addon_method( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 ); STDLIB_NAPI_ARGV_STRIDED_FLOAT64ARRAY( env, X, N, strideX, argv, 2 ); - API_SUFFIX(c_dapx_ndarray)( N, alpha, X, strideX, offsetX ); + API_SUFFIX(stdlib_strided_dapx_ndarray)( N, alpha, X, strideX, offsetX ); return NULL; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapx/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/dapx/src/main.c index 193e2fec7590..dd9e1084271b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapx/src/main.c +++ b/lib/node_modules/@stdlib/blas/ext/base/dapx/src/main.c @@ -28,9 +28,9 @@ * @param X input array * @param strideX stride length */ -void API_SUFFIX(c_dapx)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX ) { +void API_SUFFIX(stdlib_strided_dapx)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX ) { CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); - API_SUFFIX(c_dapx_ndarray)( N, alpha, X, strideX, ox ); + API_SUFFIX(stdlib_strided_dapx_ndarray)( N, alpha, X, strideX, ox ); } /** @@ -42,7 +42,7 @@ void API_SUFFIX(c_dapx)( const CBLAS_INT N, const double alpha, double *X, const * @param strideX stride length * @param offsetX starting index */ -void API_SUFFIX(c_dapx_ndarray)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { +void API_SUFFIX(stdlib_strided_dapx_ndarray)( const CBLAS_INT N, const double alpha, double *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { CBLAS_INT ix; CBLAS_INT m; CBLAS_INT i; diff --git a/lib/node_modules/@stdlib/blas/ext/base/dapxsumkbn/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/dapxsumkbn/manifest.json index f587df37ca10..e5e1db92b9e7 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/dapxsumkbn/manifest.json +++ b/lib/node_modules/@stdlib/blas/ext/base/dapxsumkbn/manifest.json @@ -79,7 +79,8 @@ "libpath": [], "dependencies": [ "@stdlib/blas/base/shared", - "@stdlib/strided/base/stride2offset" + "@stdlib/strided/base/stride2offset", + "@stdlib/math/base/special/abs" ] }, { diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/README.md b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/README.md index 804f80b75920..0abb5a891097 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/README.md +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/README.md @@ -20,7 +20,7 @@ limitations under the License. # sapxsumors -> Add a constant to each single-precision floating-point strided array element and compute the sum using ordinary recursive summation. +> Add a scalar constant to each single-precision floating-point strided array element and compute the sum using ordinary recursive summation.
@@ -36,9 +36,9 @@ limitations under the License. var sapxsumors = require( '@stdlib/blas/ext/base/sapxsumors' ); ``` -#### sapxsumors( N, alpha, x, stride ) +#### sapxsumors( N, alpha, x, strideX ) -Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. ```javascript var Float32Array = require( '@stdlib/array/float32' ); @@ -52,10 +52,11 @@ var v = sapxsumors( x.length, 5.0, x, 1 ); The function has the following parameters: - **N**: number of indexed elements. +- **alpha**: scalar constant. - **x**: input [`Float32Array`][@stdlib/array/float32]. -- **stride**: index increment for `x`. +- **strideX**: stride length. -The `N` and `stride` parameters determine which elements in the strided array are accessed at runtime. For example, to access every other element in `x`, +The `N` and stride parameters determine which elements in the strided array are accessed at runtime. For example, to access every other element: ```javascript var Float32Array = require( '@stdlib/array/float32' ); @@ -80,9 +81,9 @@ var v = sapxsumors( 4, 5.0, x1, 2 ); // returns 25.0 ``` -#### sapxsumors.ndarray( N, alpha, x, stride, offset ) +#### sapxsumors.ndarray( N, alpha, x, strideX, offsetX ) -Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. +Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. ```javascript var Float32Array = require( '@stdlib/array/float32' ); @@ -95,9 +96,9 @@ var v = sapxsumors.ndarray( x.length, 5.0, x, 1, 0 ); The function has the following additional parameters: -- **offset**: starting index for `x`. +- **offsetX**: starting index. -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 access 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 parameter supports indexing semantics based on a starting index. For example, to access every other element starting from the second element: ```javascript var Float32Array = require( '@stdlib/array/float32' ); @@ -130,12 +131,13 @@ var v = sapxsumors.ndarray( 4, 5.0, x, 2, 1 ); ```javascript -var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ).factory; -var filledarrayBy = require( '@stdlib/array/filled-by' ); +var discreteUniform = require( '@stdlib/random/array/discrete-uniform' ); var Float32Array = require( '@stdlib/array/float32' ); var sapxsumors = require( '@stdlib/blas/ext/base/sapxsumors' ); -var x = filledarrayBy( 10, 'float32', discreteUniform( 0, 100 ) ); +var x = discreteUniform( 10, -100, 100, { + 'dtype': 'float32' +}); console.log( x ); var v = sapxsumors( x.length, 5.0, x, 1 ); @@ -146,6 +148,125 @@ console.log( v ); + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/blas/ext/base/sapxsumors.h" +``` + +#### stdlib_strided_sapxsumors( N, alpha, \*X, strideX ) + +Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. + +```c +const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f }; + +float v = stdlib_strided_sapxsumors( 4, 5.0f, x, 1 ); +// returns 30.0 +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **alpha**: `[in] float` scalar constant. +- **X**: `[in] float*` input array. +- **strideX**: `[in] CBLAS_INT` stride length. + +```c +float stdlib_strided_sapxsumors( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX ); +``` + +#### stdlib_strided_sapxsumors_ndarray( N, alpha, \*X, strideX, offsetX ) + +Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. + +```c +const float x[] = { 1.0f, 2.0f, 3.0f, 4.0f }; + +float v = stdlib_strided_sapxsumors_ndarray( 4, 5.0f, x, 1, 0 ); +// returns 30.0 +``` + +The function accepts the following arguments: + +- **N**: `[in] CBLAS_INT` number of indexed elements. +- **alpha**: `[in] float` scalar constant. +- **X**: `[in] float*` input array. +- **strideX**: `[in] CBLAS_INT` stride length. +- **offsetX**: `[in] CBLAS_INT` starting index. + +```c +float stdlib_strided_sapxsumors_ndarray( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/blas/ext/base/sapxsumors.h" +#include + +int main( void ) { + // Create a strided array: + const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f }; + + // Specify the number of indexed elements: + const int N = 8; + + // Specify a stride: + const int strideX = 1; + + // Compute the sum: + float v = stdlib_strided_sapxsumors( N, 5.0f, x, strideX ); + + // Print the result: + printf( "Sum: %f\n", v ); +} +``` + +
+ + + +
+ + +
diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.js index e2e536267536..98d520cbdc14 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.js @@ -21,8 +21,7 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var uniform = require( '@stdlib/random/base/uniform' ).factory; -var filledarrayBy = require( '@stdlib/array/filled-by' ); +var uniform = require( '@stdlib/random/array/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var pkg = require( './../package.json' ).name; @@ -31,7 +30,9 @@ var sapxsumors = require( './../lib/sapxsumors.js' ); // VARIABLES // -var rand = uniform( -100.0, 100.0 ); +var options = { + 'dtype': 'float32' +}; // FUNCTIONS // @@ -44,7 +45,7 @@ var rand = uniform( -100.0, 100.0 ); * @returns {Function} benchmark function */ function createBenchmark( len ) { - var x = filledarrayBy( len, 'float32', rand ); + var x = uniform( len, -100, 100, options ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.native.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.native.js index b234d22157b4..b5f2ca916a59 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.native.js @@ -22,8 +22,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); -var uniform = require( '@stdlib/random/base/uniform' ).factory; -var filledarrayBy = require( '@stdlib/array/filled-by' ); +var uniform = require( '@stdlib/random/array/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var tryRequire = require( '@stdlib/utils/try-require' ); @@ -36,7 +35,9 @@ var sapxsumors = tryRequire( resolve( __dirname, './../lib/sapxsumors.native.js' var opts = { 'skip': ( sapxsumors instanceof Error ) }; -var rand = uniform( -100.0, 100.0 ); +var options = { + 'dtype': 'float32' +}; // FUNCTIONS // @@ -49,7 +50,7 @@ var rand = uniform( -100.0, 100.0 ); * @returns {Function} benchmark function */ function createBenchmark( len ) { - var x = filledarrayBy( len, 'float32', rand ); + var x = uniform( len, -100, 100, options ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.js index ded8d1ca0f61..650bd0aaf506 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.js @@ -21,8 +21,7 @@ // MODULES // var bench = require( '@stdlib/bench' ); -var uniform = require( '@stdlib/random/base/uniform' ).factory; -var filledarrayBy = require( '@stdlib/array/filled-by' ); +var uniform = require( '@stdlib/random/array/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var pkg = require( './../package.json' ).name; @@ -31,7 +30,9 @@ var sapxsumors = require( './../lib/ndarray.js' ); // VARIABLES // -var rand = uniform( -100.0, 100.0 ); +var options = { + 'dtype': 'float32' +}; // FUNCTIONS // @@ -44,7 +45,7 @@ var rand = uniform( -100.0, 100.0 ); * @returns {Function} benchmark function */ function createBenchmark( len ) { - var x = filledarrayBy( len, 'float32', rand ); + var x = uniform( len, -100, 100, options ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.native.js index 6a005e1a62cd..f46632ba9e49 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/benchmark.ndarray.native.js @@ -22,8 +22,7 @@ var resolve = require( 'path' ).resolve; var bench = require( '@stdlib/bench' ); -var uniform = require( '@stdlib/random/base/uniform' ).factory; -var filledarrayBy = require( '@stdlib/array/filled-by' ); +var uniform = require( '@stdlib/random/array/uniform' ); var isnan = require( '@stdlib/math/base/assert/is-nan' ); var pow = require( '@stdlib/math/base/special/pow' ); var tryRequire = require( '@stdlib/utils/try-require' ); @@ -36,7 +35,9 @@ var sapxsumors = tryRequire( resolve( __dirname, './../lib/ndarray.native.js' ) var opts = { 'skip': ( sapxsumors instanceof Error ) }; -var rand = uniform( -100.0, 100.0 ); +var options = { + 'dtype': 'float32' +}; // FUNCTIONS // @@ -49,7 +50,7 @@ var rand = uniform( -100.0, 100.0 ); * @returns {Function} benchmark function */ function createBenchmark( len ) { - var x = filledarrayBy( len, 'float32', rand ); + var x = uniform( len, -100, 100, options ); return benchmark; function benchmark( b ) { diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/c/benchmark.length.c b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/c/benchmark.length.c index b0e31cd80473..4514ab8bdbdf 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/c/benchmark.length.c +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/benchmark/c/benchmark.length.c @@ -94,7 +94,7 @@ static float rand_float( void ) { * @param len array length * @return elapsed time in seconds */ -static double benchmark( int iterations, int len ) { +static double benchmark1( int iterations, int len ) { double elapsed; float x[ len ]; float v; @@ -107,6 +107,7 @@ static double benchmark( int iterations, int len ) { v = 0.0f; t = tic(); for ( i = 0; i < iterations; i++ ) { + // cppcheck-suppress uninitvar v = stdlib_strided_sapxsumors( len, 5.0f, x, 1 ); if ( v != v ) { printf( "should not return NaN\n" ); @@ -120,6 +121,40 @@ 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 ) { + double elapsed; + float x[ len ]; + float v; + double t; + int i; + + for ( i = 0; i < len; i++ ) { + x[ i ] = ( rand_float() * 20000.0f ) - 10000.0f; + } + v = 0.0f; + t = tic(); + for ( i = 0; i < iterations; i++ ) { + // cppcheck-suppress uninitvar + v = stdlib_strided_sapxsumors_ndarray( len, 5.0f, x, 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. */ @@ -142,7 +177,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/blas/ext/base/sapxsumors/docs/repl.txt b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/repl.txt index 1936cee0559d..85989a9bdfcf 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/repl.txt +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/repl.txt @@ -1,10 +1,10 @@ -{{alias}}( N, alpha, x, stride ) - Adds a constant to each single-precision floating-point strided array - element and computes the sum using ordinary recursive summation. +{{alias}}( N, alpha, x, strideX ) + Adds a scalar constant to each single-precision floating-point strided + array element and computes the sum using ordinary recursive summation. - The `N` and `stride` parameters determine which elements in the strided - array are accessed at runtime. + The `N` and stride parameters determine which elements in the strided array + are accessed at runtime. Indexing is relative to the first index. To introduce an offset, use a typed array view. @@ -17,13 +17,13 @@ Number of indexed elements. alpha: number - Constant. + Scalar constant. x: Float32Array Input array. - stride: integer - Index increment. + strideX: integer + Stride length. Returns ------- @@ -37,29 +37,25 @@ > {{alias}}( x.length, 5.0, x, 1 ) 16.0 - // Using `N` and `stride` parameters: + // Using `N` and stride parameters: > x = new {{alias:@stdlib/array/float32}}( [ -2.0, 1.0, 1.0, -5.0, 2.0, -1.0 ] ); - > var N = {{alias:@stdlib/math/base/special/floor}}( x.length / 2 ); - > var stride = 2; - > {{alias}}( N, 5.0, x, stride ) + > {{alias}}( 3, 5.0, x, 2 ) 16.0 // Using view offsets: > var x0 = new {{alias:@stdlib/array/float32}}( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0 ] ); > var x1 = new {{alias:@stdlib/array/float32}}( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); - > N = {{alias:@stdlib/math/base/special/floor}}( x0.length / 2 ); - > stride = 2; - > {{alias}}( N, 5.0, x1, stride ) + > {{alias}}( 3, 5.0, x1, 2 ) 14.0 -{{alias}}.ndarray( N, alpha, x, stride, offset ) - Adds a constant to each single-precision floating-point strided array - element and computes the sum using ordinary recursive summation and +{{alias}}.ndarray( N, alpha, x, strideX, offsetX ) + Adds a scalar constant to each single-precision floating-point strided + array element and computes the sum using ordinary recursive summation and alternative indexing semantics. While typed array views mandate a view offset based on the underlying - buffer, the `offset` parameter supports indexing semantics based on a + buffer, the offset parameter supports indexing semantics based on a starting index. Parameters @@ -68,15 +64,15 @@ Number of indexed elements. alpha: number - Constant. + Scalar constant. x: Float32Array Input array. - stride: integer - Index increment. + strideX: integer + Stride length. - offset: integer + offsetX: integer Starting index. Returns @@ -93,8 +89,7 @@ // Using offset parameter: > var x = new {{alias:@stdlib/array/float32}}( [ 1.0, -2.0, 3.0, 2.0, 5.0, -1.0 ] ); - > var N = {{alias:@stdlib/math/base/special/floor}}( x.length / 2 ); - > {{alias}}.ndarray( N, 5.0, x, 2, 1 ) + > {{alias}}.ndarray( 3, 5.0, x, 2, 1 ) 14.0 See Also diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/types/index.d.ts b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/types/index.d.ts index f321ad409cb5..18860b460cb9 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/docs/types/index.d.ts @@ -23,12 +23,12 @@ */ interface Routine { /** - * Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. + * Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. * * @param N - number of indexed elements - * @param alpha - constant + * @param alpha - scalar constant * @param x - input array - * @param stride - stride length + * @param strideX - stride length * @returns sum * * @example @@ -39,16 +39,16 @@ interface Routine { * var v = sapxsumors( x.length, 5.0, x, 1 ); * // returns 16.0 */ - ( N: number, alpha: number, x: Float32Array, stride: number ): number; + ( N: number, alpha: number, x: Float32Array, strideX: number ): number; /** - * Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. + * Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. * * @param N - number of indexed elements - * @param alpha - constant + * @param alpha - scalar constant * @param x - input array - * @param stride - stride length - * @param offset - starting index + * @param strideX - stride length + * @param offsetX - starting index * @returns sum * * @example @@ -59,16 +59,16 @@ interface Routine { * var v = sapxsumors.ndarray( x.length, 5.0, x, 1, 0 ); * // returns 16.0 */ - ndarray( N: number, alpha: number, x: Float32Array, stride: number, offset: number ): number; + ndarray( N: number, alpha: number, x: Float32Array, strideX: number, offsetX: number ): number; } /** -* Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. * * @param N - number of indexed elements -* @param alpha - constant +* @param alpha - scalar constant * @param x - input array -* @param stride - stride length +* @param strideX - stride length * @returns sum * * @example diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/c/example.c b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/c/example.c index 78cc54d0dd41..2723281a207f 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/c/example.c +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/c/example.c @@ -17,22 +17,21 @@ */ #include "stdlib/blas/ext/base/sapxsumors.h" -#include #include int main( void ) { // Create a strided array: - const float x[] = { 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0 }; + const float x[] = { 1.0f, -2.0f, 3.0f, -4.0f, 5.0f, -6.0f, 7.0f, -8.0f }; - // Specify the number of elements: - const int64_t N = 4; + // Specify the number of indexed elements: + const int N = 8; - // Specify the stride length: - const int64_t stride = 2; + // Specify a stride: + const int strideX = 1; // Compute the sum: - float v = stdlib_strided_sapxsumors( N, 5.0f, x, stride ); + float v = stdlib_strided_sapxsumors( N, 5.0f, x, strideX ); // Print the result: - printf( "sum: %f\n", v ); + printf( "Sum: %f\n", v ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/index.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/index.js index d79569cc0357..7d89fed3d6b5 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/examples/index.js @@ -18,11 +18,12 @@ 'use strict'; -var discreteUniform = require( '@stdlib/random/base/discrete-uniform' ).factory; -var filledarrayBy = require( '@stdlib/array/filled-by' ); +var discreteUniform = require( '@stdlib/random/array/discrete-uniform' ); var sapxsumors = require( './../lib' ); -var x = filledarrayBy( 10, 'float32', discreteUniform( 0, 100 ) ); +var x = discreteUniform( 10, -100, 100, { + 'dtype': 'float32' +}); console.log( x ); var v = sapxsumors( x.length, 5.0, x, 1 ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/include/stdlib/blas/ext/base/sapxsumors.h b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/include/stdlib/blas/ext/base/sapxsumors.h index 6203f1570b44..ed5f9464f17a 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/include/stdlib/blas/ext/base/sapxsumors.h +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/include/stdlib/blas/ext/base/sapxsumors.h @@ -19,7 +19,7 @@ #ifndef STDLIB_BLAS_EXT_BASE_SAPXSUMORS_H #define STDLIB_BLAS_EXT_BASE_SAPXSUMORS_H -#include +#include "stdlib/blas/base/shared.h" /* * If C++, prevent name mangling so that the compiler emits a binary file having undecorated names, thus mirroring the behavior of a C compiler. @@ -31,7 +31,12 @@ extern "C" { /** * Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. */ -float stdlib_strided_sapxsumors( const int64_t N, const float alpha, const float *X, const int64_t stride ); +float API_SUFFIX(stdlib_strided_sapxsumors)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX ); + +/** +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. +*/ +float API_SUFFIX(stdlib_strided_sapxsumors_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ); #ifdef __cplusplus } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/index.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/index.js index 8337dd181fda..0d73bea58abd 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/index.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/index.js @@ -19,7 +19,7 @@ 'use strict'; /** -* Add a constant to each single-precision floating-point strided array element and compute the sum using ordinary recursive summation. +* Add a scalar constant to each single-precision floating-point strided array element and compute the sum using ordinary recursive summation. * * @module @stdlib/blas/ext/base/sapxsumors * @@ -34,7 +34,6 @@ * * @example * var Float32Array = require( '@stdlib/array/float32' ); -* var floor = require( '@stdlib/math/base/special/floor' ); * var sapxsumors = require( '@stdlib/blas/ext/base/sapxsumors' ); * * var x = new Float32Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.js index 0dc1c6dbc875..21ad1f288eea 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.js @@ -26,25 +26,24 @@ var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' ); // MAIN // /** -* Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. * * @param {PositiveInteger} N - number of indexed elements -* @param {number} alpha - constant +* @param {number} alpha - scalar constant * @param {Float32Array} x - input array -* @param {integer} stride - stride length -* @param {NonNegativeInteger} offset - starting index +* @param {integer} strideX - stride length +* @param {NonNegativeInteger} offsetX - starting index * @returns {number} sum * * @example * var Float32Array = require( '@stdlib/array/float32' ); -* var floor = require( '@stdlib/math/base/special/floor' ); * * var x = new Float32Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] ); * * var v = sapxsumors( 4, 5.0, x, 2, 1 ); * // returns 25.0 */ -function sapxsumors( N, alpha, x, stride, offset ) { +function sapxsumors( N, alpha, x, strideX, offsetX ) { var sum; var ix; var i; @@ -52,14 +51,14 @@ function sapxsumors( N, alpha, x, stride, offset ) { if ( N <= 0 ) { return 0.0; } - if ( N === 1 || stride === 0 ) { - return float64ToFloat32( alpha + x[ 0 ] ); + ix = offsetX; + if ( strideX === 0 ) { + return float64ToFloat32( N * float64ToFloat32( alpha + x[ ix ] ) ); } - ix = offset; sum = 0.0; for ( i = 0; i < N; i++ ) { sum = float64ToFloat32( sum + float64ToFloat32( alpha+x[ ix ] ) ); - ix += stride; + ix += strideX; } return sum; } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.native.js index 4d73c87d4a67..a91e66406049 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/ndarray.native.js @@ -20,38 +20,31 @@ // MODULES // -var minViewBufferIndex = require( '@stdlib/strided/base/min-view-buffer-index' ); -var offsetView = require( '@stdlib/strided/base/offset-view' ); -var addon = require( './sapxsumors.native.js' ); +var addon = require( './../src/addon.node' ); // MAIN // /** -* Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. * * @param {PositiveInteger} N - number of indexed elements -* @param {number} alpha - constant +* @param {number} alpha - scalar constant * @param {Float32Array} x - input array -* @param {integer} stride - stride length -* @param {NonNegativeInteger} offset - starting index +* @param {integer} strideX - stride length +* @param {NonNegativeInteger} offsetX - starting index * @returns {number} sum * * @example * var Float32Array = require( '@stdlib/array/float32' ); -* var floor = require( '@stdlib/math/base/special/floor' ); * * var x = new Float32Array( [ 2.0, 1.0, 2.0, -2.0, -2.0, 2.0, 3.0, 4.0 ] ); -* var N = floor( x.length / 2 ); * -* var v = sapxsumors( N, 5.0, x, 2, 1 ); +* var v = sapxsumors( 4, 5.0, x, 2, 1 ); * // returns 25.0 */ -function sapxsumors( N, alpha, x, stride, offset ) { - var view; - offset = minViewBufferIndex( N, stride, offset ); - view = offsetView( x, offset ); - return addon( N, alpha, view, stride ); +function sapxsumors( N, alpha, x, strideX, offsetX ) { + return addon.ndarray( N, alpha, x, strideX, offsetX ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.js index 5fcfa36f4c99..fab6b509a1bb 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.js @@ -20,18 +20,19 @@ // MODULES // -var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' ); +var stride2offset = require( '@stdlib/strided/base/stride2offset' ); +var ndarray = require( './ndarray.js' ); // MAIN // /** -* Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. * * @param {PositiveInteger} N - number of indexed elements -* @param {number} alpha - constant +* @param {number} alpha - scalar constant * @param {Float32Array} x - input array -* @param {integer} stride - stride length +* @param {integer} strideX - stride length * @returns {number} sum * * @example @@ -42,28 +43,8 @@ var float64ToFloat32 = require( '@stdlib/number/float64/base/to-float32' ); * var v = sapxsumors( x.length, 5.0, x, 1 ); * // returns 16.0 */ -function sapxsumors( N, alpha, x, stride ) { - var sum; - var ix; - var i; - - if ( N <= 0 ) { - return 0.0; - } - if ( N === 1 || stride === 0 ) { - return float64ToFloat32( alpha + x[ 0 ] ); - } - if ( stride < 0 ) { - ix = (1-N) * stride; - } else { - ix = 0; - } - sum = 0.0; - for ( i = 0; i < N; i++ ) { - sum = float64ToFloat32( sum + float64ToFloat32( alpha+x[ ix ] ) ); - ix += stride; - } - return sum; +function sapxsumors( N, alpha, x, strideX ) { + return ndarray( N, alpha, x, strideX, stride2offset( N, strideX ) ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.native.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.native.js index 72233f89ce94..bd1e6ea6b7f4 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/lib/sapxsumors.native.js @@ -26,12 +26,12 @@ var addon = require( './../src/addon.node' ); // MAIN // /** -* Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. * * @param {PositiveInteger} N - number of indexed elements -* @param {number} alpha - constant +* @param {number} alpha - scalar constant * @param {Float32Array} x - input array -* @param {integer} stride - stride length +* @param {integer} strideX - stride length * @returns {number} sum * * @example @@ -42,8 +42,8 @@ var addon = require( './../src/addon.node' ); * var v = sapxsumors( x.length, 5.0, x, 1 ); * // returns 16.0 */ -function sapxsumors( N, alpha, x, stride ) { - return addon( N, alpha, x, stride ); +function sapxsumors( N, alpha, x, strideX ) { + return addon( N, alpha, x, strideX ); } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/manifest.json b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/manifest.json index 4cbcc4fe6950..b9ce4f08365b 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/manifest.json +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/manifest.json @@ -28,50 +28,53 @@ { "task": "build", "src": [ - "./src/sapxsumors.c" + "./src/main.c" ], "include": [ "./include" ], - "libraries": [ - "-lm" - ], + "libraries": [], "libpath": [], "dependencies": [ "@stdlib/napi/export", "@stdlib/napi/argv", "@stdlib/napi/argv-int64", "@stdlib/napi/argv-float", - "@stdlib/napi/argv-strided-float32array" + "@stdlib/napi/argv-strided-float32array", + "@stdlib/napi/create-double", + "@stdlib/strided/base/stride2offset", + "@stdlib/blas/base/shared" ] }, { "task": "benchmark", "src": [ - "./src/sapxsumors.c" + "./src/main.c" ], "include": [ "./include" ], - "libraries": [ - "-lm" - ], + "libraries": [], "libpath": [], - "dependencies": [] + "dependencies": [ + "@stdlib/strided/base/stride2offset", + "@stdlib/blas/base/shared" + ] }, { "task": "examples", "src": [ - "./src/sapxsumors.c" + "./src/main.c" ], "include": [ "./include" ], - "libraries": [ - "-lm" - ], + "libraries": [], "libpath": [], - "dependencies": [] + "dependencies": [ + "@stdlib/strided/base/stride2offset", + "@stdlib/blas/base/shared" + ] } ] } diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/package.json b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/package.json index baa5f58af75e..8bee47e776c8 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/package.json +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/package.json @@ -1,7 +1,7 @@ { "name": "@stdlib/blas/ext/base/sapxsumors", "version": "0.0.0", - "description": "Add a constant to each single-precision floating-point strided array element and compute the sum using ordinary recursive summation.", + "description": "Add a scalar constant to each single-precision floating-point strided array element and compute the sum using ordinary recursive summation.", "license": "Apache-2.0", "author": { "name": "The Stdlib Authors", diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/addon.c b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/addon.c index c7df8e36096b..ca870573be01 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/addon.c +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/addon.c @@ -17,13 +17,14 @@ */ #include "stdlib/blas/ext/base/sapxsumors.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_float.h" #include "stdlib/napi/argv_strided_float32array.h" +#include "stdlib/napi/create_double.h" #include -#include /** * Receives JavaScript callback invocation data. @@ -36,14 +37,28 @@ static napi_value addon( napi_env env, napi_callback_info info ) { STDLIB_NAPI_ARGV( env, info, argv, argc, 4 ); STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); STDLIB_NAPI_ARGV_FLOAT( env, alpha, argv, 1 ); - STDLIB_NAPI_ARGV_INT64( env, stride, argv, 3 ); - STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, stride, argv, 2 ); - - napi_value v; - napi_status status = napi_create_double( env, (double)stdlib_strided_sapxsumors( N, alpha, X, stride ), &v ); - assert( status == napi_ok ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 ); + STDLIB_NAPI_CREATE_DOUBLE( env, (double)API_SUFFIX(stdlib_strided_sapxsumors)( N, alpha, X, strideX ), 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, 5 ); + STDLIB_NAPI_ARGV_INT64( env, N, argv, 0 ); + STDLIB_NAPI_ARGV_FLOAT( env, alpha, argv, 1 ); + STDLIB_NAPI_ARGV_INT64( env, strideX, argv, 3 ); + STDLIB_NAPI_ARGV_INT64( env, offsetX, argv, 4 ); + STDLIB_NAPI_ARGV_STRIDED_FLOAT32ARRAY( env, X, N, strideX, argv, 2 ); + STDLIB_NAPI_CREATE_DOUBLE( env, (double)API_SUFFIX(stdlib_strided_sapxsumors_ndarray)( N, alpha, X, strideX, offsetX ), v ); return v; } -STDLIB_NAPI_MODULE_EXPORT_FCN( addon ) +STDLIB_NAPI_MODULE_EXPORT_FCN_WITH_METHOD( addon, "ndarray", addon_method ); diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/main.c b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/main.c new file mode 100644 index 000000000000..b213d45e8ff5 --- /dev/null +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/main.c @@ -0,0 +1,65 @@ +/** +* @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/blas/ext/base/sapxsumors.h" +#include "stdlib/blas/base/shared.h" +#include "stdlib/strided/base/stride2offset.h" + +/** +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. +* +* @param N number of indexed elements +* @param alpha scalar constant +* @param X input array +* @param strideX stride length +* @return output value +*/ +float API_SUFFIX(stdlib_strided_sapxsumors)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX ) { + CBLAS_INT ox = stdlib_strided_stride2offset( N, strideX ); + return API_SUFFIX(stdlib_strided_sapxsumors_ndarray)( N, alpha, X, strideX, ox ); +} + +/** +* Adds a scalar constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation and alternative indexing semantics. +* +* @param N number of indexed elements +* @param alpha scalar constant +* @param X input array +* @param strideX stride length +* @param offsetX starting index +* @return output value +*/ +float API_SUFFIX(stdlib_strided_sapxsumors_ndarray)( const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT strideX, const CBLAS_INT offsetX ) { + CBLAS_INT ix; + CBLAS_INT i; + float sum; + + if ( N <= 0 ) { + return 0.0f; + } + ix = offsetX; + if ( strideX == 0 ) { + return N * ( alpha + X[ ix ] ); + } + sum = 0.0f; + for ( i = 0; i < N; i++ ) { + sum += alpha + X[ ix ]; + ix += strideX; + } + return sum; +} diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/sapxsumors.c b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/sapxsumors.c deleted file mode 100644 index f1a8cb27559b..000000000000 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/src/sapxsumors.c +++ /dev/null @@ -1,53 +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/blas/ext/base/sapxsumors.h" -#include - -/** -* Adds a constant to each single-precision floating-point strided array element and computes the sum using ordinary recursive summation. -* -* @param N number of indexed elements -* @param alpha constant -* @param X input array -* @param stride stride length -* @return output value -*/ -float stdlib_strided_sapxsumors( const int64_t N, const float alpha, const float *X, const int64_t stride ) { - int64_t ix; - int64_t i; - float sum; - - if ( N <= 0 ) { - return 0.0f; - } - if ( N == 1 || stride == 0 ) { - return alpha + X[ 0 ]; - } - if ( stride < 0 ) { - ix = (1-N) * stride; - } else { - ix = 0; - } - sum = 0.0f; - for ( i = 0; i < N; i++ ) { - sum += alpha + X[ ix ]; - ix += stride; - } - return sum; -} diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.js index 1c74fa313a42..739cb722f1f9 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.js @@ -135,14 +135,14 @@ tape( 'the function supports a negative `stride` parameter', function test( t ) t.end(); }); -tape( 'if provided a `stride` parameter equal to `0`, the function returns the first indexed element plus a constant', function test( t ) { +tape( 'if provided a `stride` parameter equal to `0`, the function returns the first element plus a constant repeated N times', function test( t ) { var x; var v; x = new Float32Array( [ 1.0, -2.0, -4.0, 5.0, 3.0 ] ); v = sapxsumors( x.length, 5.0, x, 0, 0 ); - t.strictEqual( v, 6.0, 'returns expected value' ); + t.strictEqual( v, x.length * (x[0]+5.0), 'returns expected value' ); t.end(); }); diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.native.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.native.js index 0f84db8210cd..3df9cfbc3f74 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.ndarray.native.js @@ -144,14 +144,14 @@ tape( 'the function supports a negative `stride` parameter', opts, function test t.end(); }); -tape( 'if provided a `stride` parameter equal to `0`, the function returns the first indexed element plus a constant', opts, function test( t ) { +tape( 'if provided a `stride` parameter equal to `0`, the function returns the first element plus a constant repeated N times', opts, function test( t ) { var x; var v; x = new Float32Array( [ 1.0, -2.0, -4.0, 5.0, 3.0 ] ); v = sapxsumors( x.length, 5.0, x, 0, 0 ); - t.strictEqual( v, 6.0, 'returns expected value' ); + t.strictEqual( v, x.length * (x[0]+5.0), 'returns expected value' ); t.end(); }); diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.js index a86cf7172505..1ca0310af316 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.js @@ -135,14 +135,14 @@ tape( 'the function supports a negative `stride` parameter', function test( t ) t.end(); }); -tape( 'if provided a `stride` parameter equal to `0`, the function returns the first element plus a constant', function test( t ) { +tape( 'if provided a `stride` parameter equal to `0`, the function returns the first element plus a constant repeated N times', function test( t ) { var x; var v; x = new Float32Array( [ 1.0, -2.0, -4.0, 5.0, 3.0 ] ); v = sapxsumors( x.length, 5.0, x, 0 ); - t.strictEqual( v, 6.0, 'returns expected value' ); + t.strictEqual( v, x.length * (x[0]+5.0), 'returns expected value' ); t.end(); }); diff --git a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.native.js b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.native.js index da0f4b3a91e9..a6d8585410f6 100644 --- a/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.native.js +++ b/lib/node_modules/@stdlib/blas/ext/base/sapxsumors/test/test.sapxsumors.native.js @@ -253,14 +253,14 @@ tape( 'the function supports a negative `stride` parameter', opts, function test t.end(); }); -tape( 'if provided a `stride` parameter equal to `0`, the function returns the first element plus a constant', opts, function test( t ) { +tape( 'if provided a `stride` parameter equal to `0`, the function returns the first element plus a constant repeated N times', opts, function test( t ) { var x; var v; x = new Float32Array( [ 1.0, -2.0, -4.0, 5.0, 3.0 ] ); v = sapxsumors( x.length, 5.0, x, 0 ); - t.strictEqual( v, 6.0, 'returns expected value' ); + t.strictEqual( v, x.length * (x[0]+5.0), 'returns expected value' ); t.end(); }); diff --git a/lib/node_modules/@stdlib/constants/float32/ln-two/README.md b/lib/node_modules/@stdlib/constants/float32/ln-two/README.md index 26f54aeb49d5..f9871b04b97b 100644 --- a/lib/node_modules/@stdlib/constants/float32/ln-two/README.md +++ b/lib/node_modules/@stdlib/constants/float32/ln-two/README.md @@ -118,6 +118,13 @@ Macro for the [natural logarithm][@stdlib/math/base/special/lnf] of `2` as a sin @@ -130,6 +137,10 @@ Macro for the [natural logarithm][@stdlib/math/base/special/lnf] of `2` as a sin +[@stdlib/constants/float32/ln-ten]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/constants/float32/ln-ten + +[@stdlib/constants/float64/ln-two]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/constants/float64/ln-two +
diff --git a/lib/node_modules/@stdlib/constants/float32/sqrt-two/README.md b/lib/node_modules/@stdlib/constants/float32/sqrt-two/README.md index 1850d2f11598..720c319fe40c 100644 --- a/lib/node_modules/@stdlib/constants/float32/sqrt-two/README.md +++ b/lib/node_modules/@stdlib/constants/float32/sqrt-two/README.md @@ -120,6 +120,13 @@ Macro for the [square root][@stdlib/math/base/special/sqrtf] of `2` as a single- @@ -132,6 +139,10 @@ Macro for the [square root][@stdlib/math/base/special/sqrtf] of `2` as a single- +[@stdlib/constants/float32/sqrt-three]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/constants/float32/sqrt-three + +[@stdlib/constants/float64/sqrt-two]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/constants/float64/sqrt-two + diff --git a/lib/node_modules/@stdlib/lapack/base/docs/types/index.d.ts b/lib/node_modules/@stdlib/lapack/base/docs/types/index.d.ts index 78ba02be4ecd..e9b598da6965 100644 --- a/lib/node_modules/@stdlib/lapack/base/docs/types/index.d.ts +++ b/lib/node_modules/@stdlib/lapack/base/docs/types/index.d.ts @@ -56,11 +56,11 @@ interface Namespace { * @example * var Float64Array = require( '@stdlib/array/float64' ); * - * var A = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0 ] ); - * var B = new Float64Array( [ 0.0, 0.0, 11.0, 312.0, 53.0, 412.0 ] ); + * var A = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + * var B = new Float64Array( 6 ); * * ns.dlacpy.ndarray( 'all', 2, 2, A, 2, 1, 1, B, 2, 1, 2 ); - * // B => [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0 ] + * // B => [ 0.0, 0.0, 2.0, 3.0, 4.0, 5.0 ] */ dlacpy: typeof dlacpy; @@ -183,11 +183,11 @@ interface Namespace { * @example * var Float32Array = require( '@stdlib/array/float32' ); * - * var A = new Float32Array( [ 0.0, 1.0, 2.0, 3.0, 4.0 ] ); - * var B = new Float32Array( [ 0.0, 0.0, 11.0, 312.0, 53.0, 412.0 ] ); + * var A = new Float32Array( [ 1.0, 2.0, 3.0, 4.0, 5.0 ] ); + * var B = new Float32Array( 6 ); * * ns.slacpy.ndarray( 'all', 2, 2, A, 2, 1, 1, B, 2, 1, 2 ); - * // B => [ 0.0, 0.0, 1.0, 2.0, 3.0, 4.0 ] + * // B => [ 0.0, 0.0, 2.0, 3.0, 4.0, 5.0 ] */ slacpy: typeof slacpy; diff --git a/lib/node_modules/@stdlib/math/base/special/cosd/README.md b/lib/node_modules/@stdlib/math/base/special/cosd/README.md index 806287504ec9..d2663660456f 100644 --- a/lib/node_modules/@stdlib/math/base/special/cosd/README.md +++ b/lib/node_modules/@stdlib/math/base/special/cosd/README.md @@ -170,6 +170,12 @@ int main( void ) { @@ -182,6 +188,8 @@ int main( void ) { +[@stdlib/math/base/special/tand]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/math/base/special/tand + diff --git a/lib/node_modules/@stdlib/math/base/special/cround/README.md b/lib/node_modules/@stdlib/math/base/special/cround/README.md index 64a8f3f122a6..911092aadfe7 100644 --- a/lib/node_modules/@stdlib/math/base/special/cround/README.md +++ b/lib/node_modules/@stdlib/math/base/special/cround/README.md @@ -224,7 +224,7 @@ int main( void ) { ## See Also -- [`@stdlib/math/base/special/cceil`][@stdlib/math/base/special/cceil]: round a double-precision complex floating-point number toward positive infinity. +- [`@stdlib/math/base/special/cceil`][@stdlib/math/base/special/cceil]: round each component of a double-precision complex floating-point number toward positive infinity. - [`@stdlib/math/base/special/cfloor`][@stdlib/math/base/special/cfloor]: round a double-precision complex floating-point number toward negative infinity. - [`@stdlib/math/base/special/croundn`][@stdlib/math/base/special/croundn]: round each component of a double-precision complex floating-point number to the nearest multiple of 10^n. diff --git a/lib/node_modules/@stdlib/number/float32/base/signbit/README.md b/lib/node_modules/@stdlib/number/float32/base/signbit/README.md index 0f70925481c0..0176e61d0c06 100644 --- a/lib/node_modules/@stdlib/number/float32/base/signbit/README.md +++ b/lib/node_modules/@stdlib/number/float32/base/signbit/README.md @@ -18,7 +18,7 @@ limitations under the License. --> -# signbit +# signbitf > Return a boolean indicating if the sign bit for a [single-precision floating-point number][ieee754] is on (true) or off (false). @@ -32,7 +32,7 @@ var signbitf = require( '@stdlib/number/float32/base/signbit' ); #### signbitf( x ) -Returns a `boolean` indicating if the sign bit for a [single-precision floating-point number][ieee754] is on (`true`) or off (`false`). +Returns a boolean indicating if the sign bit for a [single-precision floating-point number][ieee754] is on (`true`) or off (`false`). ```javascript var toFloat32 = require( '@stdlib/number/float64/base/to-float32' ); @@ -82,6 +82,88 @@ for ( i = 0; i < 100; i++ ) { + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +#include "stdlib/number/float32/base/signbit.h" +``` + +#### stdlib_base_float32_signbit( x ) + +Returns an integer indicating whether the sign bit for a single-precision floating-point number is on (`1`) or off (`0`). + +```c +#include + +int8_t out = stdlib_base_float32_signbit( 3.14f ); +``` + +The function accepts the following arguments: + +- **x**: `[in] float` input value. + +```c +int8_t stdlib_base_float32_signbit( const float x ); +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +#include "stdlib/number/float32/base/signbit.h" +#include +#include +#include + +int main( void ) { + float x[] = { 3.14f, -3.14f, 0.0f, -0.0f, 4.0f, 1.0f, -1.0f, 1.0e38f, -1.0e38f }; + + int8_t out; + int i; + for ( i = 0; i < 9; i++ ) { + stdlib_base_float32_signbit( x[ i ], &out ); + printf( "%f => signbit: %" PRId8 "\n", x[ i ], out ); + } +} +``` + +
+