El solucionador no admite restricciones del tipo MathOptInterface.VariableIndex-in-MathOptInterface.ZeroOne. Cuando se utiliza el solucionador Ipopt en Julia.

Necesito ayuda para resolver el problema de las restricciones del tipo MathOptInterface.VariableIndex-in-MathOptInterface.ZeroOne no son compatible con el solucionador. Al usar el solucionador Ipopt en Julia para el problema de BinPacking.

<h1>declaración de variables</h1>

@variable(model, x[1:length(packages[!,:PackagingID]), 1:length(items[!,:ItemNumber]), 1:length(bins[!,:binID])], Bin)
@variable(model, c[1:length(packages[!,:PackagingID]), 1:length(bins[!,:binID])], Bin)
function item<em>rule(item)
    return sum(x[i,item,k] for i in 1:length(packages[!,:PackagingID]) for k in 1:length(bins[!,:binID])) 
end
function pack</em>rule(pack,bin)
    return sum(x[pack,j,bin]<em>items.Weight[j] for j in 1:length(items[!,:ItemNumber]))
end
function pack_weight(pack,bin)
    return c[pack,bin]</em>packages.MaxWeight[pack]
end
function binpack_rule(bin)
    return sum(c[i,bin] for i in 1:length(packages[!,:PackagingID]))
end

<h1>restricciones</h1>

for j in 1:length(items[!,:ItemNumber])
    A = @NLconstraint(model, item<em>rule(j) == 1)
end
for i in 1:length(packages[!,:PackagingID]), k in 1:length(bins[!,:binID])
    B = @NLconstraint(model, pack</em>rule(i,k) <= pack<em>weight(i,k))
end
for k in 1:length(bins[!,:binID])
    C = @NLconstraint(model,binpack</em>rule(k) <= 1)
end

<h1>objetivos</h1>

@NLobjective(
    model,
    Min,
    (sum(c[i,k]<em>packages.PackageCost[i] for i in 1:length(packages[!,:PackagingID]) for k in 1:length(bins[!,:binID])) + sum(c[i,k]</em>bins.binNo[k] for i in 1:length(packages[!,:PackagingID]) for k in 1:length(bins[!,:binID])) ),
)

Aquí está la captura de pantalla del error que obtuve: