与 gurobipy 的每个元素向量乘法

Question

此代码的最后一行失败：

A = np.random.rand(d, d)*5.0
b = np.random.rand(d)*10.0
m = gp.Model()
x = m.addMVar(d, name='x')
y = m.addMVar(d, name='y', vtype=gp.GRB.BINARY)
m.addConstr(A@x <= b)
m.addConstr(A@x >= y*b) // error occurs here on y*b

对于 gurobipy 版本 9.5.1 给出了这个错误：

File "src/gurobipy/mvar.pxi", line 108, in gurobipy.MVar.__mul__
File "src/gurobipy/mvar.pxi", line 184, in gurobipy.MVar._scalar_mult
TypeError: only size-1 arrays can be converted to Python scalars

我期待每个元素与星号的乘法。如何以与 gurobipy API?

兼容的方式声明每个元素的乘法

Answer 1

据我所知，Gurobi 的矩阵界面目前不支持此功能。也许值得提出一个功能请求 here。在此期间，您可以将 element-wise 产品 y*b 表示为常规 matrix-vector 产品：

m.addConstr(A @ x >= np.diag(b) @ y)

或者，您可以编写类似于 term-based 接口：

from gurobipy import quicksum as qsum

for i in range(d):
    m.addConstr(qsum(A[i, j] * x[j] for j in range(d)) >= y[i] * b[i])

与 gurobipy 的每个元素向量乘法

per-element vector multiplication with gurobipy

python

gurobi

elementwise-operations