coin-or / qpOASES

Open-source C++ implementation of the recently proposed online active set strategy
GNU Lesser General Public License v2.1
365 stars 127 forks source link

Stuck with stepsize is 0 warning #48

Closed svigerske closed 4 years ago

svigerske commented 4 years ago

Issue created by migration from Trac.

Original creator: tristanc

Original creation time: 2016-04-14 14:29:09

Assignee: ferreau

Version: 3.2.0

I am trying to solve a quadratic problem but each time I am facing the same warning : WARNING: Stepsize is 0.000000000000000e+00 and then the following error occurs ERROR: Premature homotopy termination because QP is infeasible Any ideas would be very appreciated. thanks

svigerske commented 4 years ago

Comment by ferreau created at 2016-04-18 12:37:28

Getting a "zero stepsize" warning once in a while is typically not an issue at all. It just indicates that certain constraints exhibit some redundancy and qpOASES could not make progress along the homotopy path at the corresponding itertation while determining a suitable working set (of linearly independent active constraints).

If, in addition, QP solution fails with an error saying QP problem is infeasible, please try the following:

1) If possible, try to make sure that the QP problem you are solving is actually feasible.

2) If you are sure the QP problem you are solving is feasible, try relaxing the upper/lower limits on variables (bounds) and constraints just a little bit to make sure your QP is not simply found to be infeasible due to round-off errors.

3) Finally, make sure you construct and pass your constraint matrix A correctly. In particular, qpOASES expects all matrices in row-major format (C style) and passing them in FORTRAN-style often leads to an infeasible, but actually very different QP problem.

If all this does not help, please send your QP data (ideally including a main function/script to solve the QP using qpOASES) to the developers, such that we can have a detailled look at the problem.

svigerske commented 4 years ago

Comment by ferreau created at 2016-04-18 12:37:28

Changing status from new to assigned.

svigerske commented 4 years ago

Comment by tristanc created at 2016-04-18 12:58:06

Many thanks for the feedbacks, I'll try yours suggestions and let you know if the QP problem will be feasible.

svigerske commented 4 years ago

Comment by ferreau created at 2017-03-29 07:37:32

Resolution: fixed