Merge pull request #5057 from martin-frbg/issue5050

Replace while loop in generic C/ZGEMM_BETA to avoid going out of bounds
OpenMathLib · Jan 11, 2025 · d74eb02 · d74eb02
2 parents 30f7a41 + d91d4fa
commit d74eb02
Showing 1 changed file with 6 additions and 8 deletions.
diff --git a/kernel/generic/zgemm_beta.c b/kernel/generic/zgemm_beta.c
@@ -58,8 +58,8 @@ int CNAME(BLASLONG m, BLASLONG n, BLASLONG dummy1,
   c_offset = c;
 
   if (beta_r == 0. && beta_i == 0.) {
-    j = n;
-    do {
+
+    for (j=n;j>0;j--) {
       c_offset1 = c_offset;
       c_offset += ldc;
 
@@ -88,13 +88,12 @@ int CNAME(BLASLONG m, BLASLONG n, BLASLONG dummy1,
 	  i--;
 	} while (i > 0);
       }
-      j --;
-    } while (j > 0);
+    }
 
   } else {
 
-    j = n;
-    do {
+
+    for (j=n;j>0;j--) {
       c_offset1 = c_offset;
       c_offset += ldc;
 
@@ -151,8 +150,7 @@ int CNAME(BLASLONG m, BLASLONG n, BLASLONG dummy1,
 	  i --;
 	} while (i > 0);
       }
-      j --;
-    } while (j > 0);
+    }
   }
   return 0;
 }