Automatic Scheduling of Hypermedia Documents with Elastic Times

Similar documents
MVE165/MMG630, Applied Optimization Lecture 8 Integer linear programming algorithms. Ann-Brith Strömberg

MVE165/MMG631 Linear and integer optimization with applications Lecture 9 Discrete optimization: theory and algorithms

15.082J and 6.855J. Lagrangian Relaxation 2 Algorithms Application to LPs

Introduction to Mathematical Programming IE406. Lecture 20. Dr. Ted Ralphs

Fundamentals of Integer Programming

Assignment 3b: The traveling salesman problem

1. Lecture notes on bipartite matching February 4th,

Lagrangean relaxation - exercises

Final Exam Spring 2003

Applied Lagrange Duality for Constrained Optimization

An SDP Approach to Multi-level Crossing Minimization

Column Generation Method for an Agent Scheduling Problem

Algorithms for Integer Programming

Linear Programming. Course review MS-E2140. v. 1.1

Minimum Weight Constrained Forest Problems. Problem Definition

Math Models of OR: The Simplex Algorithm: Practical Considerations

Effective probabilistic stopping rules for randomized metaheuristics: GRASP implementations

A Computational Study of Conflict Graphs and Aggressive Cut Separation in Integer Programming

Modeling, authoring and formatting hypermedia documents in the HyperProp system

Advanced Operations Research Techniques IE316. Quiz 2 Review. Dr. Ted Ralphs

56:272 Integer Programming & Network Flows Final Examination -- December 14, 1998

Local search with perturbations for the prize collecting Steiner tree problem in graphs

Discrete Optimization with Decision Diagrams

Column Generation Based Primal Heuristics

Mathematical and Algorithmic Foundations Linear Programming and Matchings

Link Dimensioning and LSP Optimization for MPLS Networks Supporting DiffServ EF and BE Classes

Introduction to Mathematical Programming IE496. Final Review. Dr. Ted Ralphs

Lagrangian Relaxation: An overview

A Branch-and-Cut Algorithm for the Partition Coloring Problem

ACO Comprehensive Exam October 12 and 13, Computability, Complexity and Algorithms

Integer and Combinatorial Optimization

Department of Mathematics Oleg Burdakov of 30 October Consider the following linear programming problem (LP):

to the Traveling Salesman Problem 1 Susanne Timsj Applied Optimization and Modeling Group (TOM) Department of Mathematics and Physics

Methods and Models for Combinatorial Optimization Exact methods for the Traveling Salesman Problem

56:272 Integer Programming & Network Flows Final Exam -- December 16, 1997

15.083J Integer Programming and Combinatorial Optimization Fall Enumerative Methods

Integer Programming ISE 418. Lecture 7. Dr. Ted Ralphs

Integer Programming Theory

3 INTEGER LINEAR PROGRAMMING

The MIP-Solving-Framework SCIP

/ Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang

A subgradient-based branch-and-bound algorithm for the capacitated facility location problem

Heuristics in MILP. Group 1 D. Assouline, N. Molyneaux, B. Morén. Supervisors: Michel Bierlaire, Andrea Lodi. Zinal 2017 Winter School

Combinatorial Auctions: A Survey by de Vries and Vohra

9.4 SOME CHARACTERISTICS OF INTEGER PROGRAMS A SAMPLE PROBLEM

Selected Topics in Column Generation

Combinatorial Optimization

Approximation Algorithms: The Primal-Dual Method. My T. Thai

The Heuristic (Dark) Side of MIP Solvers. Asja Derviskadic, EPFL Vit Prochazka, NHH Christoph Schaefer, EPFL

Algorithms for Decision Support. Integer linear programming models

College of Computer & Information Science Fall 2007 Northeastern University 14 September 2007

Maximum flow problem CE 377K. March 3, 2015

Graph Coloring via Constraint Programming-based Column Generation

An Introduction to Dual Ascent Heuristics

and 6.855J Lagrangian Relaxation I never missed the opportunity to remove obstacles in the way of unity. Mohandas Gandhi

DETERMINISTIC OPERATIONS RESEARCH

Lecture 7: Asymmetric K-Center

lpsymphony - Integer Linear Programming in R

Recursive column generation for the Tactical Berth Allocation Problem

Primal Heuristics in SCIP

Lecture 14: Linear Programming II

Crew Scheduling Problem: A Column Generation Approach Improved by a Genetic Algorithm. Santos and Mateus (2007)

