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!