python如何调用gurobi,在Python中使用多进程构建Gurobi表达式

首先，我正在使用jupyter notebook、python3.5版本和gurobi7.0.2及其python接口，都在redhat上。在

这就是我问题的背景：我想解决一个二次问题，它有大量的变量。建立目标函数需要1-2个小时以上。在

我考虑过使用NumpyGPU加速，但是表达式有点棘手，所以，这不是一个解决方案。在

因此，我尝试使用多个线程来构建目标函数。但是，我得到了一个错误，我不知道如何处理它。在

我简化了我的代码，因此它可以更具可读性(错误仍然相同)。在from gurobipy import *

import multiprocessing as mp

import queue

mp.set_start_method('fork')

def function(obj,q):

print('We enter')

obj = x*x + x*y + y*y + y*z + z*z + 2*x

q.put(obj)

print('We end')

m = Model("qp")

obj = QuadExpr()

x = m.addVar(ub=1.0, name="x")

y = m.addVar(ub=1.0, name="y")

z = m.addVar(ub=1.0, name="z")

q = mp.Queue()

if __name__ == '__main__':

for k in range (4):

p = mp.Process(target=function, args=(obj,q,))

p.start()

obj+=q.get()

p.join()

m.setObjective(obj)

# Add constraint: x + 2 y + 3 z <= 4

m.addConstr(x + 2 * y + 3 * z >= 4, "c0")

# Add constraint: x + y >= 1

m.addConstr(x + y >= 1, "c1")

m.optimize()

for v in m.getVars():

print('%s %g' % (v.varName, v.x))

print('Obj: %g' % obj.getValue())

输出：

^{pr2}$

你知道我怎么解决这个问题吗？有没有其他方法可以“返回”这个过程的结果？。我希望避免在磁盘中写入，因为它很慢(可能最后一个资源：S)。在

提前感谢：)。在

p.D.我代码中最慢的部分是这一部分，我试图将其分为几个过程：# var is an array of GRB.BINARY

# D=edge_costs

def penalty_function(var,obj,D):

num_nodes = len(var)

for i,fil in enumerate(D):

for j,val in enumerate(fil):

# -x*D*x

if val!=0:

obj+=var[i]*var[j]*val

# -x(i)x(j)*min(Ds)

if(j>i):

for k in range(num_nodes):

if(not j==i):

minval=min(D[j][k],D[i][k])

if (minval!=0 ):

obj+=var[i]*var[j]*minval

return obj