Error in line search computation of the mu factor
In sqp.hpp line 333 (experimental branch):
mu = cost_gradient.dot(p) / ((1 - _settings.rho) * constr_l1);
mu is negative and should be positive
mu = abs(cost_gradient.dot(p)) / ((1 - _settings.rho) * constr_l1);
minus sign should work as well