diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/README.md b/lib/node_modules/@stdlib/lapack/base/dptts2/README.md
new file mode 100644
index 000000000000..5af8f11738aa
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/README.md
@@ -0,0 +1,251 @@
+
+
+# dptts2
+
+> Solve a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by [`dpttrf`][@stdlib/lapack/base/dpttrf].
+
+
+
+## Usage
+
+```javascript
+var dptts2 = require( '@stdlib/lapack/base/dptts2' );
+```
+
+#### dptts2( order, N, nrhs, D, E, B, LDB )
+
+Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by [`dpttrf`][@stdlib/lapack/base/dpttrf].
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+dptts2( 'column-major', 3, 1, D, E, B, 3 );
+// B => [ ~0.204, ~0.184, ~0.439 ]
+```
+
+The function has the following parameters:
+
+- **order**: storage layout of `B`.
+- **N**: order of the tridiagonal matrix `A`.
+- **nrhs**: number of right-hand sides (i.e., the number of columns of the matrix `B`).
+- **D**: the `N` diagonal elements of the diagonal matrix `D` from the `L * D * L^T` factorization of `A` as a [`Float64Array`][mdn-float64array].
+- **E**: the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the `L * D * L^T` factorization of `A` as a [`Float64Array`][mdn-float64array].
+- **B**: input matrix as a [`Float64Array`][mdn-float64array]. On entry, the right-hand side vectors `B`. On exit, the solution vectors `X`.
+- **LDB**: leading dimension of `B`.
+
+The input matrix `D` and `E` are typically computed by [`dpttrf`][@stdlib/lapack/base/dpttrf].
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+var dpttrf = require( '@stdlib/lapack/base/dpttrf' );
+
+var D = new Float64Array( [ 4.0, 5.0, 6.0 ] );
+var E = new Float64Array( [ 1.0, 2.0 ] );
+
+dpttrf( 3, D, E );
+
+var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+dptts2( 'column-major', 3, 1, D, E, B, 3 );
+// B => [ ~0.204, ~0.184, ~0.439 ]
+```
+
+#### dptts2.ndarray( N, nrhs, D, sd, od, E, se, oe, B, sb1, sb2, ob )
+
+Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by [`dpttrf`][@stdlib/lapack/base/dpttrf] and alternative indexing semantics.
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 );
+// B => [ ~0.204, ~0.184, ~0.439 ]
+```
+
+The function has the following parameters:
+
+- **sd**: stride length for `D`.
+- **od**: starting index for `D`.
+- **se**: stride length for `E`.
+- **oe**: starting index for `E`.
+- **sb1**: stride of the first dimension of `B`.
+- **sb2**: stride of the second dimension of `B`.
+- **ob**: starting index for `B`.
+
+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.
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+
+var D = new Float64Array( [ 0.0, 4.0, 4.75, 5.157894736842105 ] );
+var E = new Float64Array( [ 0.0, 0.25, 0.42105263157894735 ] );
+var B = new Float64Array( [ 0.0, 1.0, 2.0, 3.0 ] );
+
+dptts2.ndarray( 3, 1, D, 1, 1, E, 1, 1, B, 1, 3, 1 );
+// B => [ 0.0, ~0.204, ~0.184, ~0.439 ]
+```
+
+
+
+
+
+
+
+## Notes
+
+- `dptts2()` corresponds to the [LAPACK][LAPACK] routine [`dptts2`][lapack-dptts2].
+
+
+
+
+
+
+
+## Examples
+
+
+
+```javascript
+var Float64Array = require( '@stdlib/array/float64' );
+var dpttrf = require( '@stdlib/lapack/base/dpttrf' );
+var dptts2 = require( '@stdlib/lapack/base/dptts2' );
+
+// Define a symmetric positive definite matrix `A` via its diagonal and subdiagonal:
+var D = new Float64Array( [ 4.0, 5.0, 6.0, 7.0 ] );
+var E = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+// Compute the `L * D * L^T` factorization of `A` (overwrites `D` and `E`):
+dpttrf( D.length, D, E );
+
+// Define the right-hand side matrix `B` (column-major, two right-hand sides):
+var B = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 2.0, 1.0, 0.0, 1.0 ] );
+
+// Solve `A * X = B`, overwriting `B` with the solution `X`:
+dptts2( 'column-major', 4, 2, D, E, B, 4 );
+console.log( B );
+```
+
+
+
+
+
+
+
+* * *
+
+
+
+## C APIs
+
+
+
+
+
+
+
+
+
+
+
+### Usage
+
+```c
+TODO
+```
+
+#### TODO
+
+TODO.
+
+```c
+TODO
+```
+
+TODO
+
+```c
+TODO
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+### Examples
+
+```c
+TODO
+```
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+[lapack]: https://www.netlib.org/lapack/explore-html/
+
+[lapack-dptts2]: https://netlib.org/lapack/explore-html/d5/d5d/group__ptts2_ga35fdfa6109e8f9cbfbde271814bf0b27.html#ga35fdfa6109e8f9cbfbde271814bf0b27
+
+[mdn-float64array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/Float64Array
+
+[mdn-typed-array]: https://developer.mozilla.org/en-US/docs/Web/JavaScript/Reference/Global_Objects/TypedArray
+
+[@stdlib/lapack/base/dpttrf]: https://github.com/stdlib-js/stdlib/tree/develop/lib/node_modules/%40stdlib/lapack/base/dpttrf
+
+
+
+
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/benchmark/benchmark.js b/lib/node_modules/@stdlib/lapack/base/dptts2/benchmark/benchmark.js
new file mode 100644
index 000000000000..b008322af9a0
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/benchmark/benchmark.js
@@ -0,0 +1,105 @@
+/**
+* @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 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 dptts2 = require( './../lib/dptts2.js' );
+
+
+// VARIABLES //
+
+var options = {
+ 'dtype': 'float64'
+};
+
+
+// FUNCTIONS //
+
+/**
+* Creates a benchmark function.
+*
+* @private
+* @param {PositiveInteger} N - number of elements along each dimension
+* @returns {Function} benchmark function
+*/
+function createBenchmark( N ) {
+ var D = uniform( N, 1.0, 5.0, options );
+ var E = uniform( N-1, 0.0, 1.0, options );
+ var B = uniform( N, -1.0, 1.0, options );
+ return benchmark;
+
+ /**
+ * Benchmark function.
+ *
+ * @private
+ * @param {Benchmark} b - benchmark instance
+ */
+ function benchmark( b ) {
+ var out;
+ var i;
+
+ b.tic();
+ for ( i = 0; i < b.iterations; i++ ) {
+ out = dptts2( 'column-major', N, 1, D, E, B, N );
+ if ( isnan( out[ 0 ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ }
+ b.toc();
+ if ( isnan( out[ 0 ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ b.pass( 'benchmark finished' );
+ b.end();
+ }
+}
+
+
+// MAIN //
+
+/**
+* Main execution sequence.
+*
+* @private
+*/
+function main() {
+ var min;
+ var max;
+ var N;
+ var f;
+ var i;
+
+ min = 1; // 10^min
+ max = 6; // 10^max
+
+ for ( i = min; i <= max; i++ ) {
+ N = pow( 10, i );
+ f = createBenchmark( N );
+ bench( format( '%s:size=%d', pkg, N ), f );
+ }
+}
+
+main();
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/benchmark/benchmark.ndarray.js b/lib/node_modules/@stdlib/lapack/base/dptts2/benchmark/benchmark.ndarray.js
new file mode 100644
index 000000000000..2fc1031a948e
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/benchmark/benchmark.ndarray.js
@@ -0,0 +1,105 @@
+/**
+* @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 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 dptts2 = require( './../lib/ndarray.js' );
+
+
+// VARIABLES //
+
+var options = {
+ 'dtype': 'float64'
+};
+
+
+// FUNCTIONS //
+
+/**
+* Creates a benchmark function.
+*
+* @private
+* @param {PositiveInteger} N - number of elements along each dimension
+* @returns {Function} benchmark function
+*/
+function createBenchmark( N ) {
+ var D = uniform( N, 1.0, 5.0, options );
+ var E = uniform( N-1, 0.0, 1.0, options );
+ var B = uniform( N, -1.0, 1.0, options );
+ return benchmark;
+
+ /**
+ * Benchmark function.
+ *
+ * @private
+ * @param {Benchmark} b - benchmark instance
+ */
+ function benchmark( b ) {
+ var out;
+ var i;
+
+ b.tic();
+ for ( i = 0; i < b.iterations; i++ ) {
+ out = dptts2( N, 1, D, 1, 0, E, 1, 0, B, 1, N, 0 );
+ if ( isnan( out[ 0 ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ }
+ b.toc();
+ if ( isnan( out[ 0 ] ) ) {
+ b.fail( 'should not return NaN' );
+ }
+ b.pass( 'benchmark finished' );
+ b.end();
+ }
+}
+
+
+// MAIN //
+
+/**
+* Main execution sequence.
+*
+* @private
+*/
+function main() {
+ var min;
+ var max;
+ var N;
+ var f;
+ var i;
+
+ min = 1; // 10^min
+ max = 6; // 10^max
+
+ for ( i = min; i <= max; i++ ) {
+ N = pow( 10, i );
+ f = createBenchmark( N );
+ bench( format( '%s:ndarray:size=%d', pkg, N ), f );
+ }
+}
+
+main();
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/docs/repl.txt b/lib/node_modules/@stdlib/lapack/base/dptts2/docs/repl.txt
new file mode 100644
index 000000000000..68d03402f347
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/docs/repl.txt
@@ -0,0 +1,124 @@
+
+{{alias}}( order, N, nrhs, D, E, B, LDB )
+ Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T`
+ factorization of `A` computed by `dpttrf`.
+
+ `D` is the diagonal matrix specified in the vector `D`, `L` is the unit
+ bidiagonal matrix whose subdiagonal is specified in the vector `E`, and `X`
+ and `B` are `N`-by-`nrhs` matrices.
+
+ Indexing is relative to the first index. To introduce an offset, use typed
+ array views.
+
+ The function overwrites the input matrix `B` with the solution matrix `X`.
+
+ Parameters
+ ----------
+ order: string
+ Row-major (C-style) or column-major (Fortran-style) order of `B`. Must
+ be either 'row-major' or 'column-major'.
+
+ N: integer
+ Order of the tridiagonal matrix `A`.
+
+ nrhs: integer
+ Number of right-hand sides (i.e., the number of columns of `B`).
+
+ D: Float64Array
+ The `N` diagonal elements of the diagonal matrix `D` from the
+ `L * D * L^T` factorization of `A`.
+
+ E: Float64Array
+ The `N-1` subdiagonal elements of the unit bidiagonal factor `L` from
+ the `L * D * L^T` factorization of `A`.
+
+ B: Float64Array
+ Input matrix. On entry, the right-hand side vectors `B` for the system
+ of linear equations. On exit, the solution vectors `X`.
+
+ LDB: integer
+ Stride of the first dimension of `B` (a.k.a., leading dimension of the
+ matrix `B`).
+
+ Returns
+ -------
+ B: Float64Array
+ The solution matrix `X` (stored in `B`).
+
+ Examples
+ --------
+ > var D = new {{alias:@stdlib/array/float64}}( [ 4.0, 4.75, 5.157894736842105 ] );
+ > var E = new {{alias:@stdlib/array/float64}}( [ 0.25, 0.42105263157894735 ] );
+ > var B = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 3.0 ] );
+ > {{alias}}( 'column-major', 3, 1, D, E, B, 3 )
+ [ ~0.204, ~0.184, ~0.439 ]
+
+
+{{alias}}.ndarray( N, nrhs, D, sd, od, E, se, oe, B, sb1, sb2, ob )
+ Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T`
+ factorization of `A` computed by `dpttrf` and alternative indexing
+ semantics.
+
+ While typed array views mandate a view offset based on the underlying
+ buffer, the offset parameters support indexing semantics based on starting
+ indices.
+
+ The function overwrites the input matrix `B` with the solution matrix `X`.
+
+ Parameters
+ ----------
+ N: integer
+ Order of the tridiagonal matrix `A`.
+
+ nrhs: integer
+ Number of right-hand sides (i.e., the number of columns of `B`).
+
+ D: Float64Array
+ The `N` diagonal elements of the diagonal matrix `D` from the
+ `L * D * L^T` factorization of `A`.
+
+ sd: integer
+ Stride length for `D`.
+
+ od: integer
+ Starting index for `D`.
+
+ E: Float64Array
+ The `N-1` subdiagonal elements of the unit bidiagonal factor `L` from
+ the `L * D * L^T` factorization of `A`.
+
+ se: integer
+ Stride length for `E`.
+
+ oe: integer
+ Starting index for `E`.
+
+ B: Float64Array
+ Input matrix. On entry, the right-hand side vectors `B` for the system
+ of linear equations. On exit, the solution vectors `X`.
+
+ sb1: integer
+ Stride of the first dimension of `B`.
+
+ sb2: integer
+ Stride of the second dimension of `B`.
+
+ ob: integer
+ Starting index for `B`.
+
+ Returns
+ -------
+ B: Float64Array
+ The solution matrix `X` (stored in `B`).
+
+ Examples
+ --------
+ > var D = new {{alias:@stdlib/array/float64}}( [ 4.0, 4.75, 5.157894736842105 ] );
+ > var E = new {{alias:@stdlib/array/float64}}( [ 0.25, 0.42105263157894735 ] );
+ > var B = new {{alias:@stdlib/array/float64}}( [ 1.0, 2.0, 3.0 ] );
+ > {{alias}}.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 )
+ [ ~0.204, ~0.184, ~0.439 ]
+
+ See Also
+ --------
+
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/docs/types/index.d.ts b/lib/node_modules/@stdlib/lapack/base/dptts2/docs/types/index.d.ts
new file mode 100644
index 000000000000..64e72519335f
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/docs/types/index.d.ts
@@ -0,0 +1,120 @@
+/*
+* @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
+
+///
+
+import { Layout } from '@stdlib/types/blas';
+
+/**
+* Interface describing `dptts2`.
+*/
+interface Routine {
+ /**
+ * Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf`.
+ *
+ * @param order - storage layout of `B`
+ * @param N - order of the tridiagonal matrix `A`
+ * @param nrhs - number of right-hand sides (i.e., the number of columns of the matrix `B`)
+ * @param D - the `N` diagonal elements of the diagonal matrix `D` from the factorization of `A`
+ * @param E - the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the factorization of `A`
+ * @param B - input matrix (overwritten by the solution matrix `X`)
+ * @param LDB - leading dimension of `B`
+ * @returns solution matrix `X` (stored in `B`)
+ *
+ * @example
+ * var Float64Array = require( '@stdlib/array/float64' );
+ *
+ * var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ * var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ * var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+ *
+ * dptts2( 'column-major', 3, 1, D, E, B, 3 );
+ * // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+ */
+ ( order: Layout, N: number, nrhs: number, D: Float64Array, E: Float64Array, B: Float64Array, LDB: number ): Float64Array;
+
+ /**
+ * Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf` and alternative indexing semantics.
+ *
+ * @param N - order of the tridiagonal matrix `A`
+ * @param nrhs - number of right-hand sides (i.e., the number of columns of the matrix `B`)
+ * @param D - the `N` diagonal elements of the diagonal matrix `D` from the factorization of `A`
+ * @param sd - stride length for `D`
+ * @param od - starting index for `D`
+ * @param E - the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the factorization of `A`
+ * @param se - stride length for `E`
+ * @param oe - starting index for `E`
+ * @param B - input matrix (overwritten by the solution matrix `X`)
+ * @param sb1 - stride of the first dimension of `B`
+ * @param sb2 - stride of the second dimension of `B`
+ * @param ob - starting index for `B`
+ * @returns solution matrix `X` (stored in `B`)
+ *
+ * @example
+ * var Float64Array = require( '@stdlib/array/float64' );
+ *
+ * var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ * var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ * var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+ *
+ * dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 );
+ * // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+ */
+ ndarray( N: number, nrhs: number, D: Float64Array, sd: number, od: number, E: Float64Array, se: number, oe: number, B: Float64Array, sb1: number, sb2: number, ob: number ): Float64Array;
+}
+
+/**
+* Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf`.
+*
+* @param order - storage layout of `B`
+* @param N - order of the tridiagonal matrix `A`
+* @param nrhs - number of right-hand sides (i.e., the number of columns of the matrix `B`)
+* @param D - the `N` diagonal elements of the diagonal matrix `D` from the factorization of `A`
+* @param E - the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the factorization of `A`
+* @param B - input matrix (overwritten by the solution matrix `X`)
+* @param LDB - leading dimension of `B`
+* @returns solution matrix `X` (stored in `B`)
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2( 'column-major', 3, 1, D, E, B, 3 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*/
+declare var dptts2: Routine;
+
+
+// EXPORTS //
+
+export = dptts2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/docs/types/test.ts b/lib/node_modules/@stdlib/lapack/base/dptts2/docs/types/test.ts
new file mode 100644
index 000000000000..fdb843a454fb
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/docs/types/test.ts
@@ -0,0 +1,375 @@
+/*
+* @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 dptts2 = require( './index' );
+
+
+// TESTS //
+
+// The function returns a Float64Array...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 'column-major', 3, 1, D, E, B, 3 ); // $ExpectType Float64Array
+}
+
+// The compiler throws an error if the function is provided a first argument which is not a valid order...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 5, 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( true, 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( false, 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( null, 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( void 0, 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( [], 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( {}, 3, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( ( x: number ): number => x, 3, 1, D, E, B, 3 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a second argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 'column-major', '5', 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', true, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', false, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', null, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', void 0, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', [], 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', {}, 1, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', ( x: number ): number => x, 1, D, E, B, 3 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a third argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 'column-major', 3, '1', D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, true, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, false, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, null, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, void 0, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, [], D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, {}, D, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, ( x: number ): number => x, D, E, B, 3 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fourth argument which is not a Float64Array...
+{
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 'column-major', 3, 1, '5', E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, 5, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, true, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, false, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, null, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, void 0, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, {}, E, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, ( x: number ): number => x, E, B, 3 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a fifth argument which is not a Float64Array...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 'column-major', 3, 1, D, '5', B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, 5, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, true, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, false, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, null, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, void 0, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, {}, B, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, ( x: number ): number => x, B, 3 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a sixth argument which is not a Float64Array...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+
+ dptts2( 'column-major', 3, 1, D, E, '5', 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, 5, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, true, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, false, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, null, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, void 0, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, {}, 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, ( x: number ): number => x, 3 ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided a seventh argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2( 'column-major', 3, 1, D, E, B, '3' ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, true ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, false ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, null ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, void 0 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, [] ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, {} ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the function is provided an unsupported number of arguments...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2(); // $ExpectError
+ dptts2( 'column-major' ); // $ExpectError
+ dptts2( 'column-major', 3 ); // $ExpectError
+ dptts2( 'column-major', 3, 1 ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B ); // $ExpectError
+ dptts2( 'column-major', 3, 1, D, E, B, 3, 10 ); // $ExpectError
+}
+
+// Attached to main export is an `ndarray` method which returns a Float64Array...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectType Float64Array
+}
+
+// The compiler throws an error if the `ndarray` method is provided a first argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( '3', 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( true, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( false, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( null, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( void 0, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( [], 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( {}, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( ( x: number ): number => x, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a second argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, '1', D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, true, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, false, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, null, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, void 0, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, [], D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, {}, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, ( x: number ): number => x, D, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a third argument which is not a Float64Array...
+{
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, '5', 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, 5, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, true, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, false, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, null, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, void 0, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, {}, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, ( x: number ): number => x, 1, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a fourth argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, '1', 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, true, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, false, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, null, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, void 0, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, [], 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, {}, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, ( x: number ): number => x, 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a fifth argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, '0', E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, true, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, false, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, null, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, void 0, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, [], E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, {}, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, ( x: number ): number => x, E, 1, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a sixth argument which is not a Float64Array...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, '5', 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, 5, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, true, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, false, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, null, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, void 0, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, {}, 1, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, ( x: number ): number => x, 1, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a seventh argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, '1', 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, true, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, false, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, null, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, void 0, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, [], 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, {}, 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, ( x: number ): number => x, 0, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided an eighth argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, '0', B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, true, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, false, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, null, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, void 0, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, [], B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, {}, B, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, ( x: number ): number => x, B, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a ninth argument which is not a Float64Array...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, '5', 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, 5, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, true, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, false, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, null, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, void 0, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, {}, 1, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, ( x: number ): number => x, 1, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a tenth argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, '1', 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, true, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, false, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, null, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, void 0, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, [], 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, {}, 3, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, ( x: number ): number => x, 3, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided an eleventh argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, '3', 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, true, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, false, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, null, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, void 0, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, [], 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, {}, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, ( x: number ): number => x, 0 ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided a twelfth argument which is not a number...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, '0' ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, true ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, false ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, null ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, void 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, [] ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, {} ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, ( x: number ): number => x ); // $ExpectError
+}
+
+// The compiler throws an error if the `ndarray` method is provided an unsupported number of arguments...
+{
+ const D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+ const E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+ const B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+ dptts2.ndarray(); // $ExpectError
+ dptts2.ndarray( 3 ); // $ExpectError
+ dptts2.ndarray( 3, 1 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3 ); // $ExpectError
+ dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0, 10 ); // $ExpectError
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/examples/index.js b/lib/node_modules/@stdlib/lapack/base/dptts2/examples/index.js
new file mode 100644
index 000000000000..2ec56d107f6f
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/examples/index.js
@@ -0,0 +1,37 @@
+/**
+* @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 dpttrf = require( '@stdlib/lapack/base/dpttrf' );
+var dptts2 = require( './../lib' );
+
+// Define a symmetric positive definite matrix `A` via its diagonal and subdiagonal:
+var D = new Float64Array( [ 4.0, 5.0, 6.0, 7.0 ] );
+var E = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+
+// Compute the `L * D * L^T` factorization of `A` (overwrites `D` and `E`):
+dpttrf( D.length, D, E );
+
+// Define the right-hand side matrix `B` (column-major, two right-hand sides):
+var B = new Float64Array( [ 1.0, 2.0, 3.0, 4.0, 2.0, 1.0, 0.0, 1.0 ] );
+
+// Solve `A * X = B`, overwriting `B` with the solution `X`:
+dptts2( 'column-major', 4, 2, D, E, B, 4 );
+console.log( B );
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/lib/base.js b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/base.js
new file mode 100644
index 000000000000..4afe4695c06c
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/base.js
@@ -0,0 +1,106 @@
+/**
+* @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, max-params */
+
+'use strict';
+
+// MAIN //
+
+/**
+* Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf`.
+*
+* @private
+* @param {NonNegativeInteger} N - order of the tridiagonal matrix `A`
+* @param {NonNegativeInteger} nrhs - number of right-hand sides (i.e., the number of columns of the matrix `B`)
+* @param {Float64Array} D - the `N` diagonal elements of the diagonal matrix `D` from the factorization of `A`
+* @param {integer} strideD - stride length for `D`
+* @param {NonNegativeInteger} offsetD - starting index for `D`
+* @param {Float64Array} E - the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the factorization of `A`
+* @param {integer} strideE - stride length for `E`
+* @param {NonNegativeInteger} offsetE - starting index for `E`
+* @param {Float64Array} B - input matrix (overwritten by the solution matrix `X`)
+* @param {integer} strideB1 - stride of the first dimension of `B`
+* @param {integer} strideB2 - stride of the second dimension of `B`
+* @param {NonNegativeInteger} offsetB - starting index for `B`
+* @returns {Float64Array} solution matrix `X` (stored in `B`)
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*/
+function dptts2( N, nrhs, D, strideD, offsetD, E, strideE, offsetE, B, strideB1, strideB2, offsetB ) {
+ var da;
+ var ib;
+ var id;
+ var ie;
+ var oj;
+ var i;
+ var j;
+
+ if ( N <= 1 ) {
+ if ( N === 1 ) {
+ // Scale each right-hand side by `1 / D[0]` (multiplying by the reciprocal, matching `DSCAL`)...
+ da = 1.0 / D[ offsetD ];
+ ib = offsetB;
+ for ( j = 0; j < nrhs; j++ ) {
+ B[ ib ] *= da;
+ ib += strideB2;
+ }
+ }
+ return B;
+ }
+ // Solve `A * X = B` using the factorization `A = L * D * L^T`, overwriting each right-hand side vector with its solution...
+ for ( j = 0; j < nrhs; j++ ) {
+ oj = offsetB + ( j * strideB2 );
+
+ // Solve `L * x = b`...
+ ib = oj;
+ ie = offsetE;
+ for ( i = 1; i < N; i++ ) {
+ // B[i,j] = B[i,j] - ( B[i-1,j] * E[i-1] )
+ B[ ib+strideB1 ] -= B[ ib ] * E[ ie ];
+ ib += strideB1;
+ ie += strideE;
+ }
+ // Solve `D * L^T * x = b`. On entry, `ib` points to `B[N-1,j]`...
+ id = offsetD + ( (N-1)*strideD );
+ B[ ib ] /= D[ id ];
+ ie = offsetE + ( (N-2)*strideE );
+ for ( i = N-2; i >= 0; i-- ) {
+ id -= strideD;
+
+ // B[i,j] = ( B[i,j] / D[i] ) - ( B[i+1,j] * E[i] )
+ B[ ib-strideB1 ] = ( B[ ib-strideB1 ] / D[ id ] ) - ( B[ ib ] * E[ ie ] );
+ ib -= strideB1;
+ ie -= strideE;
+ }
+ }
+ return B;
+}
+
+
+// EXPORTS //
+
+module.exports = dptts2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/lib/dptts2.js b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/dptts2.js
new file mode 100644
index 000000000000..cc5fc8fbd5dd
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/dptts2.js
@@ -0,0 +1,77 @@
+/**
+* @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 isLayout = require( '@stdlib/blas/base/assert/is-layout' );
+var format = require( '@stdlib/string/format' );
+var base = require( './base.js' );
+
+
+// MAIN //
+
+/**
+* Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf`.
+*
+* @param {string} order - storage layout of `B`
+* @param {NonNegativeInteger} N - order of the tridiagonal matrix `A`
+* @param {NonNegativeInteger} nrhs - number of right-hand sides (i.e., the number of columns of the matrix `B`)
+* @param {Float64Array} D - the `N` diagonal elements of the diagonal matrix `D` from the factorization of `A`
+* @param {Float64Array} E - the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the factorization of `A`
+* @param {Float64Array} B - input matrix (overwritten by the solution matrix `X`)
+* @param {NonNegativeInteger} LDB - leading dimension of `B`
+* @throws {TypeError} first argument must be a valid order
+* @throws {RangeError} seventh argument must be greater than or equal to nrhs
+* @returns {Float64Array} solution matrix `X` (stored in `B`)
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2( 'column-major', 3, 1, D, E, B, 3 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*/
+function dptts2( order, N, nrhs, D, E, B, LDB ) {
+ var sb1;
+ var sb2;
+
+ if ( !isLayout( order ) ) {
+ throw new TypeError( format( 'invalid argument. First argument must be a valid order. Value: `%s`.', order ) );
+ }
+ if ( order === 'column-major' ) {
+ sb1 = 1;
+ sb2 = LDB;
+ } else { // order === 'row-major'
+ if ( LDB < nrhs ) {
+ throw new RangeError( format( 'invalid argument. Seventh argument must be greater than or equal to %d. Value: `%d`.', nrhs, LDB ) );
+ }
+ sb1 = LDB;
+ sb2 = 1;
+ }
+ return base( N, nrhs, D, 1, 0, E, 1, 0, B, sb1, sb2, 0 );
+}
+
+
+// EXPORTS //
+
+module.exports = dptts2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/lib/index.js b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/index.js
new file mode 100644
index 000000000000..f4a6e3ce2c11
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/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 solve a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf`.
+*
+* @module @stdlib/lapack/base/dptts2
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+* var dptts2 = require( '@stdlib/lapack/base/dptts2' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2( 'column-major', 3, 1, D, E, B, 3 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+* var dptts2 = require( '@stdlib/lapack/base/dptts2' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2.ndarray( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*/
+
+// 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 dptts2;
+var tmp = tryRequire( join( __dirname, './native.js' ) );
+if ( isError( tmp ) ) {
+ dptts2 = main;
+} else {
+ dptts2 = tmp;
+}
+
+
+// EXPORTS //
+
+module.exports = dptts2;
+
+// exports: { "ndarray": "dptts2.ndarray" }
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/lib/main.js b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/main.js
new file mode 100644
index 000000000000..1eb1dbc4618c
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/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 dptts2 = require( './dptts2.js' );
+var ndarray = require( './ndarray.js' );
+
+
+// MAIN //
+
+setReadOnly( dptts2, 'ndarray', ndarray );
+
+
+// EXPORTS //
+
+module.exports = dptts2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/lib/ndarray.js b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/ndarray.js
new file mode 100644
index 000000000000..9e12f0ff2682
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/lib/ndarray.js
@@ -0,0 +1,62 @@
+/**
+* @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 //
+
+/**
+* Solves a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf` and alternative indexing semantics.
+*
+* @param {NonNegativeInteger} N - order of the tridiagonal matrix `A`
+* @param {NonNegativeInteger} nrhs - number of right-hand sides (i.e., the number of columns of the matrix `B`)
+* @param {Float64Array} D - the `N` diagonal elements of the diagonal matrix `D` from the factorization of `A`
+* @param {integer} sd - stride length for `D`
+* @param {NonNegativeInteger} od - starting index for `D`
+* @param {Float64Array} E - the `N-1` subdiagonal elements of the unit bidiagonal factor `L` from the factorization of `A`
+* @param {integer} se - stride length for `E`
+* @param {NonNegativeInteger} oe - starting index for `E`
+* @param {Float64Array} B - input matrix (overwritten by the solution matrix `X`)
+* @param {integer} sb1 - stride of the first dimension of `B`
+* @param {integer} sb2 - stride of the second dimension of `B`
+* @param {NonNegativeInteger} ob - starting index for `B`
+* @returns {Float64Array} solution matrix `X` (stored in `B`)
+*
+* @example
+* var Float64Array = require( '@stdlib/array/float64' );
+*
+* var D = new Float64Array( [ 4.0, 4.75, 5.157894736842105 ] );
+* var E = new Float64Array( [ 0.25, 0.42105263157894735 ] );
+* var B = new Float64Array( [ 1.0, 2.0, 3.0 ] );
+*
+* dptts2( 3, 1, D, 1, 0, E, 1, 0, B, 1, 3, 0 );
+* // B => [ ~0.2041, ~0.1837, ~0.4388 ]
+*/
+function dptts2( N, nrhs, D, sd, od, E, se, oe, B, sb1, sb2, ob ) { // eslint-disable-line max-params
+ return base( N, nrhs, D, sd, od, E, se, oe, B, sb1, sb2, ob );
+}
+
+
+// EXPORTS //
+
+module.exports = dptts2;
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/package.json b/lib/node_modules/@stdlib/lapack/base/dptts2/package.json
new file mode 100644
index 000000000000..0da63eda392a
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/package.json
@@ -0,0 +1,69 @@
+{
+ "name": "@stdlib/lapack/base/dptts2",
+ "version": "0.0.0",
+ "description": "Solve a tridiagonal system of the form `A * X = B` using the `L * D * L^T` factorization of `A` computed by `dpttrf`.",
+ "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",
+ "dptts2",
+ "tridiagonal",
+ "solve",
+ "linear",
+ "algebra",
+ "subroutines",
+ "array",
+ "ndarray",
+ "float64",
+ "double",
+ "float64array"
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_n_eq_one.json
new file mode 100644
index 000000000000..b38c3903610b
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "column-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 5 ],
+ "sd": 1,
+ "od": 0,
+
+ "E": [],
+ "se": 1,
+ "oe": 0,
+
+ "B": [ 10, 20, 5 ],
+ "sb1": 1,
+ "sb2": 1,
+ "ob": 0,
+ "LDB": 1,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 2, 4, 1 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..f83f924a3f81
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 0,
+
+ "E": [ 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 0,
+
+ "B": [ 1, 2, 3, 4 ],
+ "sb1": 1,
+ "sb2": 4,
+ "ob": 0,
+ "LDB": 4,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 0.3087378640776699, 0.14174757281553402, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..b3fc9c727335
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/col_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 0,
+
+ "E": [ 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 0,
+
+ "B": [ 1, 2, 3, 4, 2, 1, 0, 1, 7, 7, 7, 7 ],
+ "sb1": 1,
+ "sb2": 4,
+ "ob": 0,
+ "LDB": 4,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 0.3087378640776699, 0.14174757281553402, 0.5106796116504854, 0.4563106796116505, 0.17475728155339804, -0.16504854368932037, 0.2135922330097087, 1.5223300970873785, 0.9106796116504855, 0.4621359223300972, 0.8019417475728154 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_n_eq_one.json
new file mode 100644
index 000000000000..ac367a20c5e2
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "column-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 9999, 5 ],
+ "sd": 2,
+ "od": 1,
+
+ "E": [],
+ "se": 2,
+ "oe": 1,
+
+ "B": [ 5, 9999, 9999, 20, 9999, 9999, 10 ],
+ "sb1": 2,
+ "sb2": -3,
+ "ob": 6,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 1, 9999, 9999, 4, 9999, 9999, 2 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..5c2d5b6e6eaf
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 9999, 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 1,
+
+ "E": [ 9999, 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 1,
+
+ "B": [ 1, 9999, 2, 9999, 3, 9999, 4 ],
+ "sb1": 2,
+ "sb2": -9,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 9999, 0.3087378640776699, 9999, 0.14174757281553402, 9999, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..f1944c66ee22
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_col_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 9999, 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 1,
+
+ "E": [ 9999, 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 1,
+
+ "B": [ 7, 9999, 7, 9999, 7, 9999, 7, 9999, 9999, 2, 9999, 1, 9999, 0, 9999, 1, 9999, 9999, 1, 9999, 2, 9999, 3, 9999, 4 ],
+ "sb1": 2,
+ "sb2": -9,
+ "ob": 18,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 1.5223300970873785, 9999, 0.9106796116504855, 9999, 0.4621359223300972, 9999, 0.8019417475728154, 9999, 9999, 0.4563106796116505, 9999, 0.17475728155339804, 9999, -0.16504854368932037, 9999, 0.2135922330097087, 9999, 9999, 0.17281553398058253, 9999, 0.3087378640776699, 9999, 0.14174757281553402, 9999, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_n_eq_one.json
new file mode 100644
index 000000000000..f0b3e37fd686
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "row-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 9999, 5 ],
+ "sd": 2,
+ "od": 1,
+
+ "E": [],
+ "se": 2,
+ "oe": 1,
+
+ "B": [ 5, 9999, 20, 9999, 10 ],
+ "sb1": 7,
+ "sb2": -2,
+ "ob": 4,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 1, 9999, 4, 9999, 2 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..5682b7d5c18f
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 9999, 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 1,
+
+ "E": [ 9999, 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 1,
+
+ "B": [ 1, 9999, 9999, 2, 9999, 9999, 3, 9999, 9999, 4 ],
+ "sb1": 3,
+ "sb2": -2,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 9999, 9999, 0.3087378640776699, 9999, 9999, 0.14174757281553402, 9999, 9999, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..af50ec9db46f
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/complex-access/complex_access_row_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 9999, 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 1,
+
+ "E": [ 9999, 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 1,
+
+ "B": [ 7, 9999, 2, 9999, 1, 9999, 9999, 7, 9999, 1, 9999, 2, 9999, 9999, 7, 9999, 0, 9999, 3, 9999, 9999, 7, 9999, 1, 9999, 4 ],
+ "sb1": 7,
+ "sb2": -2,
+ "ob": 4,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 1.5223300970873785, 9999, 0.4563106796116505, 9999, 0.17281553398058253, 9999, 9999, 0.9106796116504855, 9999, 0.17475728155339804, 9999, 0.3087378640776699, 9999, 9999, 0.4621359223300972, 9999, -0.16504854368932037, 9999, 0.14174757281553402, 9999, 9999, 0.8019417475728154, 9999, 0.2135922330097087, 9999, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_n_eq_one.json
new file mode 100644
index 000000000000..458ae6e9b3b7
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "column-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 5 ],
+ "sd": -2,
+ "od": 0,
+
+ "E": [],
+ "se": -2,
+ "oe": 0,
+
+ "B": [ 5, 9999, 9999, 20, 9999, 9999, 10 ],
+ "sb1": -2,
+ "sb2": -3,
+ "ob": 6,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 1, 9999, 9999, 4, 9999, 9999, 2 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..2484e746fcbd
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 5.255102040816327, 9999, 5.157894736842105, 9999, 4.75, 9999, 4 ],
+ "sd": -2,
+ "od": 6,
+
+ "E": [ 0.5816326530612245, 9999, 0.42105263157894735, 9999, 0.25 ],
+ "se": -2,
+ "oe": 4,
+
+ "B": [ 4, 9999, 3, 9999, 2, 9999, 1 ],
+ "sb1": -2,
+ "sb2": -9,
+ "ob": 6,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.5106796116504854, 9999, 0.14174757281553402, 9999, 0.3087378640776699, 9999, 0.17281553398058253 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..0e0862db1afb
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_col_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 5.255102040816327, 9999, 5.157894736842105, 9999, 4.75, 9999, 4 ],
+ "sd": -2,
+ "od": 6,
+
+ "E": [ 0.5816326530612245, 9999, 0.42105263157894735, 9999, 0.25 ],
+ "se": -2,
+ "oe": 4,
+
+ "B": [ 7, 9999, 7, 9999, 7, 9999, 7, 9999, 9999, 1, 9999, 0, 9999, 1, 9999, 2, 9999, 9999, 4, 9999, 3, 9999, 2, 9999, 1 ],
+ "sb1": -2,
+ "sb2": -9,
+ "ob": 24,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 0.8019417475728154, 9999, 0.4621359223300972, 9999, 0.9106796116504855, 9999, 1.5223300970873785, 9999, 9999, 0.2135922330097087, 9999, -0.16504854368932037, 9999, 0.17475728155339804, 9999, 0.4563106796116505, 9999, 9999, 0.5106796116504854, 9999, 0.14174757281553402, 9999, 0.3087378640776699, 9999, 0.17281553398058253 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_n_eq_one.json
new file mode 100644
index 000000000000..11ec070b6d7d
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "row-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 5 ],
+ "sd": -2,
+ "od": 0,
+
+ "E": [],
+ "se": -2,
+ "oe": 0,
+
+ "B": [ 5, 9999, 20, 9999, 10 ],
+ "sb1": -7,
+ "sb2": -2,
+ "ob": 4,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 1, 9999, 4, 9999, 2 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..96b8e65a34b5
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 5.255102040816327, 9999, 5.157894736842105, 9999, 4.75, 9999, 4 ],
+ "sd": -2,
+ "od": 6,
+
+ "E": [ 0.5816326530612245, 9999, 0.42105263157894735, 9999, 0.25 ],
+ "se": -2,
+ "oe": 4,
+
+ "B": [ 4, 9999, 9999, 3, 9999, 9999, 2, 9999, 9999, 1 ],
+ "sb1": -3,
+ "sb2": -2,
+ "ob": 9,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.5106796116504854, 9999, 9999, 0.14174757281553402, 9999, 9999, 0.3087378640776699, 9999, 9999, 0.17281553398058253 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..c981327a9b8b
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/negative-strides/neg_stride_row_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 5.255102040816327, 9999, 5.157894736842105, 9999, 4.75, 9999, 4 ],
+ "sd": -2,
+ "od": 6,
+
+ "E": [ 0.5816326530612245, 9999, 0.42105263157894735, 9999, 0.25 ],
+ "se": -2,
+ "oe": 4,
+
+ "B": [ 7, 9999, 1, 9999, 4, 9999, 9999, 7, 9999, 0, 9999, 3, 9999, 9999, 7, 9999, 1, 9999, 2, 9999, 9999, 7, 9999, 2, 9999, 1 ],
+ "sb1": -7,
+ "sb2": -2,
+ "ob": 25,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 0.8019417475728154, 9999, 0.2135922330097087, 9999, 0.5106796116504854, 9999, 9999, 0.4621359223300972, 9999, -0.16504854368932037, 9999, 0.14174757281553402, 9999, 9999, 0.9106796116504855, 9999, 0.17475728155339804, 9999, 0.3087378640776699, 9999, 9999, 1.5223300970873785, 9999, 0.4563106796116505, 9999, 0.17281553398058253 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_n_eq_one.json
new file mode 100644
index 000000000000..eb319d80bcf8
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "column-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 9999, 9999, 9999, 5 ],
+ "sd": 1,
+ "od": 3,
+
+ "E": [],
+ "se": 1,
+ "oe": 3,
+
+ "B": [ 9999, 9999, 9999, 10, 20, 5 ],
+ "sb1": 1,
+ "sb2": 1,
+ "ob": 3,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 9999, 9999, 9999, 2, 4, 1 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..59ba83caf4ef
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 9999, 9999, 9999, 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 3,
+
+ "E": [ 9999, 9999, 9999, 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 3,
+
+ "B": [ 9999, 9999, 9999, 1, 2, 3, 4 ],
+ "sb1": 1,
+ "sb2": 4,
+ "ob": 3,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 9999, 9999, 9999, 0.17281553398058253, 0.3087378640776699, 0.14174757281553402, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..bcd3e36d5660
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_col_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 9999, 9999, 9999, 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 3,
+
+ "E": [ 9999, 9999, 9999, 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 3,
+
+ "B": [ 9999, 9999, 9999, 1, 2, 3, 4, 2, 1, 0, 1, 7, 7, 7, 7 ],
+ "sb1": 1,
+ "sb2": 4,
+ "ob": 3,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 9999, 9999, 9999, 0.17281553398058253, 0.3087378640776699, 0.14174757281553402, 0.5106796116504854, 0.4563106796116505, 0.17475728155339804, -0.16504854368932037, 0.2135922330097087, 1.5223300970873785, 0.9106796116504855, 0.4621359223300972, 0.8019417475728154 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_n_eq_one.json
new file mode 100644
index 000000000000..e4e3ab1b8092
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "row-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 9999, 9999, 9999, 5 ],
+ "sd": 1,
+ "od": 3,
+
+ "E": [],
+ "se": 1,
+ "oe": 3,
+
+ "B": [ 9999, 9999, 9999, 10, 20, 5 ],
+ "sb1": 3,
+ "sb2": 1,
+ "ob": 3,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 9999, 9999, 9999, 2, 4, 1 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..3752e92c5286
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 9999, 9999, 9999, 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 3,
+
+ "E": [ 9999, 9999, 9999, 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 3,
+
+ "B": [ 9999, 9999, 9999, 1, 2, 3, 4 ],
+ "sb1": 1,
+ "sb2": 1,
+ "ob": 3,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 9999, 9999, 9999, 0.17281553398058253, 0.3087378640776699, 0.14174757281553402, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..7274afd595b6
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/offsets/offset_row_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 9999, 9999, 9999, 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 3,
+
+ "E": [ 9999, 9999, 9999, 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 3,
+
+ "B": [ 9999, 9999, 9999, 1, 2, 7, 2, 1, 7, 3, 0, 7, 4, 1, 7 ],
+ "sb1": 3,
+ "sb2": 1,
+ "ob": 3,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 9999, 9999, 9999, 0.17281553398058253, 0.4563106796116505, 1.5223300970873785, 0.3087378640776699, 0.17475728155339804, 0.9106796116504855, 0.14174757281553402, -0.16504854368932037, 0.4621359223300972, 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_n_eq_one.json
new file mode 100644
index 000000000000..18ffafe89bad
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "column-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 5 ],
+ "sd": 2,
+ "od": 0,
+
+ "E": [],
+ "se": 2,
+ "oe": 0,
+
+ "B": [ 10, 9999, 9999, 20, 9999, 9999, 5 ],
+ "sb1": 2,
+ "sb2": 3,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 2, 9999, 9999, 4, 9999, 9999, 1 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..7367b29f2a9b
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 0,
+
+ "E": [ 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 0,
+
+ "B": [ 1, 9999, 2, 9999, 3, 9999, 4 ],
+ "sb1": 2,
+ "sb2": 9,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 9999, 0.3087378640776699, 9999, 0.14174757281553402, 9999, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..85f21f568649
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_col_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "column-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 0,
+
+ "E": [ 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 0,
+
+ "B": [ 1, 9999, 2, 9999, 3, 9999, 4, 9999, 9999, 2, 9999, 1, 9999, 0, 9999, 1, 9999, 9999, 7, 9999, 7, 9999, 7, 9999, 7 ],
+ "sb1": 2,
+ "sb2": 9,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 9999, 0.3087378640776699, 9999, 0.14174757281553402, 9999, 0.5106796116504854, 9999, 9999, 0.4563106796116505, 9999, 0.17475728155339804, 9999, -0.16504854368932037, 9999, 0.2135922330097087, 9999, 9999, 1.5223300970873785, 9999, 0.9106796116504855, 9999, 0.4621359223300972, 9999, 0.8019417475728154 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_n_eq_one.json
new file mode 100644
index 000000000000..ce33fbebedbc
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "row-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 5 ],
+ "sd": 2,
+ "od": 0,
+
+ "E": [],
+ "se": 2,
+ "oe": 0,
+
+ "B": [ 10, 9999, 20, 9999, 5 ],
+ "sb1": 7,
+ "sb2": 2,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 2, 9999, 4, 9999, 1 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..4faa073bde80
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 0,
+
+ "E": [ 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 0,
+
+ "B": [ 1, 9999, 9999, 2, 9999, 9999, 3, 9999, 9999, 4 ],
+ "sb1": 3,
+ "sb2": 2,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 9999, 9999, 0.3087378640776699, 9999, 9999, 0.14174757281553402, 9999, 9999, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..504759514ef6
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/positive-strides/pos_stride_row_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 4, 9999, 4.75, 9999, 5.157894736842105, 9999, 5.255102040816327 ],
+ "sd": 2,
+ "od": 0,
+
+ "E": [ 0.25, 9999, 0.42105263157894735, 9999, 0.5816326530612245 ],
+ "se": 2,
+ "oe": 0,
+
+ "B": [ 1, 9999, 2, 9999, 7, 9999, 9999, 2, 9999, 1, 9999, 7, 9999, 9999, 3, 9999, 0, 9999, 7, 9999, 9999, 4, 9999, 1, 9999, 7 ],
+ "sb1": 7,
+ "sb2": 2,
+ "ob": 0,
+ "LDB": null,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 9999, 0.4563106796116505, 9999, 1.5223300970873785, 9999, 9999, 0.3087378640776699, 9999, 0.17475728155339804, 9999, 0.9106796116504855, 9999, 9999, 0.14174757281553402, 9999, -0.16504854368932037, 9999, 0.4621359223300972, 9999, 9999, 0.5106796116504854, 9999, 0.2135922330097087, 9999, 0.8019417475728154 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_n_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_n_eq_one.json
new file mode 100644
index 000000000000..6d1d083674f4
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_n_eq_one.json
@@ -0,0 +1,30 @@
+{
+ "order": "row-major",
+
+ "N": 1,
+ "nrhs": 3,
+
+ "D": [ 5 ],
+ "sd": 1,
+ "od": 0,
+
+ "E": [],
+ "se": 1,
+ "oe": 0,
+
+ "B": [ 10, 20, 5 ],
+ "sb1": 3,
+ "sb2": 1,
+ "ob": 0,
+ "LDB": 3,
+
+ "B_mat": [
+ [ 10, 20, 5 ]
+ ],
+
+ "expected": [ 2, 4, 1 ],
+
+ "expected_mat": [
+ [ 2, 4, 1 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_nrhs_eq_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_nrhs_eq_one.json
new file mode 100644
index 000000000000..537ee5d63897
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_nrhs_eq_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 1,
+
+ "D": [ 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 0,
+
+ "E": [ 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 0,
+
+ "B": [ 1, 2, 3, 4 ],
+ "sb1": 1,
+ "sb2": 1,
+ "ob": 0,
+ "LDB": 1,
+
+ "B_mat": [
+ [ 1 ],
+ [ 2 ],
+ [ 3 ],
+ [ 4 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 0.3087378640776699, 0.14174757281553402, 0.5106796116504854 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253 ],
+ [ 0.3087378640776699 ],
+ [ 0.14174757281553402 ],
+ [ 0.5106796116504854 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_nrhs_gt_one.json b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_nrhs_gt_one.json
new file mode 100644
index 000000000000..05e700bbb25f
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/fixtures/row_major_nrhs_gt_one.json
@@ -0,0 +1,36 @@
+{
+ "order": "row-major",
+
+ "N": 4,
+ "nrhs": 3,
+
+ "D": [ 4, 4.75, 5.157894736842105, 5.255102040816327 ],
+ "sd": 1,
+ "od": 0,
+
+ "E": [ 0.25, 0.42105263157894735, 0.5816326530612245 ],
+ "se": 1,
+ "oe": 0,
+
+ "B": [ 1, 2, 7, 2, 1, 7, 3, 0, 7, 4, 1, 7 ],
+ "sb1": 3,
+ "sb2": 1,
+ "ob": 0,
+ "LDB": 3,
+
+ "B_mat": [
+ [ 1, 2, 7 ],
+ [ 2, 1, 7 ],
+ [ 3, 0, 7 ],
+ [ 4, 1, 7 ]
+ ],
+
+ "expected": [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785, 0.3087378640776699, 0.17475728155339804, 0.9106796116504855, 0.14174757281553402, -0.16504854368932037, 0.4621359223300972, 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ],
+
+ "expected_mat": [
+ [ 0.17281553398058253, 0.4563106796116505, 1.5223300970873785 ],
+ [ 0.3087378640776699, 0.17475728155339804, 0.9106796116504855 ],
+ [ 0.14174757281553402, -0.16504854368932037, 0.4621359223300972 ],
+ [ 0.5106796116504854, 0.2135922330097087, 0.8019417475728154 ]
+ ]
+}
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.dptts2.js b/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.dptts2.js
new file mode 100644
index 000000000000..74aaaf3a239c
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.dptts2.js
@@ -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.
+*/
+
+'use strict';
+
+// MODULES //
+
+var tape = require( 'tape' );
+var Float64Array = require( '@stdlib/array/float64' );
+var dptts2 = require( './../lib/dptts2.js' );
+
+
+// FIXTURES //
+
+var COL_MAJOR_NRHS_EQ_ONE = require( './fixtures/col_major_nrhs_eq_one.json' );
+var COL_MAJOR_NRHS_GT_ONE = require( './fixtures/col_major_nrhs_gt_one.json' );
+var COL_MAJOR_N_EQ_ONE = require( './fixtures/col_major_n_eq_one.json' );
+var ROW_MAJOR_NRHS_EQ_ONE = require( './fixtures/row_major_nrhs_eq_one.json' );
+var ROW_MAJOR_NRHS_GT_ONE = require( './fixtures/row_major_nrhs_gt_one.json' );
+var ROW_MAJOR_N_EQ_ONE = require( './fixtures/row_major_n_eq_one.json' );
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+ t.ok( true, __filename );
+ t.strictEqual( typeof dptts2, 'function', 'main export is a function' );
+ t.end();
+});
+
+tape( 'the function has an arity of 7', function test( t ) {
+ t.strictEqual( dptts2.length, 7, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function throws an error if provided a first argument which is not a valid order', function test( t ) {
+ var values;
+ var data;
+ var i;
+
+ data = COL_MAJOR_NRHS_GT_ONE;
+
+ values = [
+ 'foo',
+ 'bar',
+ 'beep',
+ 'boop',
+ -5,
+ NaN,
+ true,
+ false,
+ null,
+ void 0,
+ [],
+ {},
+ function noop() {}
+ ];
+
+ for ( i = 0; i < values.length; i++ ) {
+ t.throws( badValue( values[ i ] ), TypeError, 'throws an error when provided ' + values[ i ] );
+ }
+ t.end();
+
+ function badValue( value ) {
+ return function badValue() {
+ dptts2( value, data.N, data.nrhs, new Float64Array( data.D ), new Float64Array( data.E ), new Float64Array( data.B ), data.LDB ); // eslint-disable-line max-len
+ };
+ }
+});
+
+tape( 'the function throws an error if provided a seventh argument which is not a valid `LDB` value (row-major)', function test( t ) {
+ var values;
+ var data;
+ var i;
+
+ data = ROW_MAJOR_NRHS_GT_ONE;
+
+ values = [
+ 0,
+ 1,
+ 2
+ ];
+
+ for ( i = 0; i < values.length; i++ ) {
+ t.throws( badValue( values[ i ] ), RangeError, 'throws an error when provided ' + values[ i ] );
+ }
+ t.end();
+
+ function badValue( value ) {
+ return function badValue() {
+ dptts2( data.order, data.N, data.nrhs, new Float64Array( data.D ), new Float64Array( data.E ), new Float64Array( data.B ), value ); // eslint-disable-line max-len
+ };
+ }
+});
+
+tape( 'the function returns the input array unchanged if `N` is zero', function test( t ) {
+ var expected;
+ var out;
+ var B;
+
+ B = new Float64Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+ expected = new Float64Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+ out = dptts2( 'column-major', 0, 2, new Float64Array( [] ), new Float64Array( [] ), B, 2 );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function solves `A * X = B` (column-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = COL_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.order, data.N, data.nrhs, D, E, B, data.LDB );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COL_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.order, data.N, data.nrhs, D, E, B, data.LDB );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COL_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.order, data.N, data.nrhs, D, E, B, data.LDB );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function solves `A * X = B` (row-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = ROW_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.order, data.N, data.nrhs, D, E, B, data.LDB );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = ROW_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.order, data.N, data.nrhs, D, E, B, data.LDB );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = ROW_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.order, data.N, data.nrhs, D, E, B, data.LDB );
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.js b/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.js
new file mode 100644
index 000000000000..29a922da9a04
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/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 dptts2 = 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 dptts2, '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 dptts2.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 dptts2 = proxyquire( './../lib', {
+ '@stdlib/utils/try-require': tryRequire
+ });
+
+ t.strictEqual( dptts2, 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 dptts2;
+ var main;
+
+ main = require( './../lib/dptts2.js' );
+
+ dptts2 = proxyquire( './../lib', {
+ '@stdlib/utils/try-require': tryRequire
+ });
+
+ t.strictEqual( dptts2, main, 'returns expected value' );
+ t.end();
+
+ function tryRequire() {
+ return new Error( 'Cannot find module' );
+ }
+});
diff --git a/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.ndarray.js b/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.ndarray.js
new file mode 100644
index 000000000000..b66d06ef16e1
--- /dev/null
+++ b/lib/node_modules/@stdlib/lapack/base/dptts2/test/test.ndarray.js
@@ -0,0 +1,468 @@
+/**
+* @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 dptts2 = require( './../lib/ndarray.js' );
+
+
+// FIXTURES //
+
+var COL_MAJOR_NRHS_EQ_ONE = require( './fixtures/col_major_nrhs_eq_one.json' );
+var COL_MAJOR_NRHS_GT_ONE = require( './fixtures/col_major_nrhs_gt_one.json' );
+var COL_MAJOR_N_EQ_ONE = require( './fixtures/col_major_n_eq_one.json' );
+var ROW_MAJOR_NRHS_EQ_ONE = require( './fixtures/row_major_nrhs_eq_one.json' );
+var ROW_MAJOR_NRHS_GT_ONE = require( './fixtures/row_major_nrhs_gt_one.json' );
+var ROW_MAJOR_N_EQ_ONE = require( './fixtures/row_major_n_eq_one.json' );
+var POS_STRIDE_COL_MAJOR_NRHS_EQ_ONE = require( './fixtures/positive-strides/pos_stride_col_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var POS_STRIDE_COL_MAJOR_NRHS_GT_ONE = require( './fixtures/positive-strides/pos_stride_col_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var POS_STRIDE_COL_MAJOR_N_EQ_ONE = require( './fixtures/positive-strides/pos_stride_col_major_n_eq_one.json' ); // eslint-disable-line id-length
+var POS_STRIDE_ROW_MAJOR_NRHS_EQ_ONE = require( './fixtures/positive-strides/pos_stride_row_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var POS_STRIDE_ROW_MAJOR_NRHS_GT_ONE = require( './fixtures/positive-strides/pos_stride_row_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var POS_STRIDE_ROW_MAJOR_N_EQ_ONE = require( './fixtures/positive-strides/pos_stride_row_major_n_eq_one.json' ); // eslint-disable-line id-length
+var NEG_STRIDE_COL_MAJOR_NRHS_EQ_ONE = require( './fixtures/negative-strides/neg_stride_col_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var NEG_STRIDE_COL_MAJOR_NRHS_GT_ONE = require( './fixtures/negative-strides/neg_stride_col_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var NEG_STRIDE_COL_MAJOR_N_EQ_ONE = require( './fixtures/negative-strides/neg_stride_col_major_n_eq_one.json' ); // eslint-disable-line id-length
+var NEG_STRIDE_ROW_MAJOR_NRHS_EQ_ONE = require( './fixtures/negative-strides/neg_stride_row_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var NEG_STRIDE_ROW_MAJOR_NRHS_GT_ONE = require( './fixtures/negative-strides/neg_stride_row_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var NEG_STRIDE_ROW_MAJOR_N_EQ_ONE = require( './fixtures/negative-strides/neg_stride_row_major_n_eq_one.json' ); // eslint-disable-line id-length
+var OFFSET_COL_MAJOR_NRHS_EQ_ONE = require( './fixtures/offsets/offset_col_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var OFFSET_COL_MAJOR_NRHS_GT_ONE = require( './fixtures/offsets/offset_col_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var OFFSET_COL_MAJOR_N_EQ_ONE = require( './fixtures/offsets/offset_col_major_n_eq_one.json' );
+var OFFSET_ROW_MAJOR_NRHS_EQ_ONE = require( './fixtures/offsets/offset_row_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var OFFSET_ROW_MAJOR_NRHS_GT_ONE = require( './fixtures/offsets/offset_row_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var OFFSET_ROW_MAJOR_N_EQ_ONE = require( './fixtures/offsets/offset_row_major_n_eq_one.json' );
+var COMPLEX_ACCESS_COL_MAJOR_NRHS_EQ_ONE = require( './fixtures/complex-access/complex_access_col_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var COMPLEX_ACCESS_COL_MAJOR_NRHS_GT_ONE = require( './fixtures/complex-access/complex_access_col_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var COMPLEX_ACCESS_COL_MAJOR_N_EQ_ONE = require( './fixtures/complex-access/complex_access_col_major_n_eq_one.json' ); // eslint-disable-line id-length
+var COMPLEX_ACCESS_ROW_MAJOR_NRHS_EQ_ONE = require( './fixtures/complex-access/complex_access_row_major_nrhs_eq_one.json' ); // eslint-disable-line id-length
+var COMPLEX_ACCESS_ROW_MAJOR_NRHS_GT_ONE = require( './fixtures/complex-access/complex_access_row_major_nrhs_gt_one.json' ); // eslint-disable-line id-length
+var COMPLEX_ACCESS_ROW_MAJOR_N_EQ_ONE = require( './fixtures/complex-access/complex_access_row_major_n_eq_one.json' ); // eslint-disable-line id-length
+
+
+// TESTS //
+
+tape( 'main export is a function', function test( t ) {
+ t.ok( true, __filename );
+ t.strictEqual( typeof dptts2, 'function', 'main export is a function' );
+ t.end();
+});
+
+tape( 'the function has an arity of 12', function test( t ) {
+ t.strictEqual( dptts2.length, 12, 'returns expected value' );
+ t.end();
+});
+
+tape( 'the function returns the input array unchanged if `N` is zero', function test( t ) {
+ var expected;
+ var out;
+ var B;
+
+ B = new Float64Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+ expected = new Float64Array( [ 1.0, 2.0, 3.0, 4.0 ] );
+
+ out = dptts2( 0, 2, new Float64Array( [] ), 1, 0, new Float64Array( [] ), 1, 0, B, 1, 2, 0 ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function solves `A * X = B` (column-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = COL_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COL_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COL_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function solves `A * X = B` (row-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = ROW_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = ROW_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = ROW_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports specifying positive strides (column-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = POS_STRIDE_COL_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = POS_STRIDE_COL_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = POS_STRIDE_COL_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports specifying positive strides (row-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = POS_STRIDE_ROW_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = POS_STRIDE_ROW_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = POS_STRIDE_ROW_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports specifying negative strides (column-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = NEG_STRIDE_COL_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = NEG_STRIDE_COL_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = NEG_STRIDE_COL_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports specifying negative strides (row-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = NEG_STRIDE_ROW_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = NEG_STRIDE_ROW_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = NEG_STRIDE_ROW_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports specifying offsets (column-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = OFFSET_COL_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = OFFSET_COL_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = OFFSET_COL_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports specifying offsets (row-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = OFFSET_ROW_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = OFFSET_ROW_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = OFFSET_ROW_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports complex access patterns (mixed sign strides and offsets, column-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = COMPLEX_ACCESS_COL_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COMPLEX_ACCESS_COL_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COMPLEX_ACCESS_COL_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});
+
+tape( 'the function supports complex access patterns (mixed sign strides and offsets, row-major)', function test( t ) {
+ var expected;
+ var data;
+ var out;
+ var D;
+ var E;
+ var B;
+
+ data = COMPLEX_ACCESS_ROW_MAJOR_NRHS_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COMPLEX_ACCESS_ROW_MAJOR_NRHS_GT_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ data = COMPLEX_ACCESS_ROW_MAJOR_N_EQ_ONE;
+ D = new Float64Array( data.D );
+ E = new Float64Array( data.E );
+ B = new Float64Array( data.B );
+ expected = new Float64Array( data.expected );
+ out = dptts2( data.N, data.nrhs, D, data.sd, data.od, E, data.se, data.oe, B, data.sb1, data.sb2, data.ob ); // eslint-disable-line max-len
+ t.strictEqual( out, B, 'returns a reference to `B`' );
+ t.deepEqual( out, expected, 'returns expected value' );
+
+ t.end();
+});