Re-solving stochastic programming models for airline revenue management

We study some mathematical programming formulations for the origin-destination model in airline revenue management. In particular, we focus on the traditional probabilistic model proposed in the literature. The approach we study consists of solving a sequence of two-stage stochastic programs with simple recourse, which can be viewed as an approximation to a multi-stage stochastic programming formulation to the seat allocation problem. Our theoretical results show that the proposed approximation is robust, in the sense that solving more successive two-stage programs can never worsen the expected revenue obtained with the corresponding allocation policy. Although intuitive, such a property is known not to hold for the traditional deterministic linear programming model found in the literature. We also show that this property does not hold for some bid-price policies. In addition, we propose a heuristic method to choose the re-solving points, rather than re-solving at equally-spaced times as customary. Numerical results are presented to illustrate the effectiveness of the proposed approach.     

Minimising total average cycle stock subject to practical constraints

Silver and Moon (J Opl Res Soc 50(8) (1999) 789–796) address the problem of minimising total average cycle stock subject to two practical constraints. They provide a dynamic programming formulation for obtaining an optimal solution and propose a simple and efficient heuristic algorithm. Hsieh (J Opl Res Soc 52(4) (2001) 463–470) proposes a 0–1 linear programming approach to the problem and a simple heuristic based on the relaxed 0–1 programming formulation. We show in this paper that the formulation of Hsieh can be improved for solving very large size instances of this inventory problem. So the mathematical approach is interesting for several reasons: the definition of the model is simple, its implementation is immediate by using a mathematical programming language together with a mixed integer programming software and the performance of the approach is excellent. Computational experiments carried out on the set of realistic examples considered in the above references are reported. We also show that the general framework for modelling given by mixed integer programming allows the initial model to be extended in several interesting directions.     

Outcome-Space Cutting-Plane Algorithm for Linear Multiplicative Programming

This article presents an outcome-space pure cutting-plane algorithm for globally solving the linear multiplicative programming problem. The framework of the algorithm is taken from a pure cutting-plane decision set-based method developed by Horst and Tuy for solving concave minimization problems. By adapting this method to an outcome-space reformulation of the linear multiplicative programming problem, rather than applying directly the method to the original decision-set formulation, it is expected that considerable computational savings can be obtained. Also, we show how additional computational benefits might be obtained by implementing the new algorithm appropriately. To illustrate the new algorithm, we apply it to the solution of a sample problem.     

Elzaki Transform Based Accelerated Homotopy Perturbation Method for Multi-dimensional Smoluchowski"s Coagulation and Coupled Coagulation-fragmentation Equations

This article aims to establish a semi-analytical approach based on the homotopy perturbation method (HPM) to find the closed form or approximated solutions for the population balance equations such as Smoluchowski"s coagulation, fragmentation, coupled coagulation-fragmentation and bivariate coagulation equations. An accelerated form of the HPM is combined with the Elzaki transformation to improve the accuracy and efficiency of the method. One of the significant advantages of the technique lies over the classic numerical methods as it allows solving the linear and non-linear differential equations without discretization. Further, it has benefits over the existing semi-analytical techniques such as Adomian decomposition method (ADM), optimized decomposition method (ODM), and homotopy analysis method (HAM) in the sense that computation of Adomian polynomials and convergence parameters are not required. The novelty of the scheme is shown by comparing the numerical findings with the existing results obtained via ADM, HPM, HAM and ODM for non-linear coagulation equation. This motivates us to extend the scheme for solving the other models mentioned above. The supremacy of the proposed scheme is demonstrated by taking several numerical examples for each problem. The error between exact and series solutions provided in graphs and tables show the accuracy and applicability of the method. In addition to this, convergence of the series solution is also the key attraction of the work.     

Differential-Algebraic Approach to Linear Programming

This paper presents a differential-algebraic approach for solving linear programming problems. The paper shows that the differential-algebraic approach is guaranteed to generate optimal solutions to linear programming problems with a superexponential convergence rate. The paper also shows that the path-following interior-point methods for solving linear programming problems can be viewed as a special case of the differential-algebraic approach. The results in this paper demonstrate that the proposed approach provides a promising alternative for solving linear programming problems.     

Parameter identification in financial market models with a feasible point SQP algorithm

The quickly moving market data in the finance industry requires a frequent parameter identification of the corresponding financial market models. In this paper we apply a special sequential quadratic programming algorithm to the calibration of typical equity market models. As it turns out, the projection of the iterates onto the feasible set can be efficiently computed by solving a semidefinite programming problem. Combining this approach with a Gauss-Newton framework leads to an efficient algorithm which allows to calibrate e.g. Heston’s stochastic volatility model in less than a half second on a usual 3 GHz desktop PC. Furthermore we present an appropriate regularization technique that stabilizes and significantly speeds up computations if the model parameters are chosen to be time-dependent.     

