From 02343289ca0935aeb15f658e7952a5126884f3c8 Mon Sep 17 00:00:00 2001 From: iampratik13 Date: Wed, 24 Jun 2026 19:28:09 +0530 Subject: [PATCH] feat: add lapack/base/dlarnv --- .../@stdlib/lapack/base/dlarnv/README.md | 255 +++++++ .../lapack/base/dlarnv/benchmark/benchmark.js | 97 +++ .../dlarnv/benchmark/benchmark.ndarray.js | 97 +++ .../@stdlib/lapack/base/dlarnv/docs/repl.txt | 97 +++ .../lapack/base/dlarnv/docs/types/index.d.ts | 108 +++ .../lapack/base/dlarnv/docs/types/test.ts | 246 +++++++ .../lapack/base/dlarnv/examples/index.js | 44 ++ .../@stdlib/lapack/base/dlarnv/lib/dlarnv.js | 54 ++ .../@stdlib/lapack/base/dlarnv/lib/index.js | 72 ++ .../@stdlib/lapack/base/dlarnv/lib/main.js | 35 + .../@stdlib/lapack/base/dlarnv/lib/ndarray.js | 117 ++++ .../@stdlib/lapack/base/dlarnv/package.json | 71 ++ .../dlarnv/test/fixtures/idist1_medium.json | 44 ++ .../dlarnv/test/fixtures/idist1_small.json | 34 + .../dlarnv/test/fixtures/idist2_medium.json | 44 ++ .../dlarnv/test/fixtures/idist2_small.json | 34 + .../dlarnv/test/fixtures/idist3_medium.json | 44 ++ .../dlarnv/test/fixtures/idist3_small.json | 34 + .../fixtures/large_strides/idist1_medium.json | 72 ++ .../fixtures/large_strides/idist1_small.json | 52 ++ .../fixtures/large_strides/idist2_medium.json | 72 ++ .../fixtures/large_strides/idist2_small.json | 52 ++ .../fixtures/large_strides/idist3_medium.json | 72 ++ .../fixtures/large_strides/idist3_small.json | 52 ++ .../fixtures/mixed_strides/idist1_medium.json | 72 ++ .../fixtures/mixed_strides/idist1_small.json | 52 ++ .../fixtures/mixed_strides/idist2_medium.json | 72 ++ .../fixtures/mixed_strides/idist2_small.json | 52 ++ .../fixtures/mixed_strides/idist3_medium.json | 72 ++ .../fixtures/mixed_strides/idist3_small.json | 52 ++ .../negative_strides/idist1_medium.json | 44 ++ .../negative_strides/idist1_small.json | 34 + .../negative_strides/idist2_medium.json | 44 ++ .../negative_strides/idist2_small.json | 34 + .../negative_strides/idist3_medium.json | 44 ++ .../negative_strides/idist3_small.json | 34 + .../test/fixtures/offsets/idist1_medium.json | 48 ++ .../test/fixtures/offsets/idist1_small.json | 38 ++ .../test/fixtures/offsets/idist2_medium.json | 48 ++ .../test/fixtures/offsets/idist2_small.json | 38 ++ .../test/fixtures/offsets/idist3_medium.json | 48 ++ .../test/fixtures/offsets/idist3_small.json | 38 ++ .../lapack/base/dlarnv/test/test.dlarnv.js | 223 +++++++ .../@stdlib/lapack/base/dlarnv/test/test.js | 82 +++ .../lapack/base/dlarnv/test/test.ndarray.js | 629 ++++++++++++++++++ 45 files changed, 3697 insertions(+) create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/README.md create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.ndarray.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/docs/repl.txt create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/index.d.ts create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/test.ts create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/examples/index.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/lib/dlarnv.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/lib/index.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/lib/main.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/lib/ndarray.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/package.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_medium.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_small.json create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.dlarnv.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.js create mode 100644 lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.ndarray.js diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/README.md b/lib/node_modules/@stdlib/lapack/base/dlarnv/README.md new file mode 100644 index 000000000000..5ddfdc483f1a --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/README.md @@ -0,0 +1,255 @@ + + +# dlarnv + +> Return a vector of `n` random real numbers from a uniform or normal distribution. + +
+ +## Usage + +```javascript +var dlarnv = require( '@stdlib/lapack/base/dlarnv' ); +``` + +#### dlarnv( idist, iseed, n, x ) + +Returns a vector of `n` random real numbers from a uniform or normal distribution. + +```javascript +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); + +var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +var x = new Float64Array( 5 ); + +dlarnv( 1, iseed, 5, x ); +// x => +``` + +The function has the following parameters: + +- **idist**: distribution type. Must be `1` (uniform (0,1)), `2` (uniform (-1,1)), or `3` (normal (0,1)). +- **iseed**: [`Int32Array`][@stdlib/array/int32] seed array of four integers. Each element must be between `0` and `4095`, and `iseed[3]` must be odd. On exit, the seed is updated. +- **n**: number of random numbers to generate. +- **x**: output [`Float64Array`][@stdlib/array/float64]. + +Note that indexing is relative to the first index. To introduce an offset, use [`typed array`][mdn-typed-array] views. + + + +```javascript +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); + +// Initial arrays... +var iseed0 = new Int32Array( [ 0, 0, 1, 2, 3 ] ); +var x0 = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +// Create offset views... +var iseed1 = new Int32Array( iseed0.buffer, iseed0.BYTES_PER_ELEMENT*1 ); // start at 2nd element +var x1 = new Float64Array( x0.buffer, x0.BYTES_PER_ELEMENT*1 ); // start at 2nd element + +dlarnv( 1, iseed1, 5, x1 ); +// x0 => +``` + +#### dlarnv.ndarray( idist, iseed, strideIseed, offsetIseed, n, x, strideX, offsetX ) + +Returns a vector of `n` random real numbers from a uniform or normal distribution using alternative indexing semantics. + +```javascript +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); + +var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +var x = new Float64Array( 5 ); + +dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, 0 ); +// x => +``` + +The function has the following additional parameters: + +- **strideIseed**: stride length for `iseed`. +- **offsetIseed**: starting index for `iseed`. +- **strideX**: stride length for `x`. +- **offsetX**: starting index for `x`. + +While [`typed array`][mdn-typed-array] views mandate a view offset based on the underlying buffer, the offset parameters support indexing semantics based on starting indices. For example, + + + +```javascript +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); + +var iseed = new Int32Array( [ 0, 0, 1, 2, 3 ] ); +var x = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] ); + +dlarnv.ndarray( 1, iseed, 1, 1, 5, x, 1, 1 ); +// x => +``` + +
+ + + +
+ +## Notes + +- `dlarnv()` corresponds to the [LAPACK][LAPACK] function [`dlarnv`][lapack-dlarnv]. +- Unlike `dlaruv`, which generates at most 128 numbers per call, `dlarnv` has no limit on `N` as it calls `dlaruv` internally in batches. + +
+ + + +
+ +## Examples + + + +```javascript +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); +var dlarnv = require( '@stdlib/lapack/base/dlarnv' ); + +var iseed = new Int32Array( [ 1, 23, 456, 3795 ] ); +var x; + +// Uniform (0,1): +x = new Float64Array( 10 ); +dlarnv( 1, iseed, x.length, x ); +console.log( 'Uniform(0,1):', x ); + +// Uniform (-1,1): +x = new Float64Array( 10 ); +dlarnv( 2, iseed, x.length, x ); +console.log( 'Uniform(-1,1):', x ); + +// Normal (0,1): +x = new Float64Array( 10 ); +dlarnv( 3, iseed, x.length, x ); +console.log( 'Normal(0,1):', x ); +``` + +
+ + + + + +* * * + +
+ +## C APIs + + + +
+ +
+ + + + + +
+ +### Usage + +```c +TODO +``` + +#### TODO + +TODO. + +```c +TODO +``` + +TODO + +```c +TODO +``` + +
+ + + + + +
+ +
+ + + + + +
+ +### Examples + +```c +TODO +``` + +
+ + + +
+ + + + + +
+ +
+ + + + + +
+ +[lapack]: https://www.netlib.org/lapack/explore-html/ + +[lapack-dlarnv]: https://netlib.org/lapack/explore-html/d5/dd2/group__larnv.html + +[@stdlib/array/int32]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/int32 + +[@stdlib/array/float64]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/array/float64 + +[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray + +
+ + diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.js new file mode 100644 index 000000000000..c2bea63478db --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.js @@ -0,0 +1,97 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var dlarnv = require( './../lib/dlarnv.js' ); + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var iseed = new Int32Array( [ 1619, 1381, 1011, 2459 ] ); + var x = new Float64Array( len ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + dlarnv( 1, iseed, len, x ); + if ( isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var max; + var min; + var f; + var i; + + min = 0; // 2^min + max = 7; // 2^max + + for ( i = min; i <= max; i++ ) { + len = pow( 2, i ); + f = createBenchmark( len ); + bench( format( '%s:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.ndarray.js new file mode 100644 index 000000000000..a8165d53a8b1 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/benchmark/benchmark.ndarray.js @@ -0,0 +1,97 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var bench = require( '@stdlib/bench' ); +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); +var isnan = require( '@stdlib/math/base/assert/is-nan' ); +var pow = require( '@stdlib/math/base/special/pow' ); +var format = require( '@stdlib/string/format' ); +var pkg = require( './../package.json' ).name; +var dlarnv = require( './../lib/ndarray.js' ); + + +// FUNCTIONS // + +/** +* Creates a benchmark function. +* +* @private +* @param {PositiveInteger} len - array length +* @returns {Function} benchmark function +*/ +function createBenchmark( len ) { + var iseed = new Int32Array( [ 1619, 1381, 1011, 2459 ] ); + var x = new Float64Array( len ); + return benchmark; + + /** + * Benchmark function. + * + * @private + * @param {Benchmark} b - benchmark instance + */ + function benchmark( b ) { + var i; + + b.tic(); + for ( i = 0; i < b.iterations; i++ ) { + dlarnv( 1, iseed, 1, 0, len, x, 1, 0 ); + if ( isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + } + b.toc(); + if ( isnan( x[ 0 ] ) ) { + b.fail( 'should not return NaN' ); + } + b.pass( 'benchmark finished' ); + b.end(); + } +} + + +// MAIN // + +/** +* Main execution sequence. +* +* @private +*/ +function main() { + var len; + var max; + var min; + var f; + var i; + + min = 0; // 2^min + max = 7; // 2^max (capped at 128) + + for ( i = min; i <= max; i++ ) { + len = pow( 2, i ); + f = createBenchmark( len ); + bench( format( '%s::ndarray:len=%d', pkg, len ), f ); + } +} + +main(); diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/repl.txt b/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/repl.txt new file mode 100644 index 000000000000..53bd9820b7bd --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/repl.txt @@ -0,0 +1,97 @@ + +{{alias}}( idist, ISEED, N, X ) + Returns a vector of `N` random real numbers from a uniform or normal + distribution. + + Indexing is relative to the first index. To introduce an offset, use typed + array views. + + The function mutates the input seed array `ISEED`. + + Parameters + ---------- + idist: integer + Distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1). + + ISEED: Int32Array + Seed array of four integers. Each element must be between `0` and + `4095`, and `ISEED[3]` must be odd. + + N: integer + Number of random numbers to generate. + + X: Float64Array + Output array. + + Returns + ------- + X: Float64Array + Output array `X` filled with random numbers. + + Examples + -------- + > var ISEED = new {{alias:@stdlib/array/int32}}( [ 1, 23, 456, 3795 ] ); + > var X = new {{alias:@stdlib/array/float64}}( 10 ); + > {{alias}}( 1, ISEED, X.length, X ) + + + // Using typed array views: + > var ISEED0 = new {{alias:@stdlib/array/int32}}( [ 0, 1, 2, 3, 3 ] ); + > var X0 = new {{alias:@stdlib/array/float64}}( [ 0.0, 0.0, 0.0, 0.0 ] ); + > var ISEED1 = new Int32Array( ISEED0.buffer, ISEED0.BYTES_PER_ELEMENT*1 ); + > var X1 = new Float64Array( X0.buffer, X0.BYTES_PER_ELEMENT*1 ); + > {{alias}}( 1, ISEED1, 3, X1 ) + + + +{{alias}}.ndarray( idist, ISEED, strideISEED, offsetISEED, N, X, strideX, offsetX ) + Returns a vector of `N` random real numbers from a uniform or normal + distribution using alternative indexing semantics. + + While typed array views mandate a view offset based on the underlying + buffer, the offset parameters support indexing semantics based on starting + indices. + + The function mutates the input seed array `ISEED`. + + Parameters + ---------- + idist: integer + Distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1). + + ISEED: Int32Array + Seed array of four integers. Each element must be between `0` and + `4095`, and `ISEED[3]` must be odd. + + strideISEED: integer + Stride length for `ISEED`. + + offsetISEED: integer + Starting index for `ISEED`. + + N: integer + Number of random numbers to generate. + + X: Float64Array + Output array. + + strideX: integer + Stride length for `X`. + + offsetX: integer + Starting index for `X`. + + Returns + ------- + X: Float64Array + Output array `X` filled with random numbers. + + Examples + -------- + > var ISEED = new {{alias:@stdlib/array/int32}}( [ 1, 23, 456, 3795 ] ); + > var X = new {{alias:@stdlib/array/float64}}( 10 ); + > {{alias}}.ndarray( 1, ISEED, 1, 0, X.length, X, 1, 0 ) + + + See Also + -------- diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/index.d.ts b/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/index.d.ts new file mode 100644 index 000000000000..10963ae03d85 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/index.d.ts @@ -0,0 +1,108 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +// TypeScript Version: 4.1 + +/// + +/** +* Interface describing `dlarnv`. +*/ +interface Routine { + /** + * Returns a vector of `n` random real numbers from a uniform or normal distribution. + * + * @param idist - distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1) + * @param iseed - seed array of four integers + * @param n - number of random numbers to generate + * @param x - output array + * @returns output array + * + * @example + * var Int32Array = require( '@stdlib/array/int32' ); + * var Float64Array = require( '@stdlib/array/float64' ); + * + * var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + * var x = new Float64Array( 5 ); + * + * dlarnv( 1, iseed, 5, x ); + * // x => + */ + ( idist: number, iseed: Int32Array, n: number, x: Float64Array ): Float64Array; + + /** + * Returns a vector of `n` random real numbers from a uniform or normal distribution using alternative indexing semantics. + * + * @param idist - distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1) + * @param iseed - seed array of four integers + * @param strideIseed - stride length for `iseed` + * @param offsetIseed - starting index for `iseed` + * @param n - number of random numbers to generate + * @param x - output array + * @param strideX - stride length for `x` + * @param offsetX - starting index for `x` + * @returns output array + * + * @example + * var Int32Array = require( '@stdlib/array/int32' ); + * var Float64Array = require( '@stdlib/array/float64' ); + * + * var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + * var x = new Float64Array( 5 ); + * + * dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, 0 ); + * // x => + */ + ndarray( idist: number, iseed: Int32Array, strideIseed: number, offsetIseed: number, n: number, x: Float64Array, strideX: number, offsetX: number ): Float64Array; +} + +/** +* Returns a vector of `n` random real numbers from a uniform or normal distribution. +* +* @param idist - distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1) +* @param iseed - seed array of four integers +* @param n - number of random numbers to generate +* @param x - output array +* @returns output array +* +* @example +* var Int32Array = require( '@stdlib/array/int32' ); +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +* var x = new Float64Array( 5 ); +* +* dlarnv( 1, iseed, 5, x ); +* // x => +* +* @example +* var Int32Array = require( '@stdlib/array/int32' ); +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +* var x = new Float64Array( 5 ); +* +* dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, 0 ); +* // x => +*/ +declare var dlarnv: Routine; + + +// EXPORTS // + +export = dlarnv; diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/test.ts b/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/test.ts new file mode 100644 index 000000000000..ed07172a233b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/docs/types/test.ts @@ -0,0 +1,246 @@ +/* +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +import dlarnv = require( './index' ); + + +// TESTS // + +// The function returns a Float64Array... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv( 1, iseed, 5, x ); // $ExpectType Float64Array +} + +// The compiler throws an error if the function is provided a first argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv( '5', iseed, 5, x ); // $ExpectError + dlarnv( true, iseed, 5, x ); // $ExpectError + dlarnv( false, iseed, 5, x ); // $ExpectError + dlarnv( null, iseed, 5, x ); // $ExpectError + dlarnv( void 0, iseed, 5, x ); // $ExpectError + dlarnv( [], iseed, 5, x ); // $ExpectError + dlarnv( {}, iseed, 5, x ); // $ExpectError + dlarnv( ( x: number ): number => x, iseed, 5, x ); // $ExpectError +} + +// The compiler throws an error if the function is provided a second argument which is not an Int32Array... +{ + const x = new Float64Array( 5 ); + + dlarnv( 1, '5', 5, x ); // $ExpectError + dlarnv( 1, 5, 5, x ); // $ExpectError + dlarnv( 1, true, 5, x ); // $ExpectError + dlarnv( 1, false, 5, x ); // $ExpectError + dlarnv( 1, null, 5, x ); // $ExpectError + dlarnv( 1, void 0, 5, x ); // $ExpectError + dlarnv( 1, [], 5, x ); // $ExpectError + dlarnv( 1, {}, 5, x ); // $ExpectError + dlarnv( 1, ( x: number ): number => x, 5, x ); // $ExpectError +} + +// The compiler throws an error if the function is provided a third argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv( 1, iseed, '5', x ); // $ExpectError + dlarnv( 1, iseed, true, x ); // $ExpectError + dlarnv( 1, iseed, false, x ); // $ExpectError + dlarnv( 1, iseed, null, x ); // $ExpectError + dlarnv( 1, iseed, void 0, x ); // $ExpectError + dlarnv( 1, iseed, [], x ); // $ExpectError + dlarnv( 1, iseed, {}, x ); // $ExpectError + dlarnv( 1, iseed, ( x: number ): number => x, x ); // $ExpectError +} + +// The compiler throws an error if the function is provided a fourth argument which is not a Float64Array... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + + dlarnv( 1, iseed, 5, '5' ); // $ExpectError + dlarnv( 1, iseed, 5, 5 ); // $ExpectError + dlarnv( 1, iseed, 5, true ); // $ExpectError + dlarnv( 1, iseed, 5, false ); // $ExpectError + dlarnv( 1, iseed, 5, null ); // $ExpectError + dlarnv( 1, iseed, 5, void 0 ); // $ExpectError + dlarnv( 1, iseed, 5, [] ); // $ExpectError + dlarnv( 1, iseed, 5, {} ); // $ExpectError + dlarnv( 1, iseed, 5, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the function is provided an unsupported number of arguments... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv(); // $ExpectError + dlarnv( 1 ); // $ExpectError + dlarnv( 1, iseed ); // $ExpectError + dlarnv( 1, iseed, 5 ); // $ExpectError + dlarnv( 1, iseed, 5, x, 10 ); // $ExpectError +} + +// Attached to main export is an `ndarray` method which returns a Float64Array... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectType Float64Array +} + +// The compiler throws an error if the `ndarray` method is provided a first argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( '5', iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( true, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( false, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( null, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( void 0, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( [], iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( {}, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( ( x: number ): number => x, iseed, 1, 0, 5, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a second argument which is not an Int32Array... +{ + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, '5', 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, 5, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, true, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, false, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, null, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, void 0, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, [], 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, {}, 1, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, ( x: number ): number => x, 1, 0, 5, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a third argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, iseed, '5', 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, true, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, false, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, null, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, void 0, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, [], 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, {}, 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, ( x: number ): number => x, 0, 5, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fourth argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, iseed, 1, '5', 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, true, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, false, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, null, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, void 0, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, [], 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, {}, 5, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, ( x: number ): number => x, 5, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, iseed, 1, 0, '5', x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, true, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, false, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, null, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, void 0, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, [], x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, {}, x, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, ( x: number ): number => x, x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a Float64Array... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + + dlarnv.ndarray( 1, iseed, 1, 0, 5, '5', 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, 5, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, true, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, false, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, null, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, void 0, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, [], 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, {}, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, ( x: number ): number => x, 1, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, '5', 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, true, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, false, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, null, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, void 0, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, [], 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, {}, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, ( x: number ): number => x, 0 ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an eighth argument which is not a number... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, '5' ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, true ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, false ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, null ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, void 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, [] ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, {} ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, ( x: number ): number => x ); // $ExpectError +} + +// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments... +{ + const iseed = new Int32Array( [ 0, 1, 2, 3 ] ); + const x = new Float64Array( 5 ); + + dlarnv.ndarray(); // $ExpectError + dlarnv.ndarray( 1 ); // $ExpectError + dlarnv.ndarray( 1, iseed ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1 ); // $ExpectError + dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, 0, 10 ); // $ExpectError +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/examples/index.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/examples/index.js new file mode 100644 index 000000000000..ef0e624cb436 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/examples/index.js @@ -0,0 +1,44 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +var Int32Array = require( '@stdlib/array/int32' ); +var Float64Array = require( '@stdlib/array/float64' ); +var dlarnv = require( './../lib' ); + +var iseed = new Int32Array( [ 1, 23, 456, 3795 ] ); +var x; + +// Uniform (0,1): +x = new Float64Array( 10 ); +dlarnv( 1, iseed, x.length, x ); +console.log( 'Uniform(0,1):' ); +console.log( x ); + +// Uniform (-1,1): +x = new Float64Array( 10 ); +dlarnv( 2, iseed, x.length, x ); +console.log( 'Uniform(-1,1):' ); +console.log( x ); + +// Normal (0,1): +x = new Float64Array( 10 ); +dlarnv( 3, iseed, x.length, x ); +console.log( 'Normal(0,1):' ); +console.log( x ); diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/dlarnv.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/dlarnv.js new file mode 100644 index 000000000000..5d0a6894d72b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/dlarnv.js @@ -0,0 +1,54 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +/** +* Returns a vector of `n` random real numbers from a uniform or normal distribution. +* +* @param {integer} idist - distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1) +* @param {Int32Array} iseed - seed array of four integers +* @param {NonNegativeInteger} n - number of random numbers to generate +* @param {Float64Array} x - output array +* @returns {Float64Array} output array +* +* @example +* var Int32Array = require( '@stdlib/array/int32' ); +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +* var x = new Float64Array( 5 ); +* +* dlarnv( 1, iseed, 5, x ); +* // x => +*/ +function dlarnv( idist, iseed, n, x ) { + return ndarray( idist, iseed, 1, 0, n, x, 1, 0 ); +} + + +// EXPORTS // + +module.exports = dlarnv; diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/index.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/index.js new file mode 100644 index 000000000000..ef33d5f9aebc --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/index.js @@ -0,0 +1,72 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +/** +* LAPACK routine to return a vector of `n` random real numbers from a uniform or normal distribution. +* +* @module @stdlib/lapack/base/dlarnv +* +* @example +* var Int32Array = require( '@stdlib/array/int32' ); +* var Float64Array = require( '@stdlib/array/float64' ); +* var dlarnv = require( '@stdlib/lapack/base/dlarnv' ); +* +* var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +* var x = new Float64Array( 5 ); +* +* dlarnv( 1, iseed, 5, x ); +* // x => +* +* @example +* var Int32Array = require( '@stdlib/array/int32' ); +* var Float64Array = require( '@stdlib/array/float64' ); +* var dlarnv = require( '@stdlib/lapack/base/dlarnv' ); +* +* var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +* var x = new Float64Array( 5 ); +* +* dlarnv.ndarray( 1, iseed, 1, 0, 5, x, 1, 0 ); +* // x => +*/ + +// MODULES // + +var join = require( 'path' ).join; +var tryRequire = require( '@stdlib/utils/try-require' ); +var isError = require( '@stdlib/assert/is-error' ); +var main = require( './main.js' ); + + +// MAIN // + +var dlarnv; +var tmp = tryRequire( join( __dirname, './native.js' ) ); +if ( isError( tmp ) ) { + dlarnv = main; +} else { + dlarnv = tmp; +} + + +// EXPORTS // + +module.exports = dlarnv; + +// exports: { "ndarray": "dlarnv.ndarray" } diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/main.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/main.js new file mode 100644 index 000000000000..d860d2bdcbba --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/main.js @@ -0,0 +1,35 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var setReadOnly = require( '@stdlib/utils/define-nonenumerable-read-only-property' ); +var dlarnv = require( './dlarnv.js' ); +var ndarray = require( './ndarray.js' ); + + +// MAIN // + +setReadOnly( dlarnv, 'ndarray', ndarray ); + + +// EXPORTS // + +module.exports = dlarnv; diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/ndarray.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/ndarray.js new file mode 100644 index 000000000000..bfdde6bbe36e --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/lib/ndarray.js @@ -0,0 +1,117 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len */ + +'use strict'; + +// MODULES // + +var Float64Array = require( '@stdlib/array/float64' ); +var min = require( '@stdlib/math/base/special/min' ); +var ln = require( '@stdlib/math/base/special/ln' ); +var sqrt = require( '@stdlib/math/base/special/sqrt' ); +var cos = require( '@stdlib/math/base/special/cos' ); +var dlaruv = require( '@stdlib/lapack/base/dlaruv' ).ndarray; + + +// VARIABLES // + +var LV = 128; +var TWOPI = 6.28318530717958647692528676655900576839; + + +// MAIN // + +/** +* Returns a vector of `n` random real numbers from a uniform or normal distribution using alternative indexing semantics. +* +* @param {integer} idist - distribution type: 1=uniform(0,1), 2=uniform(-1,1), 3=normal(0,1) +* @param {Int32Array} iseed - seed array of four integers +* @param {integer} strideIseed - stride length for `iseed` +* @param {NonNegativeInteger} offsetIseed - starting index for `iseed` +* @param {NonNegativeInteger} n - number of random numbers to generate +* @param {Float64Array} x - output array +* @param {integer} strideX - stride length for `x` +* @param {NonNegativeInteger} offsetX - starting index for `x` +* @returns {Float64Array} output array +* +* @example +* var Int32Array = require( '@stdlib/array/int32' ); +* var Float64Array = require( '@stdlib/array/float64' ); +* +* var iseed = new Int32Array( [ 0, 1, 2, 3 ] ); +* var x = new Float64Array( 5 ); +* +* dlarnv( 1, iseed, 1, 0, 5, x, 1, 0 ); +* // x => +*/ +function dlarnv( idist, iseed, strideIseed, offsetIseed, n, x, strideX, offsetX ) { + var il2; + var il; + var iv; + var ox; + var u; + var i; + + if ( n <= 0 ) { + return x; + } + + u = new Float64Array( LV ); + ox = offsetX; + for ( iv = 0; iv < n; iv += LV / 2 ) { + il = min( LV / 2, n - iv ); + + // For normal distribution, we need twice as many uniform numbers (Box-Muller): + if ( idist === 3 ) { + il2 = 2 * il; + } else { + il2 = il; + } + + // Generate il2 uniform (0,1) numbers: + dlaruv( il2, iseed, strideIseed, offsetIseed, u, 1, 0 ); + + if ( idist === 1 ) { + // Uniform (0,1): just copy + for ( i = 0; i < il; i++ ) { + x[ ox ] = u[ i ]; + ox += strideX; + } + } else if ( idist === 2 ) { + // Uniform (-1,1): transform + for ( i = 0; i < il; i++ ) { + x[ ox ] = ( 2.0 * u[ i ] ) - 1.0; + ox += strideX; + } + } else if ( idist === 3 ) { + // Normal (0,1): Box-Muller transform + for ( i = 0; i < il; i++ ) { + x[ ox ] = sqrt( -2.0 * ln( u[ 2 * i ] ) ) * cos( TWOPI * u[ ( 2 * i ) + 1 ] ); + ox += strideX; + } + } + } + return x; +} + + +// EXPORTS // + +module.exports = dlarnv; diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/package.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/package.json new file mode 100644 index 000000000000..eaeaec65b0ab --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/package.json @@ -0,0 +1,71 @@ +{ + "name": "@stdlib/lapack/base/dlarnv", + "version": "0.0.0", + "description": "Return a vector of n random real numbers from a uniform or normal distribution.", + "license": "Apache-2.0", + "author": { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + }, + "contributors": [ + { + "name": "The Stdlib Authors", + "url": "https://github.com/stdlib-js/stdlib/graphs/contributors" + } + ], + "main": "./lib", + "directories": { + "benchmark": "./benchmark", + "doc": "./docs", + "example": "./examples", + "lib": "./lib", + "test": "./test" + }, + "types": "./docs/types", + "scripts": {}, + "homepage": "https://github.com/stdlib-js/stdlib", + "repository": { + "type": "git", + "url": "git://github.com/stdlib-js/stdlib.git" + }, + "bugs": { + "url": "https://github.com/stdlib-js/stdlib/issues" + }, + "dependencies": {}, + "devDependencies": {}, + "engines": { + "node": ">=0.10.0", + "npm": ">2.7.0" + }, + "os": [ + "aix", + "darwin", + "freebsd", + "linux", + "macos", + "openbsd", + "sunos", + "win32", + "windows" + ], + "keywords": [ + "stdlib", + "stdmath", + "mathematics", + "math", + "lapack", + "dlarnv", + "random", + "uniform", + "normal", + "prng", + "linear", + "algebra", + "subroutines", + "array", + "ndarray", + "float64", + "double", + "float64array" + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_medium.json new file mode 100644 index 000000000000..491c6f7e34c9 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_medium.json @@ -0,0 +1,44 @@ +{ + "idist": 1, + "N": 10, + "ISEED": [ + 100, + 200, + 300, + 401 + ], + "strideISEED": 1, + "offsetISEED": 0, + "X": [ + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 0, + "X_out": [ + 0.7221151762867244, + 0.3059602026783388, + 0.7615758157127068, + 0.6079652337045935, + 0.0786368892901983, + 0.4000379493120114, + 0.10082534071584703, + 0.23506144863762657, + 0.17051886842240904, + 0.009087995793297665 + ], + "ISEED_out": [ + 37, + 919, + 1099, + 2025 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_small.json new file mode 100644 index 000000000000..21fefea173a2 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist1_small.json @@ -0,0 +1,34 @@ +{ + "idist": 1, + "N": 5, + "ISEED": [ + 0, + 1, + 2, + 3 + ], + "strideISEED": 1, + "offsetISEED": 0, + "X": [ + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 0, + "X_out": [ + 0.13185632717602047, + 0.23381133728018355, + 0.321647369501445, + 0.7105395653957878, + 0.021072610098311628 + ], + "ISEED_out": [ + 86, + 1283, + 2995, + 1711 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_medium.json new file mode 100644 index 000000000000..6af524213874 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_medium.json @@ -0,0 +1,44 @@ +{ + "idist": 2, + "N": 10, + "ISEED": [ + 100, + 200, + 300, + 401 + ], + "strideISEED": 1, + "offsetISEED": 0, + "X": [ + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 0, + "X_out": [ + 0.4442303525734488, + -0.3880795946433224, + 0.5231516314254137, + 0.2159304674091871, + -0.8427262214196034, + -0.19992410137597716, + -0.7983493185683059, + -0.5298771027247469, + -0.6589622631551819, + -0.9818240084134047 + ], + "ISEED_out": [ + 37, + 919, + 1099, + 2025 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_small.json new file mode 100644 index 000000000000..6e4e5343935b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist2_small.json @@ -0,0 +1,34 @@ +{ + "idist": 2, + "N": 5, + "ISEED": [ + 0, + 1, + 2, + 3 + ], + "strideISEED": 1, + "offsetISEED": 0, + "X": [ + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 0, + "X_out": [ + -0.7362873456479591, + -0.5323773254396329, + -0.35670526099711, + 0.42107913079157555, + -0.9578547798033767 + ], + "ISEED_out": [ + 86, + 1283, + 2995, + 1711 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_medium.json new file mode 100644 index 000000000000..a523b20d2208 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_medium.json @@ -0,0 +1,44 @@ +{ + "idist": 3, + "N": 10, + "ISEED": [ + 100, + 200, + 300, + 401 + ], + "strideISEED": 1, + "offsetISEED": 0, + "X": [ + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 0, + "X_out": [ + -0.27791440803185835, + -0.5746523001829359, + -1.8247935085233213, + 0.20076905767884418, + 1.8778433928250893, + -0.6967215479708604, + 1.2867152740566736, + 0.7660468992130739, + 0.7068809898010754, + 1.3515922455097087 + ], + "ISEED_out": [ + 4085, + 368, + 2190, + 2177 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_small.json new file mode 100644 index 000000000000..9ed0690d869b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/idist3_small.json @@ -0,0 +1,34 @@ +{ + "idist": 3, + "N": 5, + "ISEED": [ + 0, + 1, + 2, + 3 + ], + "strideISEED": 1, + "offsetISEED": 0, + "X": [ + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 0, + "X_out": [ + 0.20440003275491594, + -0.3696254660077132, + -2.397244974224152, + -1.653688956647738, + -0.8840364685948497 + ], + "ISEED_out": [ + 1897, + 2065, + 2084, + 3723 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_medium.json new file mode 100644 index 000000000000..8b78287813af --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_medium.json @@ -0,0 +1,72 @@ +{ + "idist": 1, + "N": 10, + "ISEED": [ + 100, + 9999, + 200, + 9999, + 300, + 9999, + 401, + 9999 + ], + "strideISEED": 2, + "offsetISEED": 0, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.7221151762867244, + 9999, + 0.3059602026783388, + 9999, + 0.7615758157127068, + 9999, + 0.6079652337045935, + 9999, + 0.0786368892901983, + 9999, + 0.4000379493120114, + 9999, + 0.10082534071584703, + 9999, + 0.23506144863762657, + 9999, + 0.17051886842240904, + 9999, + 0.009087995793297665, + 9999 + ], + "ISEED_out": [ + 37, + 9999, + 919, + 9999, + 1099, + 9999, + 2025, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_small.json new file mode 100644 index 000000000000..611456d4b465 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist1_small.json @@ -0,0 +1,52 @@ +{ + "idist": 1, + "N": 5, + "ISEED": [ + 0, + 9999, + 1, + 9999, + 2, + 9999, + 3, + 9999 + ], + "strideISEED": 2, + "offsetISEED": 0, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.13185632717602047, + 9999, + 0.23381133728018355, + 9999, + 0.321647369501445, + 9999, + 0.7105395653957878, + 9999, + 0.021072610098311628, + 9999 + ], + "ISEED_out": [ + 86, + 9999, + 1283, + 9999, + 2995, + 9999, + 1711, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_medium.json new file mode 100644 index 000000000000..a618eec993af --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_medium.json @@ -0,0 +1,72 @@ +{ + "idist": 2, + "N": 10, + "ISEED": [ + 100, + 9999, + 200, + 9999, + 300, + 9999, + 401, + 9999 + ], + "strideISEED": 2, + "offsetISEED": 0, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.4442303525734488, + 9999, + -0.3880795946433224, + 9999, + 0.5231516314254137, + 9999, + 0.2159304674091871, + 9999, + -0.8427262214196034, + 9999, + -0.19992410137597716, + 9999, + -0.7983493185683059, + 9999, + -0.5298771027247469, + 9999, + -0.6589622631551819, + 9999, + -0.9818240084134047, + 9999 + ], + "ISEED_out": [ + 37, + 9999, + 919, + 9999, + 1099, + 9999, + 2025, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_small.json new file mode 100644 index 000000000000..53551ae1152c --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist2_small.json @@ -0,0 +1,52 @@ +{ + "idist": 2, + "N": 5, + "ISEED": [ + 0, + 9999, + 1, + 9999, + 2, + 9999, + 3, + 9999 + ], + "strideISEED": 2, + "offsetISEED": 0, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + -0.7362873456479591, + 9999, + -0.5323773254396329, + 9999, + -0.35670526099711, + 9999, + 0.42107913079157555, + 9999, + -0.9578547798033767, + 9999 + ], + "ISEED_out": [ + 86, + 9999, + 1283, + 9999, + 2995, + 9999, + 1711, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_medium.json new file mode 100644 index 000000000000..559832ccf9bf --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_medium.json @@ -0,0 +1,72 @@ +{ + "idist": 3, + "N": 10, + "ISEED": [ + 100, + 9999, + 200, + 9999, + 300, + 9999, + 401, + 9999 + ], + "strideISEED": 2, + "offsetISEED": 0, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + -0.27791440803185835, + 9999, + -0.5746523001829359, + 9999, + -1.8247935085233213, + 9999, + 0.20076905767884418, + 9999, + 1.8778433928250893, + 9999, + -0.6967215479708604, + 9999, + 1.2867152740566736, + 9999, + 0.7660468992130739, + 9999, + 0.7068809898010754, + 9999, + 1.3515922455097087, + 9999 + ], + "ISEED_out": [ + 4085, + 9999, + 368, + 9999, + 2190, + 9999, + 2177, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_small.json new file mode 100644 index 000000000000..821a9eec0cd9 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/large_strides/idist3_small.json @@ -0,0 +1,52 @@ +{ + "idist": 3, + "N": 5, + "ISEED": [ + 0, + 9999, + 1, + 9999, + 2, + 9999, + 3, + 9999 + ], + "strideISEED": 2, + "offsetISEED": 0, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.20440003275491594, + 9999, + -0.3696254660077132, + 9999, + -2.397244974224152, + 9999, + -1.653688956647738, + 9999, + -0.8840364685948497, + 9999 + ], + "ISEED_out": [ + 1897, + 9999, + 2065, + 9999, + 2084, + 9999, + 3723, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_medium.json new file mode 100644 index 000000000000..01db6f646036 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_medium.json @@ -0,0 +1,72 @@ +{ + "idist": 1, + "N": 10, + "ISEED": [ + 401, + 9999, + 300, + 9999, + 200, + 9999, + 100, + 9999 + ], + "strideISEED": -2, + "offsetISEED": 6, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.7221151762867244, + 9999, + 0.3059602026783388, + 9999, + 0.7615758157127068, + 9999, + 0.6079652337045935, + 9999, + 0.0786368892901983, + 9999, + 0.4000379493120114, + 9999, + 0.10082534071584703, + 9999, + 0.23506144863762657, + 9999, + 0.17051886842240904, + 9999, + 0.009087995793297665, + 9999 + ], + "ISEED_out": [ + 2025, + 9999, + 1099, + 9999, + 919, + 9999, + 37, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_small.json new file mode 100644 index 000000000000..2e8c921f58b5 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist1_small.json @@ -0,0 +1,52 @@ +{ + "idist": 1, + "N": 5, + "ISEED": [ + 3, + 9999, + 2, + 9999, + 1, + 9999, + 0, + 9999 + ], + "strideISEED": -2, + "offsetISEED": 6, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.13185632717602047, + 9999, + 0.23381133728018355, + 9999, + 0.321647369501445, + 9999, + 0.7105395653957878, + 9999, + 0.021072610098311628, + 9999 + ], + "ISEED_out": [ + 1711, + 9999, + 2995, + 9999, + 1283, + 9999, + 86, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_medium.json new file mode 100644 index 000000000000..ffddbae05fa8 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_medium.json @@ -0,0 +1,72 @@ +{ + "idist": 2, + "N": 10, + "ISEED": [ + 401, + 9999, + 300, + 9999, + 200, + 9999, + 100, + 9999 + ], + "strideISEED": -2, + "offsetISEED": 6, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.4442303525734488, + 9999, + -0.3880795946433224, + 9999, + 0.5231516314254137, + 9999, + 0.2159304674091871, + 9999, + -0.8427262214196034, + 9999, + -0.19992410137597716, + 9999, + -0.7983493185683059, + 9999, + -0.5298771027247469, + 9999, + -0.6589622631551819, + 9999, + -0.9818240084134047, + 9999 + ], + "ISEED_out": [ + 2025, + 9999, + 1099, + 9999, + 919, + 9999, + 37, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_small.json new file mode 100644 index 000000000000..7eed38a72f76 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist2_small.json @@ -0,0 +1,52 @@ +{ + "idist": 2, + "N": 5, + "ISEED": [ + 3, + 9999, + 2, + 9999, + 1, + 9999, + 0, + 9999 + ], + "strideISEED": -2, + "offsetISEED": 6, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + -0.7362873456479591, + 9999, + -0.5323773254396329, + 9999, + -0.35670526099711, + 9999, + 0.42107913079157555, + 9999, + -0.9578547798033767, + 9999 + ], + "ISEED_out": [ + 1711, + 9999, + 2995, + 9999, + 1283, + 9999, + 86, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_medium.json new file mode 100644 index 000000000000..f58b4cc67c67 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_medium.json @@ -0,0 +1,72 @@ +{ + "idist": 3, + "N": 10, + "ISEED": [ + 401, + 9999, + 300, + 9999, + 200, + 9999, + 100, + 9999 + ], + "strideISEED": -2, + "offsetISEED": 6, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + -0.27791440803185835, + 9999, + -0.5746523001829359, + 9999, + -1.8247935085233213, + 9999, + 0.20076905767884418, + 9999, + 1.8778433928250893, + 9999, + -0.6967215479708604, + 9999, + 1.2867152740566736, + 9999, + 0.7660468992130739, + 9999, + 0.7068809898010754, + 9999, + 1.3515922455097087, + 9999 + ], + "ISEED_out": [ + 2177, + 9999, + 2190, + 9999, + 368, + 9999, + 4085, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_small.json new file mode 100644 index 000000000000..13f1d6c2798a --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/mixed_strides/idist3_small.json @@ -0,0 +1,52 @@ +{ + "idist": 3, + "N": 5, + "ISEED": [ + 3, + 9999, + 2, + 9999, + 1, + 9999, + 0, + 9999 + ], + "strideISEED": -2, + "offsetISEED": 6, + "X": [ + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999, + 0, + 9999 + ], + "strideX": 2, + "offsetX": 0, + "X_out": [ + 0.20440003275491594, + 9999, + -0.3696254660077132, + 9999, + -2.397244974224152, + 9999, + -1.653688956647738, + 9999, + -0.8840364685948497, + 9999 + ], + "ISEED_out": [ + 3723, + 9999, + 2084, + 9999, + 2065, + 9999, + 1897, + 9999 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_medium.json new file mode 100644 index 000000000000..333de06c74cf --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_medium.json @@ -0,0 +1,44 @@ +{ + "idist": 1, + "N": 10, + "ISEED": [ + 401, + 300, + 200, + 100 + ], + "strideISEED": -1, + "offsetISEED": 3, + "X": [ + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": -1, + "offsetX": 9, + "X_out": [ + 0.009087995793297665, + 0.17051886842240904, + 0.23506144863762657, + 0.10082534071584703, + 0.4000379493120114, + 0.0786368892901983, + 0.6079652337045935, + 0.7615758157127068, + 0.3059602026783388, + 0.7221151762867244 + ], + "ISEED_out": [ + 2025, + 1099, + 919, + 37 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_small.json new file mode 100644 index 000000000000..66203e1c91b0 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist1_small.json @@ -0,0 +1,34 @@ +{ + "idist": 1, + "N": 5, + "ISEED": [ + 3, + 2, + 1, + 0 + ], + "strideISEED": -1, + "offsetISEED": 3, + "X": [ + 0, + 0, + 0, + 0, + 0 + ], + "strideX": -1, + "offsetX": 4, + "X_out": [ + 0.021072610098311628, + 0.7105395653957878, + 0.321647369501445, + 0.23381133728018355, + 0.13185632717602047 + ], + "ISEED_out": [ + 1711, + 2995, + 1283, + 86 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_medium.json new file mode 100644 index 000000000000..391157fd1f7b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_medium.json @@ -0,0 +1,44 @@ +{ + "idist": 2, + "N": 10, + "ISEED": [ + 401, + 300, + 200, + 100 + ], + "strideISEED": -1, + "offsetISEED": 3, + "X": [ + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": -1, + "offsetX": 9, + "X_out": [ + -0.9818240084134047, + -0.6589622631551819, + -0.5298771027247469, + -0.7983493185683059, + -0.19992410137597716, + -0.8427262214196034, + 0.2159304674091871, + 0.5231516314254137, + -0.3880795946433224, + 0.4442303525734488 + ], + "ISEED_out": [ + 2025, + 1099, + 919, + 37 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_small.json new file mode 100644 index 000000000000..82966932656b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist2_small.json @@ -0,0 +1,34 @@ +{ + "idist": 2, + "N": 5, + "ISEED": [ + 3, + 2, + 1, + 0 + ], + "strideISEED": -1, + "offsetISEED": 3, + "X": [ + 0, + 0, + 0, + 0, + 0 + ], + "strideX": -1, + "offsetX": 4, + "X_out": [ + -0.9578547798033767, + 0.42107913079157555, + -0.35670526099711, + -0.5323773254396329, + -0.7362873456479591 + ], + "ISEED_out": [ + 1711, + 2995, + 1283, + 86 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_medium.json new file mode 100644 index 000000000000..83bb2de3cf98 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_medium.json @@ -0,0 +1,44 @@ +{ + "idist": 3, + "N": 10, + "ISEED": [ + 401, + 300, + 200, + 100 + ], + "strideISEED": -1, + "offsetISEED": 3, + "X": [ + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": -1, + "offsetX": 9, + "X_out": [ + 1.3515922455097087, + 0.7068809898010754, + 0.7660468992130739, + 1.2867152740566736, + -0.6967215479708604, + 1.8778433928250893, + 0.20076905767884418, + -1.8247935085233213, + -0.5746523001829359, + -0.27791440803185835 + ], + "ISEED_out": [ + 2177, + 2190, + 368, + 4085 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_small.json new file mode 100644 index 000000000000..20f891ac07d8 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/negative_strides/idist3_small.json @@ -0,0 +1,34 @@ +{ + "idist": 3, + "N": 5, + "ISEED": [ + 3, + 2, + 1, + 0 + ], + "strideISEED": -1, + "offsetISEED": 3, + "X": [ + 0, + 0, + 0, + 0, + 0 + ], + "strideX": -1, + "offsetX": 4, + "X_out": [ + -0.8840364685948497, + -1.653688956647738, + -2.397244974224152, + -0.3696254660077132, + 0.20440003275491594 + ], + "ISEED_out": [ + 3723, + 2084, + 2065, + 1897 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_medium.json new file mode 100644 index 000000000000..e389cb534a76 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_medium.json @@ -0,0 +1,48 @@ +{ + "idist": 1, + "N": 10, + "ISEED": [ + 9999, + 100, + 200, + 300, + 401 + ], + "strideISEED": 1, + "offsetISEED": 1, + "X": [ + 9999, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 1, + "X_out": [ + 9999, + 0.7221151762867244, + 0.3059602026783388, + 0.7615758157127068, + 0.6079652337045935, + 0.0786368892901983, + 0.4000379493120114, + 0.10082534071584703, + 0.23506144863762657, + 0.17051886842240904, + 0.009087995793297665 + ], + "ISEED_out": [ + 9999, + 37, + 919, + 1099, + 2025 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_small.json new file mode 100644 index 000000000000..bfdc3623d47b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist1_small.json @@ -0,0 +1,38 @@ +{ + "idist": 1, + "N": 5, + "ISEED": [ + 9999, + 0, + 1, + 2, + 3 + ], + "strideISEED": 1, + "offsetISEED": 1, + "X": [ + 9999, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 1, + "X_out": [ + 9999, + 0.13185632717602047, + 0.23381133728018355, + 0.321647369501445, + 0.7105395653957878, + 0.021072610098311628 + ], + "ISEED_out": [ + 9999, + 86, + 1283, + 2995, + 1711 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_medium.json new file mode 100644 index 000000000000..4b25b03c6aa3 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_medium.json @@ -0,0 +1,48 @@ +{ + "idist": 2, + "N": 10, + "ISEED": [ + 9999, + 100, + 200, + 300, + 401 + ], + "strideISEED": 1, + "offsetISEED": 1, + "X": [ + 9999, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 1, + "X_out": [ + 9999, + 0.4442303525734488, + -0.3880795946433224, + 0.5231516314254137, + 0.2159304674091871, + -0.8427262214196034, + -0.19992410137597716, + -0.7983493185683059, + -0.5298771027247469, + -0.6589622631551819, + -0.9818240084134047 + ], + "ISEED_out": [ + 9999, + 37, + 919, + 1099, + 2025 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_small.json new file mode 100644 index 000000000000..c8a18e2dd788 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist2_small.json @@ -0,0 +1,38 @@ +{ + "idist": 2, + "N": 5, + "ISEED": [ + 9999, + 0, + 1, + 2, + 3 + ], + "strideISEED": 1, + "offsetISEED": 1, + "X": [ + 9999, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 1, + "X_out": [ + 9999, + -0.7362873456479591, + -0.5323773254396329, + -0.35670526099711, + 0.42107913079157555, + -0.9578547798033767 + ], + "ISEED_out": [ + 9999, + 86, + 1283, + 2995, + 1711 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_medium.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_medium.json new file mode 100644 index 000000000000..0673d1ca09bb --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_medium.json @@ -0,0 +1,48 @@ +{ + "idist": 3, + "N": 10, + "ISEED": [ + 9999, + 100, + 200, + 300, + 401 + ], + "strideISEED": 1, + "offsetISEED": 1, + "X": [ + 9999, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 1, + "X_out": [ + 9999, + -0.27791440803185835, + -0.5746523001829359, + -1.8247935085233213, + 0.20076905767884418, + 1.8778433928250893, + -0.6967215479708604, + 1.2867152740566736, + 0.7660468992130739, + 0.7068809898010754, + 1.3515922455097087 + ], + "ISEED_out": [ + 9999, + 4085, + 368, + 2190, + 2177 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_small.json b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_small.json new file mode 100644 index 000000000000..6ba6e6d8e12b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/fixtures/offsets/idist3_small.json @@ -0,0 +1,38 @@ +{ + "idist": 3, + "N": 5, + "ISEED": [ + 9999, + 0, + 1, + 2, + 3 + ], + "strideISEED": 1, + "offsetISEED": 1, + "X": [ + 9999, + 0, + 0, + 0, + 0, + 0 + ], + "strideX": 1, + "offsetX": 1, + "X_out": [ + 9999, + 0.20440003275491594, + -0.3696254660077132, + -2.397244974224152, + -1.653688956647738, + -0.8840364685948497 + ], + "ISEED_out": [ + 9999, + 1897, + 2065, + 2084, + 3723 + ] +} diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.dlarnv.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.dlarnv.js new file mode 100644 index 000000000000..02ac897289e4 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.dlarnv.js @@ -0,0 +1,223 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var Int32Array = require( '@stdlib/array/int32' ); +var dlarnv = require( './../lib/dlarnv.js' ); + + +// FIXTURES // + +var IDIST1_SMALL = require( './fixtures/idist1_small.json' ); +var IDIST1_MEDIUM = require( './fixtures/idist1_medium.json' ); +var IDIST2_SMALL = require( './fixtures/idist2_small.json' ); +var IDIST2_MEDIUM = require( './fixtures/idist2_medium.json' ); +var IDIST3_SMALL = require( './fixtures/idist3_small.json' ); +var IDIST3_MEDIUM = require( './fixtures/idist3_medium.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dlarnv, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 4', function test( t ) { + t.strictEqual( dlarnv.length, 4, 'returns expected value' ); + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a uniform (0,1) distribution (idist=1, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.N, x ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a uniform (0,1) distribution (idist=1, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST1_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.N, x ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a uniform (-1,1) distribution (idist=2, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST2_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.N, x ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a uniform (-1,1) distribution (idist=2, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST2_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.N, x ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a normal (0,1) distribution (idist=3, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST3_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.N, x ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a normal (0,1) distribution (idist=3, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST3_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.N, x ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function leaves output unchanged when `N` is equal to zero', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( 1, iseed, 0, x ); + t.deepEqual( x, new Float64Array( data.X ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function produces deterministic output for the same seed', function test( t ) { + var iseed1; + var iseed2; + var x1; + var x2; + + iseed1 = new Int32Array( [ 0, 1, 2, 3 ] ); + iseed2 = new Int32Array( [ 0, 1, 2, 3 ] ); + x1 = new Float64Array( 5 ); + x2 = new Float64Array( 5 ); + + dlarnv( 1, iseed1, 5, x1 ); + dlarnv( 1, iseed2, 5, x2 ); + + t.deepEqual( x1, x2, 'returns expected value' ); + t.deepEqual( iseed1, iseed2, 'returns expected value' ); + + t.end(); +}); + +tape( 'the function generates values in the open interval (0,1) for idist=1', function test( t ) { + var iseed; + var x; + var i; + + iseed = new Int32Array( [ 100, 200, 300, 401 ] ); + x = new Float64Array( 100 ); + + dlarnv( 1, iseed, 100, x ); + + for ( i = 0; i < x.length; i++ ) { + t.ok( x[ i ] > 0.0 && x[ i ] < 1.0, 'value at index ' + i + ' is in (0,1): ' + x[ i ] ); + } + + t.end(); +}); + +tape( 'the function generates values in the open interval (-1,1) for idist=2', function test( t ) { + var iseed; + var x; + var i; + + iseed = new Int32Array( [ 100, 200, 300, 401 ] ); + x = new Float64Array( 100 ); + + dlarnv( 2, iseed, 100, x ); + + for ( i = 0; i < x.length; i++ ) { + t.ok( x[ i ] > -1.0 && x[ i ] < 1.0, 'value at index ' + i + ' is in (-1,1): ' + x[ i ] ); + } + + t.end(); +}); diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.js new file mode 100644 index 000000000000..9507571da81b --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.js @@ -0,0 +1,82 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var proxyquire = require( 'proxyquire' ); +var IS_BROWSER = require( '@stdlib/assert/is-browser' ); +var dlarnv = require( './../lib' ); + + +// VARIABLES // + +var opts = { + 'skip': IS_BROWSER +}; + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dlarnv, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'attached to the main export is a method providing an ndarray interface', function test( t ) { + t.strictEqual( typeof dlarnv.ndarray, 'function', 'method is a function' ); + t.end(); +}); + +tape( 'if a native implementation is available, the main export is the native implementation', opts, function test( t ) { + var dlarnv = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( dlarnv, mock, 'returns expected value' ); + t.end(); + + function tryRequire() { + return mock; + } + + function mock() { + // Mock... + } +}); + +tape( 'if a native implementation is not available, the main export is a JavaScript implementation', opts, function test( t ) { + var dlarnv; + var main; + + main = require( './../lib/main.js' ); + + dlarnv = proxyquire( './../lib', { + '@stdlib/utils/try-require': tryRequire + }); + + t.strictEqual( dlarnv, main, 'returns expected value' ); + t.end(); + + function tryRequire() { + return new Error( 'Cannot find module' ); + } +}); diff --git a/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.ndarray.js b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.ndarray.js new file mode 100644 index 000000000000..d3f2abd780d0 --- /dev/null +++ b/lib/node_modules/@stdlib/lapack/base/dlarnv/test/test.ndarray.js @@ -0,0 +1,629 @@ +/** +* @license Apache-2.0 +* +* Copyright (c) 2026 The Stdlib Authors. +* +* Licensed under the Apache License, Version 2.0 (the "License"); +* you may not use this file except in compliance with the License. +* You may obtain a copy of the License at +* +* http://www.apache.org/licenses/LICENSE-2.0 +* +* Unless required by applicable law or agreed to in writing, software +* distributed under the License is distributed on an "AS IS" BASIS, +* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +* See the License for the specific language governing permissions and +* limitations under the License. +*/ + +/* eslint-disable max-len */ + +'use strict'; + +// MODULES // + +var tape = require( 'tape' ); +var Float64Array = require( '@stdlib/array/float64' ); +var Int32Array = require( '@stdlib/array/int32' ); +var dlarnv = require( './../lib/ndarray.js' ); + + +// FIXTURES // + +var IDIST1_SMALL = require( './fixtures/idist1_small.json' ); +var IDIST1_MEDIUM = require( './fixtures/idist1_medium.json' ); +var IDIST2_SMALL = require( './fixtures/idist2_small.json' ); +var IDIST2_MEDIUM = require( './fixtures/idist2_medium.json' ); +var IDIST3_SMALL = require( './fixtures/idist3_small.json' ); +var IDIST3_MEDIUM = require( './fixtures/idist3_medium.json' ); +var OFFSET_IDIST1_SMALL = require( './fixtures/offsets/idist1_small.json' ); +var OFFSET_IDIST1_MEDIUM = require( './fixtures/offsets/idist1_medium.json' ); +var OFFSET_IDIST2_SMALL = require( './fixtures/offsets/idist2_small.json' ); +var OFFSET_IDIST2_MEDIUM = require( './fixtures/offsets/idist2_medium.json' ); +var OFFSET_IDIST3_SMALL = require( './fixtures/offsets/idist3_small.json' ); +var OFFSET_IDIST3_MEDIUM = require( './fixtures/offsets/idist3_medium.json' ); +var NEGATIVE_STRIDES_IDIST1_SMALL = require( './fixtures/negative_strides/idist1_small.json' ); +var NEGATIVE_STRIDES_IDIST1_MEDIUM = require( './fixtures/negative_strides/idist1_medium.json' ); +var NEGATIVE_STRIDES_IDIST2_SMALL = require( './fixtures/negative_strides/idist2_small.json' ); +var NEGATIVE_STRIDES_IDIST2_MEDIUM = require( './fixtures/negative_strides/idist2_medium.json' ); +var NEGATIVE_STRIDES_IDIST3_SMALL = require( './fixtures/negative_strides/idist3_small.json' ); +var NEGATIVE_STRIDES_IDIST3_MEDIUM = require( './fixtures/negative_strides/idist3_medium.json' ); +var LARGE_STRIDES_IDIST1_SMALL = require( './fixtures/large_strides/idist1_small.json' ); +var LARGE_STRIDES_IDIST1_MEDIUM = require( './fixtures/large_strides/idist1_medium.json' ); +var LARGE_STRIDES_IDIST2_SMALL = require( './fixtures/large_strides/idist2_small.json' ); +var LARGE_STRIDES_IDIST2_MEDIUM = require( './fixtures/large_strides/idist2_medium.json' ); +var LARGE_STRIDES_IDIST3_SMALL = require( './fixtures/large_strides/idist3_small.json' ); +var LARGE_STRIDES_IDIST3_MEDIUM = require( './fixtures/large_strides/idist3_medium.json' ); +var MIXED_STRIDES_IDIST1_SMALL = require( './fixtures/mixed_strides/idist1_small.json' ); +var MIXED_STRIDES_IDIST1_MEDIUM = require( './fixtures/mixed_strides/idist1_medium.json' ); +var MIXED_STRIDES_IDIST2_SMALL = require( './fixtures/mixed_strides/idist2_small.json' ); +var MIXED_STRIDES_IDIST2_MEDIUM = require( './fixtures/mixed_strides/idist2_medium.json' ); +var MIXED_STRIDES_IDIST3_SMALL = require( './fixtures/mixed_strides/idist3_small.json' ); +var MIXED_STRIDES_IDIST3_MEDIUM = require( './fixtures/mixed_strides/idist3_medium.json' ); + + +// TESTS // + +tape( 'main export is a function', function test( t ) { + t.ok( true, __filename ); + t.strictEqual( typeof dlarnv, 'function', 'main export is a function' ); + t.end(); +}); + +tape( 'the function has an arity of 8', function test( t ) { + t.strictEqual( dlarnv.length, 8, 'returns expected value' ); + t.end(); +}); + +// --- IDIST=1 (uniform 0,1) base tests --- + +tape( 'the function returns a vector of random numbers from a uniform (0,1) distribution (idist=1, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a uniform (0,1) distribution (idist=1, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST1_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +// --- IDIST=2 (uniform -1,1) base tests --- + +tape( 'the function returns a vector of random numbers from a uniform (-1,1) distribution (idist=2, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST2_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a uniform (-1,1) distribution (idist=2, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST2_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +// --- IDIST=3 (normal 0,1) base tests --- + +tape( 'the function returns a vector of random numbers from a normal (0,1) distribution (idist=3, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST3_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function returns a vector of random numbers from a normal (0,1) distribution (idist=3, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST3_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED_out ), 'returns expected value' ); + + t.end(); +}); + +// --- Offset tests --- + +tape( 'the function supports providing index offsets (idist=1, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = OFFSET_IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing index offsets (idist=1, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = OFFSET_IDIST1_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing index offsets (idist=2, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = OFFSET_IDIST2_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing index offsets (idist=2, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = OFFSET_IDIST2_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing index offsets (idist=3, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = OFFSET_IDIST3_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing index offsets (idist=3, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = OFFSET_IDIST3_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +// --- Negative stride tests --- + +tape( 'the function supports providing negative strides (idist=1, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = NEGATIVE_STRIDES_IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing negative strides (idist=1, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = NEGATIVE_STRIDES_IDIST1_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing negative strides (idist=2, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = NEGATIVE_STRIDES_IDIST2_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing negative strides (idist=2, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = NEGATIVE_STRIDES_IDIST2_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing negative strides (idist=3, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = NEGATIVE_STRIDES_IDIST3_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing negative strides (idist=3, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = NEGATIVE_STRIDES_IDIST3_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +// --- Mixed stride tests --- + +tape( 'the function supports providing mixed strides (idist=1, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = MIXED_STRIDES_IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing mixed strides (idist=1, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = MIXED_STRIDES_IDIST1_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing mixed strides (idist=2, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = MIXED_STRIDES_IDIST2_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing mixed strides (idist=2, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = MIXED_STRIDES_IDIST2_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing mixed strides (idist=3, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = MIXED_STRIDES_IDIST3_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing mixed strides (idist=3, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = MIXED_STRIDES_IDIST3_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +// --- Large stride tests --- + +tape( 'the function supports providing positive strides (idist=1, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = LARGE_STRIDES_IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing positive strides (idist=1, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = LARGE_STRIDES_IDIST1_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing positive strides (idist=2, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = LARGE_STRIDES_IDIST2_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing positive strides (idist=2, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = LARGE_STRIDES_IDIST2_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing positive strides (idist=3, small N)', function test( t ) { + var iseed; + var data; + var x; + + data = LARGE_STRIDES_IDIST3_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function supports providing positive strides (idist=3, medium N)', function test( t ) { + var iseed; + var data; + var x; + + data = LARGE_STRIDES_IDIST3_MEDIUM; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( data.idist, iseed, data.strideISEED, data.offsetISEED, data.N, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X_out ), 'returns expected value' ); + + t.end(); +}); + +// --- Edge case tests --- + +tape( 'the function leaves output unchanged when `N` is equal to zero', function test( t ) { + var iseed; + var data; + var x; + + data = IDIST1_SMALL; + + iseed = new Int32Array( data.ISEED ); + x = new Float64Array( data.X ); + + dlarnv( 1, iseed, data.strideISEED, data.offsetISEED, 0, x, data.strideX, data.offsetX ); + t.deepEqual( x, new Float64Array( data.X ), 'returns expected value' ); + t.deepEqual( iseed, new Int32Array( data.ISEED ), 'returns expected value' ); + + t.end(); +}); + +tape( 'the function generates values in the open interval (0,1) for idist=1', function test( t ) { + var iseed; + var x; + var i; + + iseed = new Int32Array( [ 100, 200, 300, 401 ] ); + x = new Float64Array( 100 ); + + dlarnv( 1, iseed, 1, 0, 100, x, 1, 0 ); + + for ( i = 0; i < x.length; i++ ) { + t.ok( x[ i ] > 0.0 && x[ i ] < 1.0, 'value at index ' + i + ' is in (0,1): ' + x[ i ] ); + } + + t.end(); +}); + +tape( 'the function generates values in the open interval (-1,1) for idist=2', function test( t ) { + var iseed; + var x; + var i; + + iseed = new Int32Array( [ 100, 200, 300, 401 ] ); + x = new Float64Array( 100 ); + + dlarnv( 2, iseed, 1, 0, 100, x, 1, 0 ); + + for ( i = 0; i < x.length; i++ ) { + t.ok( x[ i ] > -1.0 && x[ i ] < 1.0, 'value at index ' + i + ' is in (-1,1): ' + x[ i ] ); + } + + t.end(); +});