Compute for using orthogonal diagonalization