gpt4 book ai didi

python - ArithmeticError 导致 cvxpy 出现 "Rank(A) < p or Rank([G; A]) < n"错误

转载 作者:太空宇宙 更新时间:2023-11-03 14:27:56 24 4
gpt4 key购买 nike

我正在尝试使用 cvxpy(因此使用 cvxopt)在具有 28 个节点和 37 条线路的相对简单的网络中对最佳功率流进行建模,但得到的是“Rank(A) < p or Rank([G; A] ) < n"错误。

(使用相同的代码,我可以为具有 4 个节点和 4 条线路的更简单的网络找到最佳解决方案。)

我已经非常仔细地检查了约束是否不一致(使用我在 Stack Overflow 上对这个问题的回答: Check constraints are ok in cvxpy with actual values )

这是回溯:

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

当我查看/Library/Python/2.7/site-packages/cvxopt/coneprog.py", line 687, 它实际上是在挽救一个算术错误

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

这有意义吗?为什么会出现算术错误,例如 OverflowErrorZeroDivisionErrorFloatingPointError,除非 cvxpy 扩展了 ArithmeticError 意味着问题描述不当,即有排名(A) < p 或排名([G; A]) < n"?

最佳答案

好吧,事实证明,删除一些冗余约束确实解决了问题(如 cvxopt 谷歌群组中的建议:https://groups.google.com/forum/?fromgroups=#!topic/cvxopt/qlu3CK1TdVQ。我不确定为什么会抛出 ArithmeticError

实际上,我删除的约束不仅依赖于其他约束,而且保证满足任何输入,因为我构建问题的方式。

关于python - ArithmeticError 导致 cvxpy 出现 "Rank(A) < p or Rank([G; A]) < n"错误,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/16978763/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com