Problem_ssd, maximize
  linear(matrix_sde)
Constraint: = 1
  linear(matrix_budget)
MultiConstraint: <= vector_ub
  pm_pen(vector_th, matrix_sde)
Box: >=0, <=1
Solver: precision = 9