Exact augmented Lagrangian functions for nonlinear semidefinite programming

In this paper, we study augmented Lagrangian functions for nonlinear semidefinite programming (NSDP) problems with exactness properties. The term exact is used in the sense that the penalty parameter can be taken appropriately, so a single minimization of the augmented Lagrangian recovers a solution of the original problem. This leads to reformulations of NSDP problems into unconstrained nonlinear programming ones. Here, we first establish a unified framework for constructing these exact functions, generalizing Di Pillo and Lucidi’s work from 1996, that was aimed at solving nonlinear programming problems. Then, through our framework, we propose a practical augmented Lagrangian function for NSDP, proving that it is continuously differentiable and exact under the so-called nondegeneracy condition. We also present some preliminary numerical experiments.     

Exact SDP reformulations of adjustable robust linear programs with box uncertainties under separable quadratic decision rules via SOS representations of non-negativity

In this paper we show that two-stage adjustable robust linear programs with affinely adjustable data in the face of box data uncertainties under separable quadratic decision rules admit exact semi-definite program (SDP) reformulations in the sense that they share the same optimal values and admit a one-to-one correspondence between the optimal solutions. This result allows adjustable robust solutions of these robust linear programs to be found by simply numerically solving their SDP reformulations. We achieve this result by first proving a special sum-of-squares representation of non-negativity of a separable non-convex quadratic function over box constraints. Our reformulation scheme is illustrated via numerical experiments by applying it to an inventory-production management problem with the demand uncertainty. They demonstrate that our separable quadratic decision rule method to two-stage decision-making performs better than the single-stage approach and is capable of solving the inventory production problem with a greater degree of uncertainty in the demand.

     

Analyzing the potential of a firm: an operations research approach

An approach to analyzing the potential of a firm, which is understood as the firm's ability to provide goods or (and) services to be supplied to a marketplace under restrictions imposed by a business environment in which the firm functions, is proposed. The approach is based on using linear inequalities and, generally, mixed variables in modelling this ability for a broad spectrum of industrial, transportation, agricultural, and other types of firms and allows one to formulate problems of analyzing the potential of a firm as linear programming problems or mixed programming problems with linear constraints. This approach generalizes the one proposed by the author earlier for a more narrow class of models and allows one to effectively employ a widely available software for solving practical problems of the considered kind, especially for firms described by large scale models of mathematical programming.     

Parametric simplex algorithms for solving a special class of nonconvex minimization problems

It is shown that parametric linear programming algorithms work efficiently for a class of nonconvex quadratic programming problems called generalized linear multiplicative programming problems, whose objective function is the sum of a linear function and a product of two linear functions. Also, it is shown that the global minimum of the sum of the two linear fractional functions over a polytope can be obtained by a similar algorithm. Our numerical experiments reveal that these problems can be solved in much the same computational time as that of solving associated linear programs. Furthermore, we will show that the same approach can be extended to a more general class of nonconvex quadratic programming problems.     

Solving a class of geometric programming problems by an efficient dynamic model

In this paper, a neural network model is constructed on the basis of the duality theory, optimization theory, convex analysis theory, Lyapunov stability theory and LaSalle invariance principle to solve geometric programming (GP) problems. The main idea is to convert the GP problem into an equivalent convex optimization problem. A neural network model is then constructed for solving the obtained convex programming problem. By employing Lyapunov function approach, it is also shown that the proposed neural network model is stable in the sense of Lyapunov and it is globally convergent to an exact optimal solution of the original problem. The simulation results also show that the proposed neural network is feasible and efficient.     

On the inventory model with variable lead time and price–quantity discount

In this paper, we propose a mixed integer optimization approach for solving the inventory problem with variable lead time, crashing cost, and price–quantity discount. A linear programming relaxation based on piecewise linearization techniques is derived for the problem. It first converts non-linear terms into the sum of absolute terms, which are then linearized by goal programming techniques and linearization approaches. The proposed method can eliminate the complicated multiple-step solution process used in the traditional inventory models. In addition, the proposed model allows constraints to be added by the inventory decision-maker as deemed appropriate in real-world situations.     

Generalized linear multiplicative and fractional programming

This paper addresses itself to the algorithm for minimizing the sum of a convex function and a product of two linear functions over a polytope. It is shown that this nonconvex minimization problem can be solved by solving a sequence of convex programming problems. The basic idea of this algorithm is to embed the original problem into a problem in higher dimension and apply a parametric programming (path following) approach. Also it is shown that the same idea can be applied to a generalized linear fractional programming problem whose objective function is the sum of a convex function and a linear fractional function.     

Integer programming to schedule a hierarchical workforce with variable demands

