On relaxations of the max k-cut problem formulations

A tight continuous relaxation is a crucial factor in solving mixed integer formulations of many NP-hard combinatorial optimization problems. The (weighted) max k-cut problem is a fundamental combinatorial optimization problem with multiple notorious mixed integer optimization formulations. In this paper, we explore four existing mixed integer optimization formulations of the max k-cut problem. Specifically, we show that the continuous relaxation of a binary quadratic optimization formulation of the problem is: (i) stronger than the continuous relaxation of two mixed integer linear optimization formulations and (ii) at least as strong as the continuous relaxation of a mixed integer semidefinite optimization formulation. We also conduct a set of experiments on multiple sets of instances of the max k-cut problem using state-of-the-art solvers that empirically confirm the theoretical results in item (i). Furthermore, these numerical results illustrate the advances in the efficiency of global non-convex quadratic optimization solvers and more general mixed integer nonlinear optimization solvers. As a result, these solvers provide a promising option to solve combinatorial optimization problems. Our codes and data are available on GitHub.     

An inverse problem for an immobilized enzyme model

A method for estimating unknown kinetic parameters in a mathematical model for catalysis by an immobilized enzyme is studied. The model consists of a semilinear parabolic partial differential equation modeling the reaction‐diffusion process coupled with an ordinary differential equation for the rate transport. The well posedness of the model is proven; a PDE‐constrained optimization approach is applied to the stated inverse problem; and finally, some numerical simulations are presented.     

On the adjustment problem for linear programs

We propose a generalization of the inverse problem which we will call the adjustment problem. For an optimization problem with linear objective function and its restriction defined by a given subset of feasible solutions, the adjustment problem consists in finding the least costly perturbations of the original objective function coefficients, which guarantee that an optimal solution of the perturbed problem is also feasible for the considered restriction. We describe a method of solving the adjustment problem for continuous linear programming problems when variables in the restriction are required to be binary.     

Scope of stationary multi-objective evolutionary optimization: a case study on a hydro-thermal power dispatch problem

Many engineering design and developmental activities finally resort to an optimization task which must be solved to get an efficient and often an intelligent solution. Due to various complexities involved with objective functions, constraints, and decision variables, optimization problems are often not adequately suitable to be solved using classical point-by-point methodologies. Evolutionary optimization procedures use a population of solutions and stochastic update operators in an iteration in a manner so as to constitute a flexible search procedure thereby demonstrating promise to such difficult and practical problem-solving tasks. In this paper, we illustrate the power of evolutionary optimization algorithms in handling different kinds of optimization tasks on a hydro-thermal power dispatch optimization problem: (i) dealing with non-linear, non-differentiable objectives and constraints, (ii) dealing with more than one objectives and constraints, (iii) dealing with uncertainties in decision variables and other problem parameters, and (iv) dealing with a large number (more than 1,000) variables. The results on the static power dispatch optimization problem are compared with that reported in an existing simulated annealing based optimization procedure on a 24-variable version of the problem and new solutions are found to dominate the solutions of the existing study. Importantly, solutions found by our approach are found to satisfy theoretical Kuhn–Tucker optimality conditions by using the subdifferentials to handle non-differentiable objectives. This systematic and detail study demonstrates that evolutionary optimization procedures are not only flexible and scalable to large-scale optimization problems, but are also potentially efficient in finding theoretical optimal solutions for difficult real-world optimization problems. Kalyanmoy Deb, Deva Raj Chair Professor. Currently a Finland Distinguished Professor, Department of Business Technology, Helsinki School of Economics, 00101 Helsinki, Finland.     

Algorithm robust for the bicriteria discrete optimization problem

We apply Algorithm Robust to various problems in multiple objective discrete optimization. Algorithm Robust is a general procedure that is designed to solve bicriteria optimization problems. The algorithm performs a weight space search in which the weights are utilized in min-max type subproblems. In this paper, we experiment with Algorithm Robust on the bicriteria knapsack problem, the bicriteria assignment problem, and the bicriteria minimum cost network flow problem. We look at a heuristic variation that is based on controlling the weight space search and has an indirect control on the sample of efficient solutions generated. We then study another heuristic variation which generates samples of the efficient set with quality guarantees. We report results of computational experiments.     

A new equivalent single-level problem for bilevel problems

In this article, we give necessary optimality conditions for a bilevel optimization problem (P). An intermediate single-level problem (Q), which is equivalent to the bilevel optimization problem (P), has been introduced.     

On the Karush-Kuhn-Tucker reformulation of the bilevel optimization problem

