Use mT

Author: jessegrabowski

pytensor/tensor/slinalg.py

@@ -467,10 +467,10 @@ def solve_triangular(
     b_ndim = _default_b_ndim(b, b_ndim)
 
     if trans in [1, "T", True]:
-        a = a.T
+        a = a.mT
         lower = not lower
     if trans in [2, "C"]:
-        a = a.conj().T
+        a = a.conj().mT
         lower = not lower
 
     ret = Blockwise(
@@ -586,7 +586,7 @@ def solve(
     b_ndim = _default_b_ndim(b, b_ndim)
 
     if transposed:
-        a = a.T
+        a = a.mT
         lower = not lower
 
     return Blockwise(