# SAS/OR 14.2 User s Guide: Mathematical Programming. The Quadratic Programming Solver

Save this PDF as:

Size: px
Start display at page:

## Transcription

1 SAS/OR 14.2 User s Guide: Mathematical Programming The Quadratic Programming Solver

3 Chapter 11 The Quadratic Programming Solver Contents Overview: QP Solver Getting Started: QP Solver Syntax: QP Solver Functional Summary QP Solver Options Details: QP Solver Interior Point Algorithm: Overview Parallel Processing Iteration Log Problem Statistics Irreducible Infeasible Set Macro Variable _OROPTMODEL_ Examples: QP Solver Example 11.1: Linear Least Squares Problem Example 11.2: Portfolio Optimization Example 11.3: Portfolio Selection with Transactions References Overview: QP Solver The OPTMODEL procedure provides a framework for specifying and solving quadratic programs. Mathematically, a quadratic programming (QP) problem can be stated as follows: min subject to 1 2 xt Qx C c T x Ax f; D; g b l x u where

4 546 Chapter 11: The Quadratic Programming Solver Q 2 R nn is the quadratic (also known as Hessian) matrix A 2 R mn is the constraints matrix x 2 R n is the vector of decision variables c 2 R n is the vector of linear objective function coefficients b 2 R m is the vector of constraints right-hand sides (RHS) l 2 R n is the vector of lower bounds on the decision variables u 2 R n is the vector of upper bounds on the decision variables The quadratic matrix Q is assumed to be symmetric; that is, q ij D q j i ; 8i; j D 1; : : : ; n Indeed, it is easy to show that even if Q 6D Q T, then the simple modification QQ D 1 2.Q C QT / produces an equivalent formulation x T Qx x T QQxI hence symmetry is assumed. When you specify a quadratic matrix, it suffices to list only lower triangular coefficients. In addition to being symmetric, Q is also required to be positive semidefinite for minimization type of models: x T Qx 0; 8x 2 R n Q is required to be negative semidefinite for maximization type of models. Convexity can come as a result of a matrix-matrix multiplication Q D LL T or as a consequence of physical laws, and so on. See Figure 11.1 for examples of convex, concave, and nonconvex objective functions.

5 Getting Started: QP Solver 547 Figure 11.1 Examples of Convex, Concave, and Nonconvex Objective Functions The order of constraints is insignificant. Some or all components of l or u (lower and upper bounds, respectively) can be omitted. Getting Started: QP Solver Consider a small illustrative example. Suppose you want to minimize a two-variable quadratic function f.x 1 ; x 2 / on the nonnegative quadrant, subject to two constraints: min 2x 1 C 3x 2 C x1 2 C 10x2 2 C subject to x 1 x 2 1 2:5x 1 x 2 x 1 C 2x x 1 0 x 2 0 To use the OPTMODEL procedure, it is not necessary to fit this problem into the general QP formulation mentioned in the section Overview: QP Solver on page 545 and to compute the corresponding parameters. However, since these parameters are closely related to the data set that is used by the OPTQP procedure and has a quadratic programming system (QPS) format, you can compute these parameters as follows. The linear objective function coefficients, vector of right-hand sides, and lower and upper bounds are identified immediately as c D ; b D ; l D 0 C1 ; u D C1

6 548 Chapter 11: The Quadratic Programming Solver Carefully construct the quadratic matrix Q. Observe that you can use symmetry to separate the main-diagonal and off-diagonal elements: 1 2 xt Qx 1 2 nx i;j D1 x i q ij x j D 1 2 nx id1 q ii x 2 i C X i>j x i q ij x j The first expression 1 2 nx id1 q ii x 2 i sums the main-diagonal elements. Thus, in this case you have q 11 D 2; q 22 D 20 Notice that the main-diagonal values are doubled in order to accommodate the 1/2 factor. Now the second term X x i q ij x j i>j sums the off-diagonal elements in the strict lower triangular part of the matrix. The only off-diagonal (x i x j ; i 6D j ) term in the objective function is 2:5 x 1 x 2, so you have q 21 D 2:5 Notice that you do not need to specify the upper triangular part of the quadratic matrix. Finally, the matrix of constraints is as follows: 1 1 A D 1 2 The following OPTMODEL program formulates the preceding problem in a manner that is very close to the mathematical specification of the given problem: /* getting started */ proc optmodel; var x1 >= 0; /* declare nonnegative variable x1 */ var x2 >= 0; /* declare nonnegative variable x2 */ /* objective: quadratic function f(x1, x2) */ minimize f = /* the linear objective function coefficients */ 2 * x1 + 3 * x2 + /* quadratic <x, Qx> */ x1 * x * x1 * x * x2 * x2; /* subject to the following constraints */ con r1: x1 - x2 <= 1; con r2: x1 + 2 * x2 >= 100;

7 Getting Started: QP Solver 549 /* specify iterative interior point algorithm (QP) * in the SOLVE statement */ solve with qp; /* print the optimal solution */ print x1 x2; save qps qpsdata; quit; The with qp clause in the SOLVE statement invokes the QP solver to solve the problem. The output is shown in Figure Figure 11.2 Summaries and Optimal Solution The OPTMODEL Procedure Objective Sense Objective Function Objective Type Problem Summary Minimization f Quadratic Number of Variables 2 Bounded Above 0 Bounded Below 2 Bounded Below and Above 0 Free 0 Fixed 0 Number of Constraints 2 Linear LE (<=) 1 Linear EQ (=) 0 Linear GE (>=) 1 Linear Range 0 Constraint Coefficients 4 Performance Information Execution Mode Single-Machine Number of Threads 4

8 550 Chapter 11: The Quadratic Programming Solver Figure 11.2 continued Solution Summary Solver QP Algorithm Interior Point Objective Function f Solution Status Optimal Objective Value Primal Infeasibility 0 Dual Infeasibility E-14 Bound Infeasibility 0 Duality Gap E-16 Complementarity 0 Iterations 6 Presolve Time 0.00 Solution Time 0.02 x1 x In this example, the SAVE QPS statement is used to save the QP problem in the QPS-format data set qpsdata, shown in Figure The data set is consistent with the parameters of general quadratic programming previously computed. Also, the data set can be used as input to the OPTQP procedure. Obs FIELD1 Figure 11.3 QPS-Format Data Set FIELD2 FIELD3 FIELD4 FIELD5 FIELD6 1 NAME qpsdata.. 2 ROWS.. 3 N f.. 4 L r1.. 5 G r2.. 6 COLUMNS.. 7 x1 f 2.0 r1 1 8 x1 r x2 f 3.0 r x2 r RHS.. 12.RHS. r RHS. r QSECTION.. 15 x1 x x1 x x2 x ENDATA..

