proof-of-portfolio 0.0.104__py3-none-any.whl → 0.0.106__py3-none-any.whl

proof_of_portfolio/_version.py

@@ -1,2 +1,2 @@
 # This file is auto-generated during build
-__version__ = "0.0.104"
+__version__ = "0.0.106"

proof_of_portfolio/circuits/components/src/core/omega.nr

@@ -47,16 +47,28 @@ pub fn omega(
                 omega_loss_min
             };
 
-            // Apply cross-multiplication like in PTN reference implementation
-            let positive_sum_weighted = product_sum_positive * sum_weights_negative;
-            let negative_sum_weighted = product_sum_negative * sum_weights_positive;
-
-            let effective_denominator = if (-negative_sum_weighted) >= omega_loss_min {
-                -negative_sum_weighted
+            // Apply PTN cross-multiplication with overflow-safe scaling
+            // positive_sum = product_sum_positive * sum_weights_negative
+            // negative_sum = product_sum_negative * sum_weights_positive
+            // final = (positive_sum * RATIO_SCALE_FACTOR) / max(abs(negative_sum), omega_loss_min)
+
+            // To avoid overflow, we can rearrange: (a*b*c)/d = (a*c)*(b/d) when b > d
+            let positive_cross = product_sum_positive * sum_weights_negative;
+            let negative_cross = product_sum_negative * sum_weights_positive;
+
+            let abs_negative = if negative_cross >= 0 { negative_cross } else { -negative_cross };
+            let effective_denominator = if abs_negative >= omega_loss_min {
+                abs_negative
             } else {
                 omega_loss_min
             };
-            (positive_sum_weighted * RATIO_SCALE_FACTOR) / effective_denominator
+
+            // Avoid overflow by doing division first when possible
+            if effective_denominator >= RATIO_SCALE_FACTOR {
+                positive_cross / (effective_denominator / RATIO_SCALE_FACTOR)
+            } else {
+                (positive_cross / effective_denominator) * RATIO_SCALE_FACTOR
+            }
         } else {
             let mut positive_sum: i64 = 0;
             let mut negative_sum: i64 = 0;