MLinExpr and LinExpr incompatibilities, how to avoid that?

回答済み

Daniele Gioia

2022年03月25日 13:19 編集済み

• 
• 
• 

Hi, I am working with MVars and MLinExpr and I know that those environments are currently not compatible with  LinExpr and Vars. However, Gurobi automatically switches the type of one of my variables from MVar to Var. Specifically, in Python, 

These are three of my variables: 

import gurobipy as grb
zeta = model.addMVar(
    shape=1,
    lb=-grb.GRB.INFINITY,
    ub=grb.GRB.INFINITY,
    name="zeta",
    vtype=grb.GRB.CONTINUOUS
)
Z = model.addMVar(
    shape=100,
    vtype=grb.GRB.CONTINUOUS,
    name="Z"
)
Y = model.addMVar(
    shape=(50,100),
    vtype=grb.GRB.CONTINUOUS,
    name='Y'
)

Assuming two np arrays with the correct dimensions a and b,

 

One of my constraints is

model.addConstrs(
    Z[s] >= a @ X - b @ Y[:, s] - zeta
    for s in range(100)
)

And it raises the error:

GurobiError: Unsupported type (<class 'gurobipy.MLinExpr'>) for LinExpr addition argument

I initially thought that the problem was 'zeta', but, after some checks, the real problem is Z[s]. In particular, although Z is MVar, Z[s] is Var and generates LinExpr rather than MLinExpr. 

 

How can I solve this issue by maintaining the MVars? 

0

• 正式なコメント

  

  Simranjit Kaur

  • Gurobi Staff

  2025年11月16日 23:41

  This post is more than three years old. Some information may not be up to date. For current information, please check the Gurobi Documentation or Knowledge Base. If you need more help, please create a new post in the community forum, or try Gurobot, our chatbot interface offering instant, expert-level support.
• 

  Simon Bowly

  • Gurobi Staff

  2022年03月29日 01:01

  Hi Daniele,

  We currently have a few issues with interactions between Var/LinExpr and MVar/MLinExpr. There were some fixes added to the most recent release (version 9.5.1) which should address your case. Could you please upgrade to version 9.5.1 and additionally rewrite your constraints as below?

  model.addConstrs(
      a @ X - b @ Y[:, s] - zeta - Z[s] <= 0
      for s in range(100)
  )

  Note - I'm assuming here that X is a 1D MVar (similar to Z) and a and b are 1D numpy arrays with compatible shapes for multiplying with X and Y[:, s].

  Best regards,
  Simon

  1
• 

  Daniele Gioia

  2022年03月29日 15:16

  • 
  • 
  • 

  Thank you Simon, now it works!

  0

投稿コメントは受け付けていません。