julia - JuMP 约束中的矩阵逆

Question

我试图遵循这里提到的方法：JuMP constraints involved matrix inverse。但我仍然无法让我的代码运行。

我的代码如下：

using JuMP, Ipopt, LinearAlgebra
FP = Model(solver=IpoptSolver())
@variable(FP, x[1:2,1:2] >= 0)
@objective(FP, Max, 0)

@NLconstraint(FP, inv(x) <= 0.5*I)

status = solve(FP)

我收到以下错误：

ERROR: LoadError: Unexpected object x[i,j] >= 0 for all i in {1,2}, j in {1,2} in nonlinear expression.

我不确定出了什么问题。我正在使用 JuMP 0.18.6。能否请你帮忙？谢谢。