annotate liboctave/QP.h @ 7017:a1dbe9d80eee

[project @ 2007-10-12 21:27:11 by jwe]
author jwe
date Fri, 12 Oct 2007 21:27:37 +0000
parents 93c65f2a5668
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
1 /*
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
2
7017
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
3 Copyright (C) 1993, 1994, 1995, 1996, 1997, 2005, 2007 John W. Eaton
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
4
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
5 This file is part of Octave.
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
6
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
7 Octave is free software; you can redistribute it and/or modify it
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
8 under the terms of the GNU General Public License as published by the
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
10 option) any later version.
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
11
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
12 Octave is distributed in the hope that it will be useful, but WITHOUT
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
15 for more details.
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
16
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
17 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
19 <http://www.gnu.org/licenses/>.
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
20
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
21 */
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
22
382
e02d6b664394 [project @ 1994-03-09 21:08:03 by jwe]
jwe
parents: 238
diff changeset
23 #if !defined (octave_QP_h)
e02d6b664394 [project @ 1994-03-09 21:08:03 by jwe]
jwe
parents: 238
diff changeset
24 #define octave_QP_h 1
e02d6b664394 [project @ 1994-03-09 21:08:03 by jwe]
jwe
parents: 238
diff changeset
25
463
0bfa99f634be [project @ 1994-06-06 00:59:08 by jwe]
jwe
parents: 384
diff changeset
26 #include "dMatrix.h"
0bfa99f634be [project @ 1994-06-06 00:59:08 by jwe]
jwe
parents: 384
diff changeset
27 #include "dColVector.h"
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
28 #include "Bounds.h"
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
29 #include "LinConst.h"
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
30 #include "base-min.h"
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
31
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
32 class
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
33 QP : public base_minimizer
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
34 {
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
35 public:
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
36
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
37 QP (void)
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
38 : base_minimizer (), H (), c (), bnds (), lc () { }
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
39
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
40 QP (const ColumnVector& x, const Matrix& H_arg)
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
41 : base_minimizer (x), H (H_arg), c (), bnds (), lc ()
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
42 { make_h_symmetric (); }
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
43
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
44 QP (const ColumnVector& x, const Matrix& H_arg, const ColumnVector& c_arg)
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
45 : base_minimizer (x), H (H_arg), c (c_arg), bnds (), lc ()
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
46 { make_h_symmetric (); }
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
47
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
48 QP (const ColumnVector& x, const Matrix& H_arg, const Bounds& b)
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
49 : base_minimizer (x), H (H_arg), c (), bnds (b), lc ()
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
50 { make_h_symmetric (); }
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
51
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
52 QP (const ColumnVector& x, const Matrix& H_arg, const LinConst& l)
1866
b6c48195b552 [project @ 1996-02-04 11:07:17 by jwe]
jwe
parents: 1861
diff changeset
53 : base_minimizer (x), H (H_arg), c (), bnds (), lc (l)
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
54 { make_h_symmetric (); }
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
55
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
56 QP (const ColumnVector& x, const Matrix& H_arg, const ColumnVector& c_arg,
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
57 const Bounds& b)
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
58 : base_minimizer (x), H (H_arg), c (c_arg), bnds (b), lc ()
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
59 { make_h_symmetric (); }
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
60
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
61 QP (const ColumnVector& x, const Matrix& H_arg, const ColumnVector& c_arg,
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
62 const LinConst& l)
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
63 : base_minimizer (x), H (H_arg), c (c_arg), bnds (), lc (l)
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
64 { make_h_symmetric (); }
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
65
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
66 QP (const ColumnVector& x, const Matrix& H_arg, const Bounds& b,
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
67 const LinConst& l)
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
68 : base_minimizer (x), H (H_arg), c (), bnds (b), lc (l)
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
69 { make_h_symmetric (); }
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
70
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
71 QP (const ColumnVector& x, const Matrix& H_arg, const ColumnVector& c_arg,
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
72 const Bounds& b, const LinConst& l)
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
73 : base_minimizer (x), H (H_arg), c (c_arg), bnds (b), lc (l)
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
74 { make_h_symmetric (); }
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
75
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
76 QP (const QP& qp)
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
77 : base_minimizer (qp), H (qp.H), c (qp.c), bnds (qp.bnds), lc (qp.lc) { }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
78
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
79 QP& operator = (const QP& qp)
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
80 {
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
81 if (this != &qp)
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
82 {
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
83 base_minimizer::operator = (qp);
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
84
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
85 H = qp.H;
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
86 c = qp.c;
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
87 bnds = qp.bnds;
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
88 lc = qp.lc;
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
89 }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
90 return *this;
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
91 }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
92
1866
b6c48195b552 [project @ 1996-02-04 11:07:17 by jwe]
jwe
parents: 1861
diff changeset
93 virtual ~QP (void) { }
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
94
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
95 Matrix hessian (void) const { return H; }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
96
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
97 ColumnVector linear_obj_coeff (void) const { return c; }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
98
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
99 Bounds bounds (void) const { return bnds; }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
100
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
101 LinConst linear_constraints (void) const { return lc; }
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
102
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
103 protected:
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
104
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
105 Matrix H;
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
106 ColumnVector c;
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
107 Bounds bnds;
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
108 LinConst lc;
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
109
1861
620a65533630 [project @ 1996-02-04 10:17:34 by jwe]
jwe
parents: 1532
diff changeset
110 private:
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
111
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
112 Matrix make_h_symmetric (void) { return 0.5 * (H + H.transpose ()); }
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
113 };
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
114
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
115 #endif
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
116
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
117 /*
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
118 ;;; Local Variables: ***
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
119 ;;; mode: C++ ***
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
120 ;;; End: ***
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
121 */