Tuesday, April 10. The Network Simplex Method for Solving the Minimum Cost Flow Problem

Improved Gomory Cuts for Primal Cutting Plane Algorithms

Branch-price-and-cut for vehicle routing. Guy Desaulniers

Cutting Planes for Some Nonconvex Combinatorial Optimization Problems

Solving a Challenging Quadratic 3D Assignment Problem

3 No-Wait Job Shops with Variable Processing Times

Network Flows. 7. Multicommodity Flows Problems. Fall 2010 Instructor: Dr. Masoud Yaghini

Integer Programming ISE 418. Lecture 1. Dr. Ted Ralphs

1 Linear programming relaxation

LOCAL SEARCH FOR THE MINIMUM FUNDAMENTAL CYCLE BASIS PROBLEM

February 19, Integer programming. Outline. Problem formulation. Branch-andbound

Investigating Mixed-Integer Hulls using a MIP-Solver

Branch-Cut-and-Price solver for Vehicle Routing Problems

Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras

Lagrangean Methods bounding through penalty adjustment

A Generic Separation Algorithm and Its Application to the Vehicle Routing Problem

Minimum-Energy Broadcast and Multicast in Wireless Networks: An Integer Programming Approach and Improved Heuristic Algorithms

Mathematical Tools for Engineering and Management

Branch-and-bound: an example

1. Lecture notes on bipartite matching

Exact solutions to mixed-integer linear programming problems

Characterizing Improving Directions Unconstrained Optimization

Using Decomposition Techniques for Solving Large-Scale Capacitated Hub Location Problems with Single Assignment

Last topic: Summary; Heuristics and Approximation Algorithms Topics we studied so far:

arxiv: v1 [cs.dm] 6 May 2009

ONLY AVAILABLE IN ELECTRONIC FORM

Primal and Dual Methods for Optimisation over the Non-dominated Set of a Multi-objective Programme and Computing the Nadir Point

Discrete Covering. Location. Problems. Louis. Luangkesorn. Housekeeping. Dijkstra s Shortest Path. Discrete. Covering. Models.

Two-stage column generation

x ji = s i, i N, (1.1)

A Simpler Generalized Shortest Paths Algorithm

Improving Dual Bound for Stochastic MILP Models Using Sensitivity Analysis

Lecture 4: Primal Dual Matching Algorithm and Non-Bipartite Matching. 1 Primal/Dual Algorithm for weighted matchings in Bipartite Graphs

Resource Management in Computer Networks -- Mapping from engineering problems to mathematical formulations

In this lecture, we ll look at applications of duality to three problems:

CS599: Convex and Combinatorial Optimization Fall 2013 Lecture 14: Combinatorial Problems as Linear Programs I. Instructor: Shaddin Dughmi

A Comparison of Mixed-Integer Programming Models for Non-Convex Piecewise Linear Cost Minimization Problems

Module 10. Network Simplex Method:

Transcription:

Automatic Scheduling of Hypermedia Documents with Elastic Times Celso C. Ribeiro Joint work with M.T. Medina and L.F.G. Soares Computer Science Department Catholic University of Rio de Janeiro, Brazil Workshop on New Trends in Scheduling Parallel and Distributed Systems Luminy, October 1-5, 2001

Automatic Scheduling 2 Contents Contents: Introduction Problem formulation Primal heuristic Lagrangean relaxation Dual heuristic Concluding remarks

Automatic Scheduling 3 Introduction Hypermedia systems should have some desirable features, such as: supporting a rich set of capabilities for each media segment (video, audio, text, pictures, animation,... ): alternatives for different presentations presentation duration flexibility behavior changes during presentation allowing the explicit definition of temporal relationships among media segments with non-deterministic times supporting a temporal formatting algorithm that considers non-determinism allows on-the-fly corrections whenever unpredictable events (network delays, user interactions) occur

Automatic Scheduling 4 Introduction A presentation event is a media segment presentation d min : minimum allowed presentation time d max : maximum allowed presentation time d ideal :durationwithbestpresentationquality c: cost of shrinking or stretching the presentation duration

Automatic Scheduling 5 Introduction Scheduling algorithm: determination of starting times and optimal durations d opt (which will be tried by the formatter) so as to ensure spatial and temporal consistency at compile time: build an execution plan to guide the scheduling of the object presentation at run time: online adjustments to handle unpredictable events (network delays, user interactions) Difficult of performing fast rescheduling in run time led the systems that make time adjustments to make them only at compile time. Other systems: Firefly, Isis, Madeus (linear programming, PERT, constraints,... )