9 QP Solver Options 551 Syntax: QP Solver The following statement is available in the OPTMODEL procedure: SOLVE WITH QP < / options > ; Functional Summary Table 11.1 summarizes the list of options available for the SOLVE WITH QP statement, classified by function. Table 11.1 Options for the QP Solver Description Solver Options Enables or disables IIS detection Control Options Specifies the frequency of printing solution progress Specifies the maximum number of iterations Specifies the time limit for the optimization process Specifies the type of presolve Interior Point Algorithm Options Specifies the stopping criterion based on duality gap Specifies the stopping criterion based on dual infeasibility Specifies the stopping criterion based on primal infeasibility Specifies units of CPU time or real time Option IIS= LOGFREQ= MAXITER= MAXTIME= PRESOLVER= STOP_DG= STOP_DI= STOP_PI= TIMETYPE= QP Solver Options This section describes the options recognized by the QP solver. These options can be specified after a forward slash (/) in the SOLVE statement, provided that the QP solver is explicitly specified using a WITH clause. The QP solver does not provide an intermediate solution if the solver terminates before reaching optimality. Solver Options IIS=number j string specifies whether the QP solver attempts to identify a set of constraints and variables that form an irreducible infeasible set (IIS). Table 11.2 describes the valid values of the IIS= option.

10 552 Chapter 11: The Quadratic Programming Solver Table 11.2 Values for IIS= Option number string Description 0 OFF Disables IIS detection. 1 ON Enables IIS detection. If an IIS is found, you can find information about the infeasibilities in the.status values of the constraints and variables. The default value of this option is OFF. See the section Irreducible Infeasible Set on page 557 for details about the IIS= option. See the section Suffixes on page 131 for details about the.status suffix. Control Options LOGFREQ=k PRINTFREQ=k specifies that the printing of the solution progress to the iteration log is to occur after every k iterations. The print frequency, k, is an integer between zero and the largest four-byte signed integer, which is The value k = 0 disables the printing of the progress of the solution. The default value of this option is 1. MAXITER=k specifies the maximum number of iterations. The value k can be any integer between one and the largest four-byte signed integer, which is If you do not specify this option, the procedure does not stop based on the number of iterations performed. MAXTIME=t specifies an upper limit of t units of time for the optimization process, including problem generation time and solution time. The value of the TIMETYPE= option determines the type of units used. If you do not specify the MAXTIME= option, the solver does not stop based on the amount of time elapsed. The value of t can be any positive number; the default value is the positive number that has the largest absolute value that can be represented in your operating environment. PRESOLVER=number string PRESOL=number string specifies one of the following presolve options: number string Description 0 NONE Disables presolver. 1 AUTOMATIC Applies presolver by using default setting. You can specify the PRESOLVER= value either by a character-valued option or by an integer. The default option is AUTOMATIC.

11 QP Solver Options 553 Interior Point Algorithm Options STOP_DG=ı specifies the desired relative duality gap, ı 2 [1E 9, 1E 4]. This is the relative difference between the primal and dual objective function values and is the primary solution quality parameter. The default value is 1E 6. See the section Interior Point Algorithm: Overview on page 554 for details. STOP_DI=ˇ specifies the maximum allowed relative dual constraints violation, ˇ 2 [1E 9, 1E 4]. The default value is 1E 6. See the section Interior Point Algorithm: Overview on page 554 for details. STOP_PI= specifies the maximum allowed relative bound and primal constraints violation, 2 [1E 9, 1E 4]. The default value is 1E 6. See the section Interior Point Algorithm: Overview on page 554 for details. TIMETYPE=number j string specifies the units of time used by the MAXTIME= option and reported by the PRESOLVE_TIME and SOLUTION_TIME terms in the _OROPTMODEL_ macro variable. Table 11.4 describes the valid values of the TIMETYPE= option. Table 11.4 Values for TIMETYPE= Option number string Description 0 CPU Specifies units of CPU time. 1 REAL Specifies units of real time. The Optimization Statistics table, an output of the OPTMODEL procedure if you specify PRINT- LEVEL=2 in the PROC OPTMODEL statement, also includes the same time units for Presolver Time and Solver Time. The other times (such as Problem Generation Time) in the Optimization Statistics table are also in the same units. The default value of the TIMETYPE= option depends on the value of the NTHREADS= option in the PERFORMANCE statement of the OPTMODEL procedure. Table 11.5 describes the detailed logic for determining the default; the first context in the table that applies determines the default value. For more information about the NTHREADS= option, see the section PERFORMANCE Statement on page 19 in Chapter 4, Shared Concepts and Topics. Table 11.5 Default Value for TIMETYPE= Option Context Solver is invoked in an OPTMODEL COFOR loop NTHREADS= value is greater than 1 NTHREADS= 1 Default REAL REAL CPU

12 554 Chapter 11: The Quadratic Programming Solver Details: QP Solver Interior Point Algorithm: Overview The QP solver implements an infeasible primal-dual predictor-corrector interior point algorithm. To illustrate the algorithm and the concepts of duality and dual infeasibility, consider the following QP formulation (the primal): min subject to 1 2 xt Qx C c T x Ax b x 0 The corresponding dual formulation is 1 max 2 xt Qx C b T y subject to Qx C A T y C w D c y 0 w 0 where y 2 R m refers to the vector of dual variables and w 2 R n refers to the vector of dual slack variables. The dual makes an important contribution to the certificate of optimality for the primal. The primal and dual constraints combined with complementarity conditions define the first-order optimality conditions, also known as KKT (Karush-Kuhn-Tucker) conditions, which can be stated as follows where e.1; : : : ; 1/ T of appropriate dimension and s 2 R m is the vector of primal slack variables: Ax s D b.primal feasibility/ Qx C A T y C w D c.dual feasibility/ WXe D 0.complementarity/ SYe D 0.complementarity/ x; y; w; s 0 NOTE: Slack variables (the s vector) are automatically introduced by the solver when necessary; it is therefore recommended that you not introduce any slack variables explicitly. This enables the solver to handle slack variables much more efficiently.

13 Interior Point Algorithm: Overview 555 The letters X; Y; W; and S denote matrices with corresponding x, y, w, and s on the main diagonal and zero elsewhere, as in the following example: 2 3 x x 2 0 X 6 4 : : : :: 7 : x n If.x ; y ; w ; s / is a solution of the previously defined system of equations that represent the KKT conditions, then x is also an optimal solution to the original QP model. At each iteration the interior point algorithm solves a large, sparse system of linear equations, Y 1 S A y A T Q X 1 D W x where x and y denote the vector of search directions in the primal and dual spaces, respectively, and and constitute the vector of the right-hand sides. The preceding system is known as the reduced KKT system. The QP solver uses a preconditioned quasiminimum residual algorithm to solve this system of equations efficiently. An important feature of the interior point algorithm is that it takes full advantage of the sparsity in the constraint and quadratic matrices, thereby enabling it to efficiently solve large-scale quadratic programs. The interior point algorithm works simultaneously in the primal and dual spaces. It attains optimality when both primal and dual feasibility are achieved and when complementarity conditions hold. Therefore, it is of interest to observe the following four measures where kvk 2 is the Euclidean norm of the vector v: relative primal infeasibility measure : D kax b sk 2 kbk 2 C 1 relative dual infeasibility measure ˇ: ˇ D kqx C c AT y wk 2 kck 2 C 1 relative duality gap ı: ı D jxt Qx C c T x b T yj j 1 2 xt Qx C c T xj C 1 absolute complementarity : D nx x i w i C id1 mx y i s i id1 These measures are displayed in the iteration log.

