refactor: pointer arithmetic

---
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
---

Author: aayush0325

lib/node_modules/@stdlib/lapack/base/dgebal/lib/base.js

@@ -16,7 +16,7 @@
 * limitations under the License.
 */
 
-/* eslint-disable max-len, max-params, max-statements */
+/* eslint-disable max-len, max-params, max-statements, max-lines-per-function */
 
 'use strict';
 
@@ -226,8 +226,8 @@ function dgebal( job, N, A, strideA1, strideA2, offsetA, out, strideOut, offsetO
 		}
 	}
 
-	// Initialize `scale` for non-permuted submatrix
-	is = offsetScale + (k*strideScale);
+	// Initialize scale for non-permuted submatrix
+	is = offsetScale + ( k * strideScale );
 	for ( i = k; i <= l; i++ ) {
 		scale[ is ] = 1.0;
 		is += strideScale;
@@ -247,16 +247,25 @@ function dgebal( job, N, A, strideA1, strideA2, offsetA, out, strideOut, offsetO
 
 	noconv = true;
 	while ( noconv ) {
+		is = offsetScale + ( k * strideScale ); // Follows scale
+		ia1 = offsetA + ( k * strideA1 ) + ( k * strideA2 ); // Follows A[ k, i ]
+		ia2 = offsetA + ( k * strideA1 ) + ( k * strideA2 ); // Follows A[ i, k ]
+		ia3 = offsetA + ( k * strideA2 ); // follows `i`th column of A
+
 		noconv = false;
 		for ( i = k; i <= l; i++ ) {
-			c = dnrm2( l-k+1, A, strideA1, offsetA + (k*strideA1) + (i*strideA2) );
-			r = dnrm2( l-k+1, A, strideA2, offsetA + (i*strideA1) + (k*strideA2) );
-			ica = idamax( l+1, A, strideA1, offsetA + (i*strideA2) );
+			c = dnrm2( l-k+1, A, strideA1, ia1 );
+			r = dnrm2( l-k+1, A, strideA2, ia2 );
+			ica = idamax( l+1, A, strideA1, ia3 );
 			ca = abs( A[ offsetA + (ica*strideA1) + (i*strideA2) ] );
-			ira = idamax( N-k+1, A, strideA2, offsetA + (i*strideA1) + (k*strideA2) );
+			ira = idamax( N-k+1, A, strideA2, ia2 );
 			ra = abs( A[ offsetA + (i*strideA1) + ((ira+k)*strideA2) ] );
 
 			if ( c === 0.0 || r === 0.0 ) {
+				ia1 += strideA2;
+				ia2 += strideA1;
+				is += strideScale;
+				ia3 += strideA2;
 				continue;
 			}
 
@@ -290,27 +299,44 @@ function dgebal( job, N, A, strideA1, strideA2, offsetA, out, strideOut, offsetO
 
 			// Now balance
 			if ( ( c + r ) >= factor * s ) {
+				ia1 += strideA2;
+				ia2 += strideA1;
+				is += strideScale;
+				ia3 += strideA2;
 				continue;
 			}
 
-			if ( f < 1.0 && scale[ offsetScale + (i*strideScale) ] < 1.0 ) {
-				if ( f * scale[ offsetScale + (i*strideScale) ] <= sfmin1 ) {
+			if ( f < 1.0 && scale[ is ] < 1.0 ) {
+				if ( f * scale[ is ] <= sfmin1 ) {
+					ia1 += strideA2;
+					ia2 += strideA1;
+					is += strideScale;
+					ia3 += strideA2;
 					continue;
 				}
 			}
 
-			if ( f > 1.0 && scale[ offsetScale + (i*strideScale) ] > 1.0 ) {
-				if ( scale[ offsetScale + (i*strideScale) ] >= sfmax1 / f ) {
+			if ( f > 1.0 && scale[ is ] > 1.0 ) {
+				if ( scale[ is ] >= sfmax1 / f ) {
+					ia1 += strideA2;
+					ia2 += strideA1;
+					is += strideScale;
+					ia3 += strideA2;
 					continue;
 				}
 			}
 
 			g = 1.0 / f;
-			scale[ offsetScale + (i*strideScale) ] *= f;
+			scale[ is ] *= f;
 			noconv = true;
 
-			dscal( N-k, g, A, strideA2, offsetA + (i*strideA1) + (k*strideA2) );
-			dscal( l+1, f, A, strideA1, offsetA + (i*strideA2) );
+			dscal( N-k, g, A, strideA2, ia2 );
+			dscal( l+1, f, A, strideA1, ia3 );
+
+			ia1 += strideA2;
+			ia2 += strideA1;
+			is += strideScale;
+			ia3 += strideA2;
 		}
 	}