Python “算术错误导致”；职级（甲）<；p或秩（[G；A]）<；"；cvxpy中的错误

我试图使用cvxpy（因此cvxopt）来模拟一个相对简单的网络中的最优潮流，该网络有28个节点和37条线路，但得到了“秩（a） （使用相同的代码，我可以找到4个节点和4条线路的更简单网络的最佳解决方案。）

我非常仔细地检查了约束是否不一致（使用我对堆栈溢出问题的回答： )

以下是回溯：

File "main.py", line 201, in test_simple_optimisation
    p, bids, offers = optimize_flow(bids, offers, lines, nodes, injections_from_schedule, shift_factors, admittance.T)
  File "main.py", line 143, in optimize_flow
    p.solve()
  File "build/bdist.macosx-10.8-intel/egg/cvxpy/programs.py", line 169, in solve
obj,valid = solve_prog(new_p,quiet)
  File "build/bdist.macosx-10.8-intel/egg/cvxpy/procedures/solve_prog.py", line 44, in solve_prog
    sol = call_solver(p_expanded,quiet)
  File "build/bdist.macosx-10.8-intel/egg/cvxpy/procedures/call_solver.py", line 78, in call_solver
    r =  solvers.conelp(c,G,h,dims,A,b)
  File "/Library/Python/2.7/site-packages/cvxopt/coneprog.py", line 687, in conelp
    raise ValueError("Rank(A) < p or Rank([G; A]) < n")
ValueError: Rank(A) < p or Rank([G; A]) < n

test\u simple\u优化中的文件“main.py”，第201行
p、 投标、报价=优化流程（投标、报价、线路、节点、计划中的注入、转换系数、准入.T）
文件“main.py”，第143行，在优化流程中
p、 解决（）
文件“build/bdist.macosx-10.8-intel/egg/cvxpy/programs.py”，第169行，在solve中
obj，有效=解决程序（新程序，安静）
文件“build/bdist.macosx-10.8-intel/egg/cvxpy/procedures/solve_prog.py”，第44行，在solve_prog中
sol=调用解算器（p_展开，安静）
文件“build/bdist.macosx-10.8-intel/egg/cvxpy/procedures/call_solver.py”，第78行，在call_solver中
r=解算器。conelp（c、G、h、dims、A、b）
conelp中的文件“/Library/Python/2.7/site packages/cvxopt/coneprog.py”，第687行
升值错误（“等级（A）

当我查看/Library/Python/2.7/site packages/cvxopt/coneprog.py“第687行时，它实际上是在挽救一个算术错误

 685         try: f = kktsolver(W)
 686         except ArithmeticError:..
 687             raise ValueError("Rank(A) < p or Rank([G; A]) < n")

685尝试：f=kktsolver（W）
686除算术错误外：。。
687升值错误（“等级（A）

---

这有意义吗？为什么算术错误，例如

溢出错误

，

零分错误

，

浮点错误

，除非cvxpy扩展了算术错误，否则这意味着问题描述得不好，即秩（A）那么，事实证明，删除一些冗余约束确实解决了问题（正如在cvxopt google组上所建议的：。我不确定抛出

算术错误的原因

实际上，我移除的约束不仅依赖于其他约束，而且保证满足任何输入，因为我是如何构造问题的