The approach detailed previously requires to compute
-the variable that minimizes four convex fonctions
+the variable that minimizes four convex functions
on constraint domains, one
for each primal variable $q_i$, $P_{sh}$, $R_h$, and $x_{hl}$,
Among state of the art for this optimization step, there is