Automatic Scheduling 6 Introduction Problem: find (at compile time) the otimal starting times and durations of objects to be presented, to ensure spatial and temporal consistency of the presentation while respecting limits on shrinking and stretching the ideal duration of each object. Objectives: (1) minimize the total cost of shrinking and stretching objects; (2) minimize the number of presentation events that will have their durations different than the ideal ones; (3) spread the changes in object duration proportionally to their duration; and (4) find the shortest and longest presentation durations. Objectives (1), (3), and (4) can be dealt with at compile time by linear programming models.

Automatic Scheduling 7 Introduction Although conflicting, objectives (1) minimize the total cost of shrinking and stretching objects and (2) minimize the number of presentation events that will have their durations different than the ideal ones are the most natural. Possible strategy to handle both of them:select the best solution with respect to (2) among all those which optimize (1). The combinatorial nature of the minimization of the number of objects whose ideal duration is modified makes objective (2) the most difficult to be tackled by combinatorial algorithms. Extendable mixed 0-1 LIP formulation using objective (2) Fast and effective primal heuristic based on variable fixations Lagrangean relaxation and dual heuristic

Automatic Scheduling 8 Problem formulation Hypermedia document represented by a directed graph G =[N,A]: set Ns = {1,...,n} of synchronization points where presentation events either start or stop set N = {0} {1,...,n} {n +1} of nodes, where 0 and n +1 denote respectively the initial and the final nodes set A of arcs: presentation events are arcs (i, j) with both extremities i, j {1,...,n} and durations d min ij, d ideal ij,andd max ij dummy initial arcs (0,j)withj {1,...,n} and durations d min 0j =0,d ideal 0j =0,andd max 0j = dummy final arcs (i, n +1)withi {1,...,n} and durations d min i,n+1 =0,dideal i,n+1 =0,anddmax i,n+1 =

Automatic Scheduling 9 Problem formulation Variables: continuous variable ti 0, for every i N time in which synchronization takes place at i Ns presentation starting time: t0 = 0 end of presentation: tn+1 continuous variable xij [d min ij,d max ij ], for every (i, j) A duration of each presentation (or dummy) event (i, j) A binary variable yij {0, 1}, for every event (i, j) A yij = 1, if the ideal duration of event (i, j) ismodified yij =0,otherwise artificial costs cij associated to every arc (i, j) A cij =1,if(i, j) is a presentation event cij = 0, otherwise (i.e., (i, j) isadummyevent)

Automatic Scheduling 10 Problem formulation Formulation: mixed 0-1 integer programming problem ASHD z =min with: (i,j) A cij yij ti tj + xij =0 (i, j) A (1) xij bij yij d ideal ij (i, j) A (2) xij +aij yij d ideal ij (i, j) A (3) yij {0, 1} (i, j) A (4) ti 0 i N \{0} t0 =0 = d min d ideal aij ij ij = d max d ideal bij ij ij

Automatic Scheduling 11 Problem formulation Formulation: mixed 0-1 integer programming problem ASHD z =min with: (i,j) A cij yij ti tj + xij =0 (i, j) A (1) xij bij yij d ideal ij (i, j) A (2) xij +aij yij d ideal ij (i, j) A (3) yij {0, 1} (i, j) A (4) ti 0 i N \{0} t0 =0 = d min d ideal aij ij ij = d max d ideal bij ij ij Open question: complexity of ASHD?

Automatic Scheduling 12 Primal heuristic ASHD is the LP relaxation of ASHD: yij {0, 1} 0 yij 1 Primal heuristic PH for problem ASHD: procedure PH 1 Solve ASHD and let (y, x, t) be its optimal solution; 2 while (i, j) A :0< y ij < 1do 3 Fix to one in ASHD all variables yij such that y ij =1; 4 Fix to zero in ASHD all variables yij such that y ij =0; 5 Variable selection: select and fix at either zero or one in ASHD a variable yst with a fractional value 0 < y st < 1inASHD; 6 Solve ASHD and let (y, x, t) be its optimal solution; 7 end while; 8 return y; end PH;

