1

我正在设置一个简单的 QP 来测试 ROI R 包。但是，当它不受约束时，包不能为简单的玩具问题提供错误的解决方案。

例子，

# Maximize -1/2 x^2 + x, no constraints
> x <- OP(Q_objective(as.matrix(-1), 1), maximum = TRUE, bounds = V_bound(lb = -Inf, ub = Inf, nobj = 1))
> x
> ROI Optimization Problem:

  Maximize a quadratic objective function of length 1 with
  - 1 continuous objective variable,

  subject to
  - 0 constraints
  - 1 lower and 0 upper non-standard variable bounds.

看起来不错。但是当我解决我遇到的问题时，

> sol <- ROI_solve(x, solver = 'qpoases')
> sol
  No optimal solution found.
  The solver message was: Initialisation failed! QP could not be solved!
  The objective value is: 0.000000e+00
> sol$solution
  [1] 0
> sol$objval
  [1] 0
> sol$status
  $code
  [1] 1

  $msg
   solver qpoases
    code 36
   symbol RET_INIT_FAILED_HOTSTART
   message Initialisation failed! QP could not be solved!
   roi_code 1

那很奇怪。在相关说明中，当我使用 quadprog 求解器时，我能够获得不受约束的解决方案 (= 1)，但是由于其他原因，我不得不从使用 quadprog 切换到 qpoases。

任何帮助深表感谢。

编辑：

奇怪的是，这行得通，

> ROI_solve(OP(Q_objective(as.matrix(-1), 1), maximum = TRUE), solver = 'qpoases')
  No optimal solution found.
  The solver message was: Initialisation failed! QP could not be solved!
  The objective value is: 0.000000e+00
> ROI_solve(OP(Q_objective(as.matrix(1), -1), maximum = FALSE), solver = 'qpoases')
  Optimal solution found.
  The objective value is: -5.000000e-01
4

1 回答 1

2

hessian_type如果参数hessian_type不是由用户 ROI.plugin.qpoases 选择的，则此不同结果是由参数中设置的不同值导致的hessian_type。由于在为第一个示例选择粗麻布类型时存在错误，它为第一个示例选择hessian_type = 6（未知），为第二个示例选择正确的hessian_type = 1身份。

x1 <- OP(Q_objective(as.matrix(-1), 1), maximum = TRUE, bounds = NULL)
s1 <- ROI_solve(x1, solver = 'qpoases', hessian_type = 1L)
solution(s1)

x1 <- OP(Q_objective(as.matrix(-1), 1), maximum = TRUE, bounds = NULL)
s1 <- ROI_solve(x1, solver = 'qpoases', hessian_type = 6L)
solution(s1)

这在新版本中已修复，新版本正在向 CRAN 发送。

于 2018-12-20T07:38:52.613 回答