14 556 Chapter 11: The Quadratic Programming Solver Parallel Processing The interior point algorithm can be run in single-machine mode (in single-machine mode, the computation is executed by multiple threads on a single computer). You can specify options for parallel processing in the PERFORMANCE statement, which is documented in the section PERFORMANCE Statement on page 19 in Chapter 4, Shared Concepts and Topics. Iteration Log The following information is displayed in the iteration log: Iter Complement Duality Gap Primal Infeas Bound Infeas Dual Infeas Time indicates the iteration number. indicates the (absolute) complementarity. indicates the (relative) duality gap. indicates the (relative) primal infeasibility measure. indicates the (relative) bound infeasibility measure. indicates the (relative) dual infeasibility measure. indicates the time elapsed (in seconds). If the sequence of solutions converges to an optimal solution of the problem, you should see all columns in the iteration log converge to zero or very close to zero. Nonconvergence can be the result of insufficient iterations being performed to reach optimality. In this case, you might need to increase the value that you specify in the MAXITER= or MAXTIME= option. If the complementarity or the duality gap does not converge, the problem might be infeasible or unbounded. If the infeasibility columns do not converge, the problem might be infeasible. Problem Statistics Optimizers can encounter difficulty when solving poorly formulated models. Information about data magnitude provides a simple gauge to determine how well a model is formulated. For example, a model whose constraint matrix contains one very large entry (on the order of 10 9 ) can cause difficulty when the remaining entries are single-digit numbers. The PRINTLEVEL=2 option in the OPTMODEL procedure causes the ODS table ProblemStatistics to be generated when the QP solver is called. This table provides basic data magnitude information that enables you to improve the formulation of your models. The example output in Figure 11.4 demonstrates the contents of the ODS table ProblemStatistics.

15 Irreducible Infeasible Set 557 Figure 11.4 ODS Table ProblemStatistics The OPTMODEL Procedure Problem Statistics Number of Constraint Matrix Nonzeros 4 Maximum Constraint Matrix Coefficient 2 Minimum Constraint Matrix Coefficient 1 Average Constraint Matrix Coefficient 1.25 Number of Linear Objective Nonzeros 2 Maximum Linear Objective Coefficient 3 Minimum Linear Objective Coefficient 2 Average Linear Objective Coefficient 2.5 Number of Lower Triangular Hessian Nonzeros 1 Number of Diagonal Hessian Nonzeros 2 Maximum Hessian Coefficient 20 Minimum Hessian Coefficient 2 Average Hessian Coefficient 6.75 Number of RHS Nonzeros 2 Maximum RHS 100 Minimum RHS 1 Average RHS 50.5 Maximum Number of Nonzeros per Column 2 Minimum Number of Nonzeros per Column 2 Average Number of Nonzeros per Column 2 Maximum Number of Nonzeros per Row 2 Minimum Number of Nonzeros per Row 2 Average Number of Nonzeros per Row 2 Irreducible Infeasible Set For a quadratic programming problem, an irreducible infeasible set (IIS) is an infeasible subset of constraints and variable bounds that becomes feasible if any single constraint or variable bound is removed. It is possible to have more than one IIS in an infeasible QP. Identifying an IIS can help isolate the structural infeasibility in a QP. The IIS=ON option directs the QP solver to search for an IIS in a specified QP. Whether a quadratic programming problem is feasible or infeasible is determined by its constraints and variable bounds, which have nothing to do with its objective function. When you specify the IIS=ON option, the QP solver treats this problem as a linear programming problem by ignoring its objective function. Then finding IIS is the same as what the LP solver does with the IIS=ON option. See the section Irreducible Infeasible Set on page 272 in Chapter 7, The Linear Programming Solver, for more information about the irreducible infeasible set.

16 558 Chapter 11: The Quadratic Programming Solver Macro Variable _OROPTMODEL_ The OPTMODEL procedure always creates and initializes a SAS macro called _OROPTMODEL_. This variable contains a character string. After each PROC OROPTMODEL run, you can examine this macro by specifying %put &_OROPTMODEL_; and check the execution of the most recently invoked solver from the value of the macro variable. The various terms of the variable after the QP solver is called are interpreted as follows. STATUS indicates the solver status at termination. It can take one of the following values: OK SYNTAX_ERROR DATA_ERROR OUT_OF_MEMORY IO_ERROR SEMANTIC_ERROR ERROR The solver terminated normally. Incorrect syntax was used. The input data were inconsistent. Insufficient memory was allocated to the procedure. A problem occurred in reading or writing data. An evaluation error, such as an invalid operand type, occurred. The status cannot be classified into any of the preceding categories. ALGORITHM indicates the algorithm that produced the solution data in the macro variable. This term only appears when STATUS=OK. It can take the following value: IP The interior point algorithm produced the solution data. SOLUTION_STATUS indicates the solution status at termination. It can take one of the following values: OPTIMAL CONDITIONAL_OPTIMAL INFEASIBLE UNBOUNDED INFEASIBLE_OR_UNBOUNDED BAD_PROBLEM_TYPE ITERATION_LIMIT_REACHED TIME_LIMIT_REACHED FUNCTION_CALL_LIMIT_REACHED INTERRUPTED FAILED The solution is optimal. The solution is optimal, but some infeasibilities (primal, dual or bound) exceed tolerances due to scaling or preprocessing. The problem is infeasible. The problem is unbounded. The problem is infeasible or unbounded. The problem type is unsupported by the solver. The maximum allowable number of iterations was reached. The solver reached its execution time limit. The solver reached its limit on function evaluations. The solver was interrupted externally. The solver failed to converge, possibly due to numerical issues.

17 Examples: QP Solver 559 OBJECTIVE indicates the objective value obtained by the solver at termination. PRIMAL_INFEASIBILITY indicates the (relative) infeasibility of the primal constraints at the solution. See the section Interior Point Algorithm: Overview on page 554 for details. DUAL_INFEASIBILITY indicates the (relative) infeasibility of the dual constraints at the solution. See the section Interior Point Algorithm: Overview on page 554 for details. BOUND_INFEASIBILITY indicates the (relative) violation by the solution of the lower or upper bounds (or both). See the section Interior Point Algorithm: Overview on page 554 for details. DUALITY_GAP indicates the (relative) duality gap. See the section Interior Point Algorithm: Overview on page 554 for details. COMPLEMENTARITY indicates the (absolute) complementarity at the solution. See the section Interior Point Algorithm: Overview on page 554 for details. ITERATIONS indicates the number of iterations required to solve the problem. PRESOLVE_TIME indicates the time taken for preprocessing (in seconds). SOLUTION_TIME indicates the time (in seconds) taken to solve the problem, including preprocessing time. NOTE: The time that is reported in PRESOLVE_TIME and SOLUTION_TIME is either CPU time or real time. The type is determined by the TIMETYPE= option. Examples: QP Solver This section presents examples that illustrate the use of the OPTMODEL procedure to solve quadratic programming problems. Example 11.1 illustrates how to model a linear least squares problem and solve it by using PROC OPTMODEL. Example 11.2 and Example 11.3 show in detail how to model the portfolio optimization and selection problems. Example 11.1: Linear Least Squares Problem The linear least squares problem arises in the context of determining a solution to an overdetermined set of linear equations. In practice, these equations could arise in data fitting and estimation problems. An overdetermined system of linear equations can be defined as Ax D b

