[Feature Request] Allow a diagonalize operation analogous to symeig + root decomposition

[wjmaddox]

Logging cornellius-gp/gpytorch#1395 here as well. It's probably more useful for LinearOperator.