minimize
linearize = 1
  entropyr(matrix_p)
Constraint: == vector_beq
  linearmulti(matrix_aeq)
Box: >= 1e-013
Solver: buldozer