Я пытаюсь импортировать вывод функции из MultivariatePolynomials и использовать ее в JuMp как функцию.Ниже приведен соответствующий код, который является частью большой системы уравнений NL.
using JuMP
using MultivariatePolynomials
@polyvar k s θ_1k θ_2k θ_3k θ_4k θ_5k θ_6k θ_7k θ_8k θ_9k θ_10k θ_11k θ_12k θ_13k E ρ
d=....
j_11=subs(d,k=>q11(θ_1k,θ_2k,θ_3k,θ_4k,θ_5k,θ_6k,θ_7k,θ_8k,θ_9k),s=>f(s,ρ,E))
c_11=subs(j_11, k=>-0.7071067811865476,s=>0)
function v11(θ_1k,θ_2k,θ_3k,θ_4k,θ_5k,θ_6k,θ_7k,θ_8k,θ_9k,ρ,E)
c_11
end
m=Model(solver=IpoptSolver())
JuMP.register(m, :v11, 27, v11, autodiff=true)
@NLconstraint(m,cons1, ((v11(θ_1k,θ_2k,θ_3k,θ_4k,θ_5k,θ_6k,θ_7k,θ_8k,θ_9k,θ_10k,θ_11k,θ_12k,θ_13k,ρ, E))*w_1+*β*(1-δ+...
solve(m)
status=solve(m)
Я получаю следующую ошибку:
MethodError: no method matching parseNLExpr_runtime(::JuMP.Model, ::DynamicPolynomials.PolyVar{true}, ::Array{ReverseDiffSparse.NodeData,1}, ::Int64, ::Array{Float64,1})
Closest candidates are:
parseNLExpr_runtime(::JuMP.Model, ::Number, ::Any, ::Any, ::Any) at C:\Users\s1364322\.julia\v0.6\JuMP\src\parsenlp.jl:196
parseNLExpr_runtime(::JuMP.Model, ::JuMP.Variable, ::Any, ::Any, ::Any) at C:\Users\s1364322\.julia\v0.6\JuMP\src\parsenlp.jl:202
Кажется, что вывод c_11
не был преобразован как тип многомерной функции, как я пытался определить функцию v11
, но остался как полиномиальный тип.Как я могу решить проблему?Пожалуйста, я очень ценю вашу помощь