refactor: base implementation and consequent changes

lib/node_modules/@stdlib/blas/base/dtrsm/README.md

@@ -81,14 +81,14 @@ var A0 = new Float64Array( [ 0.0, 1.0, 3.0, 0.0, 4.0 ] );
 var B0 = new Float64Array( [ 0.0, 5.0, 7.0, 0.0, 8.0 ] );
 
 // Create offset views...
-var A1 = new Float64Array( A0.buffer, A0.BYTES_PER_ELEMENT*1 ); // start at 1st element
-var B1 = new Float64Array( B0.buffer, B0.BYTES_PER_ELEMENT*1 ); // start at 1st element
+var A1 = new Float64Array( A0.buffer, A0.BYTES_PER_ELEMENT*1 ); // start at 2nd element
+var B1 = new Float64Array( B0.buffer, B0.BYTES_PER_ELEMENT*1 ); // start at 2nd element
 
 dtrsm( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A1, 2, B1, 2 );
 // B0 => <Float64Array>[ 0.0, 30.0, 6.0, 0.0, 12.0 ]
 ```
 
-#### dtrsm.ndarray( ord, side, ul, ta, diag, m, n, α, A, LDA, oa, B, LDB, ob )
+#### dtrsm.ndarray( o, s, ul, t, d, m, n, α, A, sa1, sa2, oa, B, sb1, sb2, ob )
 
 Solves matrix equation `op(A) * X = alpha * B` or `X * op(A) = alpha * B` using alternative indexing semantics where `alpha` is a scalar, `X` and `B` are m by n matrices, `A` is a unit, or non-unit, upper or lower triangular matrix  and `op( A )` is one  of `op( A ) = A` or `op( A ) = A**T`. The matrix `X` is overwritten on `B`.
 
@@ -98,13 +98,17 @@ var Float64Array = require( '@stdlib/array/float64' );
 A = new Float64Array( [ 0.0, 0.0, 1.0, 3.0, 0.0, 4.0 ] );
 B = new Float64Array( [ 0.0, 5.0, 7.0, 0.0, 8.0 ] );
 
-dtrsm.ndarray( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 2, B, 2, 1 );
+dtrsm.ndarray( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 1, 2, B, 2, 1, 1 );
 // B => <Float64Array>[ 0.0, 30.0, 6.0, 0.0, 12.0 ]
 ```
 
 The function has the following additional parameters:
 
+-   **sa1**: stride of the first dimension of `A`.
+-   **sa2**: stride of the second dimension of `A`.
 -   **oa**: starting index for `A`.
+-   **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. For example,
@@ -115,7 +119,7 @@ var Float64Array = require( '@stdlib/array/float64' );
 A = new Float64Array( [ 0.0, 0.0, 1.0, 3.0, 0.0, 4.0 ] );
 B = new Float64Array( [ 0.0, 5.0, 7.0, 0.0, 8.0 ] );
 
-dtrsm.ndarray( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 2, B, 2, 1 );
+dtrsm.ndarray( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 1, 2, B, 2, 1, 1 );
 // B => <Float64Array>[ 0.0, 30.0, 6.0, 0.0, 12.0 ]
 ```
 

lib/node_modules/@stdlib/blas/base/dtrsm/benchmark/benchmark.ndarray.js

@@ -62,7 +62,7 @@ function createBenchmark( N ) {
 
 		b.tic();
 		for ( i = 0; i < b.iterations; i++ ) {
-			z = dtrsm( 'column-major', 'left', 'upper', 'none', 'non-unit', N, N, 1.0, A, N, 0, B, N, 0 );
+			z = dtrsm( 'column-major', 'left', 'upper', 'none', 'non-unit', N, N, 1.0, A, 1, N, 0, B, 1, N, 0 );
 			if ( isnan( z[ i%z.length ] ) ) {
 				b.fail( 'should not return NaN' );
 			}

lib/node_modules/@stdlib/blas/base/dtrsm/docs/types/index.d.ts

@@ -66,23 +66,25 @@ interface Routine {
 	* @param n - number of columns in `B`
 	* @param alpha - scalar constant
 	* @param A - input matrix `A`
-	* @param LDA - stride of the first dimension of `A` (a.k.a., leading dimension of the matrix `A`)
-	* @param offsetA - starting `A` index
+	* @param strideA1 - stride of the first dimension of `A`
+	* @param strideA2 - stride of the second dimension of `A`
+	* @param offsetA - starting index for `A`
 	* @param B - input matrix `B`
-	* @param LDB - stride of the first dimension of `B` (a.k.a., leading dimension of the matrix `B`)
-	* @param offsetB - starting `B` index
+	* @param strideB1 - stride of the first dimension of `B`
+	* @param strideB2 - stride of the second dimension of `B`
+	* @param offsetB - starting index for `B`
 	* @returns `B`
 	*
 	* @example
 	* var Float64Array = require( '@stdlib/array/float64' );
 	*
-	* var A = new Float64Array( [ 1.0, 3.0, 0.0, 4.0 ] );
-	* var B = new Float64Array( [ 5.0, 7.0, 0.0, 8.0 ] );
+	* var A = new Float64Array( [ 0.0, 0.0, 1.0, 3.0, 0.0, 4.0 ] );
+	* var B = new Float64Array( [ 0.0, 5.0, 7.0, 0.0, 8.0 ] );
 	*
-	* dtrsm( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 0, B, 2, 0 );
-	* // B => <Float64Array>[ 30.0, 6.0, 0.0, 12.0 ]
+	* dtrsm.ndarray( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 1, 2, B, 2, 1, 1 );
+	* // B => <Float64Array>[ 0.0, 30.0, 6.0, 0.0, 12.0 ]
 	*/
-	ndarray( order: Layout, side: OperationSide, uplo: MatrixTriangle, transa: TransposeOperation, diag: DiagonalType, m: number, n: number, alpha: number, A: Float64Array, LDA: number, offsetA: number, B: Float64Array, LDB: number, offsetB: number ): Float64Array;
+	ndarray( order: Layout, side: OperationSide, uplo: MatrixTriangle, transa: TransposeOperation, diag: DiagonalType, m: number, n: number, alpha: number, A: Float64Array, strideA1: number, strideA2: number, offsetA: number, B: Float64Array, strideB1: number, strideB2: number, offsetB: number ): Float64Array;
 }
 
 /**
@@ -114,11 +116,11 @@ interface Routine {
 * @example
 * var Float64Array = require( '@stdlib/array/float64' );
 *
-* var A = new Float64Array( [ 1.0, 3.0, 0.0, 4.0 ] );
-* var B = new Float64Array( [ 5.0, 7.0, 0.0, 8.0 ] );
+* var A = new Float64Array( [ 0.0, 0.0, 1.0, 3.0, 0.0, 4.0 ] );
+* var B = new Float64Array( [ 0.0, 5.0, 7.0, 0.0, 8.0 ] );
 *
-* dtrsm( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 0, B, 2, 0 );
-* // B => <Float64Array>[ 30.0, 6.0, 0.0, 12.0 ]
+* dtrsm.ndarray( 'row-major', 'left', 'upper', 'none', 'non-unit', 2, 2, 6.0, A, 2, 1, 2, B, 2, 1, 1 );
+* // B => <Float64Array>[ 0.0, 30.0, 6.0, 0.0, 12.0 ]
 */
 declare var dtrsm: Routine;