Умножение вектора на элемент с помощью gurobipy

Последняя строка этого кода не работает:

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?

Почему в Python есть оператор "pass"?
Почему в Python есть оператор "pass"?
Оператор pass в Python - это простая концепция, которую могут быстро освоить даже новички без опыта программирования.
Некоторые методы, о которых вы не знали, что они существуют в Python
Некоторые методы, о которых вы не знали, что они существуют в Python
Python - самый известный и самый простой в изучении язык в наши дни. Имея широкий спектр применения в области машинного обучения, Data Science,...
Основы Python Часть I
Основы Python Часть I
Вы когда-нибудь задумывались, почему в программах на Python вы видите приведенный ниже код?
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
Алиса и Боб имеют неориентированный граф из n узлов и трех типов ребер:
Оптимизация кода с помощью тернарного оператора Python
Оптимизация кода с помощью тернарного оператора Python
И последнее, что мы хотели бы показать вам, прежде чем двигаться дальше, это
Советы по эффективной веб-разработке с помощью Python
Советы по эффективной веб-разработке с помощью Python
Как веб-разработчик, Python может стать мощным инструментом для создания эффективных и масштабируемых веб-приложений.
1
0
21
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Насколько я знаю, в настоящее время это не поддерживается матричным интерфейсом Gurobi. Возможно, стоит сделать запрос функции здесь. Тем временем вы можете выразите поэлементное произведение y*b в виде обычного матрично-векторного произведения:

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

В качестве альтернативы вы можете написать его аналогично интерфейсу на основе терминов:

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])

Другие вопросы по теме