Adaptive wavelet methods solve(semi-)linear operator equations at optimal rates.With time evolutionary problems,the spaces can be equipped with products of temporal and spatial wavelet bases.Consequently,the evolutionary problem can be solved at a complexity of solving the stationary problem.In an adaptive wavelet scheme,the residual is used as an a posteriori error estimator.A reformulation as a system of first order enables a much more efficient approximate residual evaluation.