Computing the Optical Properties of Large Systems