Gurobi模型中求极值矩阵系数的有效方法

Question

我正在寻找一种更有效的方法来识别Gurobi模型(!)中的极值矩阵系数。模型已经建立。对于“极端”，我指的是特别小或大的系数，这些系数对模型的数值性能有负面影响。以下代码非常慢，因为它查询所有系数：

    mdl.update()
    for c in mdl.getConstrs():
        c_name = c.ConstrName
        for v in mdl.getVars():
            v_name = v.VarName
            coefficient = mdl.getCoeff(c, v)
            if (coefficient != 0.0 and abs(coefficient) < 10e-4) or abs(coefficient) > 10e7:
                # Do something

我认为Model.getA()可能会有所帮助，但由于某种原因，调用该函数会引发错误。

Answer 1

如果只需要矩阵的最大值和/或最小系数值，只需分别查询MaxCoeff或MinCoeff属性。或者，如果您需要检查所有矩阵系数，则需要对所有系数进行迭代，但只需对非零系数进行迭代就可以更有效地完成。我相信这能满足你的需要：

for c in mdl.getConstrs():
  row = mdl.getRow(c)
  for i in range(row.size()):
    coefficient = row.getCoeff(i)
    if (coefficient != 0.0 and abs(coefficient) < 1e-4) or abs(coefficient) > 1e7:
      bigvar = row.getVar(i)
      # Do something