minimize
linearize = 1
  entropyr(matrix_h)
Constraint: >= vector_bl, <= vector_b
  linearmulti(matrix_a)
Constraint: == 1.0
  linearmulti(matrix_aeq)
Box: >= 0.0