Auto-generated commit

Author: stdlib-bot

CHANGELOG.md

@@ -4,7 +4,7 @@
 
 <section class="release" id="unreleased">
 
-## Unreleased (2025-07-11)
+## Unreleased (2025-07-12)
 
 <section class="features">
 
@@ -520,6 +520,7 @@ A total of 24 issues were closed in this release:
 
 <details>
 
+-   [`a41c427`](https://github.com/stdlib-js/stdlib/commit/a41c427ee9ddd87f48288a737ea9820ea48b5498) - **refactor:** reduce pointer arithmetic operations _(by Athan Reines)_
 -   [`c9f4cb9`](https://github.com/stdlib-js/stdlib/commit/c9f4cb9b593f963d25efc20a6bd3b2aeeda515c2) - **fix:** add missing checks and tests _(by Athan Reines)_
 -   [`fffbf73`](https://github.com/stdlib-js/stdlib/commit/fffbf730b1a27175ece26a82f05c4053567e4386) - **feat:** add C implementation for `blas/base/dsyr`  [(#6566)](https://github.com/stdlib-js/stdlib/pull/6566) _(by Shabareesh Shetty, Athan Reines)_
 -   [`b7ebe14`](https://github.com/stdlib-js/stdlib/commit/b7ebe14a38ddc3ccce67834db16835f420431494) - **feat:** add `blas/ext/base/wasm/dnansumkbn2` [(#5735)](https://github.com/stdlib-js/stdlib/pull/5735) _(by Shabareesh Shetty, Athan Reines, JoyBoy, stdlib-bot)_

base/dsyr/lib/base.js

@@ -59,7 +59,7 @@ function dsyr( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offse
 	var sa1;
 	var i0;
 	var i1;
-	var oa;
+	var ia;
 	var ox;
 
 	isrm = isRowMajor( [ strideA1, strideA2 ] );
@@ -81,11 +81,12 @@ function dsyr( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offse
 		for ( i1 = 0; i1 < N; i1++ ) {
 			if ( x[ ix1 ] !== 0.0 ) {
 				tmp = alpha * x[ ix1 ];
-				oa = offsetA + (sa1*i1);
+				ia = offsetA + (sa1*i1);
 				ix0 = ox;
 				for ( i0 = 0; i0 <= i1; i0++ ) {
-					A[ oa+(sa0*i0) ] += x[ ix0 ] * tmp;
+					A[ ia ] += x[ ix0 ] * tmp;
 					ix0 += strideX;
+					ia += sa0;
 				}
 			}
 			ix1 += strideX;
@@ -97,11 +98,12 @@ function dsyr( uplo, N, alpha, x, strideX, offsetX, A, strideA1, strideA2, offse
 	for ( i1 = 0; i1 < N; i1++ ) {
 		if ( x[ ix1 ] !== 0.0 ) {
 			tmp = alpha * x[ ix1 ];
-			oa = offsetA + (sa1*i1);
+			ia = offsetA + (sa1*i1) + (sa0*i1);
 			ix0 = ix1;
 			for ( i0 = i1; i0 < N; i0++ ) {
-				A[ oa+(sa0*i0) ] += x[ ix0 ] * tmp;
+				A[ ia ] += x[ ix0 ] * tmp;
 				ix0 += strideX;
+				ia += sa0;
 			}
 		}
 		ix1 += strideX;

base/dsyr/src/dsyr_ndarray.c

@@ -44,7 +44,7 @@ void API_SUFFIX(c_dsyr_ndarray)( const CBLAS_UPLO uplo, const CBLAS_INT N, const
 	CBLAS_INT sa1;
 	CBLAS_INT i0;
 	CBLAS_INT i1;
-	CBLAS_INT oa;
+	CBLAS_INT ia;
 	CBLAS_INT ox;
 	double tmp;
 
@@ -89,11 +89,12 @@ void API_SUFFIX(c_dsyr_ndarray)( const CBLAS_UPLO uplo, const CBLAS_INT N, const
 		for ( i1 = 0; i1 < N; i1++ ) {
 			if ( X[ ix1 ] != 0.0 ) {
 				tmp = alpha * X[ ix1 ];
-				oa = offsetA + (sa1*i1);
+				ia = offsetA + (sa1*i1);
 				ix0 = ox;
 				for ( i0 = 0; i0 <= i1; i0++ ) {
-					A[ oa+(sa0*i0) ] += X[ ix0 ] * tmp;
+					A[ ia ] += X[ ix0 ] * tmp;
 					ix0 += strideX;
+					ia += sa0;
 				}
 			}
 			ix1 += strideX;
@@ -105,11 +106,12 @@ void API_SUFFIX(c_dsyr_ndarray)( const CBLAS_UPLO uplo, const CBLAS_INT N, const
 	for ( i1 = 0; i1 < N; i1++ ) {
 		if ( X[ ix1 ] != 0.0 ) {
 			tmp = alpha * X[ ix1 ];
-			oa = offsetA + (sa1*i1);
+			ia = offsetA + (sa1*i1) + (sa0*i1);
 			ix0 = ix1;
 			for ( i0 = i1; i0 < N; i0++ ) {
-				A[ oa+(sa0*i0) ] += X[ ix0 ] * tmp;
+				A[ ia ] += X[ ix0 ] * tmp;
 				ix0 += strideX;
+				ia += sa0;
 			}
 		}
 		ix1 += strideX;