This paper is mainly concerned with the classical KKT reformulation and the primal KKT reformulation (also known as an optimization problem with generalized equation constraint (OPEC)) of the optimistic bilevel optimization problem. A generalization of the MFCQ to an optimization problem with operator constraint is applied to each of these reformulations, hence leading to new constraint qualifications (CQs) for the bilevel optimization problem. M- and S-type stationarity conditions tailored for the problem are derived as well. Considering the close link between the aforementioned reformulations, similarities and relationships between the corresponding CQs and optimality conditions are highlighted. In this paper, a concept of partial calmness known for the optimal value reformulation is also introduced for the primal KKT reformulation and used to recover the M-stationarity conditions.     

The inverse optimal value problem

This paper considers the following inverse optimization problem: given a linear program, a desired optimal objective value, and a set of feasible cost vectors, determine a cost vector such that the corresponding optimal objective value of the linear program is closest to the desired value. The above problem, referred here as the inverse optimal value problem, is significantly different from standard inverse optimization problems that involve determining a cost vector for a linear program such that a pre-specified solution vector is optimal. In this paper, we show that the inverse optimal value problem is NP-hard in general. We identify conditions under which the problem reduces to a concave maximization or a concave minimization problem. We provide sufficient conditions under which the associated concave minimization problem and, correspondingly, the inverse optimal value problem is polynomially solvable. For the case when the set of feasible cost vectors is polyhedral, we describe an algorithm for the inverse optimal value problem based on solving linear and bilinear programming problems. Some preliminary computational experience is reported.Mathematics Subject Classification (1999):49N45, 90C05, 90C25, 90C26, 90C31, 90C60Acknowledgement This research has been supported in part by the National Science Foundation under CAREER Award DMII-0133943. The authors thank two anonymous reviewers for valuable comments.     

Inverse problem for a class of one-dimensional wave equations with piecewise constant coefficients

We consider the problem of reconstructing the piecewise constant coefficient of a one-dimensional wave equation on the halfline from the knowledge of the displacement on the boundary caused by an impulse at time zero. This problem is formulated as a nonlinear optimization problem. The objective function of this optimization problem has several special features that have been exploited in building an ad hoc optimization method. The optimization method is based on the solution of a nonlinear system of equations by an algorithm consisting of the evaluation of the unknowns one by one.The research of the third author has been made possible through the support and sponsorship of the Italian Government through the Ministero Pubblica Istruzione under Contract M.P.I. 60% 1987 at the Università di Roma—La Sapienza.     

Approximate solution of the trust region problem by minimization over two-dimensional subspaces

The trust region problem, minimization of a quadratic function subject to a spherical trust region constraint, occurs in many optimization algorithms. In a previous paper, the authors introduced an inexpensive approximate solution technique for this problem that involves the solution of a two-dimensional trust region problem. They showed that using this approximation in an unconstrained optimization algorithm leads to the same theoretical global and local convergence properties as are obtained using the exact solution to the trust region problem. This paper reports computational results showing that the two-dimensional minimization approach gives nearly optimal reductions in then-dimension quadratic model over a wide range of test cases. We also show that there is very little difference, in efficiency and reliability, between using the approximate or exact trust region step in solving standard test problems for unconstrained optimization. These results may encourage the application of similar approximate trust region techniques in other contexts.Research supported by ARO contract DAAG 29-84-K-0140, NSF grant DCR-8403483, and NSF cooperative agreement DCR-8420944.     

Optimization reformulations of the generalized Nash equilibrium problem using Nikaido-Isoda-type functions

We consider the generalized Nash equilibrium problem which, in contrast to the standard Nash equilibrium problem, allows joint constraints of all players involved in the game. Using a regularized Nikaido-Isoda-function, we then present three optimization problems related to the generalized Nash equilibrium problem. The first optimization problem is a complete reformulation of the generalized Nash game in the sense that the global minima are precisely the solutions of the game. However, this reformulation is nonsmooth. We then modify this approach and obtain a smooth constrained optimization problem whose global minima correspond to so-called normalized Nash equilibria. The third approach uses the difference of two regularized Nikaido-Isoda-functions in order to get a smooth unconstrained optimization problem whose global minima are, once again, precisely the normalized Nash equilibria. Conditions for stationary points to be global minima of the two smooth optimization problems are also given. Some numerical results illustrate the behaviour of our approaches.     

A global optimization method for the molecular replacement problem in X-ray crystallography

The primary technique for determining the three-dimensional structure of a protein molecule is X-ray crystallography, from which the molecular replacement (MR) problem often arises as a critical step. The MR problem is a global optimization problem to locate an optimal position of a model protein so that at this position the model will produce calculated intensities closest to those observed from an X-ray crystallography experiment involving a protein with unknown but similar atomic structure. Improving the applicability and robustness of MR methods is an important research topic because commonly used traditional MR methods, though often successful, have their limitations in solving difficult problems.We introduce a new global optimization strategy that combines a coarse-grid search, using a surrogate function, with extensive multi-start local optimization. A new MR code, called SOMoRe, based on this strategy is developed and tested on four realistic problems, including two difficult problems that traditional MR codes failed to solve directly. SOMoRe was able to solve each test problem without any complication, and SOMoRe solved an MR problem using a less complete model than the models required by three other programs. These results indicate that the new method is promising and should enhance the applicability and robustness of the MR methodology.     

