fix: implementation of inside

src/evox/problems/numerical/maf.py

@@ -12,7 +12,7 @@
 @jit
 def inside(x, a, b):
     """check if x is in [a, b] or [b, a]"""
-    return (a <= x <= b) | (b <= x <= a)
+    return jnp.minimum(a, b) <= x & x <= jnp.maximum(a, b)
 
 
 @jit