时间戳

首页 > 解决方案 > 将 MutivariatePolynomials 输出作为用户定义函数导入 JuMP

julia - 将 MutivariatePolynomials 输出作为用户定义函数导入 JuMP

问题描述

我正在尝试从 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，而是保持为多项式类型。我该如何解决这个问题？拜托，非常感谢您的帮助

标签： juliapolynomialsnonlinear-optimizationjulia-jump

解决方案

以供日后参考，已在讨论中解决。

推荐阅读