Skip to content

Commit d50308b

Browse files
committed
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 ---
1 parent 5f0fdaa commit d50308b

File tree

1 file changed

+6
-1
lines changed
  • lib/node_modules/@stdlib/lapack/base/dgebal/lib

1 file changed

+6
-1
lines changed

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

Lines changed: 6 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -194,6 +194,9 @@ function dgebal( job, N, A, strideA1, strideA2, offsetA, out, strideOut, offsetO
194194
is = offsetScale + ( k * strideScale ); // Follows scale
195195
ia1 = offsetA + ( k * strideA2 ); // Follows `j`th column of A
196196
ia2 = offsetA + ( k * strideA2 ); // Follows `k`th column of A
197+
ia3 = offsetA + ( k * strideA1 ); // Follows `j`th row of A
198+
ia4 = offsetA + ( k * strideA1 ); // Follows `k`th row of A
199+
197200
for ( j = k; j <= l; j++ ) {
198201
canSwap = true;
199202
for ( i = k; i <= l; i++ ) {
@@ -207,16 +210,18 @@ function dgebal( job, N, A, strideA1, strideA2, offsetA, out, strideOut, offsetO
207210
scale[ is ] = j;
208211
if ( j !== k ) {
209212
dswap( l+1, A, strideA1, ia1, A, strideA1, ia2 );
210-
dswap( N-k, A, strideA2, offsetA + (j*strideA1), A, strideA2, offsetA + (k*strideA1) );
213+
dswap( N-k, A, strideA2, ia3, A, strideA2, ia4 );
211214
}
212215
noconv = true;
213216
k += 1;
214217

215218
is += strideScale;
216219
ia2 += strideA2;
220+
ia4 += strideA1;
217221
}
218222

219223
ia1 += strideA2;
224+
ia3 += strideA1;
220225
}
221226
}
222227
}

0 commit comments

Comments
 (0)