Unanswered Questions
862 questions with no upvoted or accepted answers
15
votes
1
answer
309
views
Integrality gap in bilevel binary linear programming problem
I have a bilevel max-min optimization problem over binary variables, with constraints expressed using linear inequalities. The inner (minimization) problem is
$$
\begin{alignat}2
\min\limits_x&\...
13
votes
0
answers
240
views
Unpopular/Obscure solvers/algorithms that you used
There are many (numeric) optimizers which solve some subset of the set of multi-objective, black-or-grey-or-white box, robust, stochastic, mixed-integer, non-linear, non-smooth, manifold-constrained ...
12
votes
0
answers
177
views
which method has been used to automatically reformulate logical constraints in a standard MIP solver?
There are many of the formulations to linearize logical constraints by introducing new auxiliary binary or any appropriate variables and adding the corresponding constraints to the model. It can be ...
11
votes
0
answers
164
views
Characterizing the solution of a (non) linear maximization program
I have the following maximization program
\begin{align}
\max\limits_{\{q_i\}}&\quad\sum\limits_{i=1}^nq_i \\
\text{s.t.}&\quad\begin{cases} k_j \geq \sum\limits_{i=1}^n q_i^{1 \over \...
10
votes
0
answers
187
views
How to monitor/debug an optimization model in a production environment professionally?
I have seen this question How to avoid having your optimization models rusting?, which is kind of related, but I am more curious about the actual monitoring and the initial design of the software in ...
10
votes
0
answers
227
views
For the solution of complex gaming are we to leave behind approximation and fixed algorithms in favor of agents?
For an evolutionary extensive form game is the use of multi-agent systems a necessity or can more traditional approximate Nash equilibria and mean-field interactions suffice?
Few games can be ...
9
votes
0
answers
105
views
Adjusting duals heuristically
I am solving a scheduling problem-to find shifts and task schedules- using column generation. In essence, it is a set covering problem with additional constraints. The problem seems to be that the ...
9
votes
0
answers
230
views
Solving large-scale stochastic mixed integer program
What are some methods or algorithms for solving a large-scale stochastic mixed-integer optimization problem that runs on an hourly dataset for a year? Do we employ some kind of decomposition? (the ...
9
votes
0
answers
203
views
Ill-conditioned LP in Benders decomposition
I have implemented a Benders decomposition for a constrained network flow but the LP solver (Gurobi) warns me of the ill-conditioning of the subproblem dual LP. As you can see below, the coefficients ...
9
votes
0
answers
182
views
Modelling resource dependency in the assignment problem
The assignment problem is well-studied and has a nice polynomial time algorithm. I'm interested in an extension of this problem where all edges are in a certain group and taking multiple edges from ...
8
votes
1
answer
141
views
Is there any way to use Lazy constraints in Pyomo?
I would like to know if there is any way to implement a Lazy constraint in the Pyomo package. (As far as I know, the only way to implement such a constraint is by ...
8
votes
0
answers
169
views
Traveling Salesman Problem: determine k-exchange feasibility
Given a current solution $S$ and a $k$-exchange move $(v_1, .., v_{2k+1})$ with $v_1 = v_{2k+1}$, $v_i \neq v_j$, $(v_i, v_{i+1}) \in E(S)$ iff $i$ odd, i.e. we remove all edges $(v_i, v_{i+1})$ for $...
8
votes
0
answers
117
views
Multiple shipments with FILO order
I am trying to solve a problem statement with the help of jsprit.
There is a depot.
Identical items need to be delivered or picked up. An item can have two states: Bad or good.
Bad items need to be ...
8
votes
0
answers
112
views
For subset selection regression as a mixed integer program, how tightly should the bounding box be set?
When solving best subset regression as a mixed integer program, how do you decide how tightly to bound the range of values of the $X$ values? When the box is tight, the solver finds a solution ...
8
votes
0
answers
102
views
Examples for a kind of "set family hitting" problem
Given a ground set, say $[n]=\{1,2,\dots,n\}$, and a collection of subset families $\mathcal F_i\subseteq 2^{[n]}$, $i=1,2,\dots,m$, I want to select $m$ sets $B_i\in\mathcal F_i$ such that the ...
8
votes
0
answers
156
views
Is a base-stock policy optimal for a serial inventory system with upstream stockout costs?
In a serial inventory system without fixed costs, an echelon base-stock policy is known to be optimal if there is no stockout cost at any stage except the last stage. (This was proved for finite-...
7
votes
0
answers
107
views
MIP formulation for graph planarity test
In this question, it was asked wether a MIP formulation exists to test for a graph's planarity. The inputs are the graph's nodes and edges, and the output would be a certificate which guarantees that ...
7
votes
0
answers
118
views
What distuingishes a model-and-run framework solver from a programmatic solver?
I am looking for examples to understand what kind of benefits a model-and-run framework solver offers that a more programmatic solver doesn't do or doesn't do as well.
What distinguishes a model-and-...
7
votes
0
answers
118
views
List of Call for papers
Is there somewhere a list of ongoing Call for papers for special issues in OR journals?
I believe that such a list would not only create a good overview for publication opportunities but also inform ...
7
votes
0
answers
122
views
Estimate lagrangian multiplier based on instance characteristics
Assume we have a simple resource allocation problem, where all players have the same cost, but a different utility $a_s$. The resources assigned to a certain player must be between $L$ and $M$. ...
7
votes
0
answers
149
views
Is this a valid strong polynomial algorithm for deciding LP feasibility?
Let
$$A \cdot X + B \preceq 0$$ be a system of linear inequalities with $X \in \mathbb{R}^n$ $A\in \mathbb{R}^{m\times n}$ and $B \in \mathbb{R}^m$ where $m \geq n$. According to Farkas lemma, exactly ...
7
votes
0
answers
166
views
Modelers and modeling languages: how do you choose, and which features do you use?
There are many modeling languages and APIs around. One or more per solver, plus many that target multiple solvers: AMPL, GAMS, PuLP, JuMP, Pyomo...
Among all these possibilities, why did you pick a ...
7
votes
0
answers
151
views
Modeling traffic in a city
I am trying to model traffic in a city, $(i,j)$ represents a road in a city.
There are $H$ vehicles in a city they have some prescheduled set of destinations to visit, $A_{j,h}$ denotes arrival time ...
7
votes
0
answers
93
views
Building the Scheurman's Model II constraints for a multi period linear program
Scheurman's paper discusses Model I and model II Formulation to solve harvesting and scheduling problems. It is a specific implementation to solve multi period linear programs. Both models are also ...
7
votes
0
answers
80
views
Which Constraint solver is extensible to CSP algorithms?
I am novice to constraint programming and I need to implement the following algorithms.
Backtacking
Conflict-based backjumping
AC partial lookahead
MAC
Forward checking
Forward checking + dom
...
7
votes
0
answers
156
views
Why does PuLP call copy for addition and how can I avoid it?
Using a for loop to append terms to an expression seems to be much faster than summing a group of terms all at once. Constructing the expression using a for loop uses __iadd__, which does not include ...
7
votes
0
answers
81
views
Modelling a simple ordering problem to have balanced delivery days
Assuming that I should buy 50 items from 25 different vendors with pre-known delivery duration between 2-7 day for each, which day of a week should I place each order so that the delivery days be even ...
7
votes
0
answers
91
views
KKT conditions validation- one dual variable equating to two values
I have the following optimization problem:
\begin{alignat}2\min &\quad A(t)\cdot x(t)-B(t)\cdot y(t)+C(t)\cdot z(t)-D(t)\cdot k(t)\\\text{s.t.}&\quad z(t)+z_1(t)-y(t)-y_1(t)+x(t) = k(t);&...
7
votes
0
answers
176
views
Relationship between two minimization problems
Let $\mathbf{A}$ be a ${n\times J}$ matrix with $A_{ij}\geq0$ (and $A_{ij}>0$ for most $ij$, there cannot be any rows or columns that consist only of $0$s), $Q=\left\{\mathbf{q}\mid \mathbf{q}\in\...
7
votes
0
answers
885
views
Having trouble with objective function in Python: "GurobiError: Variable not in model". What else could I try?
I am trying to figure out how I can write this objective function into python using Gurobi. I have to minimize the sum of the product of three dictionary's values. The reason I am confused is that ...
7
votes
0
answers
96
views
Calculating robustness of layout plans
We have tried to design a manufacturing cell which will produce specific families of products. We figure out three layout plans for implementation. For practical reasons, we need to calculate the ...
6
votes
0
answers
142
views
What are some recommended Master's degree programs for individuals interested in specializing in the development of MI(L)P solvers?
There are tons of Masters programs across North America and Europe focused on OR/Industrial Engineering.
I am interested especially in solver technology and would be very grateful if someone could ...
6
votes
0
answers
213
views
Benders decomposition for a dense MILP
I am trying to solve a large MILP, but it seems like dense problems can be very difficult for moderns solvers. I tried to solve the problem described below considering only constraints (1) and (2) ...
6
votes
0
answers
94
views
Graph coloring problem while counting cliques
Let $G$ be a graph with a set of nodes $V$ and a set of edges $E$.
Let $G'$ be a graph with the same set of nodes $V$ but a second set of edges $E'$.
For a set of nodes $X\subset V$, we denote $f(X)$ ...
6
votes
0
answers
88
views
Robust Linear Optimization for avoiding diminishing returns
My engineering problem can be formulated as an LP as shown below
\begin{align}
\max_{\mathbf{x}}~~&\mathbf{a}^T\mathbf{x} \\
\mbox{s.t.}~~~&\mathbf{b}^T\mathbf{x} \leq B~~,~~\mathbf{1}^T\...
6
votes
0
answers
111
views
Dual instability, degeneracy, tailing off effect - Which are the causes and which are the effects?
Dual instability, degeneracy, and the tailing off effect are often mentioned together in papers. However, I cannot seem to find a clear explanation on which of these cause the other and vice versa? ...
6
votes
2
answers
171
views
Complexity of the ellipsoid method in general convex problems
The ellipsoid method is often mentioned in relation to the complexity of solving linear programs.
Is the method however polynomial in the general non-linear convex cases? Preferably I would like a ...
6
votes
0
answers
91
views
Issue of Dense columns in the formulation
I'm working on a price-selection model where we need to identify price point for each time-period (could be day/week). Objective of the model is to figure out optimal price-point for each time-point ...
6
votes
0
answers
226
views
Airline revenue management re-solving problem
I am considering a bid prices (shadow price of the capacity constraint) problem (from Chen, L. and Homem-de Mello, T. (2009)., page 14) where the acceptable classes for booking requests for ...
6
votes
0
answers
95
views
Sample Average Approximation vs. Numerical Integration
In the sense of the calculation of the expected value of objective functions, we have two choices to evaluate the value; 1. Sample Average Approximation (SAA):
$$
\frac{1}{N}\sum_{i=1}^N f(x,\xi^i).
$$...
6
votes
1
answer
248
views
Which software to use for Voronoi diagrams?
I have recently been looking at some facility location problems and how they change in the presence of the motorway (or highway, freeway) etc. Tried to visualize some observations with the Voronoi ...
6
votes
0
answers
127
views
Water quality component optimization
I have an optimization problem that I'm attempting to tackle. As you can see in the image below, there's a graph network through which water flows. I've drawn out the problem in the image to explain ...
6
votes
0
answers
135
views
Characterization for total dual integrality
A problem I study reduces to whether the polyhedron $P=\{\mathbf{x}\mid A\mathbf{x}=\mathbf{1}, \mathbf{x}\geq0\}$ is integral ($A$ is a matrix with coefficients in $\{0,1\}$). I know that the ...
6
votes
0
answers
53
views
What are useful plots/statistics/metrics when analyzing the solution sensitivity in multi-objective optimization?
Consider an optimization problem with $n>3$ objectives.
For handling this there exists often two approaches:
a) some weighting of the objectives,
b) fix an order of objectives and then optimize ...
6
votes
0
answers
106
views
Estimating multistop routing costs
In many OR problems, it is sometimes a good idea (or necessary) to relax routing constraints.
An example of this occurs in the classical facility location problem, where a warehouse can send out a ...
6
votes
0
answers
237
views
Provide basic solution to CLP
I'm using Pyomo to formulate an LP with approx 500,000 constraints and 200,000 decision variables. The LP is solved using CLP. Some instances fail to return even a feasible solution after many ...
6
votes
0
answers
145
views
Cases where RLT/SDP relaxation does not work well with standard quadratic optimization
(For people who don't know what RLT is): I am maximizing an indefinite quadratic function over a standard simplex, i.e., the standard quadratic optimization problem. A well-known approach is to relax ...
6
votes
0
answers
100
views
Proof that the leaving variable cannot be selected as the entering one in the next round
Using the Dantzig's pivoting rule, can it be proven that the leaving variable of one round cannot be selected as the entering variable in the next round?
6
votes
0
answers
87
views
How to communicate number of integer combinations to a user
I'm working on a nifty little feature for our next release, i.e., to print the number of possible integer combinations left during branch and bound.
This is really handy for the user because they ...
6
votes
0
answers
757
views
Solver for Go Programming Laguage
Has anyone used MIP solvers (open source or commercial) with Go language https://golang.org/? I am looking for a solver for simple linear MIPs for network flows and set cover types of problems.