TY - GEN A1 - Schiela, Anton A1 - Weiser, Martin T1 - Superlinear Convergence of the Control Reduced Interior Point Method for PDE Constrained Optimization N2 - A thorough convergence analysis of the Control Reduced Interior Point Method in function space is performed. This recently proposed method is a primal interior point pathfollowing scheme with the special feature, that the control variable is eliminated from the optimality system. Apart from global linear convergence we show, that this method converges locally almost quadratically, if the optimal solution satisfies a function space analogue to a non-degeneracy condition. In numerical experiments we observe, that a prototype implementation of our method behaves in compliance with our theoretical results. T3 - ZIB-Report - 05-15 KW - interior point methods in function space KW - optimal control KW - superlinear convergence Y1 - 2005 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8490 ER - TY - GEN A1 - Weiser, Martin A1 - Schiela, Anton T1 - Function space interior point methods for PDE constrained optimization N2 - A primal-dual interior point method for optimal control problems with PDE constraints is considered. The algorithm is directly applied to the infinite dimensional problem. Existence and convergence of the central path are analyzed. Numerical results from an inexact continuation method applied to a model problem are shown. T3 - ZIB-Report - 04-27 KW - interior point methods in function space KW - optimal control KW - complementarity functions Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8027 ER - TY - GEN A1 - Weiser, Martin A1 - Gänzler, Tobias A1 - Schiela, Anton T1 - A Control Reduced Primal Interior Point Method for PDE Constrained Optimization N2 - A primal interior point method for control constrained optimal control problems with PDE constraints is considered. Pointwise elimination of the control leads to a homotopy in the remaining state and dual variables, which is addressed by a short step pathfollowing method. The algorithm is applied to the continuous, infinite dimensional problem, where discretization is performed only in the innermost loop when solving linear equations. The a priori elimination of the least regular control permits to obtain the required accuracy with comparable coarse meshes. Convergence of the method and discretization errors are studied, and the method is illustrated at two numerical examples. T3 - ZIB-Report - 04-38 KW - interior point methods in function space KW - optimal control KW - finite elements KW - discretization error Y1 - 2004 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8138 ER - TY - GEN A1 - Weiser, Martin A1 - Schiela, Anton A1 - Deuflhard, Peter T1 - Asymptotic Mesh Independence of Newton's Method Revisited N2 - The paper presents a new affine invariant theory on asymptotic mesh independence of Newton's method in nonlinear PDEs. Compared to earlier attempts, the new approach is both much simpler and more natural from the algorithmic point of view. The theory is exemplified at collocation methods for ODE boundary value problems and at finite element methods for elliptic PDE problems. T3 - ZIB-Report - 03-13 KW - mesh independence KW - nonlinear partial differential equations KW - Newton method KW - finite element method KW - collocation method Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7352 ER - TY - GEN A1 - Lubkoll, Lars A1 - Schiela, Anton A1 - Weiser, Martin T1 - An optimal control problem in polyconvex hyperelasticity N2 - We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented. T3 - ZIB-Report - 12-08 KW - polyconvex elasticity KW - implant design KW - optimal control Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-14745 SN - 1438-0064 ER - TY - GEN A1 - Deuflhard, Peter A1 - Schiela, Anton A1 - Weiser, Martin T1 - Mathematical Cancer Therapy Planning in Deep Regional Hyperthermia N2 - This paper surveys the required mathematics for a typical challenging problem from computational medicine, the cancer therapy planning in deep regional hyperthermia. In the course of many years of close cooperation with clinics, the medical problem gave rise to quite a number of subtle mathematical problems, part of which had been unsolved when the common project started. Efficiency of numerical algorithms, i.e. computational speed and monitored reliability, play a decisive role for the medical treatment. Off-the-shelf software had turned out to be not sufficient to meet the requirements of medicine. Rather, new mathematical theory as well as new numerical algorithms had to be developed. In order to make our algorithms useful in the clinical environment, new visualization software, a virtual lab, including 3D geometry processing of individual virtual patients had to be designed and implemented. Moreover, before the problems could be attacked by numerical algorithms, careful mathematical modelling had to be done. Finally, parameter identification and constrained optimization for the PDEs had to be newly analyzed and realized over the individual patient's geometry. Our new techniques had an impact on the specificity of the individual patients' treatment and on the construction of an improved hyperthermia applicator. T3 - ZIB-Report - 11-39 KW - hyperthermia therapy planning KW - mathematical modelling KW - Maxwell equations KW - state constraints KW - parameter identification Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-14204 SN - 1438-0064 ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Schiela, Anton T1 - Solving Optimal Control Problems with the Kaskade 7 Finite Element Toolbox N2 - This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study. T3 - ZIB-Report - 10-25 KW - partial differential equations KW - optimal control KW - finite elements KW - generic programming KW - adaptive methods Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11909 SN - 1438-0064 ER - TY - GEN A1 - Schiela, Anton A1 - Weiser, Martin T1 - Barrier methods for a control problem from hyperthermia treatment planning N2 - We consider an optimal control problem from hyperthermia treatment planning and its barrier regularization. We derive basic results, which lay the groundwork for the computation of optimal solutions via an interior point path-following method. Further, we report on a numerical implementation of such a method and its performance at an example problem. T3 - ZIB-Report - 09-36 KW - hyperthermia treatment planning KW - optimal control KW - barrier methods Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11567 SN - 1438-0064 ER - TY - JOUR A1 - Schiela, Anton A1 - Weiser, Martin T1 - Superlinear convergence of the Control Reduced Interior Point Method for PDE Constrained Optimization JF - Computational Optimization and Applications Y1 - 2008 VL - 39 IS - 3 SP - 369 EP - 393 ER - TY - GEN A1 - Deuflhard, Peter A1 - Kornhuber, Ralf A1 - Sander, Oliver A1 - Schiela, Anton A1 - Weiser, Martin ED - Deuflhard, Peter ED - Grötschel, Martin ED - Hömberg, Dietmar ED - Horst, Ulrich ED - Kramer, Jürg ED - Mehrmann, Volker ED - Polthier, Konrad ED - Schmidt, Frank ED - Schütte, Christof ED - Skutella, Martin ED - Sprekels, Jürgen T1 - Mathematics cures virtual patients T2 - MATHEON-Mathematics for Key Technologies Y1 - 2014 VL - 1 SP - 7 EP - 25 PB - European Mathematical Society ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Schiela, Anton ED - Dedner, A. ED - Flemisch, B. ED - Klöfkorn, R. T1 - Solving Optimal Control Problems with the Kaskade 7 Finite Element Toolbox T2 - Advances in DUNE N2 - This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study. Y1 - 2012 SP - 101 EP - 112 PB - Springer ER - TY - GEN A1 - Schiela, Anton A1 - Weiser, Martin ED - Diehl, M. ED - Glineur, F. ED - Jarlebring, E. ED - Michiels, W. T1 - Barrier methods for a control problem from hyperthermia treatment planning T2 - Recent Advances in Optimization and its Applications in Engineering (Proceedings of 14th Belgian-French-German Conference on Optimization 2009) N2 - We consider an optimal control problem from hyperthermia treatment planning and its barrier regularization. We derive basic results, which lay the groundwork for the computation of optimal solutions via an interior point path-following method. Further, we report on a numerical implementation of such a method and its performance at an example problem. Y1 - 2010 SP - 419 EP - 428 PB - Springer ER - TY - JOUR A1 - Deuflhard, Peter A1 - Schiela, Anton A1 - Weiser, Martin T1 - Mathematical Cancer Therapy Planning in Deep Regional Hyperthermia JF - Acta Numerica N2 - This paper surveys the required mathematics for a typical challenging problem from computational medicine, the cancer therapy planning in deep regional hyperthermia. In the course of many years of close cooperation with clinics, the medical problem gave rise to quite a number of subtle mathematical problems, part of which had been unsolved when the common project started. Efficiency of numerical algorithms, i.e. computational speed and monitored reliability, play a decisive role for the medical treatment. Off-the-shelf software had turned out to be not sufficient to meet the requirements of medicine. Rather, new mathematical theory as well as new numerical algorithms had to be developed. In order to make our algorithms useful in the clinical environment, new visualization software, a virtual lab, including 3D geometry processing of individual virtual patients had to be designed and implemented. Moreover, before the problems could be attacked by numerical algorithms, careful mathematical modelling had to be done. Finally, parameter identification and constrained optimization for the PDEs had to be newly analyzed and realized over the individual patient's geometry. Our new techniques had an impact on the specificity of the individual patients' treatment and on the construction of an improved hyperthermia applicator. Y1 - 2012 VL - 21 SP - 307 EP - 378 ER - TY - JOUR A1 - Lubkoll, Lars A1 - Schiela, Anton A1 - Weiser, Martin T1 - An optimal control problem in polyconvex hyperelasticity JF - SIAM J. Control Opt. N2 - We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented. Y1 - 2014 U6 - https://doi.org/10.1137/120876629 VL - 52 IS - 3 SP - 1403 EP - 1422 ER - TY - JOUR A1 - Lubkoll, Lars A1 - Schiela, Anton A1 - Weiser, Martin T1 - An affine covariant composite step method for optimization with PDEs as equality constraints JF - Optimization Methods and Software N2 - We propose a composite step method, designed for equality constrained optimization with partial differential equations. Focus is laid on the construction of a globalization scheme, which is based on cubic regularization of the objective and an affine covariant damped Newton method for feasibility. We show finite termination of the inner loop and fast local convergence of the algorithm. We discuss preconditioning strategies for the iterative solution of the arising linear systems with projected conjugate gradient. Numerical results are shown for optimal control problems subject to a nonlinear heat equation and subject to nonlinear elastic equations arising from an implant design problem in craniofacial surgery. Y1 - 2017 U6 - https://doi.org/10.1080/10556788.2016.1241783 VL - 32 IS - 5 SP - 1132 EP - 1161 ER - TY - JOUR A1 - Weiser, Martin A1 - Schiela, Anton T1 - Function space interior point methods for PDE constrained optimization JF - PAMM Y1 - 2004 VL - 4 IS - 1 SP - 43 EP - 46 ER - TY - JOUR A1 - Weiser, Martin A1 - Gänzler, Tobias A1 - Schiela, Anton T1 - A control reduced primal interior point method for a class of control constrained optimal control problems JF - Comput. Optim. Appl. Y1 - 2007 VL - 41 IS - 1 SP - 127 EP - 145 ER - TY - JOUR A1 - Weiser, Martin A1 - Schiela, Anton A1 - Deuflhard, Peter T1 - Asymptotic Mesh Independence of Newton’s Method Revisited JF - SIAM J. Num. Anal. Y1 - 2005 VL - 42 IS - 5 SP - 1830 EP - 1845 ER - TY - GEN A1 - Lubkoll, Lars A1 - Schiela, Anton A1 - Weiser, Martin T1 - An affine covariant composite step method for optimization with PDEs as equality constraints N2 - We propose a composite step method, designed for equality constrained optimization with partial differential equations. Focus is laid on the construction of a globalization scheme, which is based on cubic regularization of the objective and an affine covariant damped Newton method for feasibility. We show finite termination of the inner loop and fast local convergence of the algorithm. We discuss preconditioning strategies for the iterative solution of the arising linear systems with projected conjugate gradient. Numerical results are shown for optimal control problems subject to a nonlinear heat equation and subject to nonlinear elastic equations arising from an implant design problem in craniofacial surgery. T3 - ZIB-Report - 15-09 KW - composite step methods KW - cubic regularization KW - affine covariant KW - optimization with PDEs Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-53954 SN - 1438-0064 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Schiela, Anton A1 - Weiser, Martin T1 - Kaskade 7 - a Flexible Finite Element Toolbox JF - Computers and Mathematics with Applications N2 - Kaskade 7 is a finite element toolbox for the solution of stationary or transient systems of partial differential equations, aimed at supporting application-oriented research in numerical analysis and scientific computing. The library is written in C++ and is based on the \textsc{Dune} interface. The code is independent of spatial dimension and works with different grid managers. An important feature is the mix-and-match approach to discretizing systems of PDEs with different ansatz and test spaces for all variables. We describe the mathematical concepts behind the library as well as its structure, illustrating its use at several examples on the way. Y1 - 2021 U6 - https://doi.org/10.1016/j.camwa.2020.02.011 VL - 81 SP - 444 EP - 458 ER -