diff -r 221728aeeb7e -r fa8df5a14486 src/fe_model/p2_local_model.rs
--- a/src/fe_model/p2_local_model.rs	Wed Sep 03 19:55:05 2025 -0500
+++ b/src/fe_model/p2_local_model.rs	Wed Sep 03 20:19:41 2025 -0500
@@ -32,7 +32,7 @@
 impl<'a, F: Float> Set<Loc<1, F>> for RealInterval<F> {
     #[inline]
     fn contains<I: Instance<Loc<1, F>>>(&self, z: I) -> bool {
-        z.eval_ref_decompose(|&Loc([x])| {
+        z.eval_ref(|&Loc([x])| {
             let &[Loc([x0]), Loc([x1])] = &self.0;
             (x0 < x && x < x1) || (x1 < x && x < x0)
         })
@@ -318,12 +318,7 @@
         x0: &Loc<2, F>,
         x1: &Loc<2, F>, /*, v0 : F, v1 : F*/
     ) -> (Loc<2, F>, F) {
-        let &P2LocalModel {
-            a0,
-            a1: Loc([a1, a2]),
-            a2: Loc([a11, a12, a22]),
-            ..
-        } = self;
+        let &P2LocalModel { a0, a1: Loc([a1, a2]), a2: Loc([a11, a12, a22]), .. } = self;
         let &Loc([x00, x01]) = x0;
         let d @ Loc([d0, d1]) = x1 - x0;
         let b0 = a0 + a1 * x00 + a2 * x01 + a11 * x00 * x00 + a12 * x00 * x01 + a22 * x01 * x01;