Умножение вектора на элемент с помощью 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?

Анализ настроения постов в Twitter с помощью Python, Tweepy и Flair
Анализ настроения постов в Twitter с помощью Python, Tweepy и Flair
Анализ настроения текстовых сообщений может быть настолько сложным или простым, насколько вы его сделаете. Как и в любом ML-проекте, вы можете выбрать...
7 лайфхаков для начинающих Python-программистов
7 лайфхаков для начинающих Python-программистов
В этой статье мы расскажем о хитростях и советах по Python, которые должны быть известны разработчику Python.
Установка Apache Cassandra на Mac OS
Установка Apache Cassandra на Mac OS
Это краткое руководство по установке Apache Cassandra.
Сертификатная программа "Кванты Python": Бэктестер ансамблевых методов на основе ООП
Сертификатная программа "Кванты Python": Бэктестер ансамблевых методов на основе ООП
В одном из недавних постов я рассказал о том, как я использую навыки количественных исследований, которые я совершенствую в рамках программы TPQ...
Создание персонального файлового хранилища
Создание персонального файлового хранилища
Вы когда-нибудь хотели поделиться с кем-то файлом, но он содержал конфиденциальную информацию? Многие думают, что электронная почта безопасна, но это...
Создание приборной панели для анализа данных на GCP - часть I
Создание приборной панели для анализа данных на GCP - часть I
Недавно я столкнулся с интересной бизнес-задачей - визуализацией сбоев в цепочке поставок лекарств, которую могут просматривать врачи и...
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])

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