We are interested in a class of numerical schemes for the optimization of nonlinear hyperbolic partial differential equations. We present continuous and discretized relaxation schemes for scalar, one– conservation laws. We present numerical results on tracking typew problems with nonsmooth desired states and convergence results for higher–order spatial and temporal discretization schemes.