18 560 Chapter 11: The Quadratic Programming Solver where A 2 R mn, x 2 R n, b 2 R m, and m > n. Since this system usually does not have a solution, you need to be satisfied with some sort of approximate solution. The most widely used approximation is the least squares solution, which minimizes kax bk 2 2. This problem is called a least squares problem for the following reason. Let A, x, and b be defined as previously. Let k i.x/ be the kth component of the vector Ax b: k i.x/ D a i1 x 1 C a i2 x 2 C C a in x n b i ; i D 1; 2; : : : ; m By definition of the Euclidean norm, the objective function can be expressed as follows: kax m bk 2 2 D X k i.x/ 2 id1 Therefore, the function you minimize is the sum of squares of m terms k i.x/; hence the term least squares. The following example is an illustration of the linear least squares problem; that is, each of the terms k i is a linear function of x. Consider the following least squares problem defined by A D ; b D This translates to the following set of linear equations: 4x 1 D 1; x 1 C x 2 D 0; 3x 1 C 2x 2 D 1 The corresponding least squares problem is: minimize.4x 1 1/ 2 C. x 1 C x 2 / 2 C.3x 1 C 2x 2 1/ 2 The preceding objective function can be expanded to: minimize 26x 2 1 C 5x2 2 C 10x 1x 2 14x 1 4x 2 C 2 In addition, you impose the following constraint so that the equation 3x 1 C 2x 2 D 1 is satisfied within a tolerance of 0.1: 0:9 3x 1 C 2x 2 1:1 You can use the following SAS statements to solve the least squares problem: /* example 1: linear least squares problem */ proc optmodel; /* declare free (no explicit bounds) variables x[1] and x[2] */ var x {1..2}; /* objective function: minimize the sum of squares */ minimize f = 26*x[1]^2 + 5*x[2]^2 + 10*x[1]*x[2] - 14*x[1] - 4*x[2] + 2; /* subject to the following constraint */ con R: 0.9 <= 3*x[1] + 2*x[2] <= 1.1;

19 Example 11.1: Linear Least Squares Problem 561 /* call the QP solver */ solve; /* print the optimal solution */ print x; quit; The output is shown in Output Output Summaries and Optimal Solution The OPTMODEL Procedure Objective Sense Objective Function Objective Type Problem Summary Minimization f Quadratic Number of Variables 2 Bounded Above 0 Bounded Below 0 Bounded Below and Above 0 Free 2 Fixed 0 Number of Constraints 1 Linear LE (<=) 0 Linear EQ (=) 0 Linear GE (>=) 0 Linear Range 1 Constraint Coefficients 2 Performance Information Execution Mode Single-Machine Number of Threads 4

20 562 Chapter 11: The Quadratic Programming Solver Output continued Solution Summary Solver QP Algorithm Interior Point Objective Function f Solution Status Optimal Objective Value Primal Infeasibility 0 Dual Infeasibility E-7 Bound Infeasibility 0 Duality Gap E-7 Complementarity 0 Iterations 3 Presolve Time 0.00 Solution Time 0.01 [1] x Example 11.2: Portfolio Optimization Consider a portfolio optimization example. The two competing goals of investment are (1) long-term growth of capital and (2) low risk. A good portfolio grows steadily without wild fluctuations in value. The Markowitz model is an optimization model for balancing the return and risk of a portfolio. The decision variables are the amounts invested in each asset. The objective is to minimize the variance of the portfolio s total return, subject to the constraints that (1) the expected growth of the portfolio reaches at least some target level and (2) you do not invest more capital than you have. Let x 1 ; : : : ; x n be the amount invested in each asset, B be the amount of capital you have, R be the random vector of asset returns over some period, and r be the expected value of R. Let G be the minimum growth np you hope to obtain, and C be the covariance matrix of R. The objective function is Var x i R i, which can be equivalently denoted as x T Cx. Assume, for example, n = 4. Let B = 10,000, G = 1,000, r D Œ0:05; 2 3 0:08 0:05 0:05 0:05 C D 6 0:05 0:16 0:02 0: :05 0:02 0:35 0:06 5 0:05 0:02 0:06 0:35 0:2; 0:15; 0:30, and id1

21 Example 11.2: Portfolio Optimization 563 The QP formulation can be written as: min 0:08x 2 1 0:1x 1 x 2 0:1x 1 x 3 0:1x 1 x 4 C 0:16x 2 2 0:04x 2 x 3 0:04x 2 x 4 C 0:35x 2 3 C 0:12x 3x 4 C 0:35x 2 4 subject to.budget/ x 1 C x 2 C x 3 C x growth/ 0:05x 1 0:2x 2 C 0:15x 3 C 0:30x x 1 ; x 2 ; x 3 ; x 4 0 Use the following SAS statements to solve the problem: /* example 2: portfolio optimization */ proc optmodel; /* let x1, x2, x3, x4 be the amount invested in each asset */ var x{1..4} >= 0; num coeff{1..4, 1..4} = [ ]; num r{1..4}=[ ]; /* minimize the variance of the portfolio's total return */ minimize f = sum{i in 1..4, j in 1..4}coeff[i,j]*x[i]*x[j]; /* subject to the following constraints */ con BUDGET: sum{i in 1..4}x[i] <= 10000; con GROWTH: sum{i in 1..4}r[i]*x[i] >= 1000; solve with qp; /* print the optimal solution */ print x; The summaries and the optimal solution are shown in Output

22 564 Chapter 11: The Quadratic Programming Solver Output Portfolio Optimization The OPTMODEL Procedure Objective Sense Objective Function Objective Type Problem Summary Minimization f Quadratic Number of Variables 4 Bounded Above 0 Bounded Below 4 Bounded Below and Above 0 Free 0 Fixed 0 Number of Constraints 2 Linear LE (<=) 1 Linear EQ (=) 0 Linear GE (>=) 1 Linear Range 0 Constraint Coefficients 8 Performance Information Execution Mode Single-Machine Number of Threads 4 Solution Summary Solver QP Algorithm Interior Point Objective Function f Solution Status Optimal Objective Value Primal Infeasibility 0 Dual Infeasibility E-14 Bound Infeasibility 0 Duality Gap E-16 Complementarity 0 Iterations 7 Presolve Time 0.00 Solution Time 0.01 [1] x

23 Example 11.2: Portfolio Optimization 565 Thus, the minimum variance portfolio that earns an expected return of at least 10% is x 1 = 3,452, x 2 = 0, x 3 = 1,068, x 4 D 2; 223. Asset 2 gets nothing because its expected return is 20% and its covariance with the other assets is not sufficiently negative for it to bring any diversification benefits. What if you drop the nonnegativity assumption? Financially, that means you are allowed to short-sell that is, sell low-mean-return assets and use the proceeds to invest in high-mean-return assets. In other words, you put a negative portfolio weight in low-mean assets and more than 100% in high-mean assets. To solve the portfolio optimization problem with the short-sale option, continue to submit the following SAS statements: /* example 2: portfolio optimization with short-sale option */ /* dropping nonnegativity assumption */ for {i in 1..4} x[i].lb=-x[i].ub; solve with qp; /* print the optimal solution */ print x; quit; You can see in the optimal solution displayed in Output that the decision variable x 2, denoting Asset 2, is equal to 1,563.61, which means short sale of that asset. Output Portfolio Optimization with Short-Sale Option The OPTMODEL Procedure Solver Algorithm Solution Summary Objective Function Solution Status QP Interior Point f Optimal Objective Value Primal Infeasibility E-13 Dual Infeasibility E-8 Bound Infeasibility 0 Duality Gap E-11 Complementarity 0 Iterations 5 Presolve Time 0.00 Solution Time 0.01 [1] x