Automatic Scheduling 13 Primal heuristic procedure Variable selection 1 if (i, j) A, i = 0,j = n +1:0< y ij α then do 2 (s, t) argmin (i,j) A,i =0,j =n+1 {y ij :0< y ij α}; 3 Fixyst to zero in ASHD; 4 if ASHD is not feasible then fix yst to one in ASHD; 5 else if (i, j) A, i = 0,j = n +1:1 α y ij < 1thendo 6 (s, t) argmax (i,j) A,i =0,j =n+1 {y ij :1 α y ij < 1}; 7 Fix yst to one in ASHD; 8 else do 9 (s, t) argmin (i,j) A,i =0,j =n+1 min{y ij, 1 y ij }; 10 if yst 0.5 thendo 11 Fix yst to zero in ASHD; 12 if ASHD is not feasible then fix yst to one in ASHD; 13 else fix yst to one in ASHD; 14 return (s, t); end Variable selection;

Automatic Scheduling 14 Computational results PH Primal heuristic UltraSPARC station, CPLEX 5.0, 35 instances, n = 20 1000 Best value Time (s) Duality Instance Nodes Arcs PH Exact PH Exact gap (%) a1 20 40 22 20 0.09 2.98 59.6 a2 20 80 62 62 0.29 25.07 47.8 a3 20 120 103 102 0.76 2328.75 56.4 c1 40 80 45 42 0.28 11535.65 57.3 c2 40 160 125 124 1.05 > 4hs 52.3 d1 50 100 53 53 0.33 4246.48 51.5 d2 50 200 156 154 1.54 9190.92 52.3 e1 100 200 103 102 1.29 10087.15 53.1 e2 100 400 303 304 6.74 9715.30 54.5 maximum size of the enumeration tree exceeded

Automatic Scheduling 15 Computational results PH Conclusions: PH found good solutions within small computation times for real-size instances. Reasonable computation times, except for very large instances (500 or more nodes). PH found the optimal solutions for two of seven instances for which the exact optimal is known and was off by exactly one unit for two others. Average and maximum gaps are respectively 54.5% and 59.6%. Very large duality gaps even for small instances for which PH found the optimal solution:

Automatic Scheduling 16 Computational results PH Conclusions: PH found good solutions within small computation times for real-size instances. Reasonable computation times, except for very large instances (500 or more nodes). PH found the optimal solutions for two of seven instances for which the exact optimal is known and was off by exactly one unit for two others. Average and maximum gaps are 54.5% and 59.6%. Very large duality gaps even for small instances for which PH found the optimal solution: conjecture: solutions are quite good, but bounds are poor search for polyhedral cuts or valid inequalities should be pursued if one wants to develop efficient exact methods

Automatic Scheduling 17 Lagrangean relaxation Lagrangean dual of problem ASHD: Associate non-negative dual multipliers λ + ij and λ ij to constraints (2) and (3) for each arc (i, j) A: Lagrangean function: L(x, y, t, λ +, λ )= cij yij+λ + ij (x ij bij yij d ideal ij )+λ ij ( x ij+aij yij+d ideal ij ) (i,j) A Dual function: w(λ +, λ )=min L(x, y, t, λ +, λ ) with: ti tj + xij =0 (i, j) A (1) yi,j {0, 1} (i, j) A (4) d min ij xij d max ij (i, j) A (5) ti 0 i N \{0} t0 =0.

Automatic Scheduling 18 Dual heuristic Dual heuristic DH: Combines the primal heuristic and the solution of the dual problem by subgradient optimization. Uses the multipliers obtained along the solution of the dual problem to periodically generate reduced costs c ij associated with the 0-1 variables yij. Periodically applies the primal heuristic PH. Reduced costs c ij replacing the original costs c ij in the primal heuristic are computed according with two different schemes: c ij = c ij + aij λ ij b ij λ + ij (Lagrangean multiplipliers) c ij =(1 ŷ ij) cij (primal aproximation) Stopping criteria: total number of iterations of the subgradient algorithm equals 1200 or the relative difference between two consecutive dual function values becomes less than 10 4.

Automatic Scheduling 19 Dual heuristic Dual heuristic DH for problem ASHD: procedure DH 1 Apply PH to compute a feasible solution y; 2 Set the upper bound zbest (i,j) A c ij y ij and ybest y; 3 while stopping criteria not attained do 4 Perform a certain number of subgradient iterations; 5 Recompute reduced costs c using the current dual multipliers (λ +, λ ); 6 Apply PH using c to compute a new feasible solution y; 7 if (i,j) A c ij y ij <zbest 8 then set zbest (i,j) A c ij y ij and ybest y; 9 end; 10 return y; end DH;

