julia - Julia+JuMP: 関数の可変数の引数

Question

私はJuMPを使用して、変数の数がユーザーによって決定される非線形問題を解決しようとしています-つまり、コンパイル時にはわかりません。

これを実現するために、@NLobjective行は次のようになります。

@eval @JuMP.NLobjective(m, Min, $(Expr(:call, :myf, [Expr(:ref, :x, i) for i=1:n]...)))

たとえば、 ifの場合n=3、コンパイラは次の行と同じものとして解釈します。

@JuMP.NLobjective(m, Min, myf(x[1], x[2], x[3]))

問題は@eval、グローバル スコープでのみ機能し、関数に含まれているとエラーがスローされることです。

私の質問は、関数のコンパイル時に既知ではないローカルスコープ内で、この同じ機能を達成するにはどうすればよいですか?可変数の引数で@NLobjective呼び出すことができますか?myfx[1],...,x[n]

def testme(n)
    myf(a...) = sum(collect(a).^2)

    m = JuMP.Model(solver=Ipopt.IpoptSolver())

    JuMP.register(m, :myf, n, myf, autodiff=true)
    @JuMP.variable(m, x[1:n] >= 0.5)

    @eval @JuMP.NLobjective(m, Min, $(Expr(:call, :myf, [Expr(:ref, :x, i) for i=1:n]...)))
    JuMP.solve(m)
end

testme(3)

ありがとう！