chore: add comments

---
type: pre_commit_static_analysis_report
description: Results of running static analysis checks when committing changes.
report:
  - task: lint_filenames
    status: passed
  - task: lint_editorconfig
    status: passed
  - task: lint_markdown
    status: na
  - task: lint_package_json
    status: na
  - task: lint_repl_help
    status: na
  - task: lint_javascript_src
    status: passed
  - task: lint_javascript_cli
    status: na
  - task: lint_javascript_examples
    status: na
  - task: lint_javascript_tests
    status: na
  - task: lint_javascript_benchmarks
    status: na
  - task: lint_python
    status: na
  - task: lint_r
    status: na
  - task: lint_c_src
    status: na
  - task: lint_c_examples
    status: na
  - task: lint_c_benchmarks
    status: na
  - task: lint_c_tests_fixtures
    status: na
  - task: lint_shell
    status: na
  - task: lint_typescript_declarations
    status: na
  - task: lint_typescript_tests
    status: na
  - task: lint_license_headers
    status: passed
---

lib/node_modules/@stdlib/blas/base/dtrmm/lib/base.js

@@ -55,7 +55,7 @@
 * zeros( 2, 3, X, 1, 2, 0 );
 * // X => <Float64Array>[ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 ]
 */
 function zeros( M, N, X, strideX1, strideX2, offsetX ) { // TODO: consider moving to a separate package
 	var dx0;
 	var dx1;
 	var S0;
@@ -136,27 +136,31 @@
 	var j;
 	var k;
 
+	// Note on variable naming convention: sa#, ix#, i# where # corresponds to the loop number, with `0` being the innermost loop...
+
 	isrma = isRowMajor( [ strideA1, strideA2 ] );
 	nonunit = ( diag === 'non-unit' );
 
 	if ( M === 0 || N === 0 ) {
 		return B;
 	}
 	if ( isrma ) {
-		sa0 = strideA2;
-		sa1 = strideA1;
-		sb0 = strideB2;
-		sb1 = strideB1;
+		// For row-major matrices, the last dimension has the fastest changing index...
+		sa0 = strideA2; // stride for innermost loop
+		sa1 = strideA1; // stride for outermost loop
+		sb0 = strideB2; // stride for innermost loop
+		sb1 = strideB1; // stride for outermost loop
 	} else {
-		sa0 = strideA1;
-		sa1 = strideA2;
-		sb0 = strideB1;
-		sb1 = strideB2;
+		sa0 = strideA1; // stride for innermost loop
+		sa1 = strideA2; // stride for outermost loop
+		sb0 = strideB1; // stride for innermost loop
+		sb1 = strideB2; // stride for outermost loop
 	}
 	if ( alpha === 0.0 ) {
 		zeros( M, N, B, sb0, sb1, offsetB );
 		return B;
 	}
+
 	if (
 		( isrma && side === 'left' && uplo === 'upper' && transa === 'no-transpose' ) ||
 		( !isrma && side === 'right' && uplo === 'lower' && transa === 'no-transpose' )
@@ -329,7 +333,7 @@
 		}
 		return B;
 	}
 	// ( isrma && side === 'right' && uplo === 'lower' && transa !== 'no-transpose' ) || ( !isrma && side === 'left' && uplo === 'upper' && transa !== 'no-transpose' )
 	for ( i = 0; i < M; i++ ) {
 		for ( j = N - 1; j >= 0; j-- ) {
 			ob2 = offsetB + ( i * sb1 ) + ( j * sb0 );