fix: The shape property is now consistent between LinearExpression and QuadraticExpression - both exclude internal dimensions (_term for linear, _term and

  _factor for quadratic)

Author: FBumann

linopy/expressions.py

@@ -1069,10 +1069,9 @@ def nterm(self) -> int:
     @property
     def shape(self) -> tuple[int, ...]:
         """
-        Get the total shape of the linear expression.
+        Get the shape of the expression (excluding term dimension).
         """
-        assert self.vars.shape == self.coeffs.shape
-        return self.vars.shape
+        return tuple(v for k, v in self.sizes.items() if k != TERM_DIM)
 
     @property
     def size(self) -> int:

test/test_linear_expression.py

@@ -317,7 +317,7 @@ def test_linear_expression_with_constant_multiplication(
 
     obs = expr * pd.Series([1, 2, 3], index=pd.RangeIndex(3, name="new_dim"))
     assert isinstance(obs, LinearExpression)
-    assert obs.shape == (2, 3, 1)
+    assert obs.shape == (2, 3)
 
 
 def test_linear_expression_multi_indexed(u: Variable) -> None: