Scenario Reduction: Old and New Results
|
|
- Jade Booth
- 5 years ago
- Views:
Transcription
1 Scenario Reduction: Old and New Results W. Römisch Humboldt-University Berlin Institute of Mathematics Berlin, Germany (H. Heitsch, R. Henrion, C. Küchler) Workshop in honor of Jitka Dupačová Charles University Prague, November 20, 2014
2 Introduction Most approaches for solving stochastic programs of the form { } min f 0 (x, ξ)p (dξ) : x X Ξ with a probability measure P on Ξ R d and a (normal) integrand f 0, require numerical integration techniques, i.e., replacing the integral by some quadrature formula n f 0 (x, ξ)p (dξ) p i f 0 (x, ξ i ), Ξ where p i = P ({ξ i }), n i=1 p i = 1 and ξ i Ξ, i = 1,..., n. Since f 0 is often expensive to compute, the number n should be as small as possible. i=1
3 With v(p ) and S(P ) denoting the optimal value and solution set of the stochastic program, respectively, the following estimates are known v(p ) v(q) sup f 0 (x, ξ)(p Q)(dξ) x X Ξ ) S(Q) S(P ) + Ψ P (sup f 0 (x, ξ)(p Q)(dξ), x X where X is assumed to be compact, Q is a probability distribution approximating P and the function Ψ P is the inverse of the growth function of the objective near the solution set, i.e., Ψ 1 P (t) := inf { Ξ Ξ } f 0 (x, ξ)p (dξ) v(p ) : x X, d(x, S(P )) t. Hence, the distance d F with F := {f 0 (x, ) : x X} d F (P, Q) := sup f(ξ)(p Q)(dξ) f F becomes important when approximating P. Ξ
4 For given n N and for the special case p i = 1 n, i = 1,..., n, the best possible choice of elements ξ i Ξ, i = 1,..., n (scenarios), is obtained by minimizing sup x X Ξ f 0 (x, ξ)p (dξ) n 1 n f 0 (x, ξ i ), i=1 i.e., by solving the best approximation problem where min d F(P, Q) Q P n (Ξ) P n (Ξ) := {Q : Q is a uniform probability measure with n scenarios}. It may be reformulated as a semi-infinite program. and is known as optimal quantization of P with respect to the function class F. Such optimal quantization problems of probability measures are often extremely difficult to solve. Idea: Enlarging the class F to the class of all Lipschitz functions with a uniform constant. But, then min Q P n (Ξ) d F(P, Q) = O(n 1 d ).
5 If the functions f 0 (x, ) belong to mixed Sobolev spaces, then the convergence rate min Q P n (Ξ) d F(P, Q) = O(n 1+δ ) (δ (0, 0.5]) can be achieved by certain randomized Quasi-Monte Carlo methods. The convergence rate can be improved if the functions f F satisfy a higher degree of smoothness. Aim of the talk: Solving the best approximation problem for discrete probability measures P having many scenarios and for function classes F, which are relevant for two-stage stochastic programs (scenario reduction). Additional motivation: Scenario reduction methods may be important for generating scenario trees form multistage stochastic programs.
6 Comparison of N = 256 Monte Carlo Mersenne Twister points and randomly binary shifted Sobol points in dimension d = 500, projection (8,9)
7 Linear two-stage stochastic programs min { } c, x + Φ(q(ξ), h(ξ) T (ξ)x)p (dξ) : x X, Ξ where c R m, Ξ and X are polyhedral subsets of R d and R m, respectively, P is a probability measure on Ξ and the s m-matrix T ( ), the vectors q( ) R m and h( ) R s are affine functions of ξ. Furthermore, Φ and D denote the infimum function of the linear second-stage program and its dual feasibility set, i.e., Φ(u, t) := inf{ u, y :W y = t, y Y } ((u, t) R m R s ) D := {u R m : {z R s : W z u Y } }, where q(ξ) R m are the recourse costs, W is the s m recourse matrix, W the transposed of W and Y the polar cone to the polyhedral cone Y.
8 Theorem: (Walkup-Wets 69) The function Φ(, ) is finite and continuous on the polyhedral set D W (Y ). Furthermore, the function Φ(u, ) is piecewise linear convex on the polyhedral set W (Y ) for fixed u D, and Φ(, t) is piecewise linear concave on D for fixed t W (Y ). Assumptions: (A1) relatively complete recourse: for any (ξ, x) Ξ X, h(ξ) T (ξ)x W (Y ); (A2) dual feasibility: q(ξ) D holds for all ξ Ξ. (A3) existence of second moments: Ξ ξ 2 P (dξ) < +. Note that (A1) is satisfied if W (Y ) = R s (complete recourse). In general, (A1) and (A2) impose a condition on the support of P. Extensions to certain random recourse models, i.e., to W (ξ), exist.
9 Idea: Extend the class F such that it covers all two-stage models. Fortet-Mourier metrics: (as canonical distances for two-stage models) ζ r (P, Q) := sup f(ξ)(p Q)(dξ) : f F r (Ξ), Ξ where r 1 (r {1, 2} if W (ξ) W ) F r (Ξ) := {f : Ξ R : f(ξ) f( ξ) c r (ξ, ξ), ξ, ξ Ξ}, c r (ξ, ξ) := max{1, ξ r 1, ξ r 1 } ξ ξ (ξ, ξ Ξ). Proposition: (Rachev-Rüschendorf 98) If Ξ is bounded, ζ r may be reformulated as transportation problem { } ζ r (P, Q) = inf ĉ r (ξ, ξ)η(dξ, d ξ):π 1 η =P, π 2 η =Q, Ξ Ξ where ĉ r is a metric (reduced cost) with ĉ r c r and given by { n 1 } ĉ r (ξ, ξ) := inf c r (ξ li, ξ li+1 ) : n N, ξ li Ξ, ξ l1 = ξ, ξ ln = ξ. i=1
10 Scenario reduction We consider discrete distributions P with scenarios ξ i and probabilities p i, i = 1,..., N, and Q being supported by a given subset of scenarios ξ j, j J {1,..., N}, of P. Optimal reduction of a given scenario set J: The best approximation of P with respect to ζ r by such a distribution Q exists and is denoted by Q. It has the distance D J := ζ r (P, Q ) = min ζ r(p, Q) = p i min ĉr(ξ i, ξ j ) Q j J i J and the probabilities qj = p j + p i, j J, where i J j J j := {i J : j = j(i)} and j(i) arg min ĉr(ξ i, ξ j ), i J j J (optimal redistribution).
11 Determining the optimal index set J with prescribed cardinality N n is, however, a combinatorial optimization problem: (n-median problem) min {D J : J {1,..., N}, J = N n} Hence, the problem of finding the optimal set J for deleting scenarios is N P- hard and polynomial time algorithms are not available. First idea: Reformulation as linear mixed-integer program N j=1,j i N min n 1 p j x ij ĉ r (ξ i, ξ j ) x ij + y i = 1 i,j=1 (i = 1,..., N), s.t. N y i = n, i=1 x ij y i 0 x ij 1 (i, j = 1,..., N), y i {0, 1} (1,..., N). and application of standard software or of specialized algorithms. Solution: x ij = { mini J ĉ r (ξ i,ξ j ) nĉ r (ξ i,ξ j ) 0, else., i J, j J y i = { 1, i J 0, i J.
12 Fast reduction heuristics Second idea: Application of (randomized) greedy heuristics. Starting point (n = N 1): min p l min ĉ r (ξ l, ξ j ) l {1,...,N} j l Algorithm 1: (Backward reduction) Step [0]: J [0] :=. Step [i]: Step [N-n+1]: l i arg min l J [i 1] k J [i 1] {l} J [i] := J [i 1] {l i }. Optimal redistribution. p k min ĉ r (ξ k, ξ j ). j J [i 1] {l}
13 Starting point (n = 1): min u {1,...,N} N p k ĉ r (ξ k, ξ u ) k=1 Algorithm 2: (Forward selection) Step [0]: Step [i]: Step [n+1]: J [0] := {1,..., N}. u i arg min u J [i 1] J [i] := J [i 1] \ {u i }. k J [i 1] \{u} Optimal redistribution. p k min ĉ r (ξ k, ξ j ), j J [i 1] \{u}
14 Application: Optimization of gas transport in a huge transportation network including hundreds of gas delivery nodes. A stationary situation is considered; more than 8 years of hourly data available at all delivery nodes; multivariate probability distribution for the gas output in certain temperature classes is estimated; 10 3 samples based on randomized Quasi-Monte Carlo methods are generated and later reduced by scenario reduction to 50 scenarios Hourly mean daily power in kwh/h Mean daily temperature in C
15 Mixed-integer two-stage stochastic programs We consider min where Φ is given by { } c, x + Φ(q(ξ), h(ξ) T (ξ)x)p (dξ) : x X, Ξ Φ(u, t) := inf { u 1, y 1 + u 2, y 2 W 1y 1 + W 2 y 2 t y 1 R m 1 +, y 2 Z m 2 + for all pairs (u, t) R m 1+m 2 R r, and c R m, X is a closed subset of R m, Ξ a polyhedron in R s, W 1 Q r m 1, W 2 Q r m 2, and T (ξ) R r m, q(ξ) R m 1+m 2 and h(ξ) R r are affine functions of ξ, and P is a probability measure. We again assume (A1) for W = (W 1, W 2 ) (relatively complete recourse), (A2) (dual feasibility) and (A3). }
16 Example: (Schultz-Stougie-van der Vlerk 98) Stochastic multi-knapsack problem: min = max, m = 2, m 1 = 0, m 2 = 4, c = (1.5, 4), X = [ 5, 5] 2, h(ξ) = ξ, q(ξ) q = (16, 19, 23, 28), y i {0, 1}, i = 1, 2, 3, 4, P U(5, 5.5,..., 14.5, 15} (discrete) Second stage problem: MILP with variables and 882 constraints. T = ( ) W = ( )
17 The function Φ is well understood and the function class F r,b (Ξ) := {f1l B : f F r (Ξ), B B}, is relevant, where r {1, 2}, B is a class of (convex) polyhedra in Ξ and 1l B denotes the characteristic function of the set B. The class B contains all polyhedra of the form B = {ξ Ξ : h(ξ) T (ξ)x D}, where x X and D is a polyhedron in R s each of whose facets, i.e., (s 1)- dimensional faces, is parallel to a facet of the cone W 1 (R m 1 + ) or of the unit cube [0, 1] s. Hence, B is very problem-specific. Therefore, we consider alternatively the class of rectangular sets B rect = {I 1 I 2 I d : I j is a closed interval in R} covering the situation of pure integer programs and serving as heuristic for the general case.
18 Proposition: In case F = F r,brect (Ξ), the metric d F allows the estimates d F (P, Q) max{α Brect (P, Q), ζ r (P, Q)} ( ) d F (P, Q) C ζ r (P, Q) + α Brect (P, Q) s+1 1 where C is some constant only depending on Ξ and α Brect discrepancy given by is the rectangular α Brect (P, Q) := sup P (B) Q(B) B B rect If the set Ξ is bounded, even the estimate holds α Brect (P, Q) d F (P, Q) Cα Brect (P, Q) 1 s+1. Since α Brect has a stronger influence on d F than ζ r, we consider the composite distance d λ (P, Q) = λ α Brect (P, Q) + (1 λ) ζ r (P, Q) with λ [0, 1] close to 1.
19 Scenario reduction We consider again discrete distributions P with scenarios ξ i and probabilities p i, i = 1,..., N, and Q being supported by a subset of scenarios ξ j, j J {1,..., N}, of P with weights q j, j J, where J has cardinality N n. The problem of optimal scenario reduction consists in determining such a probability measure Q deviating from P as little as possible with respect to d λ. It can be written as min {d λ ( P, j J q jδ ξj ) J {1,..., N}, J = N n q j 0 j J, j J q j = 1 }. This optimization problem may be decomposed into an outer problem for determining the index set J and an inner problem for choosing the probabilities q j, j J.
20 From 10 4 Monte Carlo samples of a two-dimensional standard normal distribution 20 scenarios are selected that represent best approximations with respect to the first order Fortet-Mourier metric (left) and the Kolmogorov distance (right), i.e., the uniform distance of probability distribution functions.
21 To this end, we denote d(p, (J, q)) := d λ ( P, j J q jδ ξj ) S n := {q R n : q j 0, j J, j J q j = 1}. Then the optimal scenario reduction problem may be rewritten as min{min d(p, (J, q)) : J {1,..., N}, J = N n} J q S n with the inner problem (optimal redistribution) min{d(p, (J, q)) : q S n } for fixed index set J. The outer problem is a N P hard combinatorial optimization problem while the inner problem may be reformulated as a linear program. Again a reformulation as linear mixed-integer program is possible. An explicit formula for D J := min q Sn d(p, (J, q)) is no longer available, but the inner problem may be rewritten.
22 For B B rect we define the system of critical index sets I(B) by and write I rect := {I(B) = {i {1,..., N} : ξ i B} : B B rect } P (B) Q(B) = i I(B) p i j I(B)\J Then, the rectangular discrepancy between P and Q is α Brect (P, Q) = max I I rect p i i I Using the reduced system of critical index sets j I\J I rect(j) := {I \ J : I I rect }, every I I rect(j) is associated with a family ϕ(i ) I rect : q j q j.. ϕ(i ) := {I I rect : I = I \ J} (I I rect(j)).
23 With the quantities γ I := max I ϕ(i ) i I p i and γ I := min p i I ϕ(i ) i J we obtain D J as infimum of the linear program min λt α + (1 λ)t ζ t α, t ζ 0, q j 0, j J q j = 1, η i,j 0, i = 1,..., N, j J, t ζ i=1,...,n,j J ĉr(ξ i, ξ j )η i,j, (I I rect(j)), j J η i,j = p i, i = 1,..., N, N i=1 η i,j = q j, j J, j I q j t α γ I, I I rect(j) j I q j t α + γ I, I I rect(j) We have I rect(j) 2 n and, hence, the LP should be solvable at least for moderate values of n.
24 How to determine I rect(j), γ I and γ I? Observation: I rect(j), γ I and γ I are determined by those rectangles B R, each of whose facets contains an element of {ξ j : j J}, such that it can not be enlarged without changing its interior s intersection with {ξ j : j J}. The rectangles in R are called supporting. Non supporting rectangle (left) and supporting rectangle (right). The dots represent the remaining scenarios ξ j, j J.
25 Proposition: It holds that I rect(j)= B R{I {1,..., N}\J : j I {ξ j }={ξ j :j J} int B} and, for every I Irect(J), γ I = max {P (int B) : B R, j I {ξ j } = {ξ j : j J} int B} γ I = p i, i I where I := {i {1,..., N} : min j I ξ j,l ξ i,l max j I ξ j,l, l = 1,..., d}. Note that R ( ) n+2 d 2!
26 Numerical results Comparison of scenario reduction for a Fortet-Mourier metric and a composite distance including α rect.
27 Optimal redistribution: α Brect versus ζ scenarios chosen by Quasi Monte Carlo out of 1000 samples from the uniform distribution on [0, 1] 2 and optimal probabilities adjusted w.r.t. λα Brect + (1 λ)ζ 2 for λ = 1 (gray balls) and λ = 0.9 (black circles)
28 Optimal redistribution w.r.t. the rectangular discrepancy α Brect : d n=5 n=10 n=15 n=20 N= N= Running times [sec] of the optimal redistribution algorithm The majority of the running time is spent for determining the supporting rectangles, while the time needed to solve the linear program is insignificant.
29 Optimal scenario reduction Forward selection: Step [0]: J [0] :=. Step [i]: l i argmin l J [i 1] inf q S i d λ J [i] := J [i 1] {l i }. Step [n+1]: Minimize d λ ( P, j J [n] q jδ ξj ) (P, ) q jδ j J[i 1] {l} ξj, s.t. q S n. N=100 n=5 n=10 n=15 d = d = d = Growth of running times (in seconds) of forward selection for λ = 1 Search for more efficient heuristics
30 Alternative heuristics (for P with independent marginals): (next neighbor) Quasi Monte Carlo: The first n numbers of randomized QMC sequences provide n equally weighted points. The closest scenarios are determined and the resulting discrepancy to the initial measure is computed for fixed probability weights. (next neighbor) adjusted Quasi Monte Carlo: The probabilities of the closest scenarios are adjusted by the optimal redistribution algorithm to obtain a minimal rectangular discrepancy to P. For general distributions P with densities transformation formulas are needed (e.g. Hlawka-Mück 71).
31 Conclusion: (Next neighbor) readjusted QMC decreases significantly the approximation error. Forward selection provides good results, but is very slow due to the optimal redistribution in each step. 0.8 discrepancy time Left: The distance d λ (λ = 1) between P and uniform (next neighbor) QMC points (dashed line) and (next neighbor) readjusted QMC points (solid line), and running time in seconds of optimal redistribution. Right: Distances α rect (solid) and ζ 2 (dashed) of 10 out of 100 scenarios, resulting from forward selection for several λ [0, 1].
32 Conclusions and outlook There exist reasonably fast heuristics for scenario reduction in linear twostage stochastic programs, It may be worth to study and compare exact solution methods with heuristics, Recursive application of the heuristics apply to generating scenario trees for multistage stochastic programs, For scenario tree reduction the heuristics have to be modified. For mixed-integer two-stage stochastic programs heuristics exist, but have to be based on different arguments. They are more expensive and so far restricted to moderate dimensions. This motivates to study exact approaches. There is hope for generating scenario trees for mixed-integer multistage models, but it is not yet supported by stability results.
33 References Avella, P., Sassano, A, Vasil ev, I.: Computational study of large-scale p-median problems, Mathematical Programming 109 (2007), Dupačová, J.; Gröwe-Kuska, N.; Römisch, W.: Scenario reduction in stochastic programming: An approach using probability metrics, Mathematical Programming 95 (2003), Heitsch, H., Römisch, W.: Scenario reduction algorithms in stochastic programming, Computational Optimization and Applications 24 (2003), Heitsch, H., Römisch, W.: A note on scenario reduction for two-stage stochastic programs, Operations Research Letters 35 (2007), Heitsch, H., Römisch, W.: Scenario tree modeling for multistage stochastic programs, Mathematical Programming 118 (2009), Heitsch, H., Römisch, W.: Scenario tree reduction for multistage stochastic programs, Computational Management Science 6 (2009), Henrion, R., Küchler, C., Römisch, W.: Scenario reduction in stochastic programming with respect to discrepancy distances, Computational Optimization and Applications 43 (2009), Henrion, R., Küchler, C., Römisch, W.: Discrepancy distances and scenario reduction in two-stage stochastic mixedinteger programming, Journal of Industrial and Management Optimization 4 (2008), Li, Z., Floudas, C. A.: Optimal scenario reduction framework based on distance of uncertainty distribution and output performance:i. Single reduction via mixed integer linear optimization, Computers and Chemical Engineering 70 (2014),
34 Römisch, W.: Stability of Stochastic Programming Problems, in: Stochastic Programming (A. Ruszczyński and A. Shapiro eds.), Handbooks in Operations Research and Management Science, Volume 10, Elsevier, Amsterdam 2003, Römisch, W., Vigerske, S.: Quantitative stability of fully random mixed-integer two-stage stochastic programs, Optimization Letters 2 (2008), Römisch, W.; Wets, R. J-B: Stability of ε-approximate solutions to convex stochastic programs, SIAM Journal on Optimization 18 (2007), Koch, T., Hiller, B., Pfetsch, M. E., Schewe, L. (Eds.): Evaluating Gas Network Capacities, SIAM-MOS (in preparation).
Jitka Dupačová and scenario reduction
Jitka Dupačová and scenario reduction W. Römisch Humboldt-University Berlin Institute of Mathematics http://www.math.hu-berlin.de/~romisch Session in honor of Jitka Dupačová ICSP 2016, Buzios (Brazil),
More informationStability of Stochastic Programming Problems
Stability of Stochastic Programming Problems W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch Page 1 of 35 Spring School Stochastic
More informationScenario generation in stochastic programming with application to energy systems
Scenario generation in stochastic programming with application to energy systems W. Römisch Humboldt-University Berlin Institute of Mathematics www.math.hu-berlin.de/~romisch SESO 2017 International Thematic
More informationA note on scenario reduction for two-stage stochastic programs
A note on scenario reduction for two-stage stochastic programs Holger Heitsch a and Werner Römisch a a Humboldt-University Berlin, Institute of Mathematics, 199 Berlin, Germany We extend earlier work on
More informationScenario generation in stochastic programming
Scenario generation in stochastic programming Werner Römisch Humboldt-University Berlin, Institute of Mathematics, 10099 Berlin, Germany, romisch@math.hu-berlin.de Abstract Stability-based methods for
More informationScenario generation Contents
Scenario generation W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Page 1 of 39 Tutorial, ICSP12, Halifax, August 14, 2010 (1) Aproximation issues (2) Generation
More informationStochastic optimization, multivariate numerical integration and Quasi-Monte Carlo methods. W. Römisch
Stochastic optimization, multivariate numerical integration and Quasi-Monte Carlo methods W. Römisch Humboldt-University Berlin Institute of Mathematics www.math.hu-berlin.de/~romisch Chemnitzer Mathematisches
More informationStochastic Programming: From statistical data to optimal decisions
Stochastic Programming: From statistical data to optimal decisions W. Römisch Humboldt-University Berlin Department of Mathematics (K. Emich, H. Heitsch, A. Möller) Page 1 of 24 6th International Conference
More informationOptimal scenario generation and reduction in stochastic programming
Noname manuscript No. (will be inserted by the editor) Optimal scenario generation and reduction in stochastic programming R. Henrion W. Römisch the date of receipt and acceptance should be inserted later
More informationDISCREPANCY DISTANCES AND SCENARIO REDUCTION IN TWO-STAGE STOCHASTIC MIXED-INTEGER PROGRAMMING. René Henrion. Christian Küchler.
JOURNAL OF INDUSTRIAL AND Website: http://aimsciences.org MANAGEMENT OPTIMIZATION Volume 4, Number 2, May 2008 pp. 363 384 DISCREPANCY DISTANCES AND SCENARIO REDUCTION IN TWO-STAGE STOCHASTIC MIXED-INTEGER
More informationStochastic Programming: Tutorial Part II
Stochastic Programming: Tutorial Part II W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch Page 1 of 77 Optimization: Computation,
More informationApproximation of stochastic optimization problems and scenario generation
Approximation of stochastic optimization problems and scenario generation W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch Page
More informationStability-based generation of scenario trees for multistage stochastic programs
Stability-based generation of scenario trees for multistage stochastic programs H. Heitsch and W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch
More informationStability of optimization problems with stochastic dominance constraints
Stability of optimization problems with stochastic dominance constraints D. Dentcheva and W. Römisch Stevens Institute of Technology, Hoboken Humboldt-University Berlin www.math.hu-berlin.de/~romisch SIAM
More informationScenario Generation in Stochastic Programming with Application to Optimizing Electricity Portfolios under Uncertainty
Scenario Generation in Stochastic Programming with Application to Optimizing Electricity Portfolios under Uncertainty W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch
More informationANOVA decomposition of convex piecewise linear functions
ANOVA decomposition of convex piecewise linear functions W. Römisch Abstract Piecewise linear convex functions arise as integrands in stochastic programs. They are Lipschitz continuous on their domain,
More informationProbability metrics and the stability of stochastic programs with recourse
Probability metrics and the stability of stochastic programs with recourse Michal Houda Charles University, Faculty of Mathematics and Physics, Prague, Czech Republic. Abstract In stochastic optimization
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 Are Quasi-Monte Carlo algorithms efficient for two-stage stochastic
More informationScenario estimation and generation
October 10, 2004 The one-period case Distances of Probability Measures Tensor products of trees Tree reduction A decision problem is subject to uncertainty Uncertainty is represented by probability To
More informationQuantitative Stability of Two Stage Linear Second-order Conic Stochastic Programs with Full Random Recourse
Quantitative Stability of Two Stage Linear Second-order Conic Stochastic Programs with Full Random Recourse Qingsong Duan, Mengwei Xu, Shaoyan Guo, and Liwei Zhang October 7, 2017 Abstract In this paper,
More informationStochastic Integer Programming An Algorithmic Perspective
Stochastic Integer Programming An Algorithmic Perspective sahmed@isye.gatech.edu www.isye.gatech.edu/~sahmed School of Industrial & Systems Engineering 2 Outline Two-stage SIP Formulation Challenges Simple
More informationConditioning of linear-quadratic two-stage stochastic programming problems
Conditioning of linear-quadratic two-stage stochastic programming problems W. Römisch Humboldt-University Berlin Institute of Mathematics http://www.math.hu-berlin.de/~romisch (K. Emich, R. Henrion) (WIAS
More informationProgress in high-dimensional numerical integration and its application to stochastic optimization
Progress in high-dimensional numerical integration and its application to stochastic optimization W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Short course,
More informationGeneration and Representation of Piecewise Polyhedral Value Functions
Generation and Representation of Piecewise Polyhedral Value Functions Ted Ralphs 1 Joint work with Menal Güzelsoy 2 and Anahita Hassanzadeh 1 1 COR@L Lab, Department of Industrial and Systems Engineering,
More informationLipschitz and differentiability properties of quasi-concave and singular normal distribution functions
Ann Oper Res (2010) 177: 115 125 DOI 10.1007/s10479-009-0598-0 Lipschitz and differentiability properties of quasi-concave and singular normal distribution functions René Henrion Werner Römisch Published
More informationDistributionally Robust Discrete Optimization with Entropic Value-at-Risk
Distributionally Robust Discrete Optimization with Entropic Value-at-Risk Daniel Zhuoyu Long Department of SEEM, The Chinese University of Hong Kong, zylong@se.cuhk.edu.hk Jin Qi NUS Business School, National
More informationA Benders Algorithm for Two-Stage Stochastic Optimization Problems With Mixed Integer Recourse
A Benders Algorithm for Two-Stage Stochastic Optimization Problems With Mixed Integer Recourse Ted Ralphs 1 Joint work with Menal Güzelsoy 2 and Anahita Hassanzadeh 1 1 COR@L Lab, Department of Industrial
More informationScenario Reduction for Risk-Averse Stochastic Programs
Scenario Reduction for Risk-Averse Stochastic Programs Sebastián Arpón 1, Tito Homem-de-Mello 1, Bernardo Pagnoncelli 1 1 School of Business, Universidad Adolfo Ibañez, Santiago, Chile Abstract In this
More informationReformulation of chance constrained problems using penalty functions
Reformulation of chance constrained problems using penalty functions Martin Branda Charles University in Prague Faculty of Mathematics and Physics EURO XXIV July 11-14, 2010, Lisbon Martin Branda (MFF
More information3.7 Strong valid inequalities for structured ILP problems
3.7 Strong valid inequalities for structured ILP problems By studying the problem structure, we can derive strong valid inequalities yielding better approximations of conv(x ) and hence tighter bounds.
More informationAsymptotics of minimax stochastic programs
Asymptotics of minimax stochastic programs Alexander Shapiro Abstract. We discuss in this paper asymptotics of the sample average approximation (SAA) of the optimal value of a minimax stochastic programming
More informationOn the Scenario-Tree Optimal-Value Error for Stochastic Programming Problems
On the Scenario-Tree Optimal-Value Error for Stochastic Programming Problems Julien Keutchayan David Munger Michel Gendreau January 2017 CIRRELT-2017-05 On the Scenario-Tree Optimal-Value Error for Stochastic
More informationConvergence of Stationary Points of Sample Average Two-Stage Stochastic Programs: A Generalized Equation Approach
MATHEMATICS OF OPERATIONS RESEARCH Vol. 36, No. 3, August 2011, pp. 568 592 issn 0364-765X eissn 1526-5471 11 3603 0568 doi 10.1287/moor.1110.0506 2011 INFORMS Convergence of Stationary Points of Sample
More informationOptimization and Optimal Control in Banach Spaces
Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,
More informationOptimization Tools in an Uncertain Environment
Optimization Tools in an Uncertain Environment Michael C. Ferris University of Wisconsin, Madison Uncertainty Workshop, Chicago: July 21, 2008 Michael Ferris (University of Wisconsin) Stochastic optimization
More informationCitation for published version (APA): van der Vlerk, M. H. (1995). Stochastic programming with integer recourse [Groningen]: University of Groningen
University of Groningen Stochastic programming with integer recourse van der Vlerk, Maarten Hendrikus IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to
More informationCONVERGENCE ANALYSIS OF SAMPLING-BASED DECOMPOSITION METHODS FOR RISK-AVERSE MULTISTAGE STOCHASTIC CONVEX PROGRAMS
CONVERGENCE ANALYSIS OF SAMPLING-BASED DECOMPOSITION METHODS FOR RISK-AVERSE MULTISTAGE STOCHASTIC CONVEX PROGRAMS VINCENT GUIGUES Abstract. We consider a class of sampling-based decomposition methods
More informationBranch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems
Branch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems Yongjia Song James R. Luedtke August 9, 2012 Abstract We study solution approaches for the design of reliably
More informationComplexity of two and multi-stage stochastic programming problems
Complexity of two and multi-stage stochastic programming problems A. Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA The concept
More informationSemi-Infinite Relaxations for a Dynamic Knapsack Problem
Semi-Infinite Relaxations for a Dynamic Knapsack Problem Alejandro Toriello joint with Daniel Blado, Weihong Hu Stewart School of Industrial and Systems Engineering Georgia Institute of Technology MIT
More informationSupplementary Material
ec1 Supplementary Material EC.1. Measures in D(S, µ, Ψ) Are Less Dispersed Lemma EC.1. Given a distributional set D(S, µ, Ψ) with S convex and a random vector ξ such that its distribution F D(S, µ, Ψ),
More informationSubmodularity in Machine Learning
Saifuddin Syed MLRG Summer 2016 1 / 39 What are submodular functions Outline 1 What are submodular functions Motivation Submodularity and Concavity Examples 2 Properties of submodular functions Submodularity
More informationDecomposition in Multistage Stochastic Programming and a Constraint Integer Programming Approach to Mixed-Integer Nonlinear Programming
Decomposition in Multistage Stochastic Programming and a Constraint Integer Programming Approach to Mixed-Integer Nonlinear Programming D I S S E R T A T I O N zur Erlangung des akademischen Grades Dr.
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationNotes on Dantzig-Wolfe decomposition and column generation
Notes on Dantzig-Wolfe decomposition and column generation Mette Gamst November 11, 2010 1 Introduction This note introduces an exact solution method for mathematical programming problems. The method is
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7, D Berlin
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7, D-14195 Berlin Georg Ch. Pug Andrzej Ruszczynski Rudiger Schultz On the Glivenko-Cantelli Problem in Stochastic Programming: Mixed-Integer
More informationTutorial on quasi-monte Carlo methods
Tutorial on quasi-monte Carlo methods Josef Dick School of Mathematics and Statistics, UNSW, Sydney, Australia josef.dick@unsw.edu.au Comparison: MCMC, MC, QMC Roughly speaking: Markov chain Monte Carlo
More informationScenario Grouping and Decomposition Algorithms for Chance-constrained Programs
Scenario Grouping and Decomposition Algorithms for Chance-constrained Programs Siqian Shen Dept. of Industrial and Operations Engineering University of Michigan Joint work with Yan Deng (UMich, Google)
More informationThe Value of Randomized Solutions in Mixed-Integer Distributionally Robust Optimization Problems
The Value of Randomized Solutions in Mixed-Integer Distributionally Robust Optimization Problems Erick Delage 1 and Ahmed Saif 2 1 Department of Decision Sciences, HEC Montréal 2 Department of Industrial
More informationAirline Network Revenue Management by Multistage Stochastic Programming
Airline Network Revenue Management by Multistage Stochastic Programming K. Emich, H. Heitsch, A. Möller, W. Römisch Humboldt-University Berlin, Department of Mathematics Page 1 of 18 GOR-Arbeitsgruppe
More information3.4 Relaxations and bounds
3.4 Relaxations and bounds Consider a generic Discrete Optimization problem z = min{c(x) : x X} with an optimal solution x X. In general, the algorithms generate not only a decreasing sequence of upper
More informationBilevel Integer Linear Programming
Bilevel Integer Linear Programming TED RALPHS SCOTT DENEGRE ISE Department COR@L Lab Lehigh University ted@lehigh.edu MOPTA 2009, Lehigh University, 19 August 2009 Thanks: Work supported in part by the
More informationThe moment-lp and moment-sos approaches
The moment-lp and moment-sos approaches LAAS-CNRS and Institute of Mathematics, Toulouse, France CIRM, November 2013 Semidefinite Programming Why polynomial optimization? LP- and SDP- CERTIFICATES of POSITIVITY
More informationAn Integer Programming Approach for Linear Programs with Probabilistic Constraints
An Integer Programming Approach for Linear Programs with Probabilistic Constraints James Luedtke Shabbir Ahmed George Nemhauser Georgia Institute of Technology 765 Ferst Drive, Atlanta, GA, USA luedtke@gatech.edu
More informationRobustness in Stochastic Programs with Risk Constraints
Robustness in Stochastic Programs with Risk Constraints Dept. of Probability and Mathematical Statistics, Faculty of Mathematics and Physics Charles University, Prague, Czech Republic www.karlin.mff.cuni.cz/~kopa
More informationBCOL RESEARCH REPORT 07.04
BCOL RESEARCH REPORT 07.04 Industrial Engineering & Operations Research University of California, Berkeley, CA 94720-1777 LIFTING FOR CONIC MIXED-INTEGER PROGRAMMING ALPER ATAMTÜRK AND VISHNU NARAYANAN
More informationHW1 solutions. 1. α Ef(x) β, where Ef(x) is the expected value of f(x), i.e., Ef(x) = n. i=1 p if(a i ). (The function f : R R is given.
HW1 solutions Exercise 1 (Some sets of probability distributions.) Let x be a real-valued random variable with Prob(x = a i ) = p i, i = 1,..., n, where a 1 < a 2 < < a n. Of course p R n lies in the standard
More informationOptimization Problems with Probabilistic Constraints
Optimization Problems with Probabilistic Constraints R. Henrion Weierstrass Institute Berlin 10 th International Conference on Stochastic Programming University of Arizona, Tucson Recommended Reading A.
More informationChance constrained optimization - applications, properties and numerical issues
Chance constrained optimization - applications, properties and numerical issues Dr. Abebe Geletu Ilmenau University of Technology Department of Simulation and Optimal Processes (SOP) May 31, 2012 This
More informationAn Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse
An Adaptive Partition-based Approach for Solving Two-stage Stochastic Programs with Fixed Recourse Yongjia Song, James Luedtke Virginia Commonwealth University, Richmond, VA, ysong3@vcu.edu University
More informationMIDAS: A Mixed Integer Dynamic Approximation Scheme
MIDAS: A Mixed Integer Dynamic Approximation Scheme Andy Philpott, Faisal Wahid, Frédéric Bonnans May 7, 2016 Abstract Mixed Integer Dynamic Approximation Scheme (MIDAS) is a new sampling-based algorithm
More informationStochastic Integer Programming
IE 495 Lecture 20 Stochastic Integer Programming Prof. Jeff Linderoth April 14, 2003 April 14, 2002 Stochastic Programming Lecture 20 Slide 1 Outline Stochastic Integer Programming Integer LShaped Method
More informationA Two-Stage Moment Robust Optimization Model and its Solution Using Decomposition
A Two-Stage Moment Robust Optimization Model and its Solution Using Decomposition Sanjay Mehrotra and He Zhang July 23, 2013 Abstract Moment robust optimization models formulate a stochastic problem with
More information3.8 Strong valid inequalities
3.8 Strong valid inequalities By studying the problem structure, we can derive strong valid inequalities which lead to better approximations of the ideal formulation conv(x ) and hence to tighter bounds.
More informationA Proximal Method for Identifying Active Manifolds
A Proximal Method for Identifying Active Manifolds W.L. Hare April 18, 2006 Abstract The minimization of an objective function over a constraint set can often be simplified if the active manifold of the
More informationExtended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications
Extended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications François Vanderbeck University of Bordeaux INRIA Bordeaux-Sud-Ouest part : Defining Extended Formulations
More informationScenario generation for stochastic optimization problems via the sparse grid method
Scenario generation for stochastic optimization problems via the sparse grid method Michael Chen Sanjay Mehrotra Dávid Papp March 7, 0 Abstract We study the use of sparse grids methods for the scenario
More informationarxiv: v2 [math.oc] 10 May 2017
Conic Programming Reformulations of Two-Stage Distributionally Robust Linear Programs over Wasserstein Balls arxiv:1609.07505v2 [math.oc] 10 May 2017 Grani A. Hanasusanto 1 and Daniel Kuhn 2 1 Graduate
More information1. Introduction. In this paper we consider stochastic optimization problems of the form
O RATES OF COVERGECE FOR STOCHASTIC OPTIMIZATIO PROBLEMS UDER O-I.I.D. SAMPLIG TITO HOMEM-DE-MELLO Abstract. In this paper we discuss the issue of solving stochastic optimization problems by means of sample
More information- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs
LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs
More informationOn John type ellipsoids
On John type ellipsoids B. Klartag Tel Aviv University Abstract Given an arbitrary convex symmetric body K R n, we construct a natural and non-trivial continuous map u K which associates ellipsoids to
More information3. Linear Programming and Polyhedral Combinatorics
Massachusetts Institute of Technology 18.433: Combinatorial Optimization Michel X. Goemans February 28th, 2013 3. Linear Programming and Polyhedral Combinatorics Summary of what was seen in the introductory
More informationEMPIRICAL ESTIMATES IN STOCHASTIC OPTIMIZATION VIA DISTRIBUTION TAILS
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 3, P AGES 459 471 EMPIRICAL ESTIMATES IN STOCHASTIC OPTIMIZATION VIA DISTRIBUTION TAILS Vlasta Kaňková Classical optimization problems depending on a probability
More informationDistributionally robust simple integer recourse
Distributionally robust simple integer recourse Weijun Xie 1 and Shabbir Ahmed 2 1 Department of Industrial and Systems Engineering, Virginia Tech, Blacksburg, VA 24061 2 School of Industrial & Systems
More informationInteger Programming ISE 418. Lecture 12. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 12 Dr. Ted Ralphs ISE 418 Lecture 12 1 Reading for This Lecture Nemhauser and Wolsey Sections II.2.1 Wolsey Chapter 9 ISE 418 Lecture 12 2 Generating Stronger Valid
More informationLIGHT ROBUSTNESS. Matteo Fischetti, Michele Monaci. DEI, University of Padova. 1st ARRIVAL Annual Workshop and Review Meeting, Utrecht, April 19, 2007
LIGHT ROBUSTNESS Matteo Fischetti, Michele Monaci DEI, University of Padova 1st ARRIVAL Annual Workshop and Review Meeting, Utrecht, April 19, 2007 Work supported by the Future and Emerging Technologies
More informationOn the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times
CORE DISCUSSION PAPER 2000/52 On the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times Andrew J. Miller 1, George L. Nemhauser 2, and Martin W.P. Savelsbergh 2 November 2000
More informationLebesgue Measure on R n
CHAPTER 2 Lebesgue Measure on R n Our goal is to construct a notion of the volume, or Lebesgue measure, of rather general subsets of R n that reduces to the usual volume of elementary geometrical sets
More informationScenario tree generation for multi-stage stochastic programs
Scenario tree generation for multi-stage stochastic programs Holger Heitsch 1 and Werner Römisch 2 1 Humboldt University Berlin, Institute of Mathematics, 10099 Berlin, Germany, heitsch@math.hu-berlin.de
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7, D Berlin
Konrad-use-entrum für Informationstechnik Berlin Takustraße 7, D-14195 Berlin Rudiger Schultz A Note on Preprocessing via Fourier - Motzkin Elimination in Two-Stage Stochastic Programming Preprint SC 96{32
More informationSolving Mixed-Integer Nonlinear Programs
Solving Mixed-Integer Nonlinear Programs (with SCIP) Ambros M. Gleixner Zuse Institute Berlin MATHEON Berlin Mathematical School 5th Porto Meeting on Mathematics for Industry, April 10 11, 2014, Porto
More informationQuadratic Two-Stage Stochastic Optimization with Coherent Measures of Risk
Noname manuscript No. (will be inserted by the editor) Quadratic Two-Stage Stochastic Optimization with Coherent Measures of Risk Jie Sun Li-Zhi Liao Brian Rodrigues Received: date / Accepted: date Abstract
More informationAPPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H.
APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION I.H. Sloan, X. Wang and H. Wozniakowski Published in Journal of Complexity, Volume 20, Number 1,
More informationK-Adaptability in Two-Stage Mixed-Integer Robust Optimization
K-Adaptability in Two-Stage Mixed-Integer Robust Optimization Anirudh Subramanyam 1, Chrysanthos E. Gounaris 1, and Wolfram Wiesemann 2 asubramanyam@cmu.edu, gounaris@cmu.edu, ww@imperial.ac.uk 1 Department
More informationOn deterministic reformulations of distributionally robust joint chance constrained optimization problems
On deterministic reformulations of distributionally robust joint chance constrained optimization problems Weijun Xie and Shabbir Ahmed School of Industrial & Systems Engineering Georgia Institute of Technology,
More informationMixed-Integer Nonlinear Programming
Mixed-Integer Nonlinear Programming Claudia D Ambrosio CNRS researcher LIX, École Polytechnique, France pictures taken from slides by Leo Liberti MPRO PMA 2016-2017 Motivating Applications Nonlinear Knapsack
More informationValid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints
Valid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints Nilay Noyan Andrzej Ruszczyński March 21, 2006 Abstract Stochastic dominance relations
More informationConvex Optimization & Parsimony of L p-balls representation
Convex Optimization & Parsimony of L p -balls representation LAAS-CNRS and Institute of Mathematics, Toulouse, France IMA, January 2016 Motivation Unit balls associated with nonnegative homogeneous polynomials
More informationMixed Integer Linear Programming Formulations for Probabilistic Constraints
Mixed Integer Linear Programming Formulations for Probabilistic Constraints J. P. Vielma a,, S. Ahmed b, G. Nemhauser b a Department of Industrial Engineering, University of Pittsburgh 1048 Benedum Hall,
More informationOn the Value Function of a Mixed Integer Linear Optimization Problem and an Algorithm for its Construction
On the Value Function of a Mixed Integer Linear Optimization Problem and an Algorithm for its Construction Ted K. Ralphs and Anahita Hassanzadeh Department of Industrial and Systems Engineering, Lehigh
More informationOn Approximations and Stability in Stochastic Programming
On Approximations and Stability in Stochastic Programming Peter Kall Abstract It has become an accepted approach to attack stochastic programming problems by approximating the given probability distribution
More informationWeighted Ehrhart polynomials and intermediate sums on polyhedra
Weighted Ehrhart polynomials and intermediate sums on polyhedra Nicole BERLINE Centre de mathématiques Laurent Schwartz, Ecole polytechnique, France Inverse Moment Problems NUS-IMS, Dec 17, 2013 with Velleda
More informationScenario Generation and Sampling Methods
Scenario Generation and Sampling Methods Tito Homem-de-Mello Güzin Bayraksan SVAN 2016 IMPA 9 13 May 2106 Let s Start with a Recap Sample Average Approximation We want to solve the true problem min {g(x)
More informationGeometric problems. Chapter Projection on a set. The distance of a point x 0 R n to a closed set C R n, in the norm, is defined as
Chapter 8 Geometric problems 8.1 Projection on a set The distance of a point x 0 R n to a closed set C R n, in the norm, is defined as dist(x 0,C) = inf{ x 0 x x C}. The infimum here is always achieved.
More informationOptimization Exercise Set n. 4 :
Optimization Exercise Set n. 4 : Prepared by S. Coniglio and E. Amaldi translated by O. Jabali 2018/2019 1 4.1 Airport location In air transportation, usually there is not a direct connection between every
More informationOn the cells in a stationary Poisson hyperplane mosaic
On the cells in a stationary Poisson hyperplane mosaic Matthias Reitzner and Rolf Schneider Abstract Let X be the mosaic generated by a stationary Poisson hyperplane process X in R d. Under some mild conditions
More informationComputational Integer Programming. Lecture 2: Modeling and Formulation. Dr. Ted Ralphs
Computational Integer Programming Lecture 2: Modeling and Formulation Dr. Ted Ralphs Computational MILP Lecture 2 1 Reading for This Lecture N&W Sections I.1.1-I.1.6 Wolsey Chapter 1 CCZ Chapter 2 Computational
More informationPolyhedral Approaches to Online Bipartite Matching
Polyhedral Approaches to Online Bipartite Matching Alejandro Toriello joint with Alfredo Torrico, Shabbir Ahmed Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Industrial
More informationStrong Formulations of Robust Mixed 0 1 Programming
Math. Program., Ser. B 108, 235 250 (2006) Digital Object Identifier (DOI) 10.1007/s10107-006-0709-5 Alper Atamtürk Strong Formulations of Robust Mixed 0 1 Programming Received: January 27, 2004 / Accepted:
More informationOptimization methods NOPT048
Optimization methods NOPT048 Jirka Fink https://ktiml.mff.cuni.cz/ fink/ Department of Theoretical Computer Science and Mathematical Logic Faculty of Mathematics and Physics Charles University in Prague
More information