Automatic Scheduling 20 Computational results DH Dual heuristic (1/2) DH PH Exact Instance Nodes Arcs zbest s zbest s z a4 20 160 142 18.23 144 1.05 141 a5 20 200 178 24.88 180 1.72 176 b1 30 60 31 8.35 32 0.16 31 b3 30 180 154 21.87 155 1.46 151 b4 30 240 209 31.88 210 2.64 208 c3 40 240 202 32.83 203 2.64 205 c4 40 320 280 49.20 283 4.71 282 c5 40 400 362 67.19 363 7.54 367 d2 50 200 154 28.63 156 1.54 154 d4 50 400 352 69.73 354 7.07 358 d5 50 500 451 95.40 453 11.62 452

Automatic Scheduling 21 Computational results DH Dual heuristic (2/2) DH PH Exact Instance Nodes Arcs zbest s zbest s z e3 100 600 505 142.63 506 16.69 520 e4 100 800 699 219.34 702 30.47 709 f1 500 1000 510 395.25 515 42.01 520 f3 500 3000 2539 3602.95 2541 782.47 2551 f5 500 5000 4514 5678.42 4519 1659.88 4553 g1 1000 2000 1043 786.30 1045 208.05 1177

Automatic Scheduling 22 Computational results DH Conclusions: The dual heuristic DH consistently improves upon the solutions found by the primal heuristic PH, finding better solutions for 17 out of the 35 test problems. However, the computational times observed for the dual heuristic DH are much larger. Need for better/faster construction algorithms and local search procedures (work-in-progress) based on variable complementation.

Automatic Scheduling 23 Concluding remarks Concluding remarks and extensions: Primal heuristic: effective in finding good solutions in small processing times (real-size instances); quite promising not only in algorithmic terms, but also concerning its integration within existing document formatting systems.

Automatic Scheduling 24 Concluding remarks Concluding remarks and extensions: Primal heuristic: effective in finding good solutions in small processing times (real-size instances); quite promising not only in algorithmic terms, but also concerning its integration within existing document formatting systems. Dual heuristic: better solutions, but at the cost of much larger computation times.

Automatic Scheduling 25 Concluding remarks Concluding remarks and extensions: Primal heuristic: effective in finding good solutions in small processing times (real-size instances); quite promising not only in algorithmic terms, but also concerning its integration within existing document formatting systems. Dual heuristic: better solutions, but at the cost of much larger computation times. Establish the computational complexity of problem ASHD.

Automatic Scheduling 26 Concluding remarks Concluding remarks and extensions: Primal heuristic: effective in finding good solutions in small processing times (real-size instances); quite promising not only in algorithmic terms, but also concerning its integration within existing document formatting systems. Dual heuristic: better solutions, but at the cost of much larger computation times. Establish the computational complexity of problem ASHD. Replace CPLEX by a customized LP algorithm to speedup the solution of the linear relaxation ASHD.

Automatic Scheduling 27 Concluding remarks Concluding remarks and extensions: Primal heuristic: effective in finding good solutions in small processing times (real-size instances); quite promising not only in algorithmic terms, but also concerning its integration within existing document formatting systems. Dual heuristic: better solutions, but at the cost of much larger computation times. Establish the computational complexity of problem ASHD. Replace CPLEX by a customized LP algorithm to speedup the solution of the linear relaxation ASHD. Improved construction procedures and local search strategies.

Automatic Scheduling 28 Concluding remarks Concluding remarks and extensions: Primal heuristic: effective in finding good solutions in small processing times (real-size instances); quite promising not only in algorithmic terms, but also concerning its integration within existing document formatting systems. Dual heuristic: better solutions, but at the cost of much larger computation times. Establish the computational complexity of problem ASHD. Replace CPLEX by a customized LP algorithm to speedup the solution of the linear relaxation ASHD. Improved construction procedures and local search strategies. Integration of such heuristics and the bicriteria optimization approach within an existing document formatter to perform compile time scheduling and even run time adjustments.

Automatic Scheduling 29 Paper and slides Automatic Scheduling of Hypermedia Documents with Elastic Times Maira T. Medina Celso C. Ribeiro Luiz Fernando G. Soares Computer Science Department Catholic University of Rio de Janeiro Paper: //www-di.inf.puc-rio.br/ celso/artigos/ashd.pdf Slides: //www-di.inf.puc-rio.br/ celso/talks/luminy.pdf