Consider an elliptic self-adjoint pseudodifferential operator A acting on m-columns of half-densities a closed manifold M, whose principal symbol is assumed to have simple eigenvalues. Relying basis projections commuting with A, we construct almost-unitary that diagonalizes modulo infinitely smoothing operator. We provide invariant algorithm for the computation its full symbol, as well explicit...