diag_Omegas function

Simultaneously diagonalize two covariance matrices