A Fast Subspace Decomposition Method for Bilinear Systems