ArithmeticError вызывает ошибку "Rank(A) <p или Rank ([G; A]) <n" в cvxpy

Question

Я пытаюсь использовать cvxpy (и, следовательно, cvxopt) для моделирования оптимального потока мощности в относительно простой сети с 28 узлами и 37 линиями, но получаю "Rank(A)

(Используя тот же код, я могу найти оптимальное решение для гораздо более простой сети с 4 узлами и 4 линиями.)

Я очень тщательно проверил, что ограничения не являются противоречивыми (используя мой ответ на этот вопрос по переполнению стека: проверьте, что ограничения в cvxpy действительны с фактическими значениями)

Вот обратная связь:

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", строка 687, это на самом деле спасает арифметическую ошибку

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

Имеет ли это смысл? Почему бы арифметическая ошибка, то есть что-то вроде OverflowError, ZeroDivisionError, FloatingPointError, если cvxpy не расширил ArithmeticError, означает, что проблема плохо описана, т.е. имеет ранг (A)

Answer 1

Что ж, получается, что удаление некоторых избыточных ограничений действительно решило проблему (как это было предложено в группах Google cvxopt: https://groups.google.com/forum/?fromgroups=. Я не являюсь конечно, почему ArithmeticError был брошен.

На самом деле, снятые мною ограничения не просто зависели от других ограничений, но гарантированно были выполнены для любого ввода из-за того, как я построил проблему.