All three control problems (P), (EB), (ED) defined in the previous sections
lead after suitable discretization to nonlinear finite-dimensional
optimization problems of the form
A discretization (P) of (P) is given in section 2 while the elliptic problems are assumed to be discretized as described in detail in [16]-[19]. symbolizes the state equation and boundary conditions while denotes both pointwise control and state constraints, the only constraints of inequality type prescribed above. Thus, alternatively, it can be written as
We state the well-known SSC for (4.1), assuming , , . Let be an admissible point satisfying the first-order necessary optimality conditions with associated Lagrange multipliers and . Let further
The point is a strict local minimizer if a exists such that, see, for example, [24]
After computing a solution an AMPL stub (or ) file is written as well as a file with the computed Lagrange multipliers. This allows to check the SSC (4.3) with the help of a Fortran, alternatively, a C or Matlab, program.
This program reads the files and verifies first the necessary first-order optimality conditions, the column regularity of and the strict complementarity. For this, it utilizes routines provided by AMPL which permit evaluation of the objective and constraint gradients. Next, it computes the QR decomposition of with the help of LAPACK routine DGEQRF. This use of dense numerical linear algebra limits the problem size and will in future work be replaced by an approach exploiting sparsity. AMPL also provides a routine to multiply the Hessian of the Lagrangian times a vector. This is called with the columns of and thus can be formed. Its eigenvalues are computed with LAPACK routine DSYEV and the smallest eigenvalue is determined.