Bit decomp interpreter

executor/src/witgen/jit/affine_symbolic_expression.rs

@@ -6,12 +6,13 @@ use std::{
 
 use itertools::Itertools;
 use num_traits::Zero;
-use powdr_number::FieldElement;
+use powdr_number::{log2_exact, FieldElement};
 
 use crate::witgen::jit::effect::Assertion;
 
 use super::{
-    super::range_constraints::RangeConstraint, effect::Effect,
+    super::range_constraints::RangeConstraint,
+    effect::{BitDecomposition, BitDecompositionComponent, Effect},
     symbolic_expression::SymbolicExpression,
 };
 
@@ -215,23 +216,19 @@ impl<T: FieldElement, V: Ord + Clone + Display> AffineSymbolicExpression<T, V> {
             }
             _ => {
                 let r = self.solve_bit_decomposition()?;
+
                 if r.complete {
                     r
                 } else {
                     let negated = -self;
-                    let r = negated.solve_bit_decomposition()?;
-                    if r.complete {
-                        r
-                    } else {
-                        let effects = self
-                            .transfer_constraints()
-                            .into_iter()
-                            .chain(negated.transfer_constraints())
-                            .collect();
-                        ProcessResult {
-                            effects,
-                            complete: false,
-                        }
+                    let effects = self
+                        .transfer_constraints()
+                        .into_iter()
+                        .chain(negated.transfer_constraints())
+                        .collect();
+                    ProcessResult {
+                        effects,
+                        complete: false,
                     }
                 }
             }
@@ -257,20 +254,32 @@ impl<T: FieldElement, V: Ord + Clone + Display> AffineSymbolicExpression<T, V> {
 
         // Check if they are mutually exclusive and compute assignments.
         let mut covered_bits: <T as FieldElement>::Integer = 0.into();
-        let mut effects = vec![];
-        for (var, coeff, constraint) in constrained_coefficients {
-            let mask = *constraint.multiple(coeff).mask();
-            if !(mask & covered_bits).is_zero() {
+        let mut components = vec![];
+        for (variable, coeff, constraint) in constrained_coefficients {
+            let is_negative = !coeff.is_in_lower_half();
+            let coeff_abs = if is_negative { -coeff } else { coeff };
+            let Some(exponent) = log2_exact(coeff_abs.to_arbitrary_integer()) else {
+                // We could work with non-powers of two, but it would require
+                // division instead of shifts.
+                return Ok(ProcessResult::empty());
+            };
+            let bit_mask = *constraint.multiple(coeff_abs).mask();
+            if !(bit_mask & covered_bits).is_zero() {
                 // Overlapping range constraints.
                 return Ok(ProcessResult::empty());
             } else {
-                covered_bits |= mask;
+                covered_bits |= bit_mask;
             }
-            let masked = -&self.offset & mask;
-            effects.push(Effect::Assignment(
-                var.clone(),
-                masked.integer_div(&coeff.into()),
-            ));
+            components.push(BitDecompositionComponent {
+                variable,
+                // We negate here because we are solving
+                // c_1 * x_1 + c_2 * x_2 + ... + offset = 0,
+                // instead of
+                // c_1 * x_1 + c_2 * x_2 + ... = offset.
+                is_negative: !is_negative,
+                exponent: exponent as u64,
+                bit_mask,
+            });
         }
 
         if covered_bits >= T::modulus() {
@@ -280,18 +289,20 @@ impl<T: FieldElement, V: Ord + Clone + Display> AffineSymbolicExpression<T, V> {
         // We need to assert that the masks cover "-offset",
         // otherwise the equation is not solvable.
         // We assert -offset & !masks == 0
+
+        // TODO can this be done even with negative coefficients?
         if let Some(offset) = self.offset.try_to_number() {
-            if (-offset).to_integer() & !covered_bits != 0.into() {
+            if offset.to_integer() & !covered_bits != 0.into() {
                 return Err(Error::ConflictingRangeConstraints);
             }
-        } else {
-            effects.push(Assertion::assert_eq(
-                -&self.offset & !covered_bits,
-                T::from(0).into(),
-            ));
         }
 
-        Ok(ProcessResult::complete(effects))
+        Ok(ProcessResult::complete(vec![Effect::BitDecomposition(
+            BitDecomposition {
+                value: self.offset.clone(),
+                components,
+            },
+        )]))
     }
 
     fn transfer_constraints(&self) -> Option<Effect<T, V>> {
@@ -521,10 +532,9 @@ mod test {
         let b = Ase::from_unknown_variable("b", rc.clone());
         let c = Ase::from_unknown_variable("c", rc.clone());
         let z = Ase::from_known_symbol("Z", Default::default());
-        // a * 0x100 + b * 0x10000 + c * 0x1000000 + 10 + Z = 0
+        // a * 0x100 - b * 0x10000 + c * 0x1000000 + 10 + Z = 0
         let ten = from_number(10);
-        let constr = mul(&a, &from_number(0x100))
-            + mul(&b, &from_number(0x10000))
+        let constr = mul(&a, &from_number(0x100)) - mul(&b, &from_number(0x10000))
             + mul(&c, &from_number(0x1000000))
             + ten.clone()
             + z.clone();
@@ -533,38 +543,39 @@ mod test {
         assert!(!result.complete && result.effects.is_empty());
         // Now add the range constraint on a, it should be solvable.
         let a = Ase::from_unknown_variable("a", rc.clone());
-        let constr = mul(&a, &from_number(0x100))
-            + mul(&b, &from_number(0x10000))
+        let constr = mul(&a, &from_number(0x100)) - mul(&b, &from_number(0x10000))
             + mul(&c, &from_number(0x1000000))
             + ten.clone()
             + z;
         let result = constr.solve().unwrap();
         assert!(result.complete);
-        let effects = result
-            .effects
-            .into_iter()
-            .map(|effect| match effect {
-                Effect::Assignment(v, expr) => format!("{v} = {expr};\n"),
-                Effect::Assertion(Assertion {
-                    lhs,
-                    rhs,
-                    expected_equal,
-                }) => {
-                    format!(
-                        "assert {lhs} {} {rhs};\n",
-                        if expected_equal { "==" } else { "!=" }
-                    )
-                }
-                _ => panic!(),
+
+        let [effect] = &result.effects[..] else {
+            panic!();
+        };
+        let Effect::BitDecomposition(BitDecomposition { value, components }) = effect else {
+            panic!();
+        };
+        assert_eq!(format!("{value}"), "(10 + Z)");
+        let formatted = components
+            .iter()
+            .map(|c| {
+                format!(
+                    "{} = (({value} & {}) >> {}){};\n",
+                    c.variable,
+                    c.bit_mask,
+                    c.exponent,
+                    if c.is_negative { " [negative]" } else { "" }
+                )
             })
-            .format("")
-            .to_string();
+            .join("");
+
         assert_eq!(
-            effects,
-            "a = ((-(10 + Z) & 0xff00) // 256);
-b = ((-(10 + Z) & 0xff0000) // 65536);
-c = ((-(10 + Z) & 0xff000000) // 16777216);
-assert (-(10 + Z) & 0xffffffff000000ff) == 0;
+            formatted,
+            "\
+a = (((10 + Z) & 65280) >> 8) [negative];
+b = (((10 + Z) & 16711680) >> 16);
+c = (((10 + Z) & 4278190080) >> 24) [negative];
 "
         );
     }

executor/src/witgen/jit/block_machine_processor.rs

@@ -521,31 +521,19 @@ main_binary::operation_id_next[0] = main_binary::operation_id[1];
 call_var(9, 0, 0) = main_binary::operation_id_next[0];
 main_binary::operation_id_next[1] = main_binary::operation_id[2];
 call_var(9, 1, 0) = main_binary::operation_id_next[1];
-main_binary::A_byte[2] = ((main_binary::A[3] & 0xff000000) // 16777216);
-main_binary::A[2] = (main_binary::A[3] & 0xffffff);
-assert (main_binary::A[3] & 0xffffffff00000000) == 0;
+2**24 * main_binary::A_byte[2] + 2**0 * main_binary::A[2] := main_binary::A[3];
 call_var(9, 2, 1) = main_binary::A_byte[2];
-main_binary::A_byte[1] = ((main_binary::A[2] & 0xff0000) // 65536);
-main_binary::A[1] = (main_binary::A[2] & 0xffff);
-assert (main_binary::A[2] & 0xffffffffff000000) == 0;
+2**16 * main_binary::A_byte[1] + 2**0 * main_binary::A[1] := main_binary::A[2];
 call_var(9, 1, 1) = main_binary::A_byte[1];
-main_binary::A_byte[0] = ((main_binary::A[1] & 0xff00) // 256);
-main_binary::A[0] = (main_binary::A[1] & 0xff);
-assert (main_binary::A[1] & 0xffffffffffff0000) == 0;
+2**8 * main_binary::A_byte[0] + 2**0 * main_binary::A[0] := main_binary::A[1];
 call_var(9, 0, 1) = main_binary::A_byte[0];
 main_binary::A_byte[-1] = main_binary::A[0];
 call_var(9, -1, 1) = main_binary::A_byte[-1];
-main_binary::B_byte[2] = ((main_binary::B[3] & 0xff000000) // 16777216);
-main_binary::B[2] = (main_binary::B[3] & 0xffffff);
-assert (main_binary::B[3] & 0xffffffff00000000) == 0;
+2**24 * main_binary::B_byte[2] + 2**0 * main_binary::B[2] := main_binary::B[3];
 call_var(9, 2, 2) = main_binary::B_byte[2];
-main_binary::B_byte[1] = ((main_binary::B[2] & 0xff0000) // 65536);
-main_binary::B[1] = (main_binary::B[2] & 0xffff);
-assert (main_binary::B[2] & 0xffffffffff000000) == 0;
+2**16 * main_binary::B_byte[1] + 2**0 * main_binary::B[1] := main_binary::B[2];
 call_var(9, 1, 2) = main_binary::B_byte[1];
-main_binary::B_byte[0] = ((main_binary::B[1] & 0xff00) // 256);
-main_binary::B[0] = (main_binary::B[1] & 0xff);
-assert (main_binary::B[1] & 0xffffffffffff0000) == 0;
+2**8 * main_binary::B_byte[0] + 2**0 * main_binary::B[0] := main_binary::B[1];
 call_var(9, 0, 2) = main_binary::B_byte[0];
 main_binary::B_byte[-1] = main_binary::B[0];
 call_var(9, -1, 2) = main_binary::B_byte[-1];
@@ -620,9 +608,7 @@ params[1] = Sub::b[0];"
         assert_eq!(
             format_code(&code),
             "SubM::a[0] = params[0];
-SubM::b[0] = ((SubM::a[0] & 0xff00) // 256);
-SubM::c[0] = (SubM::a[0] & 0xff);
-assert (SubM::a[0] & 0xffffffffffff0000) == 0;
+2**8 * SubM::b[0] + 2**0 * SubM::c[0] := SubM::a[0];
 params[1] = SubM::b[0];
 params[2] = SubM::c[0];
 call_var(1, 0, 0) = SubM::c[0];
@@ -802,4 +788,33 @@ S::Z[0] = ((S::Zi[0][0] + S::Zi[1][0]) + S::Zi[2][0]);
 params[2] = S::Z[0];"
         );
     }
+
+    #[test]
+    fn bit_decomp_negative() {
+        let input = "
+        namespace Main(256);
+            col witness a, b, c;
+            [a, b, c] is [S.Y, S.Z,  S.carry];
+        namespace S(256);
+            let BYTE: col = |i| i & 0xff;
+            let X;
+            let Y;
+            let Z;
+            let carry;
+            carry * (1 - carry) = 0;
+            [ X ] in [ BYTE ];
+            [ Y ] in [ BYTE ];
+            [ Z ] in [ BYTE ];
+            X + Y = Z + 256 * carry;
+        ";
+        let code = format_code(&generate_for_block_machine(input, "S", 2, 1).unwrap().code);
+        assert_eq!(
+            code,
+            "\
+S::Y[0] = params[0];
+S::Z[0] = params[1];
+-2**0 * S::X[0] + 2**8 * S::carry[0] := (S::Y[0] + -S::Z[0]);
+params[2] = S::carry[0];"
+        );
+    }
 }

executor/src/witgen/jit/code_cleaner.rs

@@ -50,7 +50,7 @@ fn optional_vars_in_effect<T: FieldElement>(
     required: &mut HashSet<Variable>,
 ) -> HashSet<Variable> {
     let needed = match &effect {
-        Effect::Assignment(..) | Effect::ProverFunctionCall(..) => {
+        Effect::Assignment(..) | Effect::ProverFunctionCall(..) | Effect::BitDecomposition(_) => {
             effect.written_vars().any(|(v, _)| required.contains(v))
         }
         Effect::Assertion(_) => false,