Adaptive wavelet schemes for an elliptic control problem with Dirichlet boundary control

An adaptive algorithm based on wavelets is proposed for the fast numerical solution of control problems governed by elliptic boundary value problems with Dirichlet boundary control. A quadratic cost functional representing Sobolev norms of the state and a regularization in terms of the control is to be minimized subject to linear constraints in weak form. In particular, the constraints are formulated as a saddle point problem that allows to handle the varying boundary conditions explicitly. In the framework of (biorthogonal) wavelets, a representer for the functional is derived in terms of ₂-norms of wavelet expansion coefficients and the constraints are written in form of an ₂ automorphism. Standard techniques from optimization are then used to deduce the resulting first order necessary conditions as a (still infinite) system in ₂. Applying the machinery developed in [8,9] which has been extended to control problems in [14], an adaptive method is proposed which can be interpreted as an inexact gradient method for the control. In each iteration step, in turn the primal and the adjoint saddle point system are solved up to a prescribed accuracy by an adaptive iterative Uzawa algorithm for saddle point problems which has been proposed in [10]. Under these premises, it can be shown that the adaptive algorithm containing now three layers of iterations is asymptotically optimal. This means that the convergence rate achieved for computing the solution up to a desired target tolerance is asymptotically the same as the wavelet-best N-term approximation of the solution, and the total computational work is proportional to the number of computational unknowns.AMS subject classification 65K10, 65N99, 93B40Angela Kunoth: This work has been supported partly by the Deutsche Forschungsgemeinschaft (SFB 611) at the Universität Bonn and by the European Communitys Human Potential Programme under contract HPRN-CT-2002-00286 Breaking Complexity.