The FEM calculation produces reasonable results, with RMSDs of the energy of 0.2-0.3 eV. However, it tends to underestimate all excitation energies (compared to the reference), in particular for the first four triplet states. It can be noted that in the direct calculation, the two sets of triplet states are marginally split (on the sub-meV order) in the expected 1:2:1 degeneracy pattern, whereas the FEM produces two sets of four degenerate triplet states. For the two sets of singlets, the degeneracy patterns and oscillator strengths are qualitatively correctly described by the FEM.