0

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

4

1 回答 1

0

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

于 2018-08-18T17:27:51.897 回答