Skip to content

Use set_expr_bound #1077

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 2 commits into from
Aug 29, 2024
Merged

Use set_expr_bound #1077

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
7c169dc
Use set_expr_bound
manuelma Aug 29, 2024
9e9861c
Merge branch 'v0.7.2.x-release' into use_set_expr_bound
manuelma Aug 29, 2024
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
76 changes: 13 additions & 63 deletions src/variables/variable_common.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -157,10 +157,12 @@ function _finalize_expressions!(m, name, args...)
inds, exprs, bins, ints, lbs, ubs, internal_fix_values, fix_values = _collect_info(m, args...)
_set_binary.(exprs, bins)
_set_integer.(exprs, ints)
_set_lower_bound.(exprs, lbs, Symbol(name, :_lb), inds)
_set_upper_bound.(exprs, ubs, Symbol(name, :_ub), inds)
_fix.(exprs, internal_fix_values, Symbol(name, :_fix_value), inds)
_fix.(exprs, fix_values, Symbol(name, :_internal_fix_value), inds)
m.ext[:spineopt].constraints[Symbol(name, :_lb)] = Dict(zip(inds, set_expr_bound.(exprs, >=, lbs)))
m.ext[:spineopt].constraints[Symbol(name, :_ub)] = Dict(zip(inds, set_expr_bound.(exprs, <=, ubs)))
m.ext[:spineopt].constraints[Symbol(name, :_internal_fix)] = Dict(
zip(inds, set_expr_bound.(exprs, ==, internal_fix_values))
)
m.ext[:spineopt].constraints[Symbol(name, :_fix)] = Dict(zip(inds, set_expr_bound.(exprs, ==, fix_values)))
end

function _collect_info(m, d, bin, int, lb, ub, fix_value, internal_fix_value)
Expand Down Expand Up @@ -195,75 +197,23 @@ _set_binary(expr::GenericAffExpr, bin) = bin && set_binary.(keys(expr.terms))
_set_integer(var::VariableRef, int) = int && set_integer(var)
_set_integer(expr::GenericAffExpr, int) = int && set_integer.(keys(expr.terms))

_set_lower_bound(::VariableRef, ::Nothing) = nothing
_set_lower_bound(_var, ::Nothing) = nothing
_set_lower_bound(var::VariableRef, bound::Call) = set_lower_bound(var, bound)
_set_lower_bound(var::VariableRef, bound::Number) = isfinite(bound) && set_lower_bound(var, bound)
_set_lower_bound(expr::GenericAffExpr, bound, name, ind) = _set_bound(expr, >=, bound, name, ind)

_set_upper_bound(::VariableRef, ::Nothing) = nothing
_set_upper_bound(_var, ::Nothing) = nothing
_set_upper_bound(var::VariableRef, bound::Call) = set_upper_bound(var, bound)
_set_upper_bound(var::VariableRef, bound::Number) = isfinite(bound) && set_upper_bound(var, bound)
_set_upper_bound(expr::GenericAffExpr, bound, name, ind) = _set_bound(expr, <=, bound, name, ind)

_fix(::VariableRef, ::Nothing) = nothing
_fix(var::VariableRef, x::Call) = fix(var, x)
function _fix(var::VariableRef, x::Number)
if !isnan(x)
fix(var, x; force=true)
_fix(_var, ::Nothing) = nothing
_fix(var::VariableRef, value::Call) = fix(var, value)
function _fix(var::VariableRef, value::Number)
if !isnan(value)
fix(var, value; force=true)
elseif is_fixed(var)
unfix(var)
end
end
_fix(expr::GenericAffExpr, x, name, ind) = _set_bound(expr, ==, x, name, ind)

_set_bound(_expr, _sense, _bound::Nothing, _name, _ind) = nothing
function _set_bound(expr, sense, bound::Call, name, ind)
upd = _ExpressionBoundUpdate(expr, sense, bound, name, ind)
upd()
end
function _set_bound(expr, sense, bound::Number, name, ind)
m = owner_model(expr)
m === nothing && return
bounds = get!(m.ext[:spineopt].constraints, name, Dict())
existing_constraint = get(bounds, ind, nothing)
if existing_constraint !== nothing
_update_bound_constraint(existing_constraint, expr, sense, bound)
elseif isfinite(bound)
new_constraint = build_sense_constraint(expr, sense, bound)
bounds[ind] = add_constraint(m, new_constraint)
end
end

function _update_bound_constraint(existing_constraint, expr, ::typeof(==), bound)
if isnan(bound)
for var in keys(expr.terms)
set_normalized_coefficient(existing_constraint, var, 0)
end
set_normalized_rhs(existing_constraint, 0)
else
for (var, coeff) in expr.terms
set_normalized_coefficient(existing_constraint, var, realize(coeff))
end
set_normalized_rhs(existing_constraint, bound - realize(expr.constant))
end
end
function _update_bound_constraint(existing_constraint, expr, _sense, bound)
if isfinite(bound)
set_normalized_rhs(existing_constraint, bound - realize(expr.constant))
end
end

struct _ExpressionBoundUpdate
expr
sense
bound
name
ind
end

function (upd::_ExpressionBoundUpdate)()
_set_bound(upd.expr, upd.sense, realize(upd.bound, upd), upd.name, upd.ind)
end

"""
_representative_index(ind)
Expand Down
Loading