提高約束添加在Gurobi（Python的接口）的性能

Question

我得到這個決策變量：

x={} 
for j in range(10): 
    for i in range(500000): 
     x[i,j] = m.addVar(vtype=GRB.BINARY, name="x%d%d" %(i,j)) 

所以我需要添加約束爲每個x [I，J]變量是這樣的：

for p in range(10): 
    for u in range(500000): 
     m.addConstr(x[u,p-1]<=x[u,p]) 

這花了我很多時間，多了12小時，然後在我的電腦上出現了缺少內存彈出窗口。 有人能helpme改善這個約束加法問題

Answer 1

一般備註：

• 它看起來相當昂貴，一般

具體注加500萬點的限制：

方法

• 您通過使用字典
  • 儘管有恆定訪問的複雜性是在浪費時間和空間，這些常數是大
  • 他們還浪費內存
• 在一個簡單的2維的情況是這樣的：堅持陣列！

有效性

• 你的索引是缺少第一個元素的邊界情況，所以索引休息！

嘗試此（更有效的方法;使用numpy的的陣列）：

import numpy as np 
from gurobipy import * 

N = 10 
M = 500000 

m = Model("Testmodel") 
x = np.empty((N, M), dtype=object) 
for i in range(N): 
    for j in range(M): 
     x[i,j] = m.addVar(vtype=GRB.BINARY, name="x%d%d" %(i,j)) 
m.update() 

for u in range(M): # i switched the loop-order 
    for p in range(1,N): # i'm handling the border-case 
     m.addConstr(x[p-1,u] <= x[p,u]) 

結果：

• 〜2分鐘
• 〜2.5GB存儲器（完整的程序，包括Gurobi的內部）

Answer 2

很可能，您的物理內存不足並使用虛擬（交換）內存。這不會導致您的計算機報告內存不足警告或錯誤。

我重寫代碼如下：

from gurobipy import * 

m = Model() 

x={} 
for j in range(10): 
    for i in range(500000): 
    x[i,j] = m.addVar(vtype=GRB.BINARY, name="x%d%d" %(i,j)) 
m.update() 

for p in range(10): 
    for u in range(500000): 
    try: 
     m.addConstr(x[u,p-1]<=x[u,p]) 
    except: 
     pass 
m.update() 

我測試此使用Gurobi優化6.5.2的計算機上用的Intel Xeon處理器E3-1240（3.40千兆赫）和32 GB的物理內存。它能夠在1分14秒內製定變量和限制條件。您可以使用列表來保存少量內存，但我相信Gurobi Var和Constr對象需要比Python字典或列表更多的內存。