krylov subspace exponential time domain solution of maxwell ' s equations in photonic crystal modeling .
matrix computations ( 3 . ed . ) .