Krylov subspace methods for model order reduction in computational electromagnetics