WIP

hbina · hbina · commit 4f8d25b4415b · 2025-02-08T03:28:33.000+08:00
Signed-off-by: Hanif Ariffin &lt;hanif.ariffin.4326@gmail.com&gt;
diff --git a/stdlib/src/math.rs b/stdlib/src/math.rs
@@ -635,11 +635,11 @@ mod math {
             let mut x = *obj?;
 
             let xsave = x;
-            let mut j = 0;
+            let mut i = 0;
             // This inner loop applies `hi`/`lo` summation to each
             // partial so that the list of partial sums remains exact.
-            for i in 0..partials.len() {
-                let mut y: f64 = partials[i];
+            for j in 0..partials.len() {
+                let mut y: f64 = partials[j];
                 if x.abs() < y.abs() {
                     std::mem::swap(&mut x, &mut y);
                 }
@@ -648,8 +648,8 @@ mod math {
                 let hi = x + y;
                 let lo = y - (hi - x);
                 if lo != 0.0 {
-                    partials[j] = lo;
-                    j += 1;
+                    partials[i] = lo;
+                    i += 1;
                 }
                 x = hi;
             }
@@ -668,15 +668,11 @@ mod math {
                 special_sum += xsave;
                 // reset partials
                 partials.clear();
-            }
-
-            if j >= partials.len() {
-                partials.push(x);
             } else {
-                partials[j] = x;
-                partials.truncate(j + 1);
+                partials.push(x);
             }
         }
+
         if special_sum != 0.0 {
             return if inf_sum.is_nan() {
                 Err(vm.new_value_error("-inf + inf in fsum".to_owned()))
