Currently `EigenmodesSolver` requires the `C` damping matrix to be `None`. It should be extended to allow `C != None`.