Skip to content

Add tests for semi continuous and semi integer variables #1033

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 8 commits into from
Mar 5, 2020
Merged

Add tests for semi continuous and semi integer variables #1033

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
9a2cfee
Add new tests for semi continuous and semi integer variables
rafabench Feb 21, 2020
701caa6
Bug fix
rafabench Feb 21, 2020
ef3fbc2
Fix typos and remove debug line
Feb 22, 2020
a777f91
Add MockOptimizer tests for semiconttest and semiinttest
Feb 22, 2020
d2f2305
Fix indentation and remove empty lines
rafabench Feb 28, 2020
7032a87
Use an auxiliary _semitest function
rafabench Mar 3, 2020
2d943f2
Fix variable type in _semitest
rafabench Mar 3, 2020
6c2d48a
Update model description
rafabench Mar 3, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
144 changes: 144 additions & 0 deletions src/Test/intlinear.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -681,6 +681,148 @@ function indicator4_test(model::MOI.ModelLike, config::TestConfig)
end
end

function _semitest(model::MOI.ModelLike, config::TestConfig{T}, int::Bool) where T
atol = config.atol
rtol = config.rtol

@test MOIU.supports_default_copy_to(model, #=copy_names=# false)
@test MOI.supports(model, MOI.ObjectiveFunction{MOI.ScalarAffineFunction{T}}())
@test MOI.supports(model, MOI.ObjectiveSense())
if !int
@test MOI.supports_constraint(model, MOI.SingleVariable, MOI.Semicontinuous{T})
else
@test MOI.supports_constraint(model, MOI.SingleVariable, MOI.Semiinteger{T})
end

# 2 variables
#
# min x
# st x >= y
# if !int
# x ∈ {0.0} U [2.0,3.0]
# if int
# x ∈ {0.0} U {2.0} U {3.0}
# y = 0.0

MOI.empty!(model)
@test MOI.is_empty(model)

v = MOI.add_variables(model, 2)
@test MOI.get(model, MOI.NumberOfVariables()) == 2

if !int
vc1 = MOI.add_constraint(model, MOI.SingleVariable(v[1]), MOI.Semicontinuous(T(2), T(3)))
@test MOI.get(model, MOI.NumberOfConstraints{MOI.SingleVariable,MOI.Semicontinuous{T}}()) == 1
else
vc1 = MOI.add_constraint(model, MOI.SingleVariable(v[1]), MOI.Semiinteger(T(2), T(3)))
@test MOI.get(model, MOI.NumberOfConstraints{MOI.SingleVariable,MOI.Semiinteger{T}}()) == 1
end

vc2 = MOI.add_constraint(model, MOI.SingleVariable(v[2]), MOI.EqualTo(zero(T)))
@test MOI.get(model, MOI.NumberOfConstraints{MOI.SingleVariable,MOI.EqualTo{T}}()) == 1

cf = MOI.ScalarAffineFunction{T}(MOI.ScalarAffineTerm{T}.([one(T), -one(T)], v), zero(T))
c = MOI.add_constraint(model, cf, MOI.GreaterThan(zero(T)))
@test MOI.get(model, MOI.NumberOfConstraints{MOI.ScalarAffineFunction{T},MOI.GreaterThan{T}}()) == 1

objf = MOI.ScalarAffineFunction(MOI.ScalarAffineTerm.([1.0, 0.0], v), 0.0)
MOI.set(model, MOI.ObjectiveFunction{MOI.ScalarAffineFunction{Float64}}(), objf)
MOI.set(model, MOI.ObjectiveSense(), MOI.MIN_SENSE)

@test MOI.get(model, MOI.ObjectiveSense()) == MOI.MIN_SENSE

if config.solve
@test MOI.get(model, MOI.TerminationStatus()) == MOI.OPTIMIZE_NOT_CALLED

MOI.optimize!(model)

@test MOI.get(model, MOI.TerminationStatus()) == config.optimal_status
@test MOI.get(model, MOI.ResultCount()) >= 1
@test MOI.get(model, MOI.PrimalStatus()) in [ MOI.FEASIBLE_POINT, MOI.NEARLY_FEASIBLE_POINT ]
@test MOI.get(model, MOI.ObjectiveValue()) ≈ 0.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.VariablePrimal(), v) ≈ [0, 0] atol=atol rtol=rtol
@test MOI.get(model, MOI.ConstraintPrimal(), c) ≈ 0.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.ObjectiveBound()) <= 0.0
end

# Change y fixed value

MOI.set(model, MOI.ConstraintSet(), vc2, MOI.EqualTo(one(T)))

if config.solve
MOI.optimize!(model)

@test MOI.get(model, MOI.TerminationStatus()) == config.optimal_status
@test MOI.get(model, MOI.ResultCount()) >= 1
@test MOI.get(model, MOI.PrimalStatus()) == MOI.FEASIBLE_POINT
@test MOI.get(model, MOI.ObjectiveValue()) ≈ 2.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.VariablePrimal(), v) ≈ [2.0, 1.0] atol=atol rtol=rtol
@test MOI.get(model, MOI.ConstraintPrimal(), c) ≈ 1.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.ObjectiveBound()) <= 2.0
end

MOI.set(model, MOI.ConstraintSet(), vc2, MOI.EqualTo(T(2)))

if config.solve
MOI.optimize!(model)

