krylov subspace exponential time domain solution of maxwell ' s equations in photonic crystal modeling .
computing matrix functions ∗ .