@incollection{GotzesHeineckeHilleretal.2015, author = {Gotzes, Uwe and Heinecke, Nina and Hiller, Benjamin and R{\"o}vekamp, Jessica and Koch, Thorsten}, title = {Regulatory rules for gas markets in Germany and other European countries}, booktitle = {Evaluating gas network capacities}, publisher = {Society for Industrial and Applied Mathematics}, isbn = {978-1-611973-68-6}, pages = {45 -- 64}, year = {2015}, language = {en} } @misc{BreuerBussieckCaoetal.2017, author = {Breuer, Thomas and Bussieck, Michael and Cao, Karl-Kien and Cebulla, Felix and Fiand, Frederik and Gils, Hans Christian and Gleixner, Ambros and Khabi, Dmitry and Koch, Thorsten and Rehfeldt, Daniel and Wetzel, Manuel}, title = {Optimizing Large-Scale Linear Energy System Problems with Block Diagonal Structure by Using Parallel Interior-Point Methods}, issn = {1438-0064}, doi = {10.1007/978-3-319-89920-6_85}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66183}, year = {2017}, abstract = {Current linear energy system models (ESM) acquiring to provide sufficient detail and reliability frequently bring along problems of both high intricacy and increasing scale. Unfortunately, the size and complexity of these problems often prove to be intractable even for commercial state-of-the-art linear programming solvers. This article describes an interdisciplinary approach to exploit the intrinsic structure of these large-scale linear problems to be able to solve them on massively parallel high-performance computers. A key aspect are extensions to the parallel interior-point solver PIPS-IPM originally developed for stochastic optimization problems. Furthermore, a newly developed GAMS interface to the solver as well as some GAMS language extensions to model block-structured problems will be described.}, language = {en} } @misc{BorndoerferBreuerGrimmetal.2017, author = {Bornd{\"o}rfer, Ralf and Breuer, Matthias and Grimm, Boris and Reuther, Markus and Schade, Stanley and Schlechte, Thomas}, title = {Timetable Sparsification by Rolling Stock Rotation Optimization}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-65948}, year = {2017}, abstract = {Rolling stock optimization is a task that naturally arises by operating a railway system. It could be seen with different level of details. From a strategic perspective to have a rough plan which types of fleets to be bought to a more operational perspective to decide which coaches have to be maintained first. This paper presents a new approach to deal with rolling stock optimisation in case of a (long term) strike. Instead of constructing a completely new timetable for the strike period, we propose a mixed integer programming model that is able to choose appropriate trips from a given timetable to construct efficient tours of railway vehicles covering an optimized subset of trips, in terms of deadhead kilometers and importance of the trips. The decision which trip is preferred over the other is made by a simple evaluation method that is deduced from the network and trip defining data.}, language = {en} } @misc{RehfeldtKoch2018, author = {Rehfeldt, Daniel and Koch, Thorsten}, title = {SCIP-Jack—a solver for STP and variants with parallelization extensions: An update}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66416}, year = {2018}, abstract = {The Steiner tree problem in graphs is a classical problem that commonly arises in practical applications as one of many variants. Although the different Steiner tree problem variants are usually strongly related, solution approaches employed so far have been prevalently problem-specific. Against this backdrop, the solver SCIP-Jack was created as a general-purpose framework that can be used to solve the classical Steiner tree problem and 11 of its variants. This versatility is achieved by transforming various problem variants into a general form and solving them by using a state-of-the-art MIP-framework. Furthermore, SCIP-Jack includes various newly developed algorithmic components such as preprocessing routines and heuristics. The result is a high-performance solver that can be employed in massively parallel environments and is capable of solving previously unsolved instances. After the introduction of SCIP-Jack at the 2014 DIMACS Challenge on Steiner problems, the overall performance of the solver has considerably improved. This article provides an overview on the current state.}, language = {en} } @misc{WaltherHillerSaitenmacher2017, author = {Walther, Tom and Hiller, Benjamin and Saitenmacher, Ren{\´e}}, title = {Polyhedral 3D Models for compressors in gas networks}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-65864}, year = {2017}, abstract = {Compressor machines are crucial elements in a gas transmission network, required to compensate for the pressure loss caused by friction in the pipes. Modelling all physical and technical details of a compressor machine involves a large amount of nonlinearity, which makes it hard to use such models in the optimization of large-scale gas networks. In this paper, we are going to describe a modelling approach for the operating range of a compressor machine, starting from a physical reference model and resulting in a polyhedral representation in the 3D space of mass flow throughput as well as in- and outlet pressure.}, language = {en} } @phdthesis{Schweiger2017, author = {Schweiger, Jonas}, title = {Exploiting structure in non-convex quadratic optimization and gas network planning under uncertainty}, pages = {411}, year = {2017}, abstract = {The amazing success of computational mathematical optimization over the last decades has been driven more by insights into mathematical structures than by the advance of computing technology. In this vein, we address applications, where nonconvexity in the model and uncertainty in the data pose principal difficulties. The first part of the thesis deals with non-convex quadratic programs. Branch\&Bound methods for this problem class depend on tight relaxations. We contribute in several ways: First, we establish a new way to handle missing linearization variables in the well-known Reformulation-Linearization-Technique (RLT). This is implemented into the commercial software CPLEX. Second, we study the optimization of a quadratic objective over the standard simplex or a knapsack constraint. These basic structures appear as part of many complex models. Exploiting connections to the maximum clique problem and RLT, we derive new valid inequalities. Using exact and heuristic separation methods, we demonstrate the impact of the new inequalities on the relaxation and the global optimization of these problems. Third, we strengthen the state-of-the-art relaxation for the pooling problem, a well-known non-convex quadratic problem, which is, for example, relevant in the petrochemical industry. We propose a novel relaxation that captures the essential non-convex structure of the problem but is small enough for an in-depth study. We provide a complete inner description in terms of the extreme points as well as an outer description in terms of inequalities defining its convex hull (which is not a polyhedron). We show that the resulting valid convex inequalities significantly strengthen the standard relaxation of the pooling problem. The second part of this thesis focuses on a common challenge in real world applications, namely, the uncertainty entailed in the input data. We study the extension of a gas transport network, e.g., from our project partner Open Grid Europe GmbH. For a single scenario this maps to a challenging non-convex MINLP. As the future transport patterns are highly uncertain, we propose a robust model to best prepare the network operator for an array of scenarios. We develop a custom decomposition approach that makes use of the hierarchical structure of network extensions and the loose coupling between the scenarios. The algorithm used the single-scenario problem as black-box subproblem allowing the generalization of our approach to problems with the same structure. The scenario-expanded version of this problem is out of reach for today's general-purpose MINLP solvers. Yet our approach provides primal and dual bounds for instances with up to 256 scenarios and solves many of them to optimality. Extensive computational studies show the impact of our work.}, language = {en} } @misc{OPUS4-6530, title = {Open Access als strategische Leitlinie des KOBV}, organization = {Kooperativer Bibliotheksverbund Berlin-Brandenburg (KOBV)}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-65303}, pages = {2}, year = {2017}, abstract = {Open Access ist Grundlage und Voraussetzung f{\"u}r die Nachnutzbarkeit, Transparenz und freie Zug{\"a}nglichkeit zu wissenschaftlichen Informationen, Forschungsergebnissen und Daten des kulturellen Erbes. Mit Bezug auf die vom Senat Berlin 2015 verabschiedete Open-Access-Strategie empfiehlt der Kooperative Bibliotheksverbund Berlin-Brandenburg (KOBV) seinen Bibliotheken praktische Maßnahmen zur Umsetzung von Open Access in den eigenen Einrichtungen und benennt Maßnahmen der KOBV-Zentrale zur Unterst{\"u}tzung seiner Mitgliedsbibliotheken im Zuge dieses Transformationsprozesses.}, language = {de} } @incollection{CaimiFischerSchlechte2018, author = {Caimi, Gabrio and Fischer, Frank and Schlechte, Thomas}, title = {Railway Track Allocation}, volume = {268}, booktitle = {Handbook of Optimization in the Railway Industry}, publisher = {Springer International Publishing}, isbn = {978-3-319-72152-1}, doi = {https://doi.org/10.1007/978-3-319-72153-8}, pages = {141 -- 160}, year = {2018}, abstract = {This chapter addresses the classical task to decide which train runs on which track in a railway network. In this context a track allocation defines the precise routing of trains through a railway network, which usually has only a limited capacity. Moreover, the departure and arrival times at the visited stations of each train must simultaneously meet several operational and safety requirements. The problem to find the 'best possible' allocation for all trains is called the track allocation problem (TTP). Railway systems can be modeled on a very detailed scale covering the behavior of individual trains and the safety system to a large extent. However, those microscopic models are too big and not scalable to large networks, which make them inappropriate for mathematical optimization on a network wide level. Hence, most network optimization approaches consider simplified, so called macroscopic, models. In the first part we take a look at the challenge to construct a reliable and condensed macroscopic model for the associated microscopic model and to facilitate the transition between both models of different scale. In the main part we focus on the optimization problem for macroscopic models of the railway system. Based on classical graph-theoretical tools the track allocation problem is formulated to determine conflict-free paths in corresponding time-expanded graphs. We present standard integer programming model formulations for the track allocation problem that model resource or block conflicts in terms of packing constraints. In addition, we discuss the role of maximal clique inequalities and the concept of configuration networks. We will also present classical decomposition approaches like Lagrangian relaxation and bundle methods. Furthermore, we will discuss recently developed techniques, e.g., dynamic graph generation. Finally, we will discuss the status quo and show a vision of mathematical optimization to support real world track allocation, i.e. integrated train routing and scheduling, in a data-dominated and digitized railway future.}, language = {en} } @article{BertholdHendelKoch2017, author = {Berthold, Timo and Hendel, Gregor and Koch, Thorsten}, title = {From feasibility to improvement to proof: three phases of solving mixed-integer programs}, volume = {33}, journal = {Optimization Methods and Software}, number = {3}, publisher = {Taylor \& Francis}, doi = {10.1080/10556788.2017.1392519}, pages = {499 -- 517}, year = {2017}, abstract = {Modern mixed-integer programming (MIP) solvers employ dozens of auxiliary algorithmic components to support the branch-and-bound search in finding and improving primal solutions and in strengthening the dual bound. Typically, all components are tuned to minimize the average running time to prove optimality. In this article, we take a different look at the run of a MIP solver. We argue that the solution process consists of three distinct phases, namely achieving feasibility, improving the incumbent solution, and proving optimality. We first show that the entire solving process can be improved by adapting the search strategy with respect to the phase-specific aims using different control tunings. Afterwards, we provide criteria to predict the transition between the individual phases and evaluate the performance impact of altering the algorithmic behaviour of the non-commercial MIP solver Scip at the predicted phase transition points.}, language = {en} } @article{SchmidtAssmannBurlacuetal.2017, author = {Schmidt, Martin and Assmann, Denis and Burlacu, Robert and Humpola, Jesco and Joormann, Imke and Kanelakis, Nikolaos and Koch, Thorsten and Oucherif, Djamal and Pfetsch, Marc and Schewe, Lars and Schwarz, Robert and Sirvent, Matthias}, title = {GasLib - A Library of Gas Network Instances}, volume = {2}, journal = {Data}, number = {4}, doi = {10.3390/data2040040}, year = {2017}, language = {en} } @article{Koch2017, author = {Koch, Thorsten}, title = {The ZIMPL modeling language}, volume = {103}, journal = {Optima}, pages = {8 -- 9}, year = {2017}, language = {en} } @misc{RehfeldtKoch2017, author = {Rehfeldt, Daniel and Koch, Thorsten}, title = {Combining NP-Hard Reduction Techniques and Strong Heuristics in an Exact Algorithm for the Maximum-Weight Connected Subgraph Problem}, issn = {1438-0064}, doi = {10.1137/17M1145963}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64699}, year = {2017}, abstract = {Borne out of a surprising variety of practical applications, the maximum-weight connected subgraph problem has attracted considerable interest during the past years. This interest has not only led to notable research on theoretical properties, but has also brought about several (exact) solvers-with steadily increasing performance. Continuing along this path, the following article introduces several new algorithms such as reduction techniques and heuristics and describes their integration into an exact solver. The new methods are evaluated with respect to both their theoretical and practical properties. Notably, the new exact framework allows to solve common problem instances from the literature faster than all previous approaches. Moreover, one large-scale benchmark instance from the 11th DIMACS Challenge can be solved for the first time to optimality and the primal-dual gap for two other ones can be significantly reduced.}, language = {en} } @masterthesis{Saitenmacher2016, type = {Bachelor Thesis}, author = {Saitenmacher, Ren{\´e}}, title = {Combinatorial Models of Compressor Stations in Gas Networks}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63522}, year = {2016}, abstract = {I consider models of operation modes for complex compressor stations (i.e., ones with several in- or outlets) in gas networks. In particular, my aim is to identify from a set of operation modes those that are infeasible or redundant. For this, I investigate optimization and network reduction based procedures. Ultimately, these procedures may be employed for data preprocessing in the context of gas network optimization.}, language = {en} } @misc{SchadeBorndoerferBreueretal.2017, author = {Schade, Stanley and Bornd{\"o}rfer, Ralf and Breuer, Matthias and Grimm, Boris and Reuther, Markus and Schlechte, Thomas and Siebeneicher, Patrick}, title = {Pattern Detection For Large-Scale Railway Timetables}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63390}, year = {2017}, abstract = {We consider railway timetables of our industrial partner DB Fernverkehr AG that operates the ICE high speed trains in the long-distance passenger railway network of Germany. Such a timetable covers a whole year with 364 days and, typically, includes more than 45,000 trips. A rolling stock rotation plan is not created for the whole timetable at once. Instead the timetable is divided into regular invariant sections and irregular deviations (e.g. for public holidays). A separate rotation plan with a weekly period can then be provided for each of the different sections of the timetable. We present an algorithmic approach to automatically recognize these sections. Together with the supplementing visualisation of the timetable this method has shown to be very relevant for our industrial partner.}, language = {en} } @misc{GilgKlugMartienssenetal.2017, author = {Gilg, Brady and Klug, Torsten and Martienssen, Rosemarie and Paat, Joseph and Schlechte, Thomas and Schulz, Christof and Seymen, Sinan and Tesch, Alexander}, title = {Conflict-Free Railway Track Assignment at Depots}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63843}, year = {2017}, abstract = {Managing rolling stock with no passengers aboard is a critical component of railway operations. In particular, one problem is to park the rolling stock on a given set of tracks at the end of a day or service. Depending on the parking assignment, shunting may be required in order for a parked train to depart or for an incoming train to park. Given a collection of tracks M and a collection of trains T with fixed arrival-departure timetable, the train assignment problem (TAP) is to determine the maximum number of trains from T that can be parked on M according to the timetable and without the use of shunting. Hence, efficiently solving the TAP allows to quickly compute feasible parking schedules that do not require further shunting adjustments. In this paper, we present two integer programming models for solving the TAP. To our knowledge, this is the first integrated approach that considers track lengths along with the three most common types of parking tracks. We compare these models on a theoretical level. We also prove that a decision version of the TAP is NP-complete, justifying the use of integer programming techniques. Using stochastic and robust modelling techniques, both models produce parking assignments that are optimized and robust according to random train delays. We conclude with computational results for both models, observing that they perform well on real timetables.}, language = {en} } @article{RehfeldtKochMaher2019, author = {Rehfeldt, Daniel and Koch, Thorsten and Maher, Stephen J.}, title = {Reduction Techniques for the Prize-Collecting Steiner Tree Problem and the Maximum-Weight Connected Subgraph Problem}, volume = {73}, journal = {Networks}, edition = {2}, publisher = {Wiley}, doi = {10.1002/net.21857}, pages = {206 -- 233}, year = {2019}, abstract = {The concept of reduction has frequently distinguished itself as a pivotal ingredient of exact solving approaches for the Steiner tree problem in graphs. In this paper we broaden the focus and consider reduction techniques for three Steiner problem variants that have been extensively discussed in the literature and entail various practical applications: The prize-collecting Steiner tree problem, the rooted prize-collecting Steiner tree problem and the maximum-weight connected subgraph problem. By introducing and subsequently deploying numerous new reduction methods, we are able to drastically decrease the size of a large number of benchmark instances, already solving more than 90 percent of them to optimality. Furthermore, we demonstrate the impact of these techniques on exact solving, using the example of the state-of-the-art Steiner problem solver SCIP-Jack.}, language = {en} } @inproceedings{BorndoerferSagnolSwarat2012, author = {Bornd{\"o}rfer, Ralf and Sagnol, Guillaume and Swarat, Elmar}, title = {An IP Approach to Toll Enforcement Optimization on German Motorways}, booktitle = {Operations Research Proceedings 2011}, doi = {10.1007/978-3-642-29210-1_51}, pages = {317 -- 322}, year = {2012}, abstract = {This paper proposes the first model for toll enforcement optimization on German motorways. The enforcement is done by mobile control teams and our goal is to produce a schedule achieving network-wide control, proportional to spatial and time-dependent traffic distributions. Our model consists of two parts. The first plans control tours using a vehicle routing approach with profits and some side constraints. The second plans feasible rosters for the control teams. Both problems can be modeled as Multi-Commodity Flow Problems. Adding additional coupling constraints produces a large-scale integrated integer programming formulation. We show that this model can be solved to optimality for real world instances associated with a control area in East Germany.}, language = {en} } @inproceedings{BorndoerferBuwayaSagnoletal.2013, author = {Bornd{\"o}rfer, Ralf and Buwaya, Julia and Sagnol, Guillaume and Swarat, Elmar}, title = {Optimizing Toll Enforcement in Transportation Networks: a Game-Theoretic Approach}, volume = {41}, booktitle = {Proceedings of INOC'2013}, doi = {http://dx.doi.org//10.1016/j.endm.2013.05.100}, pages = {253 -- 260}, year = {2013}, abstract = {We present a game-theoretic approach to optimize the strategies of toll enforcement on a motorway network. In contrast to previous approaches, we consider a network with an arbitrary topology, and we handle the fact that users may choose their Origin-Destination path; in particular they may take a detour to avoid sections with a high control rate. We show that a Nash equilibrium can be computed with an LP (although the game is not zero-sum), and we give a MIP for the computation of a Stackelberg equilibrium. Experimental results based on an application to the enforcement of a truck toll on German motorways are presented.}, language = {en} } @inproceedings{BorndoerferOmontSagnoletal.2012, author = {Bornd{\"o}rfer, Ralf and Omont, Bertrand and Sagnol, Guillaume and Swarat, Elmar}, title = {A Stackelberg game to optimize the distribution of controls in transportation networks}, volume = {105}, booktitle = {Proceedings of the 3rd International Conference on Game Theory for Networks (GAMENETS 2012)}, doi = {http://dx.doi.org/10.1007/978-3-642-35582-0_17}, pages = {224 -- 235}, year = {2012}, abstract = {We propose a game theoretic model for the spatial distribution of inspectors on a transportation network. The problem is to spread out the controls so as to enforce the payment of a transit toll. We formulate a linear program to find the control distribution which maximizes the expected toll revenue, and a mixed integer program for the problem of minimizing the number of evaders. Furthermore, we show that the problem of finding an optimal mixed strategy for a coalition of \$N\$ inspectors can be solved efficiently by a column generation procedure. Finally, we give experimental results from an application to the truck toll on German motorways.}, language = {en} } @article{RehfeldtKoch2018, author = {Rehfeldt, Daniel and Koch, Thorsten}, title = {Transformations for the Prize-Collecting Steiner Tree Problem and the Maximum-Weight Connected Subgraph Problem to SAP}, volume = {36}, journal = {Journal of Computational Mathematics}, number = {3}, doi = {10.4208/jcm.1709-m2017-0002}, pages = {459 -- 468}, year = {2018}, abstract = {Transformations of Steiner tree problem variants have been frequently discussed in the literature. Besides allowing to easily transfer complexity results, they constitute a central pillar of exact state-of-the-art solvers for well-known variants such as the Steiner tree problem in graphs. In this paper transformations for both the prize-collecting Steiner tree problem and the maximum-weight connected subgraph problem to the Steiner arborescence problem are introduced for the first time. Furthermore, we demonstrate the considerable implications for practical solving approaches, including the computation of strong upper and lower bounds.}, language = {en} } @misc{Rehfeldt2015, type = {Master Thesis}, author = {Rehfeldt, Daniel}, title = {A Generic Approach to Solving the Steiner Tree Problem and Variants}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57817}, pages = {185}, year = {2015}, abstract = {Spawned by practical applications, numerous variations of the classical Steiner tree problem in graphs have been studied during the last decades. Despite the strong relationship between the different variants, solution approaches employed so far have been prevalently problem-specific. In contrast, we pursue a general-purpose strategy resulting in a solver able to solve both the classical Steiner tree problem and ten of its variants without modification. These variants include well-known problems such as the prize-collecting Steiner tree problem, the maximum-weight connected subgraph problem or the rectilinear minimum Steiner tree problem. Bolstered by a variety of new methods, most notably reduction techniques, our solver is not only of unprecedented versatility, but furthermore competitive or even superior to specialized state-of-the-art programs for several Steiner problem variants.}, language = {en} } @misc{Hoppmann2016, author = {Hoppmann, Heide}, title = {An Extended Formulation for the Line Planning Problem}, issn = {1438-0064}, doi = {10.1007/978-3-319-42902-1_2}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57467}, year = {2016}, abstract = {In this paper we present a novel extended formulation for the line planning problem that is based on what we call "configurations" of lines and frequencies. Configurations account for all possible options to provide a required transportation capacity on an infrastructure edge. The proposed configuration model is strong in the sense that it implies several facet-defining inequalities for the standard model: set cover, symmetric band, MIR, and multicover inequalities. These theoretical findings can be confirmed in computational results. Further, we show how this concept can be generalized to define configurations for subsets of edges; the generalized model implies additional inequalities from the line planning literature.}, language = {en} } @inproceedings{BorndoerferLamorgeseKlugetal.2015, author = {Bornd{\"o}rfer, Ralf and Lamorgese, Leonardo and Klug, Torsten and Mannino, Carlo and Reuther, Markus and Schlechte, Thomas}, title = {Recent Success Stories on Optimization of Railway Systems}, booktitle = {Proceedings of the IAROR conference RailTokyo}, year = {2015}, abstract = {Planning and operating railway transportation systems is an extremely hard task due to the combinatorial complexity of the underlying discrete optimization problems, the technical intricacies, and the immense size of the problem instances. Because of that, however, mathematical models and optimization techniques can result in large gains for both railway cus- tomers and operators, e.g., in terms of cost reductions or service quality improvements. In the last years a large and growing group of researchers in the OR community have devoted their attention to this domain devel- oping mathematical models and optimization approaches to tackle many of the relevant problems in the railway planning process. However, there is still a gap to bridge between theory and practice, with a few notable exceptions. In this paper we address three success stories, namely, long-term freight train routing (part I), mid-term rolling stock rotation planning (part II), and real-time train dispatching (part III). In each case, we describe real-life, successful implementations. We will dis- cuss the individual problem setting, survey the optimization literature, and focus on particular aspects addressed by the mathematical models. We demonstrate on concrete applications how mathematical optimization can support railway planning and operations. This gives proof that math- ematical optimization can support the planning of rolling stock resources. Thus, mathematical models and optimization can lead to a greater effi- ciency of railway operations and will serve as a powerful and innovative tool to meet recent challenges of the railway industry.}, language = {en} } @misc{Spiegel2015, type = {Master Thesis}, author = {Spiegel, Christoph}, title = {Approximating Primitive Integand Aircraft Performance}, pages = {84}, year = {2015}, abstract = {A commercial aircraft cannot freely use the available airspace but instead has to stick to a three-dimensional network of segments similar to a car in a road network. In this network it faces several variable and interdependent costs in the form of travel time, fuel consumption and overflight charges. These are also highly dependent on other factors such as weather conditions, aircraft performance, take-off time and weight as well as the changing availability of elements in the graph. This is further complicated by the distinction of separate flight phases that challenge the standard notion of a graph with predetermined nodes and arcs. Therefore when trying to find either a distance, fuel, time or cost minimal trajectory for a specific aircraft between an origin and a destination airport, one faces a very complex shortest path problem in the airway network graph that even for strong implifications is often NP-hard. This thesis will focus on exploring the costs that occur in this graph and that are associated with the aircraft performance. Here we will rely on actual performance and weather data supplied by Lufthansa Systems AG in Frankfurt and analyze whether they meet the requirements necessary for common algorithms such as the First-in, First-out property. Since it is vital for any shortest path algorithm to have a fast and accurate way of determining the costs in the graph, we will face two problems regarding the calculation of aircraft performance during cruise as well as the calculation of the so-called air distance. So far these problems have been approached by the industry with rudimentary approximative methods. We will reformulate them as initial value problems and try to find good approximations using both general Runge-Kutta methods as well as a novel approach which relies on finding piecewise linear approximations of some primitive integrals in pre-processing. Computations will show that these approaches deliver fast and accurate results.}, language = {en} } @inproceedings{GamrathMelchioriBertholdetal.2015, author = {Gamrath, Gerald and Melchiori, Anna and Berthold, Timo and Gleixner, Ambros and Salvagnin, Domenico}, title = {Branching on Multi-aggregated Variables}, volume = {9075}, booktitle = {Integration of AI and OR Techniques in Constraint Programming. CPAIOR 2015}, doi = {10.1007/978-3-319-18008-3_10}, pages = {141 -- 156}, year = {2015}, abstract = {In mixed-integer programming, the branching rule is a key component to a fast convergence of the branch-and-bound algorithm. The most common strategy is to branch on simple disjunctions that split the domain of a single integer variable into two disjoint intervals. Multi-aggregation is a presolving step that replaces variables by an affine linear sum of other variables, thereby reducing the problem size. While this simplification typically improves the performance of MIP solvers, it also restricts the degree of freedom in variable-based branching rules. We present a novel branching scheme that tries to overcome the above drawback by considering general disjunctions defined by multi-aggregated variables in addition to the standard disjunctions based on single variables. This natural idea results in a hybrid between variable- and constraint-based branching rules. Our implementation within the constraint integer programming framework SCIP incorporates this into a full strong branching rule and reduces the number of branch-and-bound nodes on a general test set of publicly available benchmark instances. For a specific class of problems, we show that the solving time decreases significantly.}, language = {en} } @inproceedings{KlugJunoszaSzaniawskiKwasiborskietal.2015, author = {Klug, Torsten and Junosza-Szaniawski, Konstanty and Kwasiborski, Slawomir and F{\"u}genschuh, Armin and Schlechte, Thomas}, title = {Fastest, Average and Quantile Schedule}, booktitle = {SOFSEM 2015: Theory and Practice of Computer Science}, publisher = {Springer Berlin Heidelberg}, doi = {10.1007/978-3-662-46078-8_17}, pages = {201 -- 216}, year = {2015}, abstract = {We consider problems concerning the scheduling of a set of trains on a single track. For every pair of trains there is a minimum headway, which every train must wait before it enters the track after another train. The speed of each train is also given. Hence for every schedule - a sequence of trains - we may compute the time that is at least needed for all trains to travel along the track in the given order. We give the solution to three problems: the fastest schedule, the average schedule, and the problem of quantile schedules. The last problem is a question about the smallest upper bound on the time of a given fraction of all possible schedules. We show how these problems are related to the travelling salesman problem. We prove NP-completeness of the fastest schedule problem, NP-hardness of quantile of schedules problem, and polynomiality of the average schedule problem. We also describe some algorithms for all three problems. In the solution of the quantile problem we give an algorithm, based on a reverse search method, generating with polynomial delay all Eulerian multigraphs with the given degree sequence and a bound on the number of such multigraphs. A better bound is left as an open question.}, language = {en} } @inproceedings{FischerSchlechte2015, author = {Fischer, Frank and Schlechte, Thomas}, title = {Comparing two dual relaxations of large scale train timetabling problems}, booktitle = {Proceedings of Conference on Advanced Systems in Public Transport 2015}, year = {2015}, abstract = {Railway transportation and in particular train timetabling is one of the basic and source application areas of combinatorial optimization and integer programming. We will discuss two well established modeling techniques for the train timetabling problem. In this paper we focus on one major ingredient - the bounding by dual relaxations. We compare two classical dual relaxations of large scale time expanded train timetabling problems - the Lagrangean Dual and Lagrangean Decomposition. We discuss the convergence behavior and show limitations of the Lagrangean Decomposition approach for a configuration based model. We introduce a third dualization approach to overcome those limitations. Finally, we present promising preliminary computational experiments that show that our new approach indeed has superior convergence properties.}, language = {en} } @article{BertholdHendel2014, author = {Berthold, Timo and Hendel, Gregor}, title = {Shift-and-Propagate}, volume = {21}, journal = {Journal of Heuristics}, number = {1}, doi = {10.1007/s10732-014-9271-0}, pages = {73 -- 106}, year = {2014}, abstract = {In recent years, there has been a growing interest in the design of general purpose primal heuristics for use inside complete mixed integer programming solvers. Many of these heuristics rely on an optimal LP solution, which may take a significant amount of time to find. In this paper, we address this issue by introducing a pre-root primal heuristic that does not require a previously found LP solution. This heuristic, named Shift-and-Propagate , applies domain propagation techniques to quickly drive a variable assignment towards feasibility. Computational experiments indicate that this heuristic is a powerful supplement to existing rounding and propagation heuristics.}, language = {en} } @article{BorndoerferHeismann2015, author = {Bornd{\"o}rfer, Ralf and Heismann, Olga}, title = {The hypergraph assignment problem}, volume = {15}, journal = {Discrete Optimization}, doi = {10.1016/j.disopt.2014.11.002}, pages = {15 -- 25}, year = {2015}, language = {en} } @inproceedings{BentzEckelLouchetetal.2011, author = {Bentz, Winfried and Eckel, Klaus and Louchet, Hadrien and Orlowski, Sebastian and Patzak, Erwin and Seibel, Julia and Werner, Axel and Wess{\"a}ly, Roland}, title = {Sch{\"a}tze heben bei der Planung von FTTx-Netzen}, volume = {227}, booktitle = {Breitbandversorgung in Deutschland}, publisher = {VDE-Verlag}, year = {2011}, language = {en} } @inproceedings{OrlowskiWernerWessaely2012, author = {Orlowski, Sebastian and Werner, Axel and Wess{\"a}ly, Roland}, title = {Estimating trenching costs in FTTx network planning}, booktitle = {Operations Research Proceedings 2011}, editor = {Klatte, Diethart and L{\"u}thi, Hans-Jakob and Schmedders, Karl}, publisher = {Springer Verlag}, pages = {89 -- 94}, year = {2012}, language = {en} } @article{BrandtBrandt1999, author = {Brandt, Andreas and Brandt, Manfred}, title = {A note on the stability of the many-queue head-of-the-line processor-sharing system with permanent customers}, volume = {32}, journal = {Queueing Syst.}, pages = {363 -- 381}, year = {1999}, language = {en} } @article{SanyalWernerZiegler2009, author = {Sanyal, Raman and Werner, Axel and Ziegler, G{\"u}nter}, title = {On Kalai's conjectures concerning centrally symmetric polytopes.}, volume = {41}, journal = {Discrete Comput. Geom.}, number = {2}, doi = {10.1007/s00454-008-9104-8}, pages = {183 -- 198}, year = {2009}, language = {en} } @inproceedings{BentzMartensOrlowskietal.2010, author = {Bentz, Winfried and Martens, Maren and Orlowski, Sebastian and Werner, Axel and Wess{\"a}ly, Roland}, title = {FTTx-PLAN}, volume = {220}, booktitle = {Breitbandversorgung in Deutschland}, publisher = {VDE-Verlag}, year = {2010}, language = {en} } @article{BrandtBrandt1999, author = {Brandt, Andreas and Brandt, Manfred}, title = {On a two-queue priority system with impatience and its application to a call center}, volume = {1}, journal = {Methodol. Comput. Appl. Probab.}, pages = {191 -- 210}, year = {1999}, language = {en} } @article{BrandtBrandt1994, author = {Brandt, Andreas and Brandt, Manfred}, title = {On the distribution of the number of packets in the fluid flow approximation of packet arrival streams}, volume = {17}, journal = {Queueing Syst.}, pages = {275 -- 315}, year = {1994}, language = {en} } @phdthesis{Werner2009, author = {Werner, Axel}, title = {Linear constraints on face numbers of polytopes}, year = {2009}, language = {en} } @phdthesis{Achterberg2009, author = {Achterberg, Tobias}, title = {Constraint Integer Programming}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11129}, year = {2009}, abstract = {This thesis introduces the novel paradigm of "constraint integer programming" (CIP), which integrates constraint programming (CP) and mixed integer programming (MIP) modeling and solving techniques. It is supplemented by the software SCIP, which is a solver and framework for constraint integer programming that also features SAT solving techniques. SCIP is freely available in source code for academic and non-commercial purposes. Our constraint integer programming approach is a generalization of MIP that allows for the inclusion of arbitrary constraints, as long as they turn into linear constraints on the continuous variables after all integer variables have been fixed. The constraints, may they be linear or more complex, are treated by any combination of CP and MIP techniques: the propagation of the domains by constraint specific algorithms, the generation of a linear relaxation and its solving by LP methods, and the strengthening of the LP by cutting plane separation. The current version of SCIP comes with all of the necessary components to solve mixed integer programs. In the thesis, we cover most of these ingredients and present extensive computational results to compare different variants for the individual building blocks of a MIP solver. We focus on the algorithms and their impact on the overall performance of the solver. In addition to mixed integer programming, the thesis deals with chip design verification, which is an important topic of electronic design automation. Chip manufacturers have to make sure that the logic design of a circuit conforms to the specification of the chip. Otherwise, the chip would show an erroneous behavior that may cause failures in the device where it is employed. An important subproblem of chip design verification is the property checking problem, which is to verify whether a circuit satisfies a specified property. We show how this problem can be modeled as constraint integer program and provide a number of problem-specific algorithms that exploit the structure of the individual constraints and the circuit as a whole. Another set of extensive computational benchmarks compares our CIP approach to the current state-of-the-art SAT methodology and documents the success of our method.}, language = {en} } @misc{BertholdHeinzVigerske2009, author = {Berthold, Timo and Heinz, Stefan and Vigerske, Stefan}, title = {Extending a CIP framework to solve MIQCPs}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11371}, number = {09-23}, year = {2009}, abstract = {This paper discusses how to build a solver for mixed integer quadratically constrained programs (MIQCPs) by extending a framework for constraint integer programming (CIP). The advantage of this approach is that we can utilize the full power of advanced MIP and CP technologies. In particular, this addresses the linear relaxation and the discrete components of the problem. For relaxation, we use an outer approximation generated by linearization of convex constraints and linear underestimation of nonconvex constraints. Further, we give an overview of the reformulation, separation, and propagation techniques that are used to handle the quadratic constraints efficiently. We implemented these methods in the branch-cut-and-price framework SCIP. Computational experiments indicates the potential of the approach.}, language = {en} } @misc{BorndoerferGroetschelJaeger2009, author = {Bornd{\"o}rfer, Ralf and Gr{\"o}tschel, Martin and Jaeger, Ulrich}, title = {Planning Problems in Public Transit}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11252}, number = {09-13}, year = {2009}, abstract = {Every day, millions of people are transported by buses, trains, and airplanes in Germany. Public transit (PT) is of major importance for the quality of life of individuals as well as the productivity of entire regions. Quality and efficiency of PT systems depend on the political framework (state-run, market oriented) and the suitability of the infrastructure (railway tracks, airport locations), the existing level of service (timetable, flight schedule), the use of adequate technologies (information, control, and booking systems), and the best possible deployment of equipment and resources (energy, vehicles, crews). The decision, planning, and optimization problems arising in this context are often gigantic and "scream" for mathematical support because of their complexity. This article sketches the state and the relevance of mathematics in planning and operating public transit, describes today's challenges, and suggests a number of innovative actions. The current contribution of mathematics to public transit is — depending on the transportation mode — of varying depth. Air traffic is already well supported by mathematics. Bus traffic made significant advances in recent years, while rail traffic still bears significant opportunities for improvements. In all areas of public transit, the existing potentials are far from being exhausted. For some PT problems, such as vehicle and crew scheduling in bus and air traffic, excellent mathematical tools are not only available, but used in many places. In other areas, such as rolling stock rostering in rail traffic, the performance of the existing mathematical algorithms is not yet sufficient. Some topics are essentially untouched from a mathematical point of view; e.g., there are (except for air traffic) no network design or fare planning models of practical relevance. PT infrastructure construction is essentially devoid of mathematics, even though enormous capital investments are made in this area. These problems lead to questions that can only be tackled by engineers, economists, politicians, and mathematicians in a joint effort. Among other things, the authors propose to investigate two specific topics, which can be addressed at short notice, are of fundamental importance not only for the area of traffic planning, should lead to a significant improvement in the collaboration of all involved parties, and, if successful, will be of real value for companies and customers: • discrete optimal control: real-time re-planning of traffic systems in case of disruptions, • model integration: service design in bus and rail traffic. Work on these topics in interdisciplinary research projects could be funded by the German ministry of research and education (BMBF), the German ministry of economics (BMWi), or the German science foundation (DFG).}, language = {en} } @misc{BrandtBrandt2009, author = {Brandt, Manfred and Brandt, Andreas}, title = {On sojourn times in M/GI systems under state-dependent processor sharing}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11366}, number = {09-22}, year = {2009}, abstract = {We consider a system with Poisson arrivals and i.i.d. service times. The requests are served according to the state-dependent processor sharing discipline, where each request receives a service capacity which depends on the actual number of requests in the system. The linear systems of PDEs describing the residual and attained sojourn times coincide for this system, which provides time reversibility including sojourn times for this system, and their minimal non negative solution gives the LST of the sojourn time \$V(\tau)\$ of a request with required service time \$\tau\$. For the case that the service time distribution is exponential in a neighborhood of zero, we derive a linear system of ODEs, whose minimal non negative solution gives the LST of \$V(\tau)\$, and which yields linear systems of ODEs for the moments of \$V(\tau)\$ in the considered neighborhood of zero. Numerical results are presented for the variance of \$V(\tau)\$. In case of an M/GI/2-PS system, the LST of \$V(\tau)\$ is given in terms of the solution of a convolution equation in the considered neighborhood of zero. For bounded from below service times, surprisingly simple expressions for the LST and variance of \$V(\tau)\$ in this neighborhood of zero are derived, which yield in particular the LST and variance of \$V(\tau)\$ in M/D/2-PS.}, language = {en} } @misc{LuceDuintjerTebbensLiesenetal.2009, author = {Luce, Robert and Duintjer Tebbens, Jurjen and Liesen, J{\"o}rg and Nabben, Robert and Gr{\"o}tschel, Martin and Koch, Thorsten and Schenk, Olaf}, title = {On the Factorization of Simplex Basis Matrices}, organization = {TU Berlin, Zuse Institute Berlin, University of Basel}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11392}, number = {09-24}, year = {2009}, abstract = {In the simplex algorithm, solving linear systems with the basis matrix and its transpose accounts for a large part of the total computation time. We investigate various methods from modern numerical linear algebra to improve the computation speed of the basis updates arising in LPs. The experiments are executed on a large real-world test set. The most widely used solution technique is sparse LU factorization, paired with an updating scheme that allows to use the factors over several iterations. Clearly, small number of fill-in elements in the LU factors is critical for the overall performance. Using a wide range of LPs we show numerically that after a simple permutation the non-triangular part of the basis matrix is so small, that the whole matrix can be factorized with (relative) fill-in close to the optimum. This permutation has been exploited by simplex practitioners for many years. But to our knowledge no systematic numerical study has been published that demonstrates the effective reduction to a surprisingly small non-triangular problem, even for large scale LPs. For the factorization of the non-triangular part most existing simplex codes use some variant of dynamic Markowitz pivoting, which originated in the late 1950s. We also show numerically that, in terms of fill-in and in the simplex context, dynamic Markowitz is quite consistently superior to other, more recently developed techniques.}, language = {en} } @misc{BorndoerferCardonha2009, author = {Bornd{\"o}rfer, Ralf and Cardonha, Carlos}, title = {A Set Partitioning Approach to Shunting}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11326}, number = {09-18}, year = {2009}, abstract = {The Vehicle Positioning Problem (VPP) is a classical combinatorial optimization problem in public transport planning. A number of models and approaches have been suggested in the literature, which work for small problems, but not for large ones. We propose in this article a novel set partitioning model and an associated column generation solution approach for the VPP. The model provides a tight linear description of the problem. The pricing problem, and hence the LP relaxation itself, can be solved in polynomial resp. pseudo-polynomial time for some versions of the problems.}, language = {en} } @misc{BrandtBrandt2009, author = {Brandt, Andreas and Brandt, Manfred}, title = {Approximations for the second moments of sojourn times in M/GI systems under state-dependent processor sharing}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11414}, number = {09-25}, year = {2009}, abstract = {We consider a system with Poisson arrivals and general service times, where the requests are served according to the State-Dependent Processor Sharing (SDPS) discipline (Cohen's generalized processor sharing discipline), where each request receives a service capacity which depends on the actual number of requests in the system. For this system, denoted by \$M/GI/SDPS\$, we derive approximations for the squared coefficients of variation of the conditional sojourn time of a request given its service time and of the unconditional sojourn time by means of two-moment fittings of the service times. The approximations are given in terms of the squared coefficients of variation of the conditional and unconditional sojourn time in related \$M/D/SDPS\$ and \$M/M/SDPS\$ systems, respectively. The numerical results presented for \$M/GI/m-PS\$ systems illustrate that the proposed approximations work well.}, language = {en} } @misc{Bley2009, author = {Bley, Andreas}, title = {On the Hardness of Finding Small Shortest Path Routing Conflicts}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11276}, number = {09-15}, year = {2009}, abstract = {Nowadays most data networks use shortest path protocols such as OSPF or IS-IS to route traffic. Given administrative routing lengths for the links of a network, all data packets are sent along shortest paths with respect to these lengths from their source to their destination. One of the most fundamental problems in planning shortest path networks is to decide whether a given set of routing paths forms a valid routing and, if this is not the case, to find a small subset of the given paths that cannot be shortest paths simultaneously for any routing lengths. In this paper we show that it is NP-hard to approximate the size of the smallest shortest path conflict by a factor less than 7/6.}, language = {en} } @misc{BertholdHeinzPfetsch2009, author = {Berthold, Timo and Heinz, Stefan and Pfetsch, Marc}, title = {Nonlinear pseudo-Boolean optimization: relaxation or propagation?}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11232}, number = {09-11}, year = {2009}, abstract = {Pseudo-Boolean problems lie on the border between satisfiability problems, constraint programming, and integer programming. In particular, nonlinear constraints in pseudo-Boolean optimization can be handled by methods arising in these different fields: One can either linearize them and work on a linear programming relaxation or one can treat them directly by propagation. In this paper, we investigate the individual strengths of these approaches and compare their computational performance. Furthermore, we integrate these techniques into a branch-and-cut-and-propagate framework, resulting in an efficient nonlinear pseudo-Boolean solver.}, language = {en} } @misc{BertholdFeydyStuckey2010, author = {Berthold, Timo and Feydy, Thibaut and Stuckey, Peter}, title = {Rapid Learning for Binary Programs}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11663}, number = {10-04}, year = {2010}, abstract = {Learning during search allows solvers for discrete optimization problems to remember parts of the search that they have already performed and avoid revisiting redundant parts. Learning approaches pioneered by the SAT and CP communities have been successfully incorporated into the SCIP constraint integer programming platform. In this paper we show that performing a heuristic constraint programming search during root node processing of a binary program can rapidly learn useful nogoods, bound changes, primal solutions, and branching statistics that improve the remaining IP search.}, language = {en} } @misc{HeinzStephanSchlechte2009, author = {Heinz, Stefan and Stephan, R{\"u}diger and Schlechte, Thomas}, title = {Solving Steel Mill Slab Problems with Branch and Price}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11260}, number = {09-14}, year = {2009}, abstract = {The steel mill slab design problem from the CSPLib is a binpacking problem that is motivated by an application of the steel industry and that has been widely studied in the constraint programming community. Recently, several people proposed new models and methods to solve this problem. A steel mill slab library was created which contains 380 instances. A closely related binpacking problem called multiple knapsack problem with color constraints, originated from the same industrial problem, were discussed in the integer programming community. In particular, a simple integer programming for this problem has been given by Forrest et al. [3]. The aim of this paper is to bring these different studies together. Moreover, we adopt the model of [3] for the steel mill slab problem. Using a state of the art integer program solver, this model is capable to solve all instances of the steel mill slab library, mostly in less than one second, to optimality. We improved, thereby, the solution value of 76 instances.}, language = {en} } @misc{BorndoerferNeumannPfetsch2009, author = {Bornd{\"o}rfer, Ralf and Neumann, Marika and Pfetsch, Marc}, title = {The Steiner Connectivity Problem}, issn = {1438-0064}, doi = {10.1007/s10107-012-0564-5}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11171}, number = {09-07}, year = {2009}, abstract = {The Steiner connectivity problem is a generalization of the Steiner tree problem. It consists in finding a minimum cost set of simple paths to connect a subset of nodes in an undirected graph. We show that polyhedral and algorithmic results on the Steiner tree problem carry over to the Steiner connectivity problem, namely, the Steiner cut and the Steiner partition inequalities, as well as the associated polynomial time separation algorithms, can be generalized. Similar to the Steiner tree case, a directed formulation, which is stronger than the natural undirected one, plays a central role.}, language = {en} } @misc{BrandtBrandt2009, author = {Brandt, Andreas and Brandt, Manfred}, title = {Insensitivity bounds for the moments of the sojourn times in M/GI systems under state-dependent processor sharing}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11092}, number = {09-02}, year = {2009}, abstract = {We consider a system with Poisson arrivals and i.i.d. service times and where the requests are served according to the state-dependent (Cohen's generalized) processor sharing discipline, where each request in the system receives a service capacity which depends on the actual number of requests in the system. For this system we derive asymptotically tight upper bounds for the moments of the conditional sojourn time of a request with given required service time. The bounds generalize corresponding results, recently given for the single-server processor sharing system by Cheung et al. and for the state-dependent processor sharing system with exponential service times by the authors. Analogous results hold for the waiting times.}, language = {en} } @misc{TorresTorresBorndoerferetal.2008, author = {Torres, Luis Miguel and Torres, Ramiro and Bornd{\"o}rfer, Ralf and Pfetsch, Marc}, title = {Line Planning on Paths and Tree Networks with Applications to the Quito Trolebus System (Extended Abstract)}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11076}, number = {08-53}, year = {2008}, abstract = {Line planning is an important step in the strategic planning process of a public transportation system. In this paper, we discuss an optimization model for this problem in order to minimize operation costs while guaranteeing a certain level of quality of service, in terms of available transport capacity. We analyze the problem for path and tree network topologies as well as several categories of line operation that are important for the Quito Trolebus system. It turns out that, from a computational complexity worst case point of view, the problem is hard in all but the most simple variants. In practice, however, instances based on real data from the Trolebus System in Quito can be solved quite well, and significant optimization potentials can be demonstrated.}, language = {en} }