causing 2.5.0__tar.gz → 2.5.1__tar.gz

{causing-2.5.0 → causing-2.5.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: causing
-Version: 2.5.0
+Version: 2.5.1
 Summary: Causing: CAUSal INterpretation using Graphs
 Home-page: https://github.com/realrate/Causing
 Author: Dr. Holger Bartel

{causing-2.5.0 → causing-2.5.1}/causing/model.py

@@ -2,13 +2,18 @@ from __future__ import annotations
 
 from dataclasses import dataclass, field
 from typing import Sequence, Callable
-from functools import cached_property
+from functools import cached_property, reduce
 import networkx
 
 import sympy
 import numpy as np
 
 
+def _vectorized_max(*args):
+    """Elementwise maximum over any number of scalars/arrays (sympy.Max)."""
+    return reduce(np.maximum, args)
+
+
 class NumericModelError(Exception):
     pass
 
@@ -196,8 +201,11 @@ class Model:
         # Define placeholder for vectorized max function
         vectorized_max = sympy.Function("vectorized_max")
 
-        # Define custom translation mapping
-        custom_modules = [{"vectorized_max": np.maximum}, "numpy"]
+        # Define custom translation mapping. sympy.Max takes any number of
+        # arguments, np.maximum only two: a third positional argument is its
+        # `out` array, so Max(a, b, c) silently became max(a, b) (and overwrote
+        # c). Reduce pairwise over all arguments instead.
+        custom_modules = [{"vectorized_max": _vectorized_max}, "numpy"]
 
         for i, eq in enumerate(self.equations):
             # Replace sympy.Max with our placeholder

{causing-2.5.0 → causing-2.5.1}/causing.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: causing
-Version: 2.5.0
+Version: 2.5.1
 Summary: Causing: CAUSal INterpretation using Graphs
 Home-page: https://github.com/realrate/Causing
 Author: Dr. Holger Bartel

{causing-2.5.0 → causing-2.5.1}/setup.py

@@ -5,7 +5,7 @@ with open("README.md", "r", encoding="utf-8") as fh:
 
 setuptools.setup(
     name="causing",
-    version="2.5.0",
+    version="2.5.1",
     author="Dr. Holger Bartel",
     author_email="holger.bartel@realrate.ai",
     description="Causing: CAUSal INterpretation using Graphs",

{causing-2.5.0 → causing-2.5.1}/tests/test_model.py

@@ -448,5 +448,34 @@ class TestEndToEndWorkflow(unittest.TestCase):
         self.assertEqual(effects["yhat"].shape, (1, 3))
 
 
+class TestMaxWithManyArguments(unittest.TestCase):
+    """sympy.Max with more than two arguments must use all of them.
+
+    np.maximum(a, b, c) treats c as the output array, so the third argument was
+    silently dropped (RealRate de_health_insurance ZÜVN, found 14.11.2026).
+    """
+
+    def test_three_argument_max_uses_every_argument(self):
+        from sympy import Max
+
+        X1, X2, Y1 = symbols(["X1", "X2", "Y1"])
+        m = Model(xvars=[X1, X2], yvars=[Y1], equations=(Max(X1, X2, 0),), final_var=Y1)
+        xdat = np.array([[1.0, -5.0, -1.0, 3.0], [2.0, -2.0, 4.0, -7.0]])
+        x_before = xdat.copy()
+
+        yhat = m.compute(xdat)
+
+        np.testing.assert_allclose(yhat[0], [2.0, 0.0, 4.0, 3.0])
+        np.testing.assert_array_equal(xdat, x_before)  # inputs not overwritten
+
+    def test_two_argument_max_unchanged(self):
+        from sympy import Max
+
+        X1, Y1 = symbols(["X1", "Y1"])
+        m = Model(xvars=[X1], yvars=[Y1], equations=(Max(X1, 0),), final_var=Y1)
+        yhat = m.compute(np.array([[-1.0, 2.0]]))
+        np.testing.assert_allclose(yhat[0], [0.0, 2.0])
+
+
 if __name__ == "__main__":
     unittest.main()