| Autor: |
Barth, Timothy J., Griebel, Michael, Nieminen, Risto M., Roose, Dirk, Schlick, Tamar, Langer, Ulrich, Discacciati, Marco, Keyes, David E., Widlund, Olof B., Zulehner, Walter, Gonçalves, Etereldes, Mathew, Tarek P., Sarkis, Markus, Schaerer, Christian E. |
| Zdroj: |
Domain Decomposition Methods in Science & Engineering XVII; 2008, p527-534, 8p |
| Abstrakt: |
We consider an elliptic optimal control problem in two dimensions, in which the control variable corresponds to the Neumann data on a boundary segment, and where the performance functional is regularized to ensure that the problem is well posed. A finite element discretization of this control problem yields a saddle point linear system, which can be reduced to a symmetric positive definite Hessian system for determining the control variables. We formulate a robust preconditioner for this reduced Hessian system, as a matrix product involving the discrete Neumann to Dirichlet map and a mass matrix, and show that it yields a condition number bound which is uniform with respect to the mesh size and regularization parameters. On a uniform grid, this preconditioner can be implemented using a fast sine transform. Numerical tests verify the theoretical bounds. [ABSTRACT FROM AUTHOR] |
| Databáze: |
Supplemental Index |
| Externí odkaz: |
|
