Bug get/setting ConstraintSet with variable bridges about mathoptinterface.jl HOT 9 CLOSED

julia> import MathOptInterface as MOI

julia> src = MOI.Utilities.UniversalFallback(MOI.Utilities.Model{Float64}())
MOIU.UniversalFallback{MOIU.Model{Float64}}
fallback for MOIU.Model{Float64}

julia> x = MOI.add_variable(src)
MOI.VariableIndex(1)

julia> MOI.add_constraint(src, x, MOI.GreaterThan(1.0))
MathOptInterface.ConstraintIndex{MathOptInterface.VariableIndex, MathOptInterface.GreaterThan{Float64}}(1)

julia> c = MOI.add_constraint(src, 2.0 * x, MOI.EqualTo(3.0))
MathOptInterface.ConstraintIndex{MathOptInterface.ScalarAffineFunction{Float64}, MathOptInterface.EqualTo{Float64}}(1)

julia> MOI.Utilities.@model(
           Model2452,
           (),
           (),
           (MOI.Nonnegatives, MOI.Zeros),
           (),
           (),
           (),
           (MOI.VectorOfVariables,),
           (MOI.VectorAffineFunction,)
       )
MathOptInterface.Utilities.GenericModel{T, MathOptInterface.Utilities.ObjectiveContainer{T}, MathOptInterface.Utilities.VariablesContainer{T}, Model2452FunctionConstraints{T}} where T

julia> function MOI.supports_constraint(
           ::Model2452{T},
           ::Type{MOI.VariableIndex},
           ::Type{
               <:Union{
                   MOI.GreaterThan{T},
                   MOI.LessThan{T},
                   MOI.EqualTo{T},
                   MOI.Interval{T},
                   MOI.ZeroOne,
                   MOI.Integer,
               },
           },
       ) where {T}
           return false
       end

julia> function MOI.supports_constraint(
           ::Model2452{T},
           ::Type{MOI.VectorOfVariables},
           ::Type{MOI.Reals},
       ) where {T}
           return false
       end

julia> dest = MOI.instantiate(Model2452{Float64}; with_bridge_type = Float64)
MOIB.LazyBridgeOptimizer{MOIU.GenericModel{Float64, MOIU.ObjectiveContainer{Float64}, MOIU.VariablesContainer{Float64}, Model2452FunctionConstraints{Float64}}}
with 0 variable bridges
with 0 constraint bridges
with 0 objective bridges
with inner model MOIU.GenericModel{Float64, MOIU.ObjectiveContainer{Float64}, MOIU.VariablesContainer{Float64}, Model2452FunctionConstraints{Float64}}

julia> index_map = MOI.copy_to(dest, src)
MathOptInterface.Utilities.IndexMap with 3 entries:
  MOI.VariableIndex(1)             => MOI.VariableIndex(-1)
  ConstraintIndex{VariableIndex, … => ConstraintIndex{VariableIndex, GreaterThan{Float64}}(-1)
  ConstraintIndex{ScalarAffineFun… => ConstraintIndex{ScalarAffineFunction{Float64}, EqualTo{Float64…

julia> MOI.get(dest, MOI.ConstraintSet(), index_map[c])
MathOptInterface.EqualTo{Float64}(1.0)

julia> using JuMP, Clarabel

julia> model = Model(Clarabel.Optimizer)
A JuMP Model
Feasibility problem with:
Variables: 0
Model mode: AUTOMATIC
CachingOptimizer state: EMPTY_OPTIMIZER
Solver name: Clarabel

julia> MOI.Bridges.remove_bridge(
           backend(model).optimizer,
           MOI.Bridges.Variable.VectorizeBridge{Float64},
       )

julia> set_silent(model)

julia> @variable(model, x >= 1)
x

julia> @constraint(model, c, 2x == 3)
c : 2 x = 3

julia> optimize!(model)

julia> value(x)
1.4999999999999996

julia> b = normalized_rhs(c)
3.0

julia> set_normalized_rhs(c, b)

julia> optimize!(model)

julia> value(x)
1.4999999999999996

julia> using JuMP, Clarabel

julia> model = Model(Clarabel.Optimizer)
A JuMP Model
Feasibility problem with:
Variables: 0
Model mode: AUTOMATIC
CachingOptimizer state: EMPTY_OPTIMIZER
Solver name: Clarabel

julia> set_silent(model)

julia> @variable(model, x >= 1)
x

julia> @constraint(model, c, 2x == 3)
c : 2 x = 3

julia> optimize!(model)

julia> value(x)
1.4999999999999996

julia> b = normalized_rhs(c)
3.0

julia> set_normalized_rhs(c, b)

julia> optimize!(model)

julia> value(x)
0.9999999995620605

(clarabel) pkg> st
Status `/private/tmp/clarabel/Project.toml`
  [61c947e1] Clarabel v0.7.1
  [4076af6c] JuMP v1.21.0
  [b8f27783] MathOptInterface v1.27.1 `https://github.com/jump-dev/MathOptInterface.jl.git#master`