iSoron
commented
1 year ago Thanks for the bug report, @mzy2240. I am currently unable to reproduce it (on Linux). The error I get when solving your provided instance file is the one below. Does this instance have any disconnected buses or islands?
[ Info: Building model...
[ Info: Computing injection shift factors...
ERROR: LinearAlgebra.SingularException(2)
Stacktrace:
[1] checknonsingular
@ ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/factorization.jl:19 [inlined]
[2] checknonsingular
@ ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/factorization.jl:22 [inlined]
[3] #lu!#170
@ ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/lu.jl:82 [inlined]
[4] lu!
@ ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/lu.jl:80 [inlined]
[5] #lu#176
@ ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/lu.jl:299 [inlined]
[6] lu (repeats 2 times)
@ ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/lu.jl:298 [inlined]
[7] inv(A::Matrix{Float64})
@ LinearAlgebra ~/Software/julia-1.9.0/share/julia/stdlib/v1.9/LinearAlgebra/src/dense.jl:917
[8] _injection_shift_factors(; buses::Vector{UnitCommitment.Bus}, lines::Vector{UnitCommitment.TransmissionLine})
@ UnitCommitment ~/.julia/packages/UnitCommitment/G53EO/src/model/formulations/base/sensitivity.jl:23
mzy2240
Description
Failed to solve a ED (single snapshot, no contingency) and optimization stuck in an infinite loop.
Steps to Reproduce
Rename the attached txt file to json, then load just like other instance json. You can use any MIP solver to solve (HiGHS, COPT, Gurobi, etc.) example.txt
(Update) Use 1/x instead of b as the susceptance. You can do that by modifying the
_susceptance_matrixfunction:System Information