Search Results

Now showing 1 - 10 of 30
  • Item
    Joint dynamic probabilistic constraints with projected linear decision rules
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Guigues, Vincent; Henrion, René
    We consider multistage stochastic linear optimization problems combining joint dynamic probabilistic constraints with hard constraints. We develop a method for projecting decision rules onto hard constraints of wait-and-see type. We establish the relation between the original (infinite dimensional) problem and approximating problems working with projections from different subclasses of decision policies. Considering the subclass of linear decision rules and a generalized linear model for the underlying stochastic process with noises that are Gaussian or truncated Gaussian, we show that the value and gradient of the objective and constraint functions of the approximating problems can be computed analytically.
  • Item
    Task assignment, sequencing and path-planning in robotic welding cells
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Landry, Chantal; Welz, Wolfgang; Henrion, René; Hömberg, Dietmar; Skutella, Martin
    A workcell composed of a workpiece and several welding robots is considered. We are interested in minimizing the makespan in the workcell. Hence, one needs i) to assign tasks between the robots, ii) to do the sequencing of the tasks for each robot and iii) to compute the fastest collisionfree paths between the tasks. Up to now, task assignment and path-planning were always handled separately, the former being a typical Vehicle Routing Problem whereas the later is modelled using an optimal control problem. In this paper, we present a complete algorithm which combines discrete optimization techniques with collision detection and optimal control problems efficiently
  • Item
    Subdifferential characterization of probability functions under Gaussian distribution
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Hantoute, Abderrahim; Henrion, René; Pérez-Aros, Pedro
    Probability functions figure prominently in optimization problems of engineering. They may be nonsmooth even if all input data are smooth. This fact motivates the consideration of subdifferentials for such typically just continuous functions. The aim of this paper is to provide subdifferential formulae of such functions in the case of Gaussian distributions for possibly infinite-dimensional decision variables and nonsmooth (locally Lipschitzian) input data. These formulae are based on the spheric-radial decomposition of Gaussian random vectors on the one hand and on a cone of directions of moderate growth on the other. By successively adding additional hypotheses, conditions are satisfied under which the probability function is locally Lipschitzian or even differentiable.
  • Item
    Some remarks on stability of generalized equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2012) Henrion, René; Kruger, Alexander; Outrata, Jiří
    The paper concerns the computation of the graphical derivative and the regular (Fréchet) coderivative of the solution map to a class of generalized equations, where the multi-valued term amounts to the regular normal cone to a (possibly nonconvex) set given by C2 inequalities. Instead of the Linear Independence qualification condition, standardly used in this context, one assumes a combination of the Mangasarian-Fromovitz and the Constant Rank qualification conditions. On the basis of the obtained generalized derivatives, new optimality conditions for a class of mathematical programs with equilibrium constrains are derived, and a workable characterization of the isolated calmness of the considered solution map is provided.
  • Item
    A gradient formula for linear chance constraints under Gaussian distribution
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2012) Henrion, René; Möller, Andris
    We provide an explicit gradient formula for linear chance constraints under a (possibly singular) multivariate Gaussian distribution. This formula allows one to reduce the calculus of gradients to the calculus of values of the same type of chance constraints (in smaller dimension and with different distribution parameters). This is an important aspect for the numerical solution of stochastic optimization problems because existing efficient codes for e.g., calculating singular Gaussian distributions or regular Gaussian probabilities of polyhedra can be employed to calculate gradients at the same time. Moreover, the precision of gradients can be controlled by that of function values which is a great advantage over using finite difference approximations. Finally, higher order derivatives are easily derived explicitly. The use of the obtained formula is illustrated for an example of a transportation network with stochastic demands.
  • Item
    Lipschitz lower semicontinuity moduli for linear inequality systems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Cánovas, Maria Josefa; Gisbert, María Jesús; Henrion, René; Parra, Juan
    The paper is focussed on the Lipschitz lower semicontinuity of the feasible set mapping for linear (finite and infinite) inequality systems in three different perturbation frameworks: full, right-hand side and left-hand side perturbations. Inspired by [14], we introduce the Lipschitz lower semicontinuity-star as an intermediate notion between the Lipschitz lower semicontinuity and the well-known Aubin property. We provide explicit point-based formulae for the moduli (best constants) of all three Lipschitz properties in all three perturbation settings.
  • Item
    Critical objective size and calmness modulus in linear programming
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Cánovas, Maria J.; Henrion, René; Parra, Juan; Toledo, F. Javier
    This paper introduces the concept of critical objective size associated with a linear program in order to provide operative point-based formulas (only involving the nominal data, and not data in a neighborhood) for computing or estimating the calmness modulus of the optimal set (argmin) mapping under uniqueness of nominal optimal solution and perturbations of all coefficients. Our starting point is an upper bound on this modulus given in [4]. In this paper we prove that this upper bound is attained if and only if the norm of the objective function coefficient vector is less than or equal to the critical objective size. This concept also allows us to obtain operative lower bounds on the calmness modulus. We analyze in detail an illustrative example in order to xplore some strategies that can improve the referred upper and lower bounds.
  • Item
    Gradient formulae for nonlinear probabilistic constraints with Gaussian and aussian-like distributions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) van Ackooij, Wim; Henrion, René
    Probabilistic constraints represent a major model of stochastic optimization. A possible approach for solving probabilistically constrained optimization problems consists in applying nonlinear programming methods. In order to do so, one has to provide sufficiently precise approximations for values and gradients of probability functions. For linear probabilistic constraints under Gaussian distribution this can be successfully done by analytically reducing these values and gradients to values of Gaussian distribution functions and computing the latter, for instance, by Genz’ code. For nonlinear models one may fall back on the spherical-radial decomposition of Gaussian random vectors and apply, for instance, Deák’s sampling scheme for the uniform distribution on the sphere in order to compute values of corresponding probability functions. The present paper demonstrates how the same sampling scheme can be used in order to simultaneously compute gradients of these probability functions. More precisely, we prove a formula representing these gradients in the Gaussian case as a certain integral over the sphere again. Later, the result is extended to alternative distributions with an emphasis on the multivariate Student (or T-) distribution.
  • Item
    (Sub-) Gradient formulae for probability functions of random inequality systems under Gaussian distribution
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Ackooij, Wim van; Henrion, René
    We consider probability functions of parameter-dependent random inequality systems under Gaussian distribution. As a main result, we provide an upper estimate for the Clarke subdifferential of such probability functions without imposing compactness conditions. A constraint qualification ensuring continuous differentiability is formulated. Explicit formulae are derived from the general result in case of linear random inequality systems. In the case of a constant coefficient matrix an upper estimate for even the smaller Mordukhovich subdifferential is proven.
  • Item
    Conditioning of linear-quadratic two-stage stochastic optimization problems
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Emich, Konstantin; Henrion, René; Römisch, Werner
    In this paper a condition number for linear-quadratic two-stage stochastic optimization problems is introduced as the Lipschitz modulus of the multifunction assigning to a (discrete) probability distribution the solution set of the problem. Being the outer norm of the Mordukhovich coderivative of this multifunction, the condition number can be estimated from above explicitly in terms of the problem data by applying appropriate calculus rules. Here, a chain rule for the extended partial second-order subdifferential recently proved by Mordukhovich and Rockafellar plays a crucial role. The obtained results are illustrated for the example of two-stage stochastic optimization problems with simple recourse.