feat: Add pt.multiply as an alias for pt.mul (NumPy compatibility)

pytensor/tensor/math.py

@@ -2936,6 +2936,9 @@ def mul(a, *other_terms):
     # see decorator for function body
 
 
+multiply = mul
+
+
 def variadic_mul(*args):
     """Mul that accepts arbitrary number of inputs, including zero or one."""
     if not args:
@@ -4330,6 +4333,7 @@ def nan_to_num(x, nan=0.0, posinf=None, neginf=None):
     "minimum",
     "mod",
     "mul",
+    "multiply",
     "nan_to_num",
     "neg",
     "neq",

tests/tensor/test_math.py

@@ -102,6 +102,7 @@
     minimum,
     mod,
     mul,
+    multiply,
     nan_to_num,
     neg,
     neq,
@@ -3972,3 +3973,13 @@ def test_median(ndim, axis):
 
     assert np.allclose(result_odd, expected_odd)
     assert np.allclose(result_even, expected_even)
+
+
+def test_multiply():
+    x = vector()
+    y = vector()
+    z = multiply(x, y)
+    f = function([x, y], z)
+    a = np.array([1.0, 2.0, 3.0])
+    b = np.array([4.0, 5.0, 6.0])
+    assert np.allclose(f(a, b), a * b)