Fast calculation of modal interaction in large power systems using spectral decompositions of Gramians