24 566 Chapter 11: The Quadratic Programming Solver Example 11.3: Portfolio Selection with Transactions Consider a portfolio selection problem with a slight modification. You are now required to take into account the current position and transaction costs associated with buying and selling assets. The objective is to find the minimum variance portfolio. In order to understand the scenario better, consider the following data. You are given three assets. The current holding of the three assets is denoted by the vector c = [200, 300, 500], the amount of asset bought and sold is denoted by b i and s i, respectively, and the net investment in each asset is denoted by x i and is defined by the following relation: x i b i C s i D c i ; i D 1; 2; 3 Suppose that you pay a transaction fee of 0.01 every time you buy or sell. Let the covariance matrix C be defined as 2 0: : :00015 C D 4 0: : : : : : Assume that you hope to obtain at least 12% growth. Let r = [ , , ] be the vector of expected return on the three assets, and let B=1000 be the available funds. Mathematically, this problem can be written in the following manner: min 0:027489x 2 1 0:01748x 1 x 2 0:0003x 1 x 3 C 0:109449x 2 2 0:00024x 2 x 3 C 0:000766x 2 3 subject to P.return/ 3 id1 r ix i 1:12B P.budget/ 3 id1 x i C P 3 id1 0:01.b i C s i / D B.balance/ x i b i C s i D c i ; i D 1; 2; 3 x i ; b i ; s i 0; i D 1; 2; 3

25 Example 11.3: Portfolio Selection with Transactions 567 The problem can be solved by the following SAS statements: /* example 3: portfolio selection with transactions */ proc optmodel; /* let x1, x2, x3 be the amount invested in each asset */ var x{1..3} >= 0; /* let b1, b2, b3 be the amount of asset bought */ var b{1..3} >= 0; /* let s1, s2, s3 be the amount of asset sold */ var s{1..3} >= 0; /* current holdings */ num c{1..3}=[ ]; /* covariance matrix */ num coeff{1..3, 1..3} = [ ]; /* returns */ num r{1..3}=[ ]; /* minimize the variance of the portfolio's total return */ minimize f = sum{i in 1..3, j in 1..3}coeff[i,j]*x[i]*x[j]; /* subject to the following constraints */ con BUDGET: sum{i in 1..3}(x[i]+.01*b[i]+.01*s[i]) <= 1000; con RETURN: sum{i in 1..3}r[i]*x[i] >= 1120; con BALANC{i in 1..3}: x[i]-b[i]+s[i]=c[i]; solve with qp; /* print the optimal solution */ print x; quit; The output is displayed in Output

26 568 Chapter 11: The Quadratic Programming Solver Output Portfolio Selection with Transactions The OPTMODEL Procedure Objective Sense Objective Function Objective Type Problem Summary Minimization f Quadratic Number of Variables 9 Bounded Above 0 Bounded Below 9 Bounded Below and Above 0 Free 0 Fixed 0 Number of Constraints 5 Linear LE (<=) 1 Linear EQ (=) 3 Linear GE (>=) 1 Linear Range 0 Constraint Coefficients 21 Performance Information Execution Mode Single-Machine Number of Threads 4 Solution Summary Solver QP Algorithm Interior Point Objective Function f Solution Status Optimal Objective Value Primal Infeasibility E-17 Dual Infeasibility E-14 Bound Infeasibility 0 Duality Gap E-15 Complementarity 0 Iterations 11 Presolve Time 0.00 Solution Time 0.01 [1] x

27 References Freund, R. W. (1991). On Polynomial Preconditioning and Asymptotic Convergence Factors for Indefinite Hermitian Matrices. Linear Algebra and Its Applications : Freund, R. W., and Jarre, F. (1997). A QMR-Based Interior Point Algorithm for Solving Linear Programs. Mathematical Programming 76: Freund, R. W., and Nachtigal, N. M. (1996). QMRPACK: A Package of QMR Algorithms. ACM Transactions on Mathematical Software 22: Vanderbei, R. J. (1999). LOQO: An Interior Point Code for Quadratic Programming. Optimization Methods and Software 11: Wright, S. J. (1997). Primal-Dual Interior-Point Methods. Philadelphia: SIAM. 569

28 Subject Index macro variable _OROPTMODEL_, 558 OPTMODEL procedure, QP solver functional summary, 551 macro variable _OROPTMODEL_, 558 _OROPTMODEL_ macro variable, 558 positive semidefinite matrix, 546 QP Solver examples, 559 interior point algorithm overview, 554 iteration log, 556 QP solver IIS, 557 problem statistics, 556 QP solver examples covariance matrix, 562 data fitting, 559 estimation, 559 linear least squares, 559 Markowitz model, 562 portfolio optimization, 562 portfolio selection with transactions, 566 short-sell, 565 quadratic programming quadratic matrix, 546 SOLVE WITH QP statement dual infeasibility, 553 duality gap, 553 primal infeasibility, 553

29

30 Syntax Index IIS= option SOLVE WITH QP statement, 551 LOGFREQ= option SOLVE WITH QP statement, 552 MAXITER= option SOLVE WITH QP statement, 552 MAXTIME= option SOLVE WITH QP statement, 552 OPTMODEL procedure, QP solver syntax, 551 PRESOLVER= option SOLVE WITH QP statement, 552 PRINTFREQ= option SOLVE WITH QP statement, 552 SOLVE WITH QP statement IIS= option, 551 LOGFREQ= option, 552 MAXITER= option, 552 MAXTIME= option, 552 PRESOLVER= option, 552 PRINTFREQ= option, 552 STOP_DG= option, 553 STOP_DI= option, 553 STOP_PI= option, 553 TIMETYPE= option, 553 STOP_DG= option SOLVE WITH QP statement, 553 STOP_DI= option SOLVE WITH QP statement, 553 STOP_PI= option SOLVE WITH QP statement, 553 TIMETYPE= option SOLVE WITH QP statement, 553

### SAS/STAT 13.1 User s Guide. The NESTED Procedure

SAS/STAT 13.1 User s Guide The NESTED Procedure This document is an individual chapter from SAS/STAT 13.1 User s Guide. The correct bibliographic citation for the complete manual is as follows: SAS Institute

### The Mixed Integer Linear Programming Solver

SAS/OR 14.3 User s Guide Mathematical Programming The Mixed Integer Linear Programming Solver This document is an individual chapter from SAS/OR 14.3 User s Guide: Mathematical Programming. The correct

### Getting Started with SAS Factory Miner 14.2

Getting Started with SAS Factory Miner 14.2 SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016. Getting Started with SAS Factory Miner 14.2. Cary,

### SAS Simulation Studio 14.1: User s Guide. Introduction to SAS Simulation Studio

SAS Simulation Studio 14.1: User s Guide Introduction to SAS Simulation Studio This document is an individual chapter from SAS Simulation Studio 14.1: User s Guide. The correct bibliographic citation for

### Lecture 25 Nonlinear Programming. November 9, 2009

Nonlinear Programming November 9, 2009 Outline Nonlinear Programming Another example of NLP problem What makes these problems complex Scalar Function Unconstrained Problem Local and global optima: definition,

### Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems

