The present project is concerned with the efficient solution of large scale Riccati and Lyapunov equations. Such problems typically appear in state-feedback control of systems governed by partial differential equations. We apply a sparse grid discretization in order to arrive at a solution method which scales essentially linear in the number of unknowns.