A finite cutting plane method for facial disjunctive programs

This paper shows that any linear disjunctive program with a finite number of constraints can be transformed into an equivalent facial program. Based upon linear programming technique, a new, finite cutting plane method is presented for the facial programs.

---

Zusammenfassung Die Arbeit zeigt, daß jedes lineare disjunktive Optimierungsproblem mit endlich vielen Restriktionen in ein äquivalentes Fazetten-Problem transformiert werden kann. Auf der Grundlage von linearer Optimierungstechnik wird für das Fazetten-Problem ein neues, endliches Schnittebenenverfahren vorgestellt.

     

A logarithmic barrier cutting plane method for convex programming

The paper presents a logarithmic barrier cutting plane algorithm for convex (possibly non-smooth, semi-infinite) programming. Most cutting plane methods, like that of Kelley, and Cheney and Goldstein, solve a linear approximation (localization) of the problem and then generate an additional cut to remove the linear program's optimal point. Other methods, like the central cutting plane methods of Elzinga-Moore and Goffin-Vial, calculate a center of the linear approximation and then adjust the level of the objective, or separate the current center from the feasible set. In contrast to these existing techniques, we develop a method which does not solve the linear relaxations to optimality, but rather stays in the interior of the feasible set. The iterates follow the central path of a linear relaxation, until the current iterate either leaves the feasible set or is too close to the boundary. When this occurs, a new cut is generated and the algorithm iterates. We use the tools developed by den Hertog, Roos and Terlaky to analyze the effect of adding and deleting constraints in long-step logarithmic barrier methods for linear programming. Finally, implementation issues and computational results are presented. The test problems come from the class of numerically difficult convex geometric and semi-infinite programming problems.This work was completed under the support of a research grant of SHELL.On leave from the Eötvös University, Budapest, and partially supported by OTKA No. 2116.     

A combinatorial interior point method for network flow problems

For solving minimum cost flow problems, we develop a combinatorial interior point method based on a variant of the algorithm of Karmarkar, described in Gonzaga [3, 4]. Gonzaga proposes search directions generated by projecting certain directions onto the nullspace ofA. By the special combinatorial structure of networks any projection onto the nullspace ofA can be interpreted as a flow in the incremental graph ofG. In particular, to evaluate the new search direction, it is sufficient to choose a negative circuit subject to costs on the arcs depending on the current solution. That approach results in an O(mn ² L) algorithm wherem denotes the number of vertices,n denotes the number of arcs, andL denotes the total length of the input data.     

Interior proximal point algorithm for linear programs

An interior proximal point algorithm for finding a solution of a linear program is presented. The distinguishing feature of this algorithm is the addition of a quadratic proximal term to the linear objective function. This perturbation has allowed us to obtain solutions with better feasibility. Implementation of this algorithm shows that the algorithms. We also establish global convergence and local linear convergence of the algorithm.This research was supported by National Science Foundation Grants DCR-85-21228 and CCR-87-23091 and by Air Force Office of Scientific Research Grants AFOSR-86-0172 and AFOSR-89-0410. It was conducted while the author was a Graduate Student at the Computer Sciences Department, University of Wisconsin, Madison, Wisconsin.     

Interior point methods for optimal control of discrete time systems

We show that recently developed interior point methods for quadratic programming and linear complementarity problems can be put to use in solving discrete-time optimal control problems, with general pointwise constraints on states and controls. We describe interior point algorithms for a discrete-time linear-quadratic regulator problem with mixed state/control constraints and show how they can be efficiently-incorporated into an inexact sequential quadratic programming algorithm for nonlinear problems. The key to the efficiency of the interior-point method is the narrow-banded structure of the coefficient matrix which is factorized at each iteration.This research was supported by the Applied Mathematical Sciences Subprogram of the Office of Energy Research, US Department of Energy, under Contract W-31-109-Eng-38.     

A cutting plane method from analytic centers for stochastic programming

