Generator: corr. bug predicates

onera · Jan 10, 2025 · 6c40cd9 · 6c40cd9
1 parent 89be293
commit 6c40cd9
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/Cassiopee/Generator/Generator/Tetgen/predicates.cxx b/Cassiopee/Generator/Generator/Tetgen/predicates.cxx
@@ -995,10 +995,12 @@ int fast_expansion_sum_zeroelim(int elen, REAL *e, int flen, REAL *f, REAL *h)
     while ((eindex < elen) && (findex < flen)) {
       if ((fnow > enow) == (fnow > -enow)) {
         Two_Sum(Q, enow, Qnew, hh);
-        enow = e[++eindex];
+        ++eindex;
+        enow = e[(eindex < elen) ? eindex : elen-1];
       } else {
         Two_Sum(Q, fnow, Qnew, hh);
-        fnow = f[++findex];
+        ++findex;
+        fnow = f[(findex < flen) ? findex : flen-1];
       }
       Q = Qnew;
       if (hh != 0.0) {
@@ -1008,15 +1010,17 @@ int fast_expansion_sum_zeroelim(int elen, REAL *e, int flen, REAL *f, REAL *h)
   }
   while (eindex < elen) {
     Two_Sum(Q, enow, Qnew, hh);
-    enow = e[++eindex];
+    ++eindex;
+    enow = e[(eindex < elen) ? eindex : elen-1];
     Q = Qnew;
     if (hh != 0.0) {
       h[hindex++] = hh;
     }
   }
   while (findex < flen) {
     Two_Sum(Q, fnow, Qnew, hh);
-    fnow = f[++findex];
+    ++findex;
+    fnow = f[(findex < flen) ? findex : flen-1];
     Q = Qnew;
     if (hh != 0.0) {
       h[hindex++] = hh;