Subspace methods for large inverse problems with multiple parameter classes