Let calc_l_r_roots() special-case eyemat matrices again.
Let excitations work for r not being a special matrix.