Paper SAS1502-2015 Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems Rob Pratt, SAS Institute Inc. ABSTRACT Mathematical optimization is a powerful paradigm for modeling and solving business

### SAS Enterprise Miner : Tutorials and Examples

SAS Enterprise Miner : Tutorials and Examples SAS Documentation February 13, 2018 The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2017. SAS Enterprise Miner : Tutorials

### SAS Federation Server 4.2: Migration Guide

SAS Federation Server 4.2: Migration Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016. SAS Federation Server 4.2: Migration Guide. Cary,

### CSDP User s Guide. Brian Borchers. August 15, 2006

CSDP User s Guide Brian Borchers August 5, 6 Introduction CSDP is a software package for solving semidefinite programming problems. The algorithm is a predictor corrector version of the primal dual barrier

### Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems. Rob Pratt, Senior R&D Manager, SAS/OR

Using the OPTMODEL Procedure in SAS/OR to Solve Complex Problems Rob Pratt, Senior R&D Manager, SAS/OR Outline 1 Recent Features in PROC OPTMODEL 2 Graph Partitioning with Connectivity Constraints 2 /

### BCN Decision and Risk Analysis. Syed M. Ahmed, Ph.D.

Linear Programming Module Outline Introduction The Linear Programming Model Examples of Linear Programming Problems Developing Linear Programming Models Graphical Solution to LP Problems The Simplex Method

### Linear Programming with Bounds

Chapter 481 Linear Programming with Bounds Introduction Linear programming maximizes (or minimizes) a linear objective function subject to one or more constraints. The technique finds broad use in operations

### SAS Clinical Data Integration 2.4

SAS Clinical Data Integration 2.4 User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2013. SAS Clinical Data Integration 2.4: User's Guide.

### Applied Lagrange Duality for Constrained Optimization

Applied Lagrange Duality for Constrained Optimization Robert M. Freund February 10, 2004 c 2004 Massachusetts Institute of Technology. 1 1 Overview The Practical Importance of Duality Review of Convexity

### SAS University Edition: Installation Guide for Windows

SAS University Edition: Installation Guide for Windows i 17 June 2014 The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2014. SAS University Edition: Installation Guide

### SAS Contextual Analysis 14.3: Administrator s Guide

SAS Contextual Analysis 14.3: Administrator s Guide SAS Documentation August 25, 2017 The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2017. SAS Contextual Analysis

### AM 121: Intro to Optimization Models and Methods Fall 2017

AM 121: Intro to Optimization Models and Methods Fall 2017 Lecture 10: Dual Simplex Yiling Chen SEAS Lesson Plan Interpret primal simplex in terms of pivots on the corresponding dual tableau Dictionaries

### MATLAB Solution of Linear Programming Problems

MATLAB Solution of Linear Programming Problems The simplex method is included in MATLAB using linprog function. All is needed is to have the problem expressed in the terms of MATLAB definitions. Appendix

### Grid Computing in SAS 9.4

Grid Computing in SAS 9.4 SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2013. Grid Computing in SAS 9.4. Cary, NC: SAS Institute Inc. Grid Computing

### Discrete Optimization. Lecture Notes 2

Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The

### SAS. Social Network Analysis Server 6.2: Installation and Configuration Guide, Third Edition. SAS Documentation

SAS Social Network Analysis Server 6.2: Installation and Configuration Guide, Third Edition SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016.

### SAS Add-In 7.1 for Microsoft Office: Getting Started in Microsoft Excel, Microsoft Word, and Microsoft PowerPoint, Second Edition

SAS Add-In 7.1 for Microsoft Office: Getting Started in Microsoft Excel, Microsoft Word, and Microsoft PowerPoint, Second Edition SAS Documentation The correct bibliographic citation for this manual is

### Programs. Introduction

16 Interior Point I: Linear Programs Lab Objective: For decades after its invention, the Simplex algorithm was the only competitive method for linear programming. The past 30 years, however, have seen

### CSDP User s Guide. Brian Borchers. November 1, 2006

CSDP User s Guide Brian Borchers November, 6 Introduction CSDP is a software package for solving semidefinite programming problems. The algorithm is a predictor corrector version of the primal dual barrier

### Appendix A MATLAB s Optimization Toolbox Algorithms

Appendix A MATLAB s Optimization Toolbox Algorithms Abstract MATLAB s Optimization Toolbox (version 7:2) includes a family of algorithms for solving optimization problems. The toolbox provides functions

### CSCI 5654 (Fall 2013): Network Simplex Algorithm for Transshipment Problems.

CSCI 5654 (Fall 21): Network Simplex Algorithm for Transshipment Problems. Sriram Sankaranarayanan November 14, 21 1 Introduction e will present the basic ideas behind the network Simplex algorithm for

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

Introduction to Mathematical Programming IE406 Lecture 20 Dr. Ted Ralphs IE406 Lecture 20 1 Reading for This Lecture Bertsimas Sections 10.1, 11.4 IE406 Lecture 20 2 Integer Linear Programming An integer

### SAS/STAT 13.1 User s Guide. The SCORE Procedure

SAS/STAT 13.1 User s Guide The SCORE Procedure This document is an individual chapter from SAS/STAT 13.1 User s Guide. The correct bibliographic citation for the complete manual is as follows: SAS Institute

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

. Tuesday, April The Network Simplex Method for Solving the Minimum Cost Flow Problem Quotes of the day I think that I shall never see A poem lovely as a tree. -- Joyce Kilmer Knowing trees, I understand

### CSDP User s Guide

CSDP 6.. User s Guide Brian Borchers July 4, 7 Introduction CSDP is a software package for solving semidefinite programming problems. The algorithm is a predictor corrector version of the primal dual barrier

### Section Notes 5. Review of Linear Programming. Applied Math / Engineering Sciences 121. Week of October 15, 2017

Section Notes 5 Review of Linear Programming Applied Math / Engineering Sciences 121 Week of October 15, 2017 The following list of topics is an overview of the material that was covered in the lectures

### The EMCLUS Procedure. The EMCLUS Procedure

The EMCLUS Procedure Overview Procedure Syntax PROC EMCLUS Statement VAR Statement INITCLUS Statement Output from PROC EMCLUS EXAMPLES-SECTION Example 1: Syntax for PROC FASTCLUS Example 2: Use of the

### SAS/ETS 13.2 User s Guide. The TIMEID Procedure

SAS/ETS 13.2 User s Guide The TIMEID Procedure This document is an individual chapter from SAS/ETS 13.2 User s Guide. The correct bibliographic citation for the complete manual is as follows: SAS Institute

### Linear Programming MT 4.0

Linear Programming MT 4.0 for GAUSS TM Mathematical and Statistical System Aptech Systems, Inc. Information in this document is subject to change without notice and does not represent a commitment on the

### DataFlux Web Studio 2.5. Installation and Configuration Guide

DataFlux Web Studio 2.5 Installation and Configuration Guide The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2014. DataFlux Web Studio 2.5: Installation and Configuration

### Chapter 1 Introduction to Optimization

Chapter 1 Introduction to Optimization Chapter Contents OVERVIEW................................... 15 DATA FLOW................................... 16 PROC LP................................... 17 PROC

### SAS. OnDemand for Academics: User s Guide. SAS Documentation

SAS OnDemand for Academics: User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016. SAS OnDemand for Academics: User s Guide. Cary, NC:

