Krylov subspace methods for model order reduction of bilinear control systems