still the problem with equality constraints

[xiaoyuanzh]

If I put CE and ce0 to 0 it says terminate called after throwing an instance of 'std::runtime_error' what(): Constraints are linearly dependent

[liuq]

Yes, this is true since the equality constraints must be linearly independent. If you do not have equality constraints set their number to zero.

[xiaoyuanzh]

Oh, thx you a lot. Yes, If I write this "m = 0; CE.resize(n, m); CE[0][0] = 0.0; CE[1][0] = 0.0; ce0.resize(m);" It works, and your code run very fast, which I appreciate it a lot. By the way, have you considered projection gradient method on a linear constrainted face. In fact, that is the reason I find your code!