### A Short SVM (Support Vector Machine) Tutorial

A Short SVM (Support Vector Machine) Tutorial j.p.lewis CGIT Lab / IMSC U. Southern California version 0.zz dec 004 This tutorial assumes you are familiar with linear algebra and equality-constrained optimization/lagrange

### 5. DUAL LP, SOLUTION INTERPRETATION, AND POST-OPTIMALITY

5. DUAL LP, SOLUTION INTERPRETATION, AND POST-OPTIMALITY 5.1 DUALITY Associated with every linear programming problem (the primal) is another linear programming problem called its dual. If the primal involves

SAS Data Loader 2.4 for Hadoop vapp Deployment Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2015. SAS Data Loader 2.4 for Hadoop: vapp Deployment

### 3 INTEGER LINEAR PROGRAMMING

3 INTEGER LINEAR PROGRAMMING PROBLEM DEFINITION Integer linear programming problem (ILP) of the decision variables x 1,..,x n : (ILP) subject to minimize c x j j n j= 1 a ij x j x j 0 x j integer n j=

### Mathematical Programming and Research Methods (Part II)

Mathematical Programming and Research Methods (Part II) 4. Convexity and Optimization Massimiliano Pontil (based on previous lecture by Andreas Argyriou) 1 Today s Plan Convex sets and functions Types

### 60 2 Convex sets. {x a T x b} {x ã T x b}