多损失WCVaR模型的等价性定理

研究了多概率分布簇下的多损失下的WCVaR(Multi Worst Conditional Value-at-Risk)模型等价性定理, 根据概率分布簇的VaR测度值, 定义了多损失下的WCVaR风险测度值和对应的多目标优化模型(MWCVaR), 证明了多目标优化模型(MWCVaR)等价另一个多目标优化模型求解. 对于有限分布簇情形, 在一定条件下, 证明了用有限个分布簇就可以近似计算多损失(MWCVaR)优化模型.     

A linear optimization problem to derive relative weights using an interval judgement matrix

In this paper, we propose a new Decision Making model, enabling to assess a finite number of alternatives according to a set of bounds on the preference ratios for the pairwise comparisons between alternatives, that is, an “interval judgement matrix”. In the case in which these bounds cannot be achieved by any assessment vector, we analyze the problem of determining of an efficient or Pareto-optimal solution from a multi-objective optimization problem. This multi-objective formulation seeks for assessment vectors that are near to simultaneously fulfil all the bound requirements imposed by the interval judgement matrix. Our new model introduces a linear optimization problem in order to define a consistency index for the interval matrix. By solving this optimization problem it can be associated a weakly efficient assessment vector to the consistency index in those cases in which the bound requirements are infeasible. Otherwise, this assessment vector fulfils all the bound requirements and has geometrical properties that make it appropriate as a representative assessment vector of the set of feasible weights.     

Bilevel aggregator-prosumers' optimization problem in real-time: A convex optimization approach

This paper proposes a Real-Time Market (RTM) platform for an aggregator and its corresponding prosumers to participate in the electricity wholesale market. The proposed energy market platform is modeled as a bilevel optimization problem where the aggregator and the prosumers are considered as self-interested agents. We present a convex optimization problem which can capture a subset of the set of global optima of the bilevel problem as its optimal solution.     

On the calculation of a membership function for the solution of a fuzzy linear optimization problem

In the present paper the fuzzy linear optimization problem (with fuzzy coefficients in the objective function) is considered. Recent concepts of fuzzy solution to the fuzzy optimization problem based on the level-cut and the set of Pareto optimal solutions of a multiobjective optimization problem are applied. Chanas and Kuchta suggested one approach to determine the membership function values of fuzzy optimal solutions of the fuzzy optimization problem, which is based on calculating the sum of lengths of certain intervals. The purpose of this paper is to determine a method for realizing this idea. We derive explicit formulas for the bounds of these intervals in the case of triangular fuzzy numbers and show that only one interval needs to be considered.     

Tikhonov regularization for a general nonlinear constrained optimization problem

The goal of this study is to analyze the Tikhonov regularization method as applied to a general nonlinear optimization problem that has been previously reduced to an unconstrained optimization problem. The stability properties of the method are examined, and its convergence is proved. The text was submitted by the author in English.     

A fuzzy goal programming approach to multi-objective optimization problem with priorities

In goal programming problem, the general equilibrium and optimization are often two conflicting factors. This paper proposes a generalized varying-domain optimization method for fuzzy goal programming (FGP) incorporating multiple priorities. According to the three possible styles of the objective function, the varying-domain optimization method and its generalization are proposed. This method can generate the results consistent with the decision-maker (DM)’s expectation, that the goal with higher priority may have higher level of satisfaction. Using this new method, it is a simple process to balance between the equilibrium and optimization, and the result is the consequence of a synthetic decision between them. In contrast to the previous method, the proposed method can make that the higher priority achieving the higher satisfactory degree. To get the global solution of the nonlinear nonconvex programming problem resulting from the original problem and the varying-domain optimization method, the co-evolutionary genetic algorithms (GAs), called GENOCOPIII, is used instead of the SQP method. In this way the DM can get the optimum of the optimization problem. We demonstrate the power of this proposed method by illustrative examples.     

Generalized-Hukuhara penalty method for optimization problem with interval-valued functions and its application in interval-valued portfolio optimization problems

In this study, a gH-penalty method is developed to obtain efficient solutions to constrained optimization problems with interval-valued functions. The algorithmic implementation of the proposed method is illustrated. In order to develop the gH-penalty method, an interval-valued penalty function is defined and the characterization of efficient solutions of a CIOP is done. As an application of the proposed method, a portfolio optimization problem with interval-valued return is solved.