The stochastic linear programming problem with recourse has a dual block-angular structure. It can thus be handled by Benders' decomposition or by Kelley's method of cutting planes; equivalently the dual problem has a primal block-angular structure and can be handled by Dantzig-Wolfe decomposition—the two approaches are in fact identical by duality. Here we shall investigate the use of the method of cutting planes from analytic centers applied to similar formulations. The only significant difference form the aforementioned methods is that new cutting planes (or columns, by duality) will be generated not from the optimum of the linear programming relaxation, but from the analytic center of the set of localization.This research has been supported by the Fonds National de la Recherche Scientifique Suisse (grant # 12-26434.89), NSERC-Canada and FCAR-Quebec.Corresponding author.     

Reduced gradient method combined with augmented Lagrangian and barrier for the optimal power flow problem

A new approach for solving the optimal power flow (OPF) problem is established by combining the reduced gradient method and the augmented Lagrangian method with barriers and exploring specific characteristics of the relations between the variables of the OPF problem. Computer simulations on IEEE 14-bus and IEEE 30-bus test systems illustrate the method.     

Local radial point interpolation method for the fully developed magnetohydrodynamic flow

In this paper, a local radial point interpolation method (LRPIM) is presented to obtain the numerical solutions of the coupled equations in velocity and magnetic field for the fully developed magnetohydrodynamic (MHD) flow through a straight duct of rectangular section with arbitrary wall conductivity and orientation of applied magnetic field. Local weak forms are developed using weighted residual method locally for the governing equations of fully developed MHD flow. The shape functions from LRPIM possess the delta function property. Therefore, essential boundary conditions can be applied as easily as that in the finite-element method. The implementation procedure of LRPIM method is node based, and it doesn’t need any “mesh” or “element”. Computations have been carried out for different Hartmann numbers, wall conductivities and orientations of applied magnetic field.     

A boundary perturbation interior point homotopy method for solving fixed point problems

In this paper, a boundary perturbation interior point homotopy method is proposed to give a constructive proof of the general Brouwer fixed point theorem and thus solve fixed point problems in a class of nonconvex sets. Compared with the previous results, by using the newly proposed method, initial points can be chosen in the whole space of Rn, which may improve greatly the computational efficiency of reduced predictor-corrector algorithms resulted from that method. Some numerical examples are given to illustrate the results of this paper.     

Generalized fractional programming and cutting plane algorithms

In this paper, we introduce a variant of a cutting plane algorithm and show that this algorithm reduces to the well-known Dinkelbach-type procedure of Crouzeix, Ferland, and Schaible if the optimization problem is a generalized fractional program. By this observation, an easy geometrical interpretation of one of the most important algorithms in generalized fractional programming is obtained. Moreover, it is shown that the convergence of the Dinkelbach-type procedure is a direct consequence of the properties of this cutting plane method. Finally, a class of generalized fractional programs is considered where the standard positivity assumption on the denominators of the ratios of the objective function has to be imposed explicitly. It is also shown that, when using a Dinkelbach-type approach for this class of programs, the constraints ensuring the positivity on the denominators can be dropped.The authors like to thank the anonymous referees and Frank Plastria for their constructive remarks on an earlier version of this paper.This research was carried out at Erasmus University, Rotterdam, The Netherlands and was supported by JNICT, Lisboa, Portugal, under Contract BD/707/90-RM.     

An active-set strategy in an interior point method for linear programming

We will present a potential reduction method for linear programming where only the constraints with relatively small dual slacks—termed active constraints—will be taken into account to form the ellipsoid constraint at each iteration of the process. The algorithm converges to the optimal feasible solution in O( L) iterations with the same polynomial bound as in the full constraints case, wheren is the number of variables andL is the data length. If a small portion of the constraints is active near the optimal solution, the computational cost to find the next direction of movement in one iteration may be considerably reduced by the proposed strategy.This research was partially done in June 1990 while the author was visiting the Department of Mathematics, University of Pisa.     

A primal interior point method for the linear semidefinite programming problem

The linear semidefinite programming problem is examined. A primal interior point method is proposed to solve this problem. It extends the barrier-projection method used for linear programs. The basic properties of the proposed method are discussed, and its local convergence is proved.     

Interior point method for long-term generation scheduling of large-scale hydrothermal systems

This paper presents an interior point method for the long-term generation scheduling of large-scale hydrothermal systems. The problem is formulated as a nonlinear programming one due to the nonlinear representation of hydropower production and thermal fuel cost functions. Sparsity exploitation techniques and an heuristic procedure for computing the interior point method search directions have been developed. Numerical tests in case studies with systems of different dimensions and inflow scenarios have been carried out in order to evaluate the proposed method. Three systems were tested, with the largest being the Brazilian hydropower system with 74 hydro plants distributed in several cascades. Results show that the proposed method is an efficient and robust tool for solving the long-term generation scheduling problem.     

An implementation of a parallel primal-dual interior point method for block-structured linear programs

An implementation of the primal-dual predictor-corrector interior point method is specialized to solve block-structured linear programs with side constraints. The block structure of the constraint matrix is exploited via parallel computation. The side constraints require the Cholesky factorization of a dense matrix, where a method that exploits parallelism for the dense Cholesky factorization is used. For testing, multicommodity flow problems were used. The resulting implementation is 65%–90% efficient, depending on the problem instance. For a problem with K commodities, an approximate speedup for the interior point method of 0.8K is realized.     

On the optimal parameters of GMSSOR method for saddle point problems

For solving saddle point problems, SOR-type methods are investigated by many researchers in the literature. In this short note, we study the GMSSOR method for solving saddle point problems and obtain the optimal parameters which minimize the spectral (or pseudo-spectral) radii of the iteration matrices.     

Immersions of the projective plane with one triple point

We consider C^∞ generic immersions of the projective plane into the 3-sphere. Pinkall has shown that every immersion of the projective plane is homotopic through immersions to Boy's immersion, or its mirror. There is another lesser-known immersion of the projective plane with self-intersection set equivalent to Boy's but whose image is not homeomorphic to Boy's. We show that any C^∞ generic immersion of the projective plane whose self-intersection set in the 3-sphere is connected and has a single triple point is ambiently isotopic to precisely one of these two models, or their mirrors. We further show that any generic immersion of the projective plane with one triple point can be obtained by a sequence of toral and spherical surgical modifications of these models. Finally we present some simple applications of the theorem regarding discrete ambient automorphism groups; image-homology of immersions with one triple point; and almost tight ambient isotopy classes.     

Combined branch-and-bound and cutting plane methods for solving a class of nonlinear programming problems

We propose unified branch-and-bound and cutting plane algorithms for global minimization of a functionf(x, y) over a certain closed set. By formulating the problem in terms of two groups of variables and two groups of constraints we obtain new relaxation bounding and adaptive branching operations. The branching operation takes place in y-space only and uses the iteration points obtained through the bounding operation. The cutting is performed in parallel with the branch-and-bound procedure. The method can be applied implementably for a certain class of nonconvex programming problems.On leave from Institute of Mathematics, Hanoi, by a grant from Alexander-von-Humboldt-Stiftung.     

一种求解不等式约束优化问题的光滑化算法

利用光滑函数建立了不等式约束优化问题KT条件的一个扰动方程组,提出了一个新的内点型算法.该算法在有限步终止时当前迭代点即为优化问题的一个精确稳定点.在一定条件下算法具有全局收敛性,数值试验表明该算法是有效的.