60 2 Convex sets Exercises Definition of convexity 21 Let C R n be a convex set, with x 1,, x k C, and let θ 1,, θ k R satisfy θ i 0, θ 1 + + θ k = 1 Show that θ 1x 1 + + θ k x k C (The definition of convexity

### PRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING. 1. Introduction

PRIMAL-DUAL INTERIOR POINT METHOD FOR LINEAR PROGRAMMING KELLER VANDEBOGERT AND CHARLES LANNING 1. Introduction Interior point methods are, put simply, a technique of optimization where, given a problem

### Time Series Studio 13.1

SAS Time Series Studio 13.1 Administrator s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2013. SAS Time Series Studio 13.1: Administrator's

### SAS Visual Analytics 7.3: Installation and Configuration Guide (Distributed SAS LASR )

SAS Visual Analytics 7.3: Installation and Configuration Guide (Distributed SAS LASR ) SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2015. SAS Visual

### SAS Theme Designer 4.7 for Flex

SAS Theme Designer 4.7 for Flex User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2015. SAS Theme Designer 4.7 for Flex: User's Guide.

### 16.410/413 Principles of Autonomy and Decision Making

16.410/413 Principles of Autonomy and Decision Making Lecture 17: The Simplex Method Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology November 10, 2010 Frazzoli (MIT)

### ADDENDUM TO THE SEDUMI USER GUIDE VERSION 1.1

ADDENDUM TO THE SEDUMI USER GUIDE VERSION 1.1 IMRE PÓLIK 1. Introduction The main goal of this reference guide is to give a summary of all the options in SeDuMi. The default value of the options is satisfactory

### MLR Institute of Technology

Course Name : Engineering Optimization Course Code : 56021 Class : III Year Branch : Aeronautical Engineering Year : 2014-15 Course Faculty : Mr Vamsi Krishna Chowduru, Assistant Professor Course Objective

### SAS/FSP 9.2. Procedures Guide

SAS/FSP 9.2 Procedures Guide The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2008. SAS/FSP 9.2 Procedures Guide. Cary, NC: SAS Institute Inc. SAS/FSP 9.2 Procedures

### Affine function. suppose f : R n R m is affine (f(x) =Ax + b with A R m n, b R m ) the image of a convex set under f is convex

Affine function suppose f : R n R m is affine (f(x) =Ax + b with A R m n, b R m ) the image of a convex set under f is convex S R n convex = f(s) ={f(x) x S} convex the inverse image f 1 (C) of a convex

### OPERATIONS RESEARCH. Linear Programming Problem

OPERATIONS RESEARCH Chapter 1 Linear Programming Problem Prof. Bibhas C. Giri Department of Mathematics Jadavpur University Kolkata, India Email: bcgiri.jumath@gmail.com 1.0 Introduction Linear programming

### 3 No-Wait Job Shops with Variable Processing Times

3 No-Wait Job Shops with Variable Processing Times In this chapter we assume that, on top of the classical no-wait job shop setting, we are given a set of processing times for each operation. We may select

### Lecture 7: Support Vector Machine

Lecture 7: Support Vector Machine Hien Van Nguyen University of Houston 9/28/2017 Separating hyperplane Red and green dots can be separated by a separating hyperplane Two classes are separable, i.e., each

### SAS 9.4 Data Quality Accelerator for Teradata: User s Guide

SAS 9.4 Data Quality Accelerator for Teradata: User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016. SAS 9.4 Data Quality Accelerator

### SAS. Contextual Analysis 13.2: User s Guide. SAS Documentation

SAS Contextual Analysis 13.2: User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2014. SAS Contextual Analysis 13.2: User's Guide. Cary,

### Demand Classification and Clustering 6.1

SAS Demand Classification and Clustering 6.1 User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2014. SAS Demand Classification and Clustering

### Generalized Network Flow Programming

Appendix C Page Generalized Network Flow Programming This chapter adapts the bounded variable primal simplex method to the generalized minimum cost flow problem. Generalized networks are far more useful

### MATH3016: OPTIMIZATION

MATH3016: OPTIMIZATION Lecturer: Dr Huifu Xu School of Mathematics University of Southampton Highfield SO17 1BJ Southampton Email: h.xu@soton.ac.uk 1 Introduction What is optimization? Optimization is

### Linear programming and duality theory

Linear programming and duality theory Complements of Operations Research Giovanni Righini Linear Programming (LP) A linear program is defined by linear constraints, a linear objective function. Its variables

### Using AMPL/MINOS. 1 Applicability

Using AMPL/MINOS MINOS is an optimization package for linear and nonlinear mathematical programs in continuous variables. This supplement to AMPL: A Modeling Language for Mathematical Programming describes

### Clinical Standards Toolkit 1.7

SAS Clinical Standards Toolkit 1.7 Migration Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2015. SAS Clinical Standards Toolkit 1.7: Migration

### Cutting Planes for Some Nonconvex Combinatorial Optimization Problems

Cutting Planes for Some Nonconvex Combinatorial Optimization Problems Ismael Regis de Farias Jr. Department of Industrial Engineering Texas Tech Summary Problem definition Solution strategy Multiple-choice

### A Nonlinear Presolve Algorithm in AIMMS

A Nonlinear Presolve Algorithm in AIMMS By Marcel Hunting marcel.hunting@aimms.com November 2011 This paper describes the AIMMS presolve algorithm for nonlinear problems. This presolve algorithm uses standard

### SAS/OR 14.1 User s Guide: Local Search Optimization. The GA Procedure

SAS/OR 14.1 User s Guide: Local Search Optimization The GA Procedure This document is an individual chapter from SAS/OR 14.1 User s Guide: Local Search Optimization. The correct bibliographic citation

### mai Installation Instructions for SAS 9.4 Electronic Software Delivery for Basic Installations on z /OS

mai Installation Instructions for SAS 9.4 Electronic Software Delivery for Basic Installations on z /OS The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2017. Installation

### The Simplex Algorithm

The Simplex Algorithm Uri Feige November 2011 1 The simplex algorithm The simplex algorithm was designed by Danzig in 1947. This write-up presents the main ideas involved. It is a slight update (mostly

### Bilinear Programming

Bilinear Programming Artyom G. Nahapetyan Center for Applied Optimization Industrial and Systems Engineering Department University of Florida Gainesville, Florida 32611-6595 Email address: artyom@ufl.edu

### Lecture 16 October 23, 2014

CS 224: Advanced Algorithms Fall 2014 Prof. Jelani Nelson Lecture 16 October 23, 2014 Scribe: Colin Lu 1 Overview In the last lecture we explored the simplex algorithm for solving linear programs. While

### Convexity: an introduction

Convexity: an introduction Geir Dahl CMA, Dept. of Mathematics and Dept. of Informatics University of Oslo 1 / 74 1. Introduction 1. Introduction what is convexity where does it arise main concepts and

SAS Data Loader 2.4 for Hadoop: User s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016. SAS Data Loader 2.4 for Hadoop: User s Guide. Cary,

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

Approximation Algorithms: The Primal-Dual Method My T. Thai 1 Overview of the Primal-Dual Method Consider the following primal program, called P: min st n c j x j j=1 n a ij x j b i j=1 x j 0 Then the

### Module 1 Lecture Notes 2. Optimization Problem and Model Formulation

Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization

### Math 5593 Linear Programming Lecture Notes

Math 5593 Linear Programming Lecture Notes Unit II: Theory & Foundations (Convex Analysis) University of Colorado Denver, Fall 2013 Topics 1 Convex Sets 1 1.1 Basic Properties (Luenberger-Ye Appendix B.1).........................

### Theoretical Concepts of Machine Learning

Theoretical Concepts of Machine Learning Part 2 Institute of Bioinformatics Johannes Kepler University, Linz, Austria Outline 1 Introduction 2 Generalization Error 3 Maximum Likelihood 4 Noise Models 5

### 5 The Theory of the Simplex Method

5 The Theory of the Simplex Method Chapter 4 introduced the basic mechanics of the simplex method. Now we shall delve a little more deeply into this algorithm by examining some of its underlying theory.

### A New Unblocking Technique to Warmstart Interior Point Methods based on Sensitivity Analysis. December 19th, 2006 revised June 1, 2007 MS

A New Unblocking Technique to Warmstart Interior Point Methods based on Sensitivity Analysis Jacek Gondzio Andreas Grothey December 19th, 2006 revised June 1, 2007 MS-06-005 1 A New Unblocking Technique

### Convex sets and convex functions

Convex sets and convex functions Convex optimization problems Convex sets and their examples Separating and supporting hyperplanes Projections on convex sets Convex functions, conjugate functions ECE 602,

### SAS Studio 3.5. Developer s Guide to Repositories. SAS Documentation

SAS Studio 3.5 Developer s Guide to Repositories SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2016. SAS Studio 3.5: Developer's Guide to Repositories.

### APPENDIX 2 Customizing SAS/ASSIST Software

241 APPENDIX 2 Customizing SAS/ASSIST Software Introduction 241 Setting User Profile Options 241 Creating an Alternate Menu Bar 243 Introduction This appendix describes how you can customize your SAS/ASSIST

### Lecture 2. Topology of Sets in R n. August 27, 2008

Lecture 2 Topology of Sets in R n August 27, 2008 Outline Vectors, Matrices, Norms, Convergence Open and Closed Sets Special Sets: Subspace, Affine Set, Cone, Convex Set Special Convex Sets: Hyperplane,

### Introduction to Constrained Optimization

Introduction to Constrained Optimization Duality and KKT Conditions Pratik Shah {pratik.shah [at] lnmiit.ac.in} The LNM Institute of Information Technology www.lnmiit.ac.in February 13, 2013 LNMIIT MLPR

### Using Analytic QP and Sparseness to Speed Training of Support Vector Machines

Using Analytic QP and Sparseness to Speed Training of Support Vector Machines John C. Platt Microsoft Research 1 Microsoft Way Redmond, WA 9805 jplatt@microsoft.com Abstract Training a Support Vector Machine

### VARIANTS OF THE SIMPLEX METHOD

C H A P T E R 6 VARIANTS OF THE SIMPLEX METHOD By a variant of the Simplex Method (in this chapter) we mean an algorithm consisting of a sequence of pivot steps in the primal system using alternative rules

### Object-Oriented Software for Quadratic Programming

Object-Oriented Software for Quadratic Programming E. MICHAEL GERTZ and STEPHEN J. WRIGHT University of Wisconsin-Madison The object-oriented software package OOQP for solving convex quadratic programming

### SAS IT Resource Management 3.3

SAS IT Resource Management 3.3 Migration Documentation SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2012. SAS IT Resource Management 3.3: Migration

SAS Business Rules Manager 3.1: Administrator s Guide SAS Documentation The correct bibliographic citation for this manual is as follows: SAS Institute Inc. 2015. SAS Business Rules Manager 3.1: Administrator

### qpoases - Online Active Set Strategy for Fast Linear MPC

qpoases - Online Active Set Strategy for Fast Linear MPC Moritz Diehl, Hans Joachim Ferreau, Lieboud Vanden Broeck, Jan Swevers Dept. ESAT and Center of Excellence for Optimization in Engineering OPTEC

### Machine Learning: Think Big and Parallel

Day 1 Inderjit S. Dhillon Dept of Computer Science UT Austin CS395T: Topics in Multicore Programming Oct 1, 2013 Outline Scikit-learn: Machine Learning in Python Supervised Learning day1 Regression: Least

### 5.4 Pure Minimal Cost Flow

Pure Minimal Cost Flow Problem. Pure Minimal Cost Flow Networks are especially convenient for modeling because of their simple nonmathematical structure that can be easily portrayed with a graph. This

### Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity

Discrete Optimization 2010 Lecture 5 Min-Cost Flows & Total Unimodularity Marc Uetz University of Twente m.uetz@utwente.nl Lecture 5: sheet 1 / 26 Marc Uetz Discrete Optimization Outline 1 Min-Cost Flows

### USER S GUIDE FOR TOMLAB /CONOPT 1

USER S GUIDE FOR TOMLAB /CONOPT 1 Kenneth Holmström 2, Anders O. Göran 3 and Marcus M. Edvall 4 April 25, 2007 1 More information available at the TOMLAB home page: http://tomopt.com. E-mail: tomlab@tomopt.com.

### Using SAS Files. Introduction CHAPTER 5

123 CHAPTER 5 Using SAS Files Introduction 123 SAS Data Libraries 124 Accessing SAS Files 124 Advantages of Using Librefs Rather than OpenVMS Logical Names 124 Assigning Librefs 124 Using the LIBNAME Statement

### DATA Step Debugger APPENDIX 3

1193 APPENDIX 3 DATA Step Debugger Introduction 1194 Definition: What is Debugging? 1194 Definition: The DATA Step Debugger 1194 Basic Usage 1195 How a Debugger Session Works 1195 Using the Windows 1195

### Support Vector Machines. James McInerney Adapted from slides by Nakul Verma

Support Vector Machines James McInerney Adapted from slides by Nakul Verma Last time Decision boundaries for classification Linear decision boundary (linear classification) The Perceptron algorithm Mistake