 vector v
 and with further es_constraints given in 'es_constraints':

 Input: dmu_min, dmu_max, sign(v), rho, dmu_limit, dmu_fix, and dmu_eqconstraint

 But DISREGARD all other entries in 'es_constraints' (about mu, c, c0, Keq)

 Compute some dmu vectors at the boundaries of the feasible region
 set by the (previously sampled or determined) flux directions 
 compute the corresponding mu vectors using the pseudoinverse of N'</pre></div>

<!-- crossreference -->
<h2><a name="_cross"></a>CROSS-REFERENCE INFOR