TY - GEN A1 - Götschel, Sebastian A1 - von Tycowicz, Christoph A1 - Polthier, Konrad A1 - Weiser, Martin T1 - Reducing Memory Requirements in Scientific Computing and Optimal Control N2 - In high accuracy numerical simulations and optimal control of time-dependent processes, often both many time steps and fine spatial discretizations are needed. Adjoint gradient computation, or post-processing of simulation results, requires the storage of the solution trajectories over the whole time, if necessary together with the adaptively refined spatial grids. In this paper we discuss various techniques to reduce the memory requirements, focusing first on the storage of the solution data, which typically are double precision floating point values. We highlight advantages and disadvantages of the different approaches. Moreover, we present an algorithm for the efficient storage of adaptively refined, hierarchic grids, and the integration with the compressed storage of solution data. T3 - ZIB-Report - 13-64 KW - optimal control KW - trajectory storage KW - mesh compression KW - compression Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42695 SN - 1438-0064 ER - TY - THES A1 - Lubkoll, Lars T1 - Optimal Control in Implant Shape Design KW - optimal control Y1 - 2011 ER - TY - GEN A1 - Schiela, Anton T1 - State constrained optimal control problems with states of low regularity N2 - We consider first order optimality conditions for state constrained optimal control problems. In particular we study the case where the state equation has not enough regularity to admit existence of a Slater point in function space. We overcome this difficulty by a special transformation. Under a density condition we show existence of Lagrange multipliers, which have a representation via measures and additional regularity properties. T3 - ZIB-Report - 08-24 KW - optimal control KW - state constraints KW - optimality conditions Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-509 SN - 1438-0064 ER - TY - GEN A1 - Weiser, Martin T1 - Interior Point Methods in Function Space N2 - A primal-dual interior point method for optimal control problems is considered. The algorithm is directly applied to the infinite dimensional problem. Existence and convergence of the central path are analyzed, and linear convergence of a short step pathfollowing method is established. T3 - ZIB-Report - 03-35 KW - interior point methods in function space KW - optimal control KW - complementarity functions Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7578 ER - 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 - THES A1 - Götschel, Sebastian T1 - Adaptive Lossy Trajectory Compression for Optimal Control of Parabolic PDEs N2 - Optimal control problems governed by nonlinear, time-dependent PDEs on three-dimensional spatial domains are an important tool in many fields, ranging from engineering applications to medicine. For the solution of such optimization problems, methods working on the reduced objective functional are often employed to avoid a full spatio-temporal discretization of the problem. The evaluation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, requiring the storage of a full 4D data set. If Newton-CG methods are used, two additional trajectories have to be stored. To get numerical results that are accurate enough, in many cases very fine discretizations in time and space are necessary, leading to a significant amount of data to be stored and transmitted to mass storage. This thesis deals with the development and analysis of methods for lossy compression of such finite element solutions. The algorithms are based on a change of basis to reduce correlations in the data, combined with quantization. This is achieved by transforming the finite element coefficient vector from the nodal to the hierarchical basis, followed by rounding the coefficients to a prescribed precision. Due to the inexact reconstruction, and thus inexact data for the adjoint equation, the error induced in the reduced gradient, and reduced Hessian, has to be controlled, to not impede convergence of the optimization. Accuracy requirements of different optimization methods are analyzed, and computable error estimates for the influence of lossy trajectory storage are derived. These tools are used to adaptively control the accuracy of the compressed data. The efficiency of the algorithms is demonstrated on several numerical examples, ranging from a simple linear, scalar equation to a semi-linear system of reaction-diffusion equations. In all examples considerable reductions in storage space and bandwidth requirements are achieved, without significantly influencing the convergence behavior of the optimization methods. Finally, to go beyond pointwise error control, the hierarchical basis transform can be replaced by more sophisticated wavelet transforms. Numerical experiments indicate that choosing suitable norms for error control allows higher compression factors. KW - optimal control KW - semi-linear parabolic PDEs KW - adjoint gradient computation KW - trajectory storage KW - lossy compression Y1 - 2015 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000098552 ER - TY - GEN A1 - Zhang, Wei A1 - Wang, Han A1 - Hartmann, Carsten A1 - Weber, Marcus A1 - Schütte, Christof T1 - Applications of the cross-entropy method to importance sampling and optimal control of diffusions N2 - We study the cross-entropy method for diffusions. One of the results is a versatile cross-entropy algorithm that can be used to design efficient importance sampling strategies for rare events or to solve optimal control problems. The approach is based on the minimization of a suitable cross-entropy functional, with a parametric family of exponentially tilted probability distributions. We illustrate the new algorithm with several numerical examples and discuss algorithmic issues and possible extensions of the method. T3 - ZIB-Report - 14-10 KW - important sampling KW - optimal control KW - cross-entropy method KW - rare events KW - change of measure Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-49720 SN - 1438-0064 ER - TY - GEN A1 - Sarich, Marco A1 - Banisch, Ralf A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Markov State Models for Rare Events in Molecular Dynamics N2 - Rare but important transition events between long lived states are a key feature of many molecular systems. In many cases the computation of rare event statistics by direct molecular dynamics (MD) simulations is infeasible even on the most powerful computers because of the immensely long simulation timescales needed. Recently a technique for spatial discretization of the molecular state space designed to help overcome such problems, so-called Markov State Models (MSMs), has attracted a lot of attention. We review the theoretical background and algorithmic realization of MSMs and illustrate their use by some numerical examples. Furthermore we introduce a novel approach to using MSMs for the efficient solution of optimal control problems that appear in applications where one desires to optimize molecular properties by means of external controls. T3 - ZIB-Report - 13-52 KW - rare events KW - Markov state models KW - long timescales KW - optimal control Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42420 SN - 1438-0064 ER - TY - GEN A1 - Schiela, Anton T1 - An Interior Point Method in Function Space for the Efficient Solution of State Constrained Optimal Control Problems N2 - We propose and analyse an interior point path-following method in function space for state constrained optimal control. Our emphasis is on proving convergence in function space and on constructing a practical path-following algorithm. In particular, the introduction of a pointwise damping step leads to a very efficient method, as verified by numerical experiments. T3 - ZIB-Report - 07-44 KW - interior point methods in function space KW - optimal control KW - state constraints Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10471 SN - 1438-0064 ER - TY - GEN A1 - Schiela, Anton T1 - Optimality Conditions for Convex State Constrained Optimal Control Problems with Discontinuous States N2 - We discuss first order optimality conditions for state constrained optimal control problems. Our concern is the treatment of problems, where the solution of the state equation is not known to be continuous, as in the case of boundary control in three space dimensions or optimal control with parabolic partial differential equations. We show existence of measure valued Lagrangian multipliers, which have just enough additional regularity to be applicable to all possibly discontinuous solutions of the state equation. T3 - ZIB-Report - 07-35 KW - optimal control KW - state constraints KW - optimality conditions Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10342 SN - 1438-0064 ER - TY - JOUR A1 - Borndörfer, Ralf A1 - Danecker, Fabian A1 - Weiser, Martin T1 - A Discrete-Continuous Algorithm for Free Flight Planning JF - Algorithms N2 - We propose a hybrid discrete-continuous algorithm for flight planning in free flight airspaces. In a first step, our DisCOptER method discrete-continuous optimization for enhanced resolution) computes a globally optimal approximate flight path on a discretization of the problem using the A* method. This route initializes a Newton method that converges rapidly to the smooth optimum in a second step. The correctness, accuracy, and complexity of the method are goverened by the choice of the crossover point that determines the coarseness of the discretization. We analyze the optimal choice of the crossover point and demonstrate the asymtotic superority of DisCOptER over a purely discrete approach. KW - shortest path KW - flight planning KW - free flight KW - discrete-continuous algorithm KW - optimal control KW - discrete optimization Y1 - 2020 U6 - https://doi.org/10.3390/a14010004 SN - 1438-0064 VL - 14 IS - 1 SP - 4 PB - MDPI ER -