diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/README.md b/lib/node_modules/@stdlib/lapack/base/dlasv2/README.md
new file mode 100644
index 000000000000..59fe5266d112
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/README.md
@@ -0,0 +1,262 @@
+
+
+# dlasv2
+
+> Compute the singular value decomposition of a 2x2 upper triangular matrix.
+
+
+
+The singular value decomposition of a 2x2 upper triangular matrix $ A = \begin{bmatrix} F & G \\ 0 & H \end{bmatrix}$ is given by
+
+
+
+```math
+\left[\begin{array}{rr}\mathrm{CSL} & \mathrm{SNL} \\
+-\mathrm{SNL} & \mathrm{CSL}\end{array}\right]
+
+\left[\begin{array}{rr}F & G \\
+0 & H\end{array}\right]
+
+\left[\begin{array}{rr}\mathrm{CSR} & -\mathrm{SNR} \\
+\mathrm{SNR} & \mathrm{CSR}\end{array}\right]
+
+=
+
+\left[\begin{array}{rr}\mathrm{SSMAX} & 0 \\
+0 & \mathrm{SSMIN}\end{array}\right]
+```
+
+
+
+where `CSL`,`SNL`,`CSR`,`SNR` define orthogonal rotations and $\mathrm{SSMAX}^2, \mathrm{SSMIN}^2$ are the eigenvalues of $A^{T}A$.
+
+
+
+
+
+
+
+## Usage
+
+```javascript
+var dlasv2 = require( '@stdlib/lapack/base/dlasv2' );
+```
+
+#### dlasv2( F, G, H, out )
+
+Computes the singular value decomposition of a 2x2 upper triangular matrix.
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+var out = new Float64Array( 6 );
+
+var v = dlasv2( 2.0, 3.0, 4.0, out );
+// returns [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+
+var bool = ( v === out );
+// returns true
+```
+
+The function has the following parameters:
+
+- **F**: the (0,0) element of matrix.
+- **G**: the (0,1) element of matrix.
+- **H**: the (1,1) element of matrix.
+- **out**: [`Float64Array`][@stdlib/array/float64] output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively.
+
+On return, `abs(SSMAX)` is the larger singular value, `abs(SSMIN)` is the smaller singular value, and (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)` as described in the decomposition.
+
+Note that indexing is relative to the first index. To introduce an offset, use [`typed array`][mdn-typed-array] views.
+
+
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+// Initial arrays...
+var OUT0 = new Float64Array( [ 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+// Create offset views...
+var OUT1 = new Float64Array( OUT0.buffer, OUT0.BYTES_PER_ELEMENT*1 ); // start at 2nd element
+
+dlasv2( 2.0, 3.0, 4.0, OUT1 );
+// OUT0 => [ 0.0, ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+```
+
+#### dlasv2.ndarray( F, G, H, out, strideOut, offsetOut )
+
+Computes the singular value decomposition of a 2x2 upper triangular matrix using alternative indexing semantics.
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+var out = new Float64Array( 6 );
+
+dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 0 );
+// out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+```
+
+The function has the following additional parameters:
+
+- **strideOut**: stride length for `Out`.
+- **offsetOut**: starting index for `Out`.
+
+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 Float64Array = require( '@stdlib/array/float64' );
+var out = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
+
+dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 1 );
+// out => [ 0.0, ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+```
+
+
+
+
+
+
+
+## Notes
+
+- `dlasv2()` corresponds to the [LAPACK][LAPACK] routine [`dlasv2`][lapack-dlasv2].
+
+
+
+
+
+
+
+## Examples
+
+
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+var dlasv2 = require( '@stdlib/lapack/base/dlasv2' );
+
+var out = new Float64Array( [ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ] );
+
+dlasv2( 1.0, 1e6, 2.0, out );
+console.log( out );
+
+dlasv2.ndarray( 1.0, 1e6, 2.0, out, 1, 0 );
+console.log( out );
+```
+
+
+
+
+
+
+
+* * *
+
+
+
+## C APIs
+
+
+
+
+
+
+
+
+
+
+
+### Usage
+
+```c
+TODO
+```
+
+#### TODO
+
+TODO.
+
+```c
+TODO
+```
+
+TODO
+
+```c
+TODO
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+### Examples
+
+```c
+TODO
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+[lapack]: https://www.netlib.org/lapack/explore-html/
+
+[lapack-dlasv2]: https://www.netlib.org/lapack/explore-html/d8/da7/group__lasv2_ga96f9f244300d82921950e2c393b4b20f.html#ga96f9f244300d82921950e2c393b4b20f
+
+[@stdlib/array/float64]: https://stdlib.io/docs/api/latest/@stdlib/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/dlasv2/benchmark/benchmark.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/benchmark/benchmark.js
new file mode 100644
index 000000000000..1ec202f960e8
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/benchmark/benchmark.js
@@ -0,0 +1,69 @@
+/**
+* @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 uniform = require( '@stdlib/random/array/uniform' );
+var Float64Array = require( '@stdlib/array/float64' );
+var format = require( '@stdlib/string/format' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var pkg = require( './../package.json' ).name;
+var dlasv2 = require( './../lib/dlasv2.js' );
+
+
+// VARIABLES //
+
+var options = {
+ 'dtype': 'float64'
+};
+
+
+// MAIN //
+
+bench( format( '%s', pkg ), function benchmark( b ) {
+ var out;
+ var N;
+ var f;
+ var g;
+ var h;
+ var i;
+
+ N = 100;
+ f = uniform( N, -500.0, 500.0, options );
+ g = uniform( N, -500.0, 500.0, options );
+ h = uniform( N, -500.0, 500.0, options );
+
+ out = new Float64Array( 6 );
+
+ b.tic();
+ for ( i = 0; i < b.iterations; i++ ) {
+ dlasv2( f[ i%f.length ], g[ i%g.length ], h[ i%h.length ], out );
+ if ( isnan( out[ i%out.length ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ }
+ b.toc();
+ if ( isnan( out[ i%out.length ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ b.pass( 'benchmark finished' );
+ b.end();
+});
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/benchmark/benchmark.ndarray.js
new file mode 100644
index 000000000000..b309bc51a4a2
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/benchmark/benchmark.ndarray.js
@@ -0,0 +1,69 @@
+/**
+* @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 uniform = require( '@stdlib/random/array/uniform' );
+var Float64Array = require( '@stdlib/array/float64' );
+var format = require( '@stdlib/string/format' );
+var isnan = require( '@stdlib/math/base/assert/is-nan' );
+var pkg = require( './../package.json' ).name;
+var dlasv2 = require( './../lib/ndarray.js' );
+
+
+// VARIABLES //
+
+var options = {
+ 'dtype': 'float64'
+};
+
+
+// MAIN //
+
+bench( format( '%s:ndarray', pkg ), function benchmark( b ) {
+ var out;
+ var N;
+ var f;
+ var g;
+ var h;
+ var i;
+
+ N = 100;
+ f = uniform( N, -500.0, 500.0, options );
+ g = uniform( N, -500.0, 500.0, options );
+ h = uniform( N, -500.0, 500.0, options );
+
+ out = new Float64Array( 6 );
+
+ b.tic();
+ for ( i = 0; i < b.iterations; i++ ) {
+ dlasv2( f[ i%f.length ], g[ i%g.length ], h[ i%h.length ], out, 1, 0 );
+ if ( isnan( out[ i%out.length ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ }
+ b.toc();
+ if ( isnan( out[ i%out.length ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ b.pass( 'benchmark finished' );
+ b.end();
+});
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/repl.txt b/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/repl.txt
new file mode 100644
index 000000000000..2d7095b4e433
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/repl.txt
@@ -0,0 +1,113 @@
+
+{{alias}}( F, G, H, out )
+ Computes the singular value decomposition of a 2x2 upper triangular matrix.
+
+ Indexing is relative to the first index. To introduce an offset, use typed
+ array views.
+
+ Parameters
+ ----------
+ F: number
+ The (0,0) element of matrix.
+
+ G: number
+ The (0,1) element of matrix.
+
+ H: number
+ The (1,1) element of matrix.
+
+ out: Float64Array
+ Output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and
+ `CSL` respectively.
+
+ Returns
+ -------
+ out: Float64Array
+ Output array.
+
+ Note that on return, `abs(SSMAX)` is the larger singular value,
+ `abs(SSMIN)` is the smaller singular value, and (`CSL`,`SNL`)
+ and (`CSR`,`SNR`) are the left and right singular vectors for
+ `abs(SSMAX)`.
+
+ Examples
+ --------
+ > var out = new {{alias:@stdlib/array/float64}}( 6 );
+ > {{alias}}( 2.0, 3.0, 4.0, out );
+ > out[ 0 ]
+ ~1.5513
+ > out[ 1 ]
+ ~5.1569
+ > out[ 2 ]
+ ~0.9665
+ > out[ 3 ]
+ ~0.2567
+ > out[ 4 ]
+ ~0.7497
+ > out[ 5 ]
+ ~0.6618
+
+
+{{alias}}.ndarray( F, G, H, out, strideOut, offsetOut )
+ Computes the singular value decomposition of a 2x2 upper triangular
+ matrix using alternative indexing semantics.
+
+ Note that on return, `abs(SSMAX)` is the larger singular value,
+ `abs(SSMIN)` is the smaller singular value, and (`CSL`,`SNL`)
+ and (`CSR`,`SNR`) are the left and right singular vectors for
+ `abs(SSMAX)`.
+
+ While typed array views mandate a view offset based on the underlying
+ buffer, the offset parameter support indexing semantics based on
+ a starting index.
+
+ Parameters
+ ----------
+ F: number
+ The (0,0) element of matrix.
+
+ G: number
+ The (0,1) element of matrix.
+
+ H: number
+ The (1,1) element of matrix.
+
+ out: Float64Array
+ Output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`,
+ and `CSL` respectively.
+
+ strideOut: integer
+ Stride length for `Out`.
+
+ offsetOut: integer
+ Starting index for `Out`.
+
+ Returns
+ -------
+ out: Float64Array
+ Output array.
+
+ Note that on return, `abs(SSMAX)` is the larger singular value,
+ `abs(SSMIN)` is the smaller singular value, and (`CSL`,`SNL`)
+ and (`CSR`,`SNR`) are the left and right singular vectors for
+ `abs(SSMAX)`.
+
+ Examples
+ --------
+ > var out = new {{alias:@stdlib/array/float64}}( 6 );
+ > {{alias}}.ndarray( 2.0, 3.0, 4.0, out, 1, 0 );
+ > out[ 0 ]
+ ~1.5513
+ > out[ 1 ]
+ ~5.1569
+ > out[ 2 ]
+ ~0.9665
+ > out[ 3 ]
+ ~0.2567
+ > out[ 4 ]
+ ~0.7497
+ > out[ 5 ]
+ ~0.6618
+
+ See Also
+ --------
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/types/index.d.ts b/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/types/index.d.ts
new file mode 100644
index 000000000000..d40b45104f42
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/types/index.d.ts
@@ -0,0 +1,116 @@
+/**
+* @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 `dlasv2`.
+*/
+interface Routine {
+ /**
+ * Computes the singular value decomposition of a 2x2 upper triangular matrix.
+ *
+ * ## Notes
+ *
+ * - `abs(SSMAX)` is the larger singular value.
+ * - `abs(SSMIN)` is the smaller singular value.
+ * - (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)`.
+ *
+ * @param F - the (0,0) element of matrix
+ * @param G - the (0,1) element of matrix
+ * @param H - the (1,1) element of matrix
+ * @param out - output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively
+ * @returns output array
+ *
+ * @example
+ * var Float64Array = require( '@stdlib/array/float64' );
+ *
+ * var out = new Float64Array( 6 );
+ *
+ * dlasv2( 2.0, 3.0, 4.0, out );
+ * // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+ */
+ ( F: number, G: number, H: number, out: Float64Array ): Float64Array;
+
+ /**
+ * Computes the singular value decomposition of a 2x2 upper triangular matrix using alternative indexing semantics.
+ *
+ * ## Notes
+ *
+ * - `abs(SSMAX)` is the larger singular value.
+ * - `abs(SSMIN)` is the smaller singular value.
+ * - (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)`.
+ *
+ * @param F - the (0,0) element of matrix
+ * @param G - the (0,1) element of matrix
+ * @param H - the (1,1) element of matrix
+ * @param out - output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively
+ * @param strideOut - stride length for `out`
+ * @param offsetOut - starting index of `out`
+ * @returns output array
+ *
+ * @example
+ * var Float64Array = require( '@stdlib/array/float64' );
+ *
+ * var out = new Float64Array( 6 );
+ *
+ * dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 0 );
+ * // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+ */
+ ndarray( F: number, G: number, H: number, out: Float64Array, strideOut: number, offsetOut: number ): Float64Array;
+}
+
+/**
+* Computes the singular value decomposition of a 2x2 upper triangular matrix.
+*
+* ## Notes
+*
+* - `abs(SSMAX)` is the larger singular value.
+* - `abs(SSMIN)` is the smaller singular value.
+* - (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)`.
+*
+* @param F - the (0,0) element of matrix
+* @param G - the (0,1) element of matrix
+* @param H - the (1,1) element of matrix
+* @param out - output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively
+* @returns output array
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var out = new Float64Array( 6 );
+*
+* dlasv2( 2.0, 3.0, 4.0, out );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var out = new Float64Array( 6 );
+*
+* dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 0 );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*/
+declare var dlasv2: Routine;
+
+
+// EXPORTS //
+
+export = dlasv2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/types/test.ts b/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/types/test.ts
new file mode 100644
index 000000000000..1adfd989e517
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/docs/types/test.ts
@@ -0,0 +1,200 @@
+/*
+* @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 dlasv2 = require( './index' );
+
+
+// TESTS //
+
+// The function returns a Float64Array...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2( 2.0, 3.0, 4.0, out ); // $ExpectType Float64Array
+}
+
+// The compiler throws an error if the function is provided a first argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2( '5', 3.0, 4.0, out ); // $ExpectError
+ dlasv2( true, 3.0, 4.0, out ); // $ExpectError
+ dlasv2( false, 3.0, 4.0, out ); // $ExpectError
+ dlasv2( null, 3.0, 4.0, out ); // $ExpectError
+ dlasv2( void 0, 3.0, 4.0, out ); // $ExpectError
+ dlasv2( [], 3.0, 4.0, out ); // $ExpectError
+ dlasv2( {}, 3.0, 4.0, out ); // $ExpectError
+ dlasv2( ( x: number ): number => x, 3.0, 4.0, out ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a second argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2( 2.0, '5', 4.0, out ); // $ExpectError
+ dlasv2( 2.0, true, 4.0, out ); // $ExpectError
+ dlasv2( 2.0, false, 4.0, out ); // $ExpectError
+ dlasv2( 2.0, null, 4.0, out ); // $ExpectError
+ dlasv2( 2.0, void 0, 4.0, out ); // $ExpectError
+ dlasv2( 2.0, [], 4.0, out ); // $ExpectError
+ dlasv2( 2.0, {}, 4.0, out ); // $ExpectError
+ dlasv2( 2.0, ( x: number ): number => x, 4.0, out ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a third argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2( 2.0, 3.0, '5', out ); // $ExpectError
+ dlasv2( 2.0, 3.0, true, out ); // $ExpectError
+ dlasv2( 2.0, 3.0, false, out ); // $ExpectError
+ dlasv2( 2.0, 3.0, null, out ); // $ExpectError
+ dlasv2( 2.0, 3.0, void 0, out ); // $ExpectError
+ dlasv2( 2.0, 3.0, [], out ); // $ExpectError
+ dlasv2( 2.0, 3.0, {}, out ); // $ExpectError
+ dlasv2( 2.0, 3.0, ( x: number ): number => x, out ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fourth argument which is not a Float64Array...
+{
+
+ dlasv2( 2.0, 3.0, 4.0, '5' ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, 5 ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, true ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, false ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, null ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, void 0 ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, [] ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, {} ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided an unsupported number of arguments...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2(); // $ExpectError
+ dlasv2( 2.0 ); // $ExpectError
+ dlasv2( 2.0, 3.0 ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0 ); // $ExpectError
+ dlasv2( 2.0, 3.0, 4.0, out, 10 ); // $ExpectError
+}
+
+// Attached to main export is an `ndarray` method which returns a Float64Array...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 0 ); // $ExpectType Float64Array
+}
+
+// The compiler throws an error if the function is provided a first argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray( '5', 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( true, 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( false, 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( null, 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( void 0, 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( [], 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( {}, 3.0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( ( x: number ): number => x, 3.0, 4.0, out, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a second argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray( 2.0, '5', 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, true, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, false, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, null, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, void 0, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, [], 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, {}, 4.0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, ( x: number ): number => x, 4.0, out, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a third argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray( 2.0, 3.0, '5', out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, true, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, false, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, null, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, void 0, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, [], out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, {}, out, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, ( x: number ): number => x, out, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fourth argument which is not a Float64Array...
+{
+
+ dlasv2.ndarray( 2.0, 3.0, 4.0, '5', 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, 5, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, true, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, false, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, null, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, void 0, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, [], 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, {}, 1, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, ( x: number ): number => x, 1, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fifth argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, '5', 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, true, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, false, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, null, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, void 0, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, [], 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, {}, 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, ( x: number ): number => x, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a sixth argument which is not a number...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, '5' ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, true ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, false ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, null ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, void 0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, [] ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, {} ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments...
+{
+ const out = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 ] );
+
+ dlasv2.ndarray(); // $ExpectError
+ dlasv2.ndarray( 2.0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1 ); // $ExpectError
+ dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 0, 10 ); // $ExpectError
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/examples/index.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/examples/index.js
new file mode 100644
index 000000000000..2cfba28f0ae6
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/examples/index.js
@@ -0,0 +1,30 @@
+/**
+* @license Apache-2.0
+*
+* Copyright (c) 2026 The Stdlib Authors.
+*
+* Licensed under the Apache License, Version 2.0 (the "License");
+* you may not use this file except in compliance with the License.
+* You may obtain a copy of the License at
+*
+* http://www.apache.org/licenses/LICENSE-2.0
+*
+* Unless required by applicable law or agreed to in writing, software
+* distributed under the License is distributed on an "AS IS" BASIS,
+* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+* See the License for the specific language governing permissions and
+* limitations under the License.
+*/
+
+'use strict';
+
+var Float64Array = require( '@stdlib/array/float64' );
+var dlasv2 = require( './../lib/' );
+
+var out = new Float64Array( 6 );
+
+dlasv2( 1.0, 1e6, 2.0, out );
+console.log( out );
+
+dlasv2.ndarray( 1.0, 1e6, 2.0, out, 1, 0 );
+console.log( out );
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/base.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/base.js
new file mode 100644
index 000000000000..e97cf89473d3
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/base.js
@@ -0,0 +1,245 @@
+/**
+* @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-statements */
+
+'use strict';
+
+// MODULES //
+
+var dlamch = require( '@stdlib/lapack/base/dlamch' );
+var sqrt = require( '@stdlib/math/base/special/sqrt' );
+var abs = require( '@stdlib/math/base/special/abs' );
+var copysign = require( '@stdlib/math/base/special/copysign' );
+
+
+// VARIABLES //
+
+var EPS = dlamch( 'E' );
+
+
+// MAIN //
+
+/**
+* Computes the singular value decomposition of a 2x2 upper triangular matrix.
+*
+* ## Notes
+*
+* - `abs(SSMAX)` is the larger singular value.
+* - `abs(SSMIN)` is the smaller singular value.
+* - (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)`.
+*
+* @private
+* @param {number} F - the (0,0) element of matrix
+* @param {number} G - the (0,1) element of matrix
+* @param {number} H - the (1,1) element of matrix
+* @param {Float64Array} out - output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively
+* @param {integer} strideOut - stride length for `out`
+* @param {NonNegativeInteger} offsetOut - starting index of `out`
+* @returns {Float64Array} output array
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var out = new Float64Array( 6 );
+*
+* dlasv2( 2.0, 3.0, 4.0, out, 1, 0 );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*/
+function dlasv2( F, G, H, out, strideOut, offsetOut ) {
+ var gasmal;
+ var ssmin;
+ var ssmax;
+ var tsign;
+ var PMAX;
+ var swap;
+ var clt;
+ var crt;
+ var csl;
+ var csr;
+ var idx;
+ var slt;
+ var snl;
+ var snr;
+ var srt;
+ var tmp;
+ var fa;
+ var ft;
+ var ga;
+ var gt;
+ var ha;
+ var ht;
+ var mm;
+ var tt;
+ var a;
+ var d;
+ var l;
+ var m;
+ var r;
+ var s;
+ var t;
+
+ ft = F;
+ fa = abs( ft );
+ ht = H;
+ ha = abs( ht );
+
+ /*
+ * PMAX points to the maximum absolute element of the matrix.
+ *
+ * PMAX = 1 if F is the largest in absolute value
+ * PMAX = 2 if G is the largest in absolute value
+ * PMAX = 3 if H is the largest in absolute value
+ */
+ PMAX = 1;
+ swap = ha > fa;
+ if ( swap ) {
+ PMAX = 3;
+ tmp = ft;
+ ft = ht;
+ ht = tmp;
+ tmp = fa;
+ fa = ha;
+ ha = tmp;
+
+ // Now FA >= HA
+ }
+ gt = G;
+ ga = abs( gt );
+ if ( ga === 0.0 ) {
+ // Diagonal matrix
+ ssmin = ha;
+ ssmax = fa;
+ clt = 1.0;
+ crt = 1.0;
+ slt = 0.0;
+ srt = 0.0;
+ } else {
+ gasmal = true;
+ if ( ga > fa ) {
+ PMAX = 2;
+ if ( ( fa / ga ) < EPS ) {
+ // Case of very large GA
+ gasmal = false;
+ ssmax = ga;
+ if ( ha > 1.0 ) {
+ ssmin = fa / ( ga / ha );
+ } else {
+ ssmin = ( fa / ga ) * ha;
+ }
+ clt = 1.0;
+ slt = ht / gt;
+ srt = 1.0;
+ crt = ft / gt;
+ }
+ }
+ if ( gasmal === true ) {
+ // Normal case
+ d = fa - ha;
+ if ( d === fa ) {
+ // Copes with infinite F or H
+ l = 1.0;
+ } else {
+ l = d / fa;
+ }
+
+ // Note that 0 < l < 1
+ m = gt / ft;
+
+ // Note that abs( m ) < 1 / macheps
+ t = 2.0 - l;
+
+ // Note that T > 1
+ mm = m * m;
+ tt = t * t;
+ s = sqrt( tt + mm );
+
+ // Note that 1 < S < 1 + 1 / macheps
+ if ( l === 0.0 ) {
+ r = abs( m );
+ } else {
+ r = sqrt( ( l * l ) + mm );
+ }
+
+ // Note that 0 < R < 1 + 1 / macheps
+ a = 0.5 * ( s + r);
+
+ // Note that 1 < A < 1 + abs( m )
+ ssmin = ha / a;
+ ssmax = fa * a;
+ if ( mm === 0.0 ) {
+ // Note that M is very tiny
+ if ( l === 0.0 ) {
+ t = copysign( 2.0, ft ) * copysign( 1.0, gt );
+ } else {
+ t = ( ( gt / copysign( d, ft ) ) + ( m / t ) );
+ }
+ } else {
+ t = ( ( m / ( s + t ) ) + ( m / ( r + l ) ) ) * ( 1.0 + a );
+ }
+ l = sqrt( ( t * t ) + 4.0 );
+ crt = 2.0 / l;
+ srt = t / l;
+ clt = ( crt + ( srt * m ) ) / a;
+ slt = ( ht / ft ) * ( srt / a );
+ }
+ }
+ if ( swap ) {
+ csl = srt;
+ snl = crt;
+ csr = slt;
+ snr = clt;
+ } else {
+ csl = clt;
+ snl = slt;
+ csr = crt;
+ snr = srt;
+ }
+
+ // Correct signs of ssmax and ssmin
+ if ( PMAX === 1 ) {
+ tsign = copysign( 1.0, csr ) * copysign( 1.0, csl ) * copysign( 1.0, F );
+ }
+ if ( PMAX === 2 ) {
+ tsign = copysign( 1.0, snr ) * copysign( 1.0, csl ) * copysign( 1.0, G );
+ }
+ if ( PMAX === 3 ) {
+ tsign = copysign( 1.0, snr ) * copysign( 1.0, snl ) * copysign( 1.0, H );
+ }
+ ssmax = copysign( ssmax, tsign );
+ ssmin = copysign( ssmin, tsign * copysign( 1.0, F ) * copysign( 1.0, H ) );
+
+ idx = offsetOut;
+ out[ idx ] = ssmin;
+ idx += strideOut;
+ out[ idx ] = ssmax;
+ idx += strideOut;
+ out[ idx ] = snr;
+ idx += strideOut;
+ out[ idx ] = csr;
+ idx += strideOut;
+ out[ idx ] = snl;
+ idx += strideOut;
+ out[ idx ] = csl;
+ return out;
+}
+
+
+// EXPORTS //
+
+module.exports = dlasv2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/dlasv2.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/dlasv2.js
new file mode 100644
index 000000000000..e19ebb9a93b8
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/dlasv2.js
@@ -0,0 +1,58 @@
+/**
+* @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 base = require( './base.js' );
+
+
+// MAIN //
+
+/**
+* Computes singular value decomposition of a 2x2 upper triangular matrix.
+*
+* ## Notes
+*
+* - `abs(SSMAX)` is the larger singular value.
+* - `abs(SSMIN)` is the smaller singular value.
+* - (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)`.
+*
+* @param {number} F - the (0,0) element of matrix
+* @param {number} G - the (0,1) element of matrix
+* @param {number} H - the (1,1) element of matrix
+* @param {Float64Array} out - output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively
+* @returns {Float64Array} output array
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var out = new Float64Array( 6 );
+*
+* dlasv2( 2.0, 3.0, 4.0, out );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*/
+function dlasv2( F, G, H, out ) {
+ return base( F, G, H, out, 1, 0 );
+}
+
+
+// EXPORTS //
+
+module.exports = dlasv2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/index.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/index.js
new file mode 100644
index 000000000000..60bd4a5a4263
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/index.js
@@ -0,0 +1,66 @@
+/**
+* @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 compute the singular value decomposition of a 2x2 upper triangular matrix.
+*
+* @module @stdlib/lapack/base/dlasv2
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+* var dlasv2 = require( '@stdlib/lapack/base/dlasv2' );
+*
+* var out = new Float64Array( 6 );
+*
+* dlasv2( 2.0, 3.0, 4.0, out );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+* var dlasv2 = require( '@stdlib/lapack/base/dlasv2' );
+*
+* var out = new Float64Array( 6 );
+*
+* dlasv2.ndarray( 2.0, 3.0, 4.0, out, 1, 0 );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*/
+
+// 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 dlasv2;
+var tmp = tryRequire( join( __dirname, './native.js' ) );
+if ( isError( tmp ) ) {
+ dlasv2 = main;
+} else {
+ dlasv2 = tmp;
+}
+
+
+// EXPORTS //
+
+module.exports = dlasv2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/main.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/main.js
new file mode 100644
index 000000000000..79357225880a
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/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 dlasv2 = require( './dlasv2.js' );
+var ndarray = require( './ndarray.js' );
+
+
+// MAIN //
+
+setReadOnly( dlasv2, 'ndarray', ndarray );
+
+
+// EXPORTS //
+
+module.exports = dlasv2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/ndarray.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/ndarray.js
new file mode 100644
index 000000000000..ec346777c7fa
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/lib/ndarray.js
@@ -0,0 +1,60 @@
+/**
+* @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 base = require( './base.js' );
+
+
+// MAIN //
+
+/**
+* Computes the singular value decomposition of a 2x2 upper triangular matrix using alternative indexing semantics.
+*
+* ## Notes
+*
+* - `abs(SSMAX)` is the larger singular value.
+* - `abs(SSMIN)` is the smaller singular value.
+* - (`CSL`,`SNL`) and (`CSR`,`SNR`) are the left and right singular vectors for `abs(SSMAX)`.
+*
+* @param {number} F - the (0,0) element of matrix
+* @param {number} G - the (0,1) element of matrix
+* @param {number} H - the (1,1) element of matrix
+* @param {Float64Array} out - output array containing `SSMIN`, `SSMAX`, `SNR`, `CSR`, `SNL`, and `CSL` respectively
+* @param {integer} strideOut - stride length for `out`
+* @param {NonNegativeInteger} offsetOut - starting index of `out`
+* @returns {Float64Array} output array
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var out = new Float64Array( 6 );
+*
+* out = dlasv2( 2.0, 3.0, 4.0, out, 1, 0 );
+* // out => [ ~1.5513, ~5.1569, ~0.9665, ~0.2567, ~0.7497, ~0.6618 ]
+*/
+function dlasv2( F, G, H, out, strideOut, offsetOut ) {
+ return base( F, G, H, out, strideOut, offsetOut );
+}
+
+
+// EXPORTS //
+
+module.exports = dlasv2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/package.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/package.json
new file mode 100644
index 000000000000..df04f78241e7
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/package.json
@@ -0,0 +1,72 @@
+{
+ "name": "@stdlib/lapack/base/dlasv2",
+ "version": "0.0.0",
+ "description": "Compute the singular value decomposition of a 2x2 upper triangular matrix.",
+ "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",
+ "svd",
+ "decomposition",
+ "dlasv2",
+ "exchange",
+ "permute",
+ "permutedims",
+ "linear",
+ "algebra",
+ "subroutines",
+ "array",
+ "ndarray",
+ "float64",
+ "double",
+ "float64array"
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/decreasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/decreasing_order_elements.json
new file mode 100644
index 000000000000..35409522867c
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/decreasing_order_elements.json
@@ -0,0 +1,15 @@
+{
+ "F": 4.0,
+ "G": 3.0,
+ "H": 2.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 1.5513263285176897,
+ 5.1568776039816795,
+ 0.6618025632357402,
+ 0.7496781758158659,
+ 0.25666793515702424,
+ 0.9664996487646696
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/diagonal_matrix.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/diagonal_matrix.json
new file mode 100644
index 000000000000..814517855e8e
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/diagonal_matrix.json
@@ -0,0 +1,15 @@
+{
+ "F": 2.0,
+ "G": 0.0,
+ "H": -12.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 2.0,
+ -12.0,
+ 1.0,
+ 0.0,
+ 1.0,
+ 0.0
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/g_gt_f_and_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/g_gt_f_and_h.json
new file mode 100644
index 000000000000..44ff548c100a
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/g_gt_f_and_h.json
@@ -0,0 +1,15 @@
+{
+ "F": 4.0,
+ "G": 5.0,
+ "H": 2.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 1.2125424343860751,
+ 6.597707241520579,
+ 0.8090383403577232,
+ 0.5877558709457702,
+ 0.2452483296822559,
+ 0.9694602914962858
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/increasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/increasing_order_elements.json
new file mode 100644
index 000000000000..9925344bda65
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/increasing_order_elements.json
@@ -0,0 +1,15 @@
+{
+ "F": 2.0,
+ "G": 3.0,
+ "H": 4.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 1.5513263285176897,
+ 5.1568776039816795,
+ 0.9664996487646696,
+ 0.25666793515702424,
+ 0.7496781758158659,
+ 0.6618025632357402
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_g_and_h_gt_1.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_g_and_h_gt_1.json
new file mode 100644
index 000000000000..0bff731a8a68
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_g_and_h_gt_1.json
@@ -0,0 +1,15 @@
+{
+ "F": 2.0,
+ "G": 1.0E+20,
+ "H": 3.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 6.0E-20,
+ 1.0E+20,
+ 1.0,
+ 2.0E-20,
+ 3.0E-20,
+ 1.0
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_g_diagonal.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_g_diagonal.json
new file mode 100644
index 000000000000..ca5f429368ff
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_g_diagonal.json
@@ -0,0 +1,15 @@
+{
+ "F": 2.0,
+ "G": 1.0E+20,
+ "H": 1.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 2.0E-20,
+ 1.0E+20,
+ 1.0,
+ 2.0E-20,
+ 1.0E-20,
+ 1.0
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/decreasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/decreasing_order_elements.json
new file mode 100644
index 000000000000..b321ae33bcf7
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/decreasing_order_elements.json
@@ -0,0 +1,21 @@
+{
+ "F": 4,
+ "G": 3,
+ "H": 2,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 1.5513263285176897,
+ 9999,
+ 5.1568776039816795,
+ 9999,
+ 0.6618025632357402,
+ 9999,
+ 0.7496781758158659,
+ 9999,
+ 0.25666793515702424,
+ 9999,
+ 0.9664996487646696,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/diagonal_matrix.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/diagonal_matrix.json
new file mode 100644
index 000000000000..f299d333f323
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/diagonal_matrix.json
@@ -0,0 +1,21 @@
+{
+ "F": 2,
+ "G": 0,
+ "H": -12,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 2,
+ 9999,
+ -12,
+ 9999,
+ 1,
+ 9999,
+ 0,
+ 9999,
+ 1,
+ 9999,
+ 0,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/g_gt_f_and_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/g_gt_f_and_h.json
new file mode 100644
index 000000000000..c3c181a0fdf9
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/g_gt_f_and_h.json
@@ -0,0 +1,21 @@
+{
+ "F": 4,
+ "G": 5,
+ "H": 2,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 1.2125424343860751,
+ 9999,
+ 6.597707241520579,
+ 9999,
+ 0.8090383403577232,
+ 9999,
+ 0.5877558709457702,
+ 9999,
+ 0.2452483296822559,
+ 9999,
+ 0.9694602914962858,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/increasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/increasing_order_elements.json
new file mode 100644
index 000000000000..f8ebcd12a5b4
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/increasing_order_elements.json
@@ -0,0 +1,21 @@
+{
+ "F": 2,
+ "G": 3,
+ "H": 4,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 1.5513263285176897,
+ 9999,
+ 5.1568776039816795,
+ 9999,
+ 0.9664996487646696,
+ 9999,
+ 0.25666793515702424,
+ 9999,
+ 0.7496781758158659,
+ 9999,
+ 0.6618025632357402,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/large_g_and_h_gt_1.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/large_g_and_h_gt_1.json
new file mode 100644
index 000000000000..2c45c3d30921
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/large_g_and_h_gt_1.json
@@ -0,0 +1,21 @@
+{
+ "F": 2,
+ "G": 100000000000000000000,
+ "H": 3,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 6e-20,
+ 9999,
+ 100000000000000000000,
+ 9999,
+ 1,
+ 9999,
+ 2e-20,
+ 9999,
+ 3e-20,
+ 9999,
+ 1,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/large_g_diagonal.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/large_g_diagonal.json
new file mode 100644
index 000000000000..192510267c14
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/large_g_diagonal.json
@@ -0,0 +1,21 @@
+{
+ "F": 2,
+ "G": 100000000000000000000,
+ "H": 1,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 2e-20,
+ 9999,
+ 100000000000000000000,
+ 9999,
+ 1,
+ 9999,
+ 2e-20,
+ 9999,
+ 1e-20,
+ 9999,
+ 1,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/small_g_and_f_eq_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/small_g_and_f_eq_h.json
new file mode 100644
index 000000000000..13e8e65a7ec3
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/large_strides/small_g_and_f_eq_h.json
@@ -0,0 +1,21 @@
+{
+ "F": 3,
+ "G": 5e-324,
+ "H": 3,
+ "strideOut": 2,
+ "offsetOut": 0,
+ "out": [
+ 3,
+ 9999,
+ 3,
+ 9999,
+ 0.7071067811865475,
+ 9999,
+ 0.7071067811865475,
+ 9999,
+ 0.7071067811865475,
+ 9999,
+ 0.7071067811865475,
+ 9999
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/decreasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/decreasing_order_elements.json
new file mode 100644
index 000000000000..51d6510a1774
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/decreasing_order_elements.json
@@ -0,0 +1,15 @@
+{
+ "F": 4,
+ "G": 3,
+ "H": 2,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 0.9664996487646696,
+ 0.25666793515702424,
+ 0.7496781758158659,
+ 0.6618025632357402,
+ 5.1568776039816795,
+ 1.5513263285176897
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/diagonal_matrix.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/diagonal_matrix.json
new file mode 100644
index 000000000000..9e7b8cbed009
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/diagonal_matrix.json
@@ -0,0 +1,15 @@
+{
+ "F": 2,
+ "G": 0,
+ "H": -12,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 0,
+ 1,
+ 0,
+ 1,
+ -12,
+ 2
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/g_gt_f_and_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/g_gt_f_and_h.json
new file mode 100644
index 000000000000..1c99e07a1c96
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/g_gt_f_and_h.json
@@ -0,0 +1,15 @@
+{
+ "F": 4.0,
+ "G": 5.0,
+ "H": 2.0,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 0.9694602914962858,
+ 0.2452483296822559,
+ 0.5877558709457702,
+ 0.8090383403577232,
+ 6.597707241520579,
+ 1.2125424343860751
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/increasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/increasing_order_elements.json
new file mode 100644
index 000000000000..b16c6019c733
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/increasing_order_elements.json
@@ -0,0 +1,15 @@
+{
+ "F": 2,
+ "G": 3,
+ "H": 4,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 0.6618025632357402,
+ 0.7496781758158659,
+ 0.25666793515702424,
+ 0.9664996487646696,
+ 5.1568776039816795,
+ 1.5513263285176897
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/large_g_and_h_gt_1.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/large_g_and_h_gt_1.json
new file mode 100644
index 000000000000..77c59e41a96d
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/large_g_and_h_gt_1.json
@@ -0,0 +1,15 @@
+{
+ "F": 2,
+ "G": 100000000000000000000,
+ "H": 3,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 1,
+ 3e-20,
+ 2e-20,
+ 1,
+ 100000000000000000000,
+ 6e-20
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/large_g_diagonal.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/large_g_diagonal.json
new file mode 100644
index 000000000000..1a1bde69d0b0
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/large_g_diagonal.json
@@ -0,0 +1,15 @@
+{
+ "F": 2,
+ "G": 100000000000000000000,
+ "H": 1,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 1,
+ 1e-20,
+ 2e-20,
+ 1,
+ 100000000000000000000,
+ 2e-20
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/small_g_and_f_eq_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/small_g_and_f_eq_h.json
new file mode 100644
index 000000000000..add225b46762
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/negative_strides/small_g_and_f_eq_h.json
@@ -0,0 +1,15 @@
+{
+ "F": 3,
+ "G": 5e-324,
+ "H": 3,
+ "strideOut": -1,
+ "offsetOut": 5,
+ "out": [
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 3,
+ 3
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/decreasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/decreasing_order_elements.json
new file mode 100644
index 000000000000..31f88d89ed1e
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/decreasing_order_elements.json
@@ -0,0 +1,16 @@
+{
+ "F": 4,
+ "G": 3,
+ "H": 2,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 1.5513263285176897,
+ 5.1568776039816795,
+ 0.6618025632357402,
+ 0.7496781758158659,
+ 0.25666793515702424,
+ 0.9664996487646696
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/diagonal_matrix.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/diagonal_matrix.json
new file mode 100644
index 000000000000..3c6fe9fe9c62
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/diagonal_matrix.json
@@ -0,0 +1,16 @@
+{
+ "F": 2,
+ "G": 0,
+ "H": -12,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 2,
+ -12,
+ 1,
+ 0,
+ 1,
+ 0
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/g_gt_f_and_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/g_gt_f_and_h.json
new file mode 100644
index 000000000000..3d432ed52fc5
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/g_gt_f_and_h.json
@@ -0,0 +1,16 @@
+{
+ "F": 4,
+ "G": 5,
+ "H": 2,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 1.2125424343860751,
+ 6.597707241520579,
+ 0.8090383403577232,
+ 0.5877558709457702,
+ 0.2452483296822559,
+ 0.9694602914962858
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/increasing_order_elements.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/increasing_order_elements.json
new file mode 100644
index 000000000000..b864705ad0f2
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/increasing_order_elements.json
@@ -0,0 +1,16 @@
+{
+ "F": 2,
+ "G": 3,
+ "H": 4,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 1.5513263285176897,
+ 5.1568776039816795,
+ 0.9664996487646696,
+ 0.25666793515702424,
+ 0.7496781758158659,
+ 0.6618025632357402
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/large_g_and_h_gt_1.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/large_g_and_h_gt_1.json
new file mode 100644
index 000000000000..c1ac1cf0d8ca
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/large_g_and_h_gt_1.json
@@ -0,0 +1,16 @@
+{
+ "F": 2,
+ "G": 100000000000000000000,
+ "H": 3,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 6e-20,
+ 100000000000000000000,
+ 1,
+ 2e-20,
+ 3e-20,
+ 1
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/large_g_diagonal.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/large_g_diagonal.json
new file mode 100644
index 000000000000..c5830d834549
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/large_g_diagonal.json
@@ -0,0 +1,16 @@
+{
+ "F": 2,
+ "G": 100000000000000000000,
+ "H": 1,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 2e-20,
+ 100000000000000000000,
+ 1,
+ 2e-20,
+ 1e-20,
+ 1
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/small_g_and_f_eq_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/small_g_and_f_eq_h.json
new file mode 100644
index 000000000000..a1cdea93f161
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/offsets/small_g_and_f_eq_h.json
@@ -0,0 +1,16 @@
+{
+ "F": 3,
+ "G": 5e-324,
+ "H": 3,
+ "strideOut": 1,
+ "offsetOut": 1,
+ "out": [
+ 9999,
+ 3,
+ 3,
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 0.7071067811865475
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/small_g_and_f_eq_h.json b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/small_g_and_f_eq_h.json
new file mode 100644
index 000000000000..1cd2c29cbca2
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/fixtures/small_g_and_f_eq_h.json
@@ -0,0 +1,15 @@
+{
+ "F": 3.0,
+ "G": 5.0E-324,
+ "H": 3.0,
+ "strideOut": 1,
+ "offsetOut": 0,
+ "out": [
+ 3.0,
+ 3.0,
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 0.7071067811865475,
+ 0.7071067811865475
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.dlasv2.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.dlasv2.js
new file mode 100644
index 000000000000..70080489d985
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.dlasv2.js
@@ -0,0 +1,169 @@
+/**
+* @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 isAlmostSameValueFloat64Array = require( '@stdlib/assert/is-almost-same-value-float64array' );
+var PINF = require( '@stdlib/constants/float64/pinf' );
+var dlasv2 = require( './../lib/dlasv2.js' );
+
+
+// FIXTURES //
+
+// NOTE: Expected outputs in fixtures were generated using the reference LAPACK Fortran implementation of DLASV2.
+var INC_ORD_ELE = require( './fixtures/increasing_order_elements.json' );
+var DEC_ORD_ELE = require( './fixtures/decreasing_order_elements.json' );
+var LARGE_G = require( './fixtures/large_g_diagonal.json' );
+var LARGE_G_AND_H_GT_1 = require( './fixtures/large_g_and_h_gt_1.json' );
+var SMALL_G_AND_F_EQ_H = require( './fixtures/small_g_and_f_eq_h.json' );
+var DIA_MAT = require( './fixtures/diagonal_matrix.json' );
+var G_GT_F_AND_H = require( './fixtures/g_gt_f_and_h.json' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+ t.ok( true, __filename );
+ t.strictEqual( typeof dlasv2, 'function', 'main export is a function' );
+ t.end();
+});
+
+tape( 'the function has an arity of 4', function test( t ) {
+ t.strictEqual( dlasv2.length, 4, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in increasing order', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = INC_ORD_ELE;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in decreasing order', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = DEC_ORD_ELE;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with parameter `G` greater than `F` and `H`', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = G_GT_F_AND_H;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = DIA_MAT;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G`', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LARGE_G;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` when `H` is greater than 1', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LARGE_G_AND_H_GT_1;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with an infinite `F`', function test( t ) {
+ var expected;
+ var out;
+
+ expected = new Float64Array( [ 1.0, PINF, 0.0, 1.0, 0.0, 1.0 ] );
+ out = new Float64Array( 6 );
+ out = dlasv2( PINF, 2.0, 1.0, out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with tiny off-diagonal entries when `F` and `H` are equal', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = SMALL_G_AND_F_EQ_H;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.js
new file mode 100644
index 000000000000..82644bcbfe4d
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/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 dlasv2 = 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 dlasv2, '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 dlasv2.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 dlasv2 = proxyquire( './../lib', {
+ '@stdlib/utils/try-require': tryRequire
+ });
+
+ t.strictEqual( dlasv2, 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 dlasv2;
+ var main;
+
+ main = require( './../lib/dlasv2.js' );
+
+ dlasv2 = proxyquire( './../lib', {
+ '@stdlib/utils/try-require': tryRequire
+ });
+
+ t.strictEqual( dlasv2, main, 'returns expected value' );
+ t.end();
+
+ function tryRequire() {
+ return new Error( 'Cannot find module' );
+ }
+});
diff --git a/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.ndarray.js b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.ndarray.js
new file mode 100644
index 000000000000..a021e663049a
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dlasv2/test/test.ndarray.js
@@ -0,0 +1,493 @@
+/**
+* @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 isAlmostSameValueFloat64Array = require( '@stdlib/assert/is-almost-same-value-float64array' );
+var dlasv2 = require( './../lib/ndarray.js' );
+
+
+// FIXTURES //
+
+// NOTE: Expected outputs in fixtures were generated using the reference LAPACK Fortran implementation of DLASV2.
+var INC_ORD_ELE = require( './fixtures/increasing_order_elements.json' );
+var DEC_ORD_ELE = require( './fixtures/decreasing_order_elements.json' );
+var LARGE_G = require( './fixtures/large_g_diagonal.json' );
+var LARGE_G_AND_H_GT_1 = require( './fixtures/large_g_and_h_gt_1.json' );
+var SMALL_G_AND_F_EQ_H = require( './fixtures/small_g_and_f_eq_h.json' );
+var DIA_MAT = require( './fixtures/diagonal_matrix.json' );
+var G_GT_F_AND_H = require( './fixtures/g_gt_f_and_h.json' );
+var LAR_STR_INC_ORD_ELE = require( './fixtures/large_strides/increasing_order_elements.json' );
+var LAR_STR_DEC_ORD_ELE = require( './fixtures/large_strides/decreasing_order_elements.json' );
+var LAR_STR_LARGE_G = require( './fixtures/large_strides/large_g_diagonal.json' );
+var LAR_STR_LARGE_G_AND_H_GT_1 = require( './fixtures/large_strides/large_g_and_h_gt_1.json' );
+var LAR_STR_SMALL_G_AND_F_EQ_H = require( './fixtures/large_strides/small_g_and_f_eq_h.json' );
+var LAR_STR_DIA_MAT = require( './fixtures/large_strides/diagonal_matrix.json' );
+var LAR_STR_G_GT_F_AND_H = require( './fixtures/large_strides/g_gt_f_and_h.json' );
+var NEG_STR_INC_ORD_ELE = require( './fixtures/negative_strides/increasing_order_elements.json' );
+var NEG_STR_DEC_ORD_ELE = require( './fixtures/negative_strides/decreasing_order_elements.json' );
+var NEG_STR_LARGE_G = require( './fixtures/negative_strides/large_g_diagonal.json' );
+var NEG_STR_LARGE_G_AND_H_GT_1 = require( './fixtures/negative_strides/large_g_and_h_gt_1.json' );
+var NEG_STR_SMALL_G_AND_F_EQ_H = require( './fixtures/negative_strides/small_g_and_f_eq_h.json' );
+var NEG_STR_DIA_MAT = require( './fixtures/negative_strides/diagonal_matrix.json' );
+var NEG_STR_G_GT_F_AND_H = require( './fixtures/negative_strides/g_gt_f_and_h.json' );
+var OFF_INC_ORD_ELE = require( './fixtures/offsets/increasing_order_elements.json' );
+var OFF_DEC_ORD_ELE = require( './fixtures/offsets/decreasing_order_elements.json' );
+var OFF_LARGE_G = require( './fixtures/offsets/large_g_diagonal.json' );
+var OFF_LARGE_G_AND_H_GT_1 = require( './fixtures/offsets/large_g_and_h_gt_1.json' );
+var OFF_SMALL_G_AND_F_EQ_H = require( './fixtures/offsets/small_g_and_f_eq_h.json' );
+var OFF_DIA_MAT = require( './fixtures/offsets/diagonal_matrix.json' );
+var OFF_G_GT_F_AND_H = require( './fixtures/offsets/g_gt_f_and_h.json' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+ t.ok( true, __filename );
+ t.strictEqual( typeof dlasv2, 'function', 'main export is a function' );
+ t.end();
+});
+
+tape( 'the function has an arity of 6', function test( t ) {
+ t.strictEqual( dlasv2.length, 6, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in increasing order', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = INC_ORD_ELE;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in decreasing order', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = DEC_ORD_ELE;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with parameter `G` greater than `F` and `H`', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = G_GT_F_AND_H;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = DIA_MAT;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G`', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LARGE_G;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` when `H` is greater than 1', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LARGE_G_AND_H_GT_1;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with tiny off-diagonal entries when `F` and `H` are equal', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = SMALL_G_AND_F_EQ_H;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in increasing order (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_INC_ORD_ELE;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in decreasing order (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_DEC_ORD_ELE;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with parameter `G` greater than `F` and `H` (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_G_GT_F_AND_H;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_DIA_MAT;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_LARGE_G;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` when `H` is greater than 1 (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_LARGE_G_AND_H_GT_1;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with tiny off-diagonal entries when `F` and `H` are equal (large stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = LAR_STR_SMALL_G_AND_F_EQ_H;
+
+ out = new Float64Array( [ 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999, 0, 9999 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in increasing order (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_INC_ORD_ELE;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in decreasing order (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_DEC_ORD_ELE;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with parameter `G` greater than `F` and `H` (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_G_GT_F_AND_H;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_DIA_MAT;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_LARGE_G;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` when `H` is greater than 1 (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_LARGE_G_AND_H_GT_1;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with tiny off-diagonal entries when `F` and `H` are equal (negative stride)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = NEG_STR_SMALL_G_AND_F_EQ_H;
+
+ out = new Float64Array( 6 );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in increasing order (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_INC_ORD_ELE;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with elements in decreasing order (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_DEC_ORD_ELE;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a triangular matrix with parameter `G` greater than `F` and `H` (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_G_GT_F_AND_H;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_DIA_MAT;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_LARGE_G;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with a very large `G` when `H` is greater than 1 (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_LARGE_G_AND_H_GT_1;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function computes the singular value decomposition of a diagonal matrix with tiny off-diagonal entries when `F` and `H` are equal (offset)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+
+ data = OFF_SMALL_G_AND_F_EQ_H;
+
+ out = new Float64Array( [ 9999, 0, 0, 0, 0, 0, 0 ] );
+
+ dlasv2( data.F, data.G, data.H, out, data.strideOut, data.offsetOut );
+ expected = new Float64Array( data.out );
+ t.strictEqual( isAlmostSameValueFloat64Array( out, expected, 1 ), true, 'returns expected value' );
+ t.end();
+});