線量JuMP.setRHSはジュリアの非線形モデルに適用されますか?
Aug 22 2020
次のコードでは、NL制約のRHSを変更する必要があります。しかし、エラーが発生します。エラー:UndefVarError:setRHSが定義されていません。このエラーが発生する理由を教えてください。あなたの助けに感謝します
using JuMP,CPLEX, Ipopt
#parameters--------------------------------------------------------
sig=0.86;
#---------------------------------------------------------------------------
ALT= Model(optimizer_with_attributes(Juniper.Optimizer, "nl_solver"=>optimizer_with_attributes(Ipopt.Optimizer, "print_level" => 0),
"mip_solver"=>optimizer_with_attributes(CPLEx.Optimizer, "logLevel" => 0),"registered_functions" =>[Juniper.register( :f, 1, f; autodiff = true)])
)
# variables-----------------------------------------------------------------
f(x) = cdf(Normal(0, 1), x);
JuMP.register(ALT, :f, 1, f; autodiff = true);
@variable(ALT, h >= 0.1);
@variable(ALT, L >= 0.0001);
@variable(ALT, n>=2, Int);
#-------------------------------------------------------------------
@NLexpression(ALT,k7,1-f(L-sig*sqrt(n))+f(-L-sig*sqrt(n)));
@NLexpression(ALT,f2,1/k7)
#constraints--------------------------------------------------------
@NLconstraint(ALT, f(-L) <= 1/400);
@NLconstraint(ALT,rf2,f2<=10000);
#-------------------------------------------------------------------
@NLobjective(ALT, Min, f2)
optimize!(ALT)
JuMP.setRHS(rf2,getvalueNLobjective(1/k7))
回答
1 PrzemyslawSzufel Aug 22 2020 at 14:24
古いJuMPバージョンの例を使用しています。今日の時点であなたは使用する必要がありますset_normalized_rhs
:
set_normalized_rhs(con_ref, some_rhs_value)
これにより、JuMPによって事前計算された後の正規化されたRHSが設定されることに注意してください。たとえば@constraint(model, 2x - 5 <= 2)
、正規化された値の場合は7
です。
も参照してください https://jump.dev/JuMP.jl/v0.21/constraints/#Constraint-modifications-1 詳細については。