Skip to content

Commit 62abab2

Browse files
committed
refactor: c implementation
--- 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: na - 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: missing_dependencies - 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 ca04fbe commit 62abab2

File tree

1 file changed

+7
-5
lines changed

1 file changed

+7
-5
lines changed

lib/node_modules/@stdlib/blas/base/ssyr/src/ssyr_ndarray.c

Lines changed: 7 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -44,7 +44,7 @@ void API_SUFFIX(c_ssyr_ndarray)( const CBLAS_UPLO uplo, const CBLAS_INT N, const
4444
CBLAS_INT sa1;
4545
CBLAS_INT i0;
4646
CBLAS_INT i1;
47-
CBLAS_INT oa;
47+
CBLAS_INT ia;
4848
CBLAS_INT ox;
4949
float tmp;
5050

@@ -89,11 +89,12 @@ void API_SUFFIX(c_ssyr_ndarray)( const CBLAS_UPLO uplo, const CBLAS_INT N, const
8989
for ( i1 = 0; i1 < N; i1++ ) {
9090
if ( X[ ix1 ] != 0.0f ) {
9191
tmp = alpha * X[ ix1 ];
92-
oa = offsetA + (sa1*i1);
92+
ia = offsetA + (sa1*i1);
9393
ix0 = ox;
9494
for ( i0 = 0; i0 <= i1; i0++ ) {
95-
A[ oa+(sa0*i0) ] += X[ ix0 ] * tmp;
95+
A[ ia ] += X[ ix0 ] * tmp;
9696
ix0 += strideX;
97+
ia += sa0;
9798
}
9899
}
99100
ix1 += strideX;
@@ -105,11 +106,12 @@ void API_SUFFIX(c_ssyr_ndarray)( const CBLAS_UPLO uplo, const CBLAS_INT N, const
105106
for ( i1 = 0; i1 < N; i1++ ) {
106107
if ( X[ ix1 ] != 0.0f ) {
107108
tmp = alpha * X[ ix1 ];
108-
oa = offsetA + (sa1*i1);
109+
ia = offsetA + (sa1*i1) + (sa0*i1);
109110
ix0 = ix1;
110111
for ( i0 = i1; i0 < N; i0++ ) {
111-
A[ oa+(sa0*i0) ] += X[ ix0 ] * tmp;
112+
A[ ia ] += X[ ix0 ] * tmp;
112113
ix0 += strideX;
114+
ia += sa0;
113115
}
114116
}
115117
ix1 += strideX;

0 commit comments

Comments
 (0)