[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem
From: |
asherikov |
Subject: |
[Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem |
Date: |
Tue, 19 Feb 2013 22:54:52 +0000 |
User-agent: |
Uzbl (Webkit 1.4.3) (FreeBSD amd64 [amd64]) |
URL:
<http://savannah.gnu.org/bugs/?38353>
Summary: qp() fails on a feasible problem
Project: GNU Octave
Submitted by: asherikov
Submitted on: Tue 19 Feb 2013 10:54:50 PM GMT
Category: None
Severity: 3 - Normal
Priority: 5 - Normal
Item Group: Incorrect Result
Status: None
Assigned to: None
Originator Name:
Originator Email:
Open/Closed: Open
Discussion Lock: Any
Release: 3.6.1
Operating System: *BSD
_______________________________________________________
Details:
All matrices, which are used below, are saved in the attached file.
When qp() is called in the following way
[X, OBJ, INFO, LAMBDA] = qp ([], H, q, Ge, ge, [], [], [], G, G_ub)
it fails with status 6 'The problem is infeasible.'
However, when the known solution is given as initial guess
[X, OBJ, INFO, LAMBDA] = qp (solution, H, q, Ge, ge, [], [], [], G, G_ub)
the function finishes successfully.
I've made a quick check and found out that the last part of this condition in
qp.m fails
if ((status == 180 || status == 181 || status == 151)
&& all (abs (P(n-n_eq+1:end)) < rtol * (1 + norm (btmp))))
I suspect that the problem is related to tolerances used in qp() and glpk().
_______________________________________________________
File Attachments:
-------------------------------------------------------
Date: Tue 19 Feb 2013 10:54:50 PM GMT Name: qp_debug.dat.gz Size: 4kB By:
asherikov
<http://savannah.gnu.org/bugs/download.php?file_id=27483>
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?38353>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem,
asherikov <=
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, Rik, 2013/02/22
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, Michael Goffioul, 2013/02/22
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, Rik, 2013/02/22
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, Michael Goffioul, 2013/02/22
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, Rik, 2013/02/26
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, Michael Goffioul, 2013/02/26
- [Octave-bug-tracker] [bug #38353] qp() fails on a feasible problem, asherikov, 2013/02/28