We consider a hierarchical workforce in which a higher qualified worker can substitute for a lower qualified one, but not vice versa. Daily labor requirements within a week may vary, but each worker must receive n off-days in the week. This problem has been considered by Hung (R. Hung, Eur. J. Oper. Res. 78(1) (1994) 49–57), who discusses a necessary and sufficient condition for a labor mix to be feasible and presents a simple one-pass method that frequently gives the least cost labor mix. We show in this paper that the integer programming approach is well suited for solving this problem: the definition of the integer programming model is simple, its implementation is immediate by using, for example, the Mathematical programming language (MPL) and the integer programming solver XA, the computation times are low (generally a few seconds on a small microcomputer) and finally the powerful of the integer programming approach allows us to extend the model in two interesting directions.     

On Applied Nonlinear and Bilevel Programming or Pursuit-Evasion Games

Motivated by the benefits of discretization in optimal control problems, we consider the possibility of discretizing pursuit-evasion games. Two approaches are introduced. In the first approach, the solution of the necessary conditions of the continuous-time game is decomposed into ordinary optimal control problems that can be solved using discretization and nonlinear programming techniques. In the second approach, the game is discretized and transformed into a bilevel programming problem, which is solved using a first-order feasible direction method. Although the starting points of the approaches are different, they lead in practice to the same solution algorithm. We demonstrate the usability of the discretization by solving some open-loop representations of feedback solutions for a complex pursuit-evasion game between a realistically modeled aircraft and a missile, with terminal time as the payoff. The solutions are compared with those obtained via an indirect method.     

An unconstrained dual approach to solving Karmarkar-type linear programs using conventional barrier functions

This paper proposes an unconstrained dual approach and an efficient algorithm for solving Karmarkar-type linear programming problems. Conventional barrier functions are incorporated as a perturbation term in the derivation of the associated duality theory. An optimal solution of the original linear program can be obtained by solving a sequence of unconstrained concave programs, or be approximated by solving one such dual program with a sufficiently small perturbation parameter. A globally convergent curved-search algorithm with a quadratic rate of convergence is designed for this purpose. Based on our testing results, we find that the computational procedure is very efficient and can be a viable approach for solving linear programming problems.     

A hybrid approach combining interior-point and branch-and-bound methods applied to the problem of sugar cane waste

This paper proposes a hybrid approach for solving the multi-objective model related to the minimisation of sugar cane waste collection costs and/or the maximisation of produced energy by this waste, with the aid of strategies for solving multi-objective problems, which transform the problem into a set of single-objective problems. This approach combines the predictor-corrector primal-dual interior-point and branch-and-bound methods in order to solve these single-objective problems. The model consists in identifying the sugar cane varieties with the lowest waste collection costs, while simultaneously it aims to obtain the greatest amount of produced energy by this waste. The hybrid methods are implemented in C++ programming language, and tests are performed to determine the efficient solutions in Pareto optimal sense of the multi-objective model and compare the performance of the hybrid method using the integrality test and without considering it. The mathematical results confirm that the proposed hybrid method for solving the aforementioned models presents good computational performance and reliable solutions.     

Mathematical programming time-based decomposition algorithm for discrete event simulation

Mathematical programming has been proposed in the literature as an alternative technique to simulating a special class of Discrete Event Systems. There are several benefits to using mathematical programs for simulation, such as the possibility of performing sensitivity analysis and the ease of better integrating the simulation and optimisation. However, applications are limited by the usually long computational times. This paper proposes a time-based decomposition algorithm that splits the mathematical programming model into a number of submodels that can be solved sequentially to make the mathematical programming approach viable for long running simulations. The number of required submodels is the solution of an optimisation problem that minimises the expected time for solving all of the submodels. In this way, the solution time becomes a linear function of the number of simulated entities.     

A Practical Approach for Robust and Flexible Vehicle Routing Using Metaheuristics and Monte Carlo Sampling

In this paper, we investigate how robust and flexible solutions of a number of stochastic variants of the capacitated vehicle routing problem can be obtained. To this end, we develop and discuss a method that combines a sampling based approach to estimate the robustness or flexibility of a solution with a metaheuristic optimization technique. This combination allows us to solve larger problems with more complex stochastic structures than traditional methods based on stochastic programming. It is also more flexible in the sense that adaptation of the approach to more complex problems can be easily done. We explicitly recognize the fact that the decision maker’s risk preference should be taken into account when choosing a robust or flexible solution and show how this can be done using our approach.     

The tolerance approach in multiobjective linear fractional programming

When solving a multiobjective programming problem by the weighted sum approach, weights represent the relative importance associated to the objectives. As these values are usually imprecise, it is important to analyze the sensitivity of the solution under possible deviations on the estimated values. In this sense, the tolerance approach provides a direct measure of how weights may vary simultaneously and independently from their estimated values while still retaining the same efficient solution. This paper provides an explicit expression to the maximum tolerance on weights in a multiobjective linear fractional programming problem when all the denominators are equal. An application is also presented to illustrate how the results may help the decision maker to choose a most satisfactory solution in a production problem.