04.02.2015 Views

Stochastic Programming - Index of

Stochastic Programming - Index of

Stochastic Programming - Index of

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

288 STOCHASTIC PROGRAMMING<br />

a solution, the column is not part <strong>of</strong> the frame, and can be removed. An<br />

important property <strong>of</strong> this procedure is that to determine if a column can be<br />

discarded, we have to use all other columns in the test. This is a major reason<br />

why procedure framebylp is so slow when the number <strong>of</strong> columns gets very<br />

large.<br />

So, a generator w ∗ <strong>of</strong> the cone pos W ∗ has the property that a right-hand<br />

side h must satisfy h T w ∗ ≤ 0 to be feasible. In the uncapacitated network<br />

case we saw that a right-hand side β had to satisfy b(Y ) T β ≤ 0torepresenta<br />

feasible problem. Therefore the index vector b(Y ) corresponds exactly to the<br />

column w ∗ . And calling procedure framebylp to remove those columns that<br />

are not in the frame <strong>of</strong> the cone pos W ∗ corresponds to using Proposition 6.5.<br />

Therefore the index vector <strong>of</strong> a node set from Proposition 6.5 corresponds to<br />

the columns in W ∗ .<br />

Computationally there are major differences, though. First, to find a<br />

candidate for W ∗ , we had to start out with W , and use procedure support,<br />

which is an iterative procedure. The network inequalities, on the other hand,<br />

are produced more directly by looking at all subsets <strong>of</strong> nodes. But the<br />

most important difference is that, while the use <strong>of</strong> procedure framebylp,<br />

as just explained, requires all columns to be available in order to determine if<br />

one should be discarded, Proposition 6.5 is totally local. We can pick up<br />

an inequality and determine if it is needed without looking at any other<br />

inequalities. With possibly millions <strong>of</strong> candidates, this difference is crucial.<br />

We did not develop the LP case for explicit bounds on variables. If such<br />

bounds exist, they can, however, be put in as explicit constraints. If so, a<br />

column w ∗ from W ∗ corresponds to the index vector<br />

(<br />

b(Y )<br />

−a(Q + )<br />

)<br />

.<br />

6.3 Generating Relatively Complete Recourse<br />

Let us now discuss how the results obtained in the previous section can help<br />

us, and how they can be used in a setting that deserves the term preprocessing.<br />

Let us first repeat some <strong>of</strong> our terminology, in order to see how this fits in<br />

with our discussions in the LP setting.<br />

A two-stage stochastic linear programming problem where the second-stage<br />

problem is a directed capacitated network flow problem can be formulated as<br />

follows:<br />

[<br />

min x c T x + Q(x) ]<br />

s.t. Ax = b, x ≥ 0,<br />

where<br />

Q(x) = ∑ Q(x, ξ j )p j

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!