@test MOI.get(model, MOI.TerminationStatus()) == config.optimal_status
@test MOI.get(model, MOI.ResultCount()) >= 1
@test MOI.get(model, MOI.PrimalStatus()) == MOI.FEASIBLE_POINT
@test MOI.get(model, MOI.ObjectiveValue()) ≈ 2.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.VariablePrimal(), v) ≈ [2.0, 2.0] atol=atol rtol=rtol
@test MOI.get(model, MOI.ConstraintPrimal(), c) ≈ 0.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.ObjectiveBound()) <= 2.0
end

MOI.set(model, MOI.ConstraintSet(), vc2, MOI.EqualTo(T(5//2)))

if config.solve
MOI.optimize!(model)

@test MOI.get(model, MOI.TerminationStatus()) == config.optimal_status
@test MOI.get(model, MOI.ResultCount()) >= 1
@test MOI.get(model, MOI.PrimalStatus()) == MOI.FEASIBLE_POINT
if !int
@test MOI.get(model, MOI.ObjectiveValue()) ≈ 2.5 atol=atol rtol=rtol
@test MOI.get(model, MOI.VariablePrimal(), v) ≈ [2.5, 2.5] atol=atol rtol=rtol
@test MOI.get(model, MOI.ConstraintPrimal(), c) ≈ 0.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.ObjectiveBound()) <= 2.5
else
@test MOI.get(model, MOI.ObjectiveValue()) ≈ 3.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.VariablePrimal(), v) ≈ [3.0, 2.5] atol=atol rtol=rtol
@test MOI.get(model, MOI.ConstraintPrimal(), c) ≈ 0.5 atol=atol rtol=rtol
@test MOI.get(model, MOI.ObjectiveBound()) <= 3.0
end
end

MOI.set(model, MOI.ConstraintSet(), vc2, MOI.EqualTo(T(3)))

if config.solve
MOI.optimize!(model)

@test MOI.get(model, MOI.TerminationStatus()) == config.optimal_status
@test MOI.get(model, MOI.ResultCount()) >= 1
@test MOI.get(model, MOI.PrimalStatus()) == MOI.FEASIBLE_POINT
@test MOI.get(model, MOI.ObjectiveValue()) ≈ 3.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.VariablePrimal(), v) ≈ [3.0, 3.0] atol=atol rtol=rtol
@test MOI.get(model, MOI.ConstraintPrimal(), c) ≈ 0.0 atol=atol rtol=rtol
@test MOI.get(model, MOI.ObjectiveBound()) <= 3.0
end

MOI.set(model, MOI.ConstraintSet(), vc2, MOI.EqualTo(T(4)))

if config.solve
MOI.optimize!(model)

@test MOI.get(model, MOI.TerminationStatus()) == MOI.INFEASIBLE ||
MOI.get(model, MOI.TerminationStatus()) == MOI.INFEASIBLE_OR_UNBOUNDED
end
end

semiconttest(model::MOI.ModelLike, config::TestConfig) = _semitest(model, config, false)
semiinttest(model::MOI.ModelLike, config::TestConfig) = _semitest(model, config, true)

const intlineartests = Dict("knapsack" => knapsacktest,
"int1" => int1test,
"int2" => int2test,
Expand All @@ -689,6 +831,8 @@ const intlineartests = Dict("knapsack" => knapsacktest,
"indicator2" => indicator2_test,
"indicator3" => indicator3_test,
"indicator4" => indicator4_test,
"semiconttest" => semiconttest,
"semiinttest" => semiinttest
)

@moitestset intlinear
48 changes: 48 additions & 0 deletions test/Test/intlinear.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -42,3 +42,51 @@ MOIU.set_mock_optimize!(mock,
(mock::MOIU.MockOptimizer) -> MOIU.mock_optimize!(mock, [1.25, 8.75, 0., 1.])
)
MOIT.indicator4_test(mock, config)
MOIU.set_mock_optimize!(mock,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 0.0)
MOIU.mock_optimize!(mock, [0.0, 0.0])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 2.0)
MOIU.mock_optimize!(mock, [2.0, 1.0])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 2.0)
MOIU.mock_optimize!(mock, [2.0, 2.0])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 2.5)
MOIU.mock_optimize!(mock, [2.5, 2.5])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 3.0)
MOIU.mock_optimize!(mock, [3.0, 3.0])
end,
(mock::MOIU.MockOptimizer) -> MOI.set(mock, MOI.TerminationStatus(), MOI.INFEASIBLE)
)
MOIT.semiconttest(mock,config)
MOIU.set_mock_optimize!(mock,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 0.0)
MOIU.mock_optimize!(mock, [0.0, 0.0])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 2.0)
MOIU.mock_optimize!(mock, [2.0, 1.0])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 2.0)
MOIU.mock_optimize!(mock, [2.0, 2.0])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 3.0)
MOIU.mock_optimize!(mock, [3.0, 2.5])
end,
(mock::MOIU.MockOptimizer) -> begin
MOI.set(mock, MOI.ObjectiveBound(), 3.0)
MOIU.mock_optimize!(mock, [3.0, 3.0])
end,
(mock::MOIU.MockOptimizer) -> MOI.set(mock, MOI.TerminationStatus(), MOI.INFEASIBLE)
)
MOIT.semiinttest(mock,config)