Why is that?
If it’s not just a feature lacking temporarily, what alternatives are out there?  I assume torch.solve works in most cases but that doesn’t leverage the positive definitiveness of a matrix:
l = torch.cholesky(x)
x_inv_u = torch.cholesky_solve(u, l)

vs.
x_inv_u = torch.solve(u, x)