Don't introduce Second in AlegrabicCanonizer because of shape specialization (only for broadcasting)

Author: ricardoV94

pytensor/tensor/rewriting/math.py

@@ -1135,10 +1135,12 @@ def same(x, y):
         if new.type.dtype != out.type.dtype:
             new = cast(new, out.type.dtype)
 
-        if new.type != out.type:
+        if new.type.broadcastable != out.type.broadcastable:
             new = fill_chain(new, node.inputs)[0]
 
-        if new.type == out.type:
+        if (new.type.dtype == out.type.dtype) and (
+            new.type.broadcastable == out.type.broadcastable
+        ):
             new.tag.values_eq_approx = values_eq_approx_remove_inf_nan
             copy_stack_trace(out, new)
             return [new]

tests/tensor/rewriting/test_math.py

@@ -30,7 +30,7 @@
 from pytensor.misc.safe_asarray import _asarray
 from pytensor.printing import debugprint
 from pytensor.tensor import inplace
-from pytensor.tensor.basic import Alloc, join, switch
+from pytensor.tensor.basic import Alloc, join, second, switch
 from pytensor.tensor.blas import Dot22, Gemv
 from pytensor.tensor.blas_c import CGemv
 from pytensor.tensor.elemwise import CAReduce, DimShuffle, Elemwise
@@ -979,6 +979,28 @@ def test_mismatching_types(self):
         # No rewrite was applied
         assert z_rewritten is z
 
+    def test_specified_shape_by_constant(self):
+        x = vector("x")
+        const = np.full(shape=(5,), fill_value=2.0).astype(config.floatX)
+        out = x * const
+
+        new_out = rewrite_graph(
+            out, custom_rewrite=in2out(local_mul_canonizer, name="test")
+        )
+        expected_out = [2.0] * specify_shape(x, (5,))
+        assert equal_computations([new_out], [expected_out])
+
+    def test_broadcasted_by_constant(self):
+        x = vector("x")
+        const = np.full(shape=(3, 5), fill_value=2.0).astype(config.floatX)
+        out = x * const
+
+        new_out = rewrite_graph(
+            out, custom_rewrite=in2out(local_mul_canonizer, name="test")
+        )
+        expected_out = second(const, [[2.0]] * x)
+        assert equal_computations([new_out], [expected_out])
+
 
 def test_local_merge_abs():
     x, y, z = matrices("xyz")