Mercurial > octave
annotate scripts/optimization/sqp.m @ 31653:972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
* sqp.m: Avoid setting lower and upper bounds to +/-realmax when
unspecified in the function call. Set bounds to +/-Inf in main
function. Update constraint subfunction to correctly return empty
distance to constraint when bounds are unspecified.
author | Julien Bect <julien.bect@supelec.fr> |
---|---|
date | Thu, 08 Dec 2022 15:23:29 -0500 |
parents | c05ef94a2bbc |
children | 5f11de0e7440 |
rev | line source |
---|---|
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
1 ######################################################################## |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
2 ## |
30564
796f54d4ddbf
update Octave Project Developers copyright for the new year
John W. Eaton <jwe@octave.org>
parents:
30330
diff
changeset
|
3 ## Copyright (C) 2005-2022 The Octave Project Developers |
27918
b442ec6dda5c
use centralized file for copyright info for individual contributors
John W. Eaton <jwe@octave.org>
parents:
26619
diff
changeset
|
4 ## |
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
5 ## See the file COPYRIGHT.md in the top-level directory of this |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
6 ## distribution or <https://octave.org/copyright/>. |
5289 | 7 ## |
8 ## This file is part of Octave. | |
9 ## | |
24534
194eb4bd202b
maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
10 ## Octave is free software: you can redistribute it and/or modify it |
5289 | 11 ## under the terms of the GNU General Public License as published by |
24534
194eb4bd202b
maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
12 ## the Free Software Foundation, either version 3 of the License, or |
22755
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22323
diff
changeset
|
13 ## (at your option) any later version. |
5289 | 14 ## |
15 ## Octave is distributed in the hope that it will be useful, but | |
16 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
22755
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22323
diff
changeset
|
17 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22323
diff
changeset
|
18 ## GNU General Public License for more details. |
5289 | 19 ## |
20 ## You should have received a copy of the GNU General Public License | |
7016 | 21 ## along with Octave; see the file COPYING. If not, see |
24534
194eb4bd202b
maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
22 ## <https://www.gnu.org/licenses/>. |
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
23 ## |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
24 ######################################################################## |
5289 | 25 |
26 ## -*- texinfo -*- | |
20852
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20735
diff
changeset
|
27 ## @deftypefn {} {[@var{x}, @var{obj}, @var{info}, @var{iter}, @var{nf}, @var{lambda}] =} sqp (@var{x0}, @var{phi}) |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20735
diff
changeset
|
28 ## @deftypefnx {} {[@dots{}] =} sqp (@var{x0}, @var{phi}, @var{g}) |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20735
diff
changeset
|
29 ## @deftypefnx {} {[@dots{}] =} sqp (@var{x0}, @var{phi}, @var{g}, @var{h}) |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20735
diff
changeset
|
30 ## @deftypefnx {} {[@dots{}] =} sqp (@var{x0}, @var{phi}, @var{g}, @var{h}, @var{lb}, @var{ub}) |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20735
diff
changeset
|
31 ## @deftypefnx {} {[@dots{}] =} sqp (@var{x0}, @var{phi}, @var{g}, @var{h}, @var{lb}, @var{ub}, @var{maxiter}) |
28948
36cd3569f702
maint: match names in documentation and input parameters in function.
Rik <rik@octave.org>
parents:
28942
diff
changeset
|
32 ## @deftypefnx {} {[@dots{}] =} sqp (@var{x0}, @var{phi}, @var{g}, @var{h}, @var{lb}, @var{ub}, @var{maxiter}, @var{tolerance}) |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
33 ## Minimize an objective function using sequential quadratic programming (SQP). |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
34 ## |
5289 | 35 ## Solve the nonlinear program |
6741 | 36 ## @tex |
37 ## $$ | |
38 ## \min_x \phi (x) | |
39 ## $$ | |
40 ## @end tex | |
41 ## @ifnottex | |
5289 | 42 ## |
43 ## @example | |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
44 ## @group |
14327
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
45 ## min phi (x) |
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
46 ## x |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
47 ## @end group |
5289 | 48 ## @end example |
49 ## | |
6741 | 50 ## @end ifnottex |
51 ## subject to | |
5289 | 52 ## @tex |
6741 | 53 ## $$ |
8167
17352ccd860e
describe additional arguments in sqp() documentation string
Ivan Sutoris <ivan.sutoris@gmail.com>
parents:
8047
diff
changeset
|
54 ## g(x) = 0 \qquad h(x) \geq 0 \qquad lb \leq x \leq ub |
6741 | 55 ## $$ |
5289 | 56 ## @end tex |
6741 | 57 ## @ifnottex |
5289 | 58 ## |
59 ## @example | |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
60 ## @group |
14327
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
61 ## g(x) = 0 |
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
62 ## h(x) >= 0 |
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
63 ## lb <= x <= ub |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
64 ## @end group |
5289 | 65 ## @end example |
10821
693e22af08ae
Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents:
10678
diff
changeset
|
66 ## |
6741 | 67 ## @end ifnottex |
5289 | 68 ## @noindent |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
69 ## using a sequential quadratic programming method. |
5289 | 70 ## |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
71 ## The first argument is the initial guess for the vector @var{x0}. |
5289 | 72 ## |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
73 ## The second argument is a function handle pointing to the objective function |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
74 ## @var{phi}. The objective function must accept one vector argument and |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
75 ## return a scalar. |
5289 | 76 ## |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
77 ## The second argument may also be a 2- or 3-element cell array of function |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
78 ## handles. The first element should point to the objective function, the |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
79 ## second should point to a function that computes the gradient of the |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
80 ## objective function, and the third should point to a function that computes |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
81 ## the Hessian of the objective function. If the gradient function is not |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
82 ## supplied, the gradient is computed by finite differences. If the Hessian |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
83 ## function is not supplied, a BFGS update formula is used to approximate the |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
84 ## Hessian. |
5289 | 85 ## |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
86 ## When supplied, the gradient function @code{@var{phi}@{2@}} must accept one |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
87 ## vector argument and return a vector. When supplied, the Hessian function |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
88 ## @code{@var{phi}@{3@}} must accept one vector argument and return a matrix. |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
89 ## |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
90 ## The third and fourth arguments @var{g} and @var{h} are function handles |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
91 ## pointing to functions that compute the equality constraints and the |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
92 ## inequality constraints, respectively. If the problem does not have |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
93 ## equality (or inequality) constraints, then use an empty matrix ([]) for |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
94 ## @var{g} (or @var{h}). When supplied, these equality and inequality |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
95 ## constraint functions must accept one vector argument and return a vector. |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
96 ## |
5289 | 97 ## The third and fourth arguments may also be 2-element cell arrays of |
98 ## function handles. The first element should point to the constraint | |
99 ## function and the second should point to a function that computes the | |
100 ## gradient of the constraint function: | |
6741 | 101 ## @tex |
102 ## $$ | |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
103 ## \Bigg( {\partial f(x) \over \partial x_1}, |
6741 | 104 ## {\partial f(x) \over \partial x_2}, \ldots, |
105 ## {\partial f(x) \over \partial x_N} \Bigg)^T | |
106 ## $$ | |
107 ## @end tex | |
108 ## @ifnottex | |
10821
693e22af08ae
Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents:
10678
diff
changeset
|
109 ## |
5289 | 110 ## @example |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
111 ## @group |
14327
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
112 ## [ d f(x) d f(x) d f(x) ] |
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
113 ## transpose ( [ ------ ----- ... ------ ] ) |
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
114 ## [ dx_1 dx_2 dx_N ] |
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
115 ## @end group |
5289 | 116 ## @end example |
10821
693e22af08ae
Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents:
10678
diff
changeset
|
117 ## |
6741 | 118 ## @end ifnottex |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
119 ## The fifth and sixth arguments, @var{lb} and @var{ub}, contain lower and |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
120 ## upper bounds on @var{x}. These must be consistent with the equality and |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
121 ## inequality constraints @var{g} and @var{h}. If the arguments are vectors |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
122 ## then @var{x}(i) is bound by @var{lb}(i) and @var{ub}(i). A bound can also |
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
123 ## be a scalar in which case all elements of @var{x} will share the same |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
124 ## bound. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
125 ## |
13207
88bd1d1d6657
Reword sqp's docstring
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
11587
diff
changeset
|
126 ## The seventh argument @var{maxiter} specifies the maximum number of |
13929
9cae456085c2
Grammarcheck of documentation before 3.6.0 release.
Rik <octave@nomad.inbox5.com>
parents:
13305
diff
changeset
|
127 ## iterations. The default value is 100. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
128 ## |
28948
36cd3569f702
maint: match names in documentation and input parameters in function.
Rik <rik@octave.org>
parents:
28942
diff
changeset
|
129 ## The eighth argument @var{tolerance} specifies the tolerance for the stopping |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19833
diff
changeset
|
130 ## criteria. The default value is @code{sqrt (eps)}. |
8167
17352ccd860e
describe additional arguments in sqp() documentation string
Ivan Sutoris <ivan.sutoris@gmail.com>
parents:
8047
diff
changeset
|
131 ## |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
132 ## The value returned in @var{info} may be one of the following: |
10821
693e22af08ae
Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents:
10678
diff
changeset
|
133 ## |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
134 ## @table @asis |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
135 ## @item 101 |
19597
db92e7e28e1f
strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents:
19596
diff
changeset
|
136 ## The algorithm terminated normally. |
18255
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
137 ## All constraints meet the specified tolerance. |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
138 ## |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
139 ## @item 102 |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
140 ## The BFGS update failed. |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
141 ## |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
142 ## @item 103 |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
143 ## The maximum number of iterations was reached. |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
144 ## |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
145 ## @item 104 |
19597
db92e7e28e1f
strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents:
19596
diff
changeset
|
146 ## The stepsize has become too small, i.e., |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
147 ## @tex |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
148 ## $\Delta x,$ |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
149 ## @end tex |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
150 ## @ifnottex |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
151 ## delta @var{x}, |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
152 ## @end ifnottex |
13207
88bd1d1d6657
Reword sqp's docstring
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents:
11587
diff
changeset
|
153 ## is less than @code{@var{tol} * norm (x)}. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
154 ## @end table |
8167
17352ccd860e
describe additional arguments in sqp() documentation string
Ivan Sutoris <ivan.sutoris@gmail.com>
parents:
8047
diff
changeset
|
155 ## |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
156 ## An example of calling @code{sqp}: |
5289 | 157 ## |
158 ## @example | |
7031 | 159 ## function r = g (x) |
160 ## r = [ sumsq(x)-10; | |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
161 ## x(2)*x(3)-5*x(4)*x(5); |
7031 | 162 ## x(1)^3+x(2)^3+1 ]; |
163 ## endfunction | |
164 ## | |
165 ## function obj = phi (x) | |
14327
4d917a6a858b
doc: Use Octave coding conventions in @example blocks of docstrings.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
166 ## obj = exp (prod (x)) - 0.5*(x(1)^3+x(2)^3+1)^2; |
7031 | 167 ## endfunction |
5289 | 168 ## |
7031 | 169 ## x0 = [-1.8; 1.7; 1.9; -0.8; -0.8]; |
170 ## | |
171 ## [x, obj, info, iter, nf, lambda] = sqp (x0, @@phi, @@g, []) | |
5289 | 172 ## |
7031 | 173 ## x = |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
174 ## |
7031 | 175 ## -1.71714 |
176 ## 1.59571 | |
177 ## 1.82725 | |
178 ## -0.76364 | |
179 ## -0.76364 | |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
180 ## |
7031 | 181 ## obj = 0.053950 |
182 ## info = 101 | |
183 ## iter = 8 | |
184 ## nf = 10 | |
185 ## lambda = | |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
186 ## |
7031 | 187 ## -0.0401627 |
188 ## 0.0379578 | |
189 ## -0.0052227 | |
5289 | 190 ## @end example |
191 ## | |
5642 | 192 ## @seealso{qp} |
5289 | 193 ## @end deftypefn |
194 | |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
195 function [x, obj, info, iter, nf, lambda] = sqp (x0, objf, cef, cif, lb, ub, maxiter, tolerance) |
5289 | 196 |
17336
b81b9d079515
Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents:
17245
diff
changeset
|
197 globals = struct (); # data and handles, needed and changed by subfunctions |
5289 | 198 |
28789
28de41192f3c
Eliminate unneeded verification of nargin, nargout in m-files.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
199 if (nargin < 2 || nargin == 5) |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
200 print_usage (); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
201 endif |
5289 | 202 |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
203 if (! isvector (x0)) |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
204 error ("sqp: X0 must be a vector"); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
205 endif |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
206 if (rows (x0) == 1) |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
207 x0 = x0'; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
208 endif |
5289 | 209 |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
210 have_grd = 0; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
211 have_hess = 0; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
212 if (iscell (objf)) |
11120
a44f979a35ce
style fixes for some .m files
John W. Eaton <jwe@octave.org>
parents:
10821
diff
changeset
|
213 switch (numel (objf)) |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
214 case 1 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
215 obj_fcn = objf{1}; |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
216 obj_grd = @(x, obj) fd_obj_grd (x, obj, obj_fcn); |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
217 case 2 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
218 obj_fcn = objf{1}; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
219 obj_grd = objf{2}; |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
220 have_grd = 1; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
221 case 3 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
222 obj_fcn = objf{1}; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
223 obj_grd = objf{2}; |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
224 obj_hess = objf{3}; |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
225 have_grd = 1; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
226 have_hess = 1; |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
227 otherwise |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
228 error ("sqp: invalid objective function specification"); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
229 endswitch |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
230 else |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
231 obj_fcn = objf; # No cell array, only obj_fcn set |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
232 obj_grd = @(x, obj) fd_obj_grd (x, obj, obj_fcn); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
233 endif |
5289 | 234 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
235 ce_fcn = @empty_cf; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
236 ce_grd = @empty_jac; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
237 if (nargin > 2) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
238 if (iscell (cef)) |
11120
a44f979a35ce
style fixes for some .m files
John W. Eaton <jwe@octave.org>
parents:
10821
diff
changeset
|
239 switch (numel (cef)) |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
240 case 1 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
241 ce_fcn = cef{1}; |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
242 ce_grd = @(x) fd_ce_jac (x, ce_fcn); |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
243 case 2 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
244 ce_fcn = cef{1}; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
245 ce_grd = cef{2}; |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
246 otherwise |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
247 error ("sqp: invalid equality constraint function specification"); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
248 endswitch |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
249 elseif (! isempty (cef)) |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
250 ce_fcn = cef; # No cell array, only constraint equality function set |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
251 ce_grd = @(x) fd_ce_jac (x, ce_fcn); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
252 endif |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
253 endif |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
254 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
255 ci_fcn = @empty_cf; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
256 ci_grd = @empty_jac; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
257 if (nargin > 3) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
258 ## constraint function given by user with possible gradient |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
259 globals.cif = cif; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
260 ## constraint function given by user without gradient |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
261 globals.cifcn = @empty_cf; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
262 if (iscell (cif)) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
263 if (length (cif) > 0) |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
264 globals.cifcn = cif{1}; |
5289 | 265 endif |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
266 elseif (! isempty (cif)) |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
267 globals.cifcn = cif; |
5289 | 268 endif |
269 | |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
270 if (nargin < 5 || (nargin > 5 && isempty (lb) && isempty (ub))) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
271 ## constraint inequality function only without any bounds |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
272 ci_grd = @(x) fd_ci_jac (x, globals.cifcn); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
273 if (iscell (cif)) |
16933
e39f00a32dc7
maint: Use parentheses around condition for switch(),while(),if() statements.
Rik <rik@octave.org>
parents:
16064
diff
changeset
|
274 switch (length (cif)) |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
275 case 1 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
276 ci_fcn = cif{1}; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
277 case 2 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
278 ci_fcn = cif{1}; |
17174
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
279 ci_grd = cif{2}; |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
280 otherwise |
c3c1ebfaa7dc
maint: Use common indentation for switch statement.
Rik <rik@octave.org>
parents:
16933
diff
changeset
|
281 error ("sqp: invalid inequality constraint function specification"); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
282 endswitch |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
283 elseif (! isempty (cif)) |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
284 ci_fcn = cif; # No cell array, only constraint inequality function set |
5289 | 285 endif |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
286 else |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
287 ## constraint inequality function with bounds present |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
288 lb_idx = ub_idx = true (size (x0)); |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
289 ub_grad = - (lb_grad = eye (rows (x0))); |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
290 |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
291 ## if unspecified set ub and lb to +/-Inf, preserving single type |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
292 if (isempty (lb)) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
293 if (isa (x0, "single")) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
294 lb = - inf (size (x0), "single"); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
295 else |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
296 lb = - inf (size (x0)); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
297 endif |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
298 endif |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
299 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
300 if (isvector (lb)) |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
301 lb = lb(:); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
302 lb_idx(:) = (lb != -Inf); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
303 globals.lb = lb(lb_idx, 1); |
11347
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
304 lb_grad = lb_grad(lb_idx, :); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
305 else |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
306 error ("sqp: invalid lower bound"); |
6768 | 307 endif |
308 | |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
309 if (isempty (ub)) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
310 if (isa (x0, "single")) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
311 ub = inf (size (x0), "single"); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
312 else |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
313 ub = inf (size (x0)); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
314 endif |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
315 endif |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
316 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
317 if (isvector (ub)) |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
318 ub = ub(:); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
319 ub_idx(:) = (ub != Inf); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
320 globals.ub = ub(ub_idx, 1); |
11347
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
321 ub_grad = ub_grad(ub_idx, :); |
6768 | 322 else |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
323 error ("sqp: invalid upper bound"); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
324 endif |
6768 | 325 |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
326 if (any (globals.lb > globals.ub)) |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
327 error ("sqp: upper bound smaller than lower bound"); |
6768 | 328 endif |
11347
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
329 bounds_grad = [lb_grad; ub_grad]; |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
330 ci_fcn = @(x) cf_ub_lb (x, lb_idx, ub_idx, globals); |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
331 ci_grd = @(x) cigrad_ub_lb (x, bounds_grad, globals); |
6768 | 332 endif |
5289 | 333 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
334 endif # if (nargin > 3) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
335 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
336 iter_max = 100; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
337 if (nargin > 6 && ! isempty (maxiter)) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
338 if (isscalar (maxiter) && maxiter > 0 && fix (maxiter) == maxiter) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
339 iter_max = maxiter; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
340 else |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
341 error ("sqp: invalid number of maximum iterations"); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
342 endif |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
343 endif |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
344 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
345 tol = sqrt (eps); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
346 if (nargin > 7 && ! isempty (tolerance)) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
347 if (isscalar (tolerance) && tolerance > 0) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
348 tol = tolerance; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
349 else |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
11371
diff
changeset
|
350 error ("sqp: invalid value for TOLERANCE"); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
351 endif |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
352 endif |
5289 | 353 |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
354 ## Initialize variables for search loop |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
355 ## Seed x with initial guess and evaluate objective function, constraints, |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
356 ## and gradients at initial value x0. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
357 ## |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
358 ## obj_fcn -- objective function |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
359 ## obj_grad -- objective gradient |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
360 ## ce_fcn -- equality constraint functions |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
361 ## ci_fcn -- inequality constraint functions |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
362 ## A == [grad_{x_1} cx_fun, grad_{x_2} cx_fun, ..., grad_{x_n} cx_fun]^T |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
363 x = x0; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
364 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
365 obj = feval (obj_fcn, x0); |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
366 globals.nfev = 1; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
367 |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
368 if (have_grd) |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
369 c = feval (obj_grd, x0); |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
370 else |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
371 c = feval (obj_grd, x0, obj); |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
372 endif |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
373 |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
374 ## Choose an initial NxN symmetric positive definite Hessian approximation B. |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
375 n = length (x0); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
376 if (have_hess) |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
377 B = feval (obj_hess, x0); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
378 else |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
379 B = eye (n, n); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
380 endif |
5289 | 381 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
382 ce = feval (ce_fcn, x0); |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
383 F = feval (ce_grd, x0); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
384 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
385 ci = feval (ci_fcn, x0); |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
386 C = feval (ci_grd, x0); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
387 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
388 A = [F; C]; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
389 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
390 ## Choose an initial lambda (x is provided by the caller). |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
391 lambda = 100 * ones (rows (A), 1); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
392 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
393 qp_iter = 1; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
394 alpha = 1; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
395 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
396 info = 0; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
397 iter = 0; |
28942
fc4bb4bd1d5e
maint: Use '##' as lead-in for full-line comments.
Rik <rik@octave.org>
parents:
28896
diff
changeset
|
398 ## report (); # Called with no arguments to initialize reporting |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
399 ## report (iter, qp_iter, alpha, __sqp_nfev__, obj); |
5289 | 400 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
401 while (++iter < iter_max) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
402 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
403 ## Check convergence. This is just a simple check on the first |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
404 ## order necessary conditions. |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
405 nr_f = rows (F); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
406 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
407 lambda_e = lambda((1:nr_f)'); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
408 lambda_i = lambda((nr_f+1:end)'); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
409 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
410 con = [ce; ci]; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
411 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
412 t0 = norm (c - A' * lambda); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
413 t1 = norm (ce); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
414 t2 = all (ci >= 0); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
415 t3 = all (lambda_i >= 0); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
416 t4 = norm (lambda .* con); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
417 |
18255
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
418 ## Normal convergence. All constraints are satisfied |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
419 ## and objective has converged. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
420 if (t2 && t3 && max ([t0; t1; t4]) < tol) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
421 info = 101; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
422 break; |
6382 | 423 endif |
424 | |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
425 ## Compute search direction p by solving QP. |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
426 g = -ce; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
427 d = -ci; |
5289 | 428 |
22224
20e0c0b8820c
Allow sequential quadratic programs with infeasible QPs (bug #36015).
Lachlan Andrew <lachlanbis@gmail.com>
parents:
21758
diff
changeset
|
429 old_lambda = lambda; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
430 [p, obj_qp, INFO, lambda] = qp (x, B, c, F, g, [], [], d, C, |
26332
d4211810202a
qp and sqp: Non-fixed tolerance for qp (bug #53506).
Maor Shutman <maorus12@gmail.com>
parents:
26237
diff
changeset
|
431 Inf (size (d)), struct ("TolX", tol)); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
432 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
433 info = INFO.info; |
5289 | 434 |
18255
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
435 ## FIXME: check QP solution and attempt to recover if it has failed. |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
436 ## For now, just warn about possible problems. |
19593
446c46af4b42
strip trailing whitespace from most source files
John W. Eaton <jwe@octave.org>
parents:
17744
diff
changeset
|
437 |
13206
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
438 id = "Octave:SQP-QP-subproblem"; |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
439 switch (info) |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
440 case 2 |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
441 warning (id, "sqp: QP subproblem is non-convex and unbounded"); |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
442 case 3 |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
443 warning (id, "sqp: QP subproblem failed to converge in %d iterations", |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
444 INFO.solveiter); |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
445 case 6 |
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
446 warning (id, "sqp: QP subproblem is infeasible"); |
22224
20e0c0b8820c
Allow sequential quadratic programs with infeasible QPs (bug #36015).
Lachlan Andrew <lachlanbis@gmail.com>
parents:
21758
diff
changeset
|
447 lambda = old_lambda; # The return value was size 0x0 in this case. |
13206
658aa9fd8665
sqp: provide warnings for QP subproblem failures
John W. Eaton <jwe@octave.org>
parents:
11587
diff
changeset
|
448 endswitch |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
449 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
450 ## Choose mu such that p is a descent direction for the chosen |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
451 ## merit function phi. |
17245
7babcdb9bc13
Use ... instead of \ for line continuation marker.
Stefan Mahr <dac922@gmx.de>
parents:
17174
diff
changeset
|
452 [x_new, alpha, obj_new, globals] = ... |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
453 linesearch_L1 (x, p, obj_fcn, obj_grd, ce_fcn, ci_fcn, lambda, ... |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
454 obj, c, globals); |
26617
98afb8bbd1f6
maint: Strip trailing whitespace from source files.
John W. Eaton <jwe@octave.org>
parents:
26376
diff
changeset
|
455 |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
456 delx = x_new - x; |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
457 |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
458 ## Check if step size has become too small (indicates lack of progress). |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
459 if (norm (delx) < tol * norm (x)) |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
460 info = 104; |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
461 break; |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
462 endif |
5289 | 463 |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
464 ## Evaluate objective function, constraints, and gradients at x_new. |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
465 if (have_grd) |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
466 c_new = feval (obj_grd, x_new); |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
467 else |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
468 c_new = feval (obj_grd, x_new, obj_new); |
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
469 endif |
5289 | 470 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
471 ce_new = feval (ce_fcn, x_new); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
472 F_new = feval (ce_grd, x_new); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
473 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
474 ci_new = feval (ci_fcn, x_new); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
475 C_new = feval (ci_grd, x_new); |
5289 | 476 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
477 A_new = [F_new; C_new]; |
5289 | 478 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
479 ## Set |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
480 ## |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
481 ## s = alpha * p |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
482 ## y = grad_x L (x_new, lambda) - grad_x L (x, lambda}) |
6527 | 483 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
484 y = c_new - c; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
485 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
486 if (! isempty (A)) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
487 t = ((A_new - A)'*lambda); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
488 y -= t; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
489 endif |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
490 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
491 if (have_hess) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
492 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
493 B = feval (obj_hess, x); |
5289 | 494 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
495 else |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
496 ## Update B using a quasi-Newton formula. |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
497 delxt = delx'; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
498 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
499 ## Damped BFGS. Or maybe we would actually want to use the Hessian |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
500 ## of the Lagrangian, computed directly? |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
501 d1 = delxt*B*delx; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
502 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
503 t1 = 0.2 * d1; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
504 t2 = delxt*y; |
5289 | 505 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
506 if (t2 < t1) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
507 theta = 0.8*d1/(d1 - t2); |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
508 else |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
509 theta = 1; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
510 endif |
5289 | 511 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
512 r = theta*y + (1-theta)*B*delx; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
513 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
514 d2 = delxt*r; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
515 |
18255
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
516 ## Check if the next BFGS update will work properly. |
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
517 ## If d1 or d2 vanish, the BFGS update will fail. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
518 if (d1 == 0 || d2 == 0) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
519 info = 102; |
10549 | 520 break; |
5289 | 521 endif |
522 | |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
523 B = B - B*delx*delxt*B/d1 + r*r'/d2; |
5289 | 524 |
525 endif | |
526 | |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
527 x = x_new; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
528 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
529 obj = obj_new; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
530 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
531 c = c_new; |
5289 | 532 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
533 ce = ce_new; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
534 F = F_new; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
535 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
536 ci = ci_new; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
537 C = C_new; |
5289 | 538 |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
539 A = A_new; |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
540 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
541 ## report (iter, qp_iter, alpha, __sqp_nfev__, obj); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
542 |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
543 endwhile |
5289 | 544 |
18255
f58a6cd3f909
sqp.m: Fix return values in sqp (bug #32008).
Arun Giridhar <arungiridhar@gmail.com>
parents:
17744
diff
changeset
|
545 ## Check if we've spent too many iterations without converging. |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
546 if (iter >= iter_max) |
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
547 info = 103; |
5289 | 548 endif |
549 | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
550 nf = globals.nfev; |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
551 |
7399 | 552 endfunction |
5289 | 553 |
554 | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
555 function [merit, obj, globals] = phi_L1 (obj, obj_fcn, ce_fcn, ci_fcn, ... |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
556 x, mu, globals) |
5289 | 557 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
558 ce = feval (ce_fcn, x); |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
559 ci = feval (ci_fcn, x); |
5289 | 560 |
561 idx = ci < 0; | |
562 | |
563 con = [ce; ci(idx)]; | |
564 | |
565 if (isempty (obj)) | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
566 obj = feval (obj_fcn, x); |
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
567 globals.nfev += 1; |
5289 | 568 endif |
569 | |
570 merit = obj; | |
571 t = norm (con, 1) / mu; | |
572 | |
573 if (! isempty (t)) | |
574 merit += t; | |
575 endif | |
576 | |
7399 | 577 endfunction |
5289 | 578 |
579 | |
17245
7babcdb9bc13
Use ... instead of \ for line continuation marker.
Stefan Mahr <dac922@gmx.de>
parents:
17174
diff
changeset
|
580 function [x_new, alpha, obj, globals] = ... |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
581 linesearch_L1 (x, p, obj_fcn, obj_grd, ce_fcn, ci_fcn, lambda, obj, c, globals) |
5289 | 582 |
583 ## Choose parameters | |
584 ## | |
585 ## eta in the range (0, 0.5) | |
586 ## tau in the range (0, 1) | |
587 | |
588 eta = 0.25; | |
589 tau = 0.5; | |
590 | |
591 delta_bar = sqrt (eps); | |
592 | |
593 if (isempty (lambda)) | |
594 mu = 1 / delta_bar; | |
595 else | |
596 mu = 1 / (norm (lambda, Inf) + delta_bar); | |
597 endif | |
598 | |
599 alpha = 1; | |
600 | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
601 ce = feval (ce_fcn, x); |
5289 | 602 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
603 [phi_x_mu, obj, globals] = phi_L1 (obj, obj_fcn, ce_fcn, ci_fcn, x, ... |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
604 mu, globals); |
5289 | 605 |
606 D_phi_x_mu = c' * p; | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
607 d = feval (ci_fcn, x); |
5289 | 608 ## only those elements of d corresponding |
609 ## to violated constraints should be included. | |
610 idx = d < 0; | |
611 t = - norm ([ce; d(idx)], 1) / mu; | |
612 if (! isempty (t)) | |
613 D_phi_x_mu += t; | |
614 endif | |
615 | |
616 while (1) | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
617 [p1, obj, globals] = phi_L1 ([], obj_fcn, ce_fcn, ci_fcn, ... |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
618 x+alpha*p, mu, globals); |
5289 | 619 p2 = phi_x_mu+eta*alpha*D_phi_x_mu; |
620 if (p1 > p2) | |
621 ## Reset alpha = tau_alpha * alpha for some tau_alpha in the | |
622 ## range (0, tau). | |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
623 tau_alpha = 0.9 * tau; # ?? |
5289 | 624 alpha = tau_alpha * alpha; |
625 else | |
626 break; | |
627 endif | |
628 endwhile | |
629 | |
630 x_new = x + alpha * p; | |
631 | |
7399 | 632 endfunction |
5289 | 633 |
634 | |
26237
34cc6edf2041
sqp.m: Eliminate redundant objective function calls (patch #8073).
Ken Marek <gm.kmarek@gmail.com>
parents:
25054
diff
changeset
|
635 function grd = fdgrd (f, x, y0) |
5289 | 636 |
637 if (! isempty (f)) | |
638 nx = length (x); | |
639 grd = zeros (nx, 1); | |
640 deltax = sqrt (eps); | |
641 for i = 1:nx | |
642 t = x(i); | |
643 x(i) += deltax; | |
644 grd(i) = (feval (f, x) - y0) / deltax; | |
645 x(i) = t; | |
646 endfor | |
647 else | |
648 grd = zeros (0, 1); | |
649 endif | |
650 | |
7399 | 651 endfunction |
5289 | 652 |
653 | |
654 function jac = fdjac (f, x) | |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
655 |
6768 | 656 nx = length (x); |
5289 | 657 if (! isempty (f)) |
658 y0 = feval (f, x); | |
659 nf = length (y0); | |
660 nx = length (x); | |
661 jac = zeros (nf, nx); | |
662 deltax = sqrt (eps); | |
663 for i = 1:nx | |
664 t = x(i); | |
665 x(i) += deltax; | |
666 jac(:,i) = (feval (f, x) - y0) / deltax; | |
667 x(i) = t; | |
668 endfor | |
669 else | |
670 jac = zeros (0, nx); | |
671 endif | |
672 | |
7399 | 673 endfunction |
5289 | 674 |
675 | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
676 function grd = fd_obj_grd (x, obj, obj_fcn) |
5289 | 677 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
678 grd = fdgrd (obj_fcn, x, obj); |
5289 | 679 |
7399 | 680 endfunction |
5289 | 681 |
682 | |
683 function res = empty_cf (x) | |
684 | |
685 res = zeros (0, 1); | |
686 | |
7399 | 687 endfunction |
5289 | 688 |
689 | |
690 function res = empty_jac (x) | |
691 | |
692 res = zeros (0, length (x)); | |
693 | |
7399 | 694 endfunction |
5289 | 695 |
696 | |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
697 function jac = fd_ce_jac (x, ce_fcn) |
5289 | 698 |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
699 jac = fdjac (ce_fcn, x); |
5289 | 700 |
7399 | 701 endfunction |
5289 | 702 |
703 | |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
704 function jac = fd_ci_jac (x, cifcn) |
5289 | 705 |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
706 ## cifcn = constraint function without gradients and lb or ub |
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
707 jac = fdjac (cifcn, x); |
6768 | 708 |
7399 | 709 endfunction |
6768 | 710 |
7017 | 711 |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
712 function res = cf_ub_lb (x, lbidx, ubidx, globals) |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
713 ## function returning constraint evaluated at x, and distance between ub/lb |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
714 ## and x, when they are specified (otherwise return empty) |
5289 | 715 |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
716 ## inequality constraints |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
717 if (! isempty (globals.cifcn)) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
718 ci = feval (globals.cifcn, x); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
719 else |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
720 ci = []; |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
721 endif |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
722 |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
723 ## lower bounds |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
724 if (! isempty (globals.lb)) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
725 lb = x(lbidx, 1) - globals.lb; |
11347
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
726 else |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
727 lb = []; |
6768 | 728 endif |
5289 | 729 |
31653
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
730 ## upper bounds |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
731 if (! isempty (globals.ub)) |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
732 ub = globals.ub - x(ubidx, 1); |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
733 else |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
734 ub = []; |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
735 end |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
736 |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
737 res = [ci; lb; ub]; |
972dcc46bb41
sqp.m: Ignore missing lower/upper bounds (bug #32008)
Julien Bect <julien.bect@supelec.fr>
parents:
31493
diff
changeset
|
738 |
7399 | 739 endfunction |
6768 | 740 |
7017 | 741 |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
742 function res = cigrad_ub_lb (x, bgrad, globals) |
6768 | 743 |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
744 cigradfcn = @(x) fd_ci_jac (x, globals.cifcn); |
6768 | 745 |
14595
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
746 if (iscell (globals.cif) && length (globals.cif) > 1) |
11a9d448fdc3
Get rid of global variables in sqp.
Olaf Till <i7tiol@t-online.de>
parents:
14552
diff
changeset
|
747 cigradfcn = globals.cif{2}; |
6768 | 748 endif |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
749 |
11347
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
750 if (isempty (cigradfcn)) |
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
751 res = bgrad; |
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
752 else |
2726132f77f6
sqp.m: Remove never violated Inf bounds from computation (bug #31742)
Rik <octave@nomad.inbox5.com>
parents:
11290
diff
changeset
|
753 res = [feval(cigradfcn,x); bgrad]; |
6768 | 754 endif |
755 | |
7399 | 756 endfunction |
7361 | 757 |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
758 ## Utility function used to debug sqp |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
759 function report (iter, qp_iter, alpha, nfev, obj) |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
760 |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
761 if (nargin == 0) |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
762 printf (" Itn ItQP Step Nfev Objective\n"); |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
763 else |
30893
e1788b1a315f
maint: Use "fcn" as preferred abbreviation for "function" in m-files.
Rik <rik@octave.org>
parents:
30564
diff
changeset
|
764 printf ("%5d %4d %8.1g %5d %13.6e\n", iter, qp_iter, alpha, nfev, obj); |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
765 endif |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
766 |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
767 endfunction |
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
768 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
769 |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
770 ################################################################################ |
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
771 ## Test Code |
11348
2ae0ca4ee36b
sqp.m: Change docstring to refer to x0 as the initial seed vector
Rik <octave@nomad.inbox5.com>
parents:
11347
diff
changeset
|
772 |
13305
63463570d9fe
Add %!endfunction block keyword to test.m
Rik <octave@nomad.inbox5.com>
parents:
13268
diff
changeset
|
773 %!function r = __g (x) |
7371 | 774 %! r = [sumsq(x)-10; |
775 %! x(2)*x(3)-5*x(4)*x(5); | |
776 %! x(1)^3+x(2)^3+1 ]; | |
13305
63463570d9fe
Add %!endfunction block keyword to test.m
Rik <octave@nomad.inbox5.com>
parents:
13268
diff
changeset
|
777 %!endfunction |
7361 | 778 %! |
13305
63463570d9fe
Add %!endfunction block keyword to test.m
Rik <octave@nomad.inbox5.com>
parents:
13268
diff
changeset
|
779 %!function obj = __phi (x) |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
780 %! obj = exp (prod (x)) - 0.5*(x(1)^3 + x(2)^3 + 1)^2; |
13305
63463570d9fe
Add %!endfunction block keyword to test.m
Rik <octave@nomad.inbox5.com>
parents:
13268
diff
changeset
|
781 %!endfunction |
7361 | 782 %! |
783 %!test | |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
784 %! |
7361 | 785 %! x0 = [-1.8; 1.7; 1.9; -0.8; -0.8]; |
786 %! | |
13305
63463570d9fe
Add %!endfunction block keyword to test.m
Rik <octave@nomad.inbox5.com>
parents:
13268
diff
changeset
|
787 %! [x, obj, info, iter, nf, lambda] = sqp (x0, @__phi, @__g, []); |
7361 | 788 %! |
789 %! x_opt = [-1.717143501952599; | |
790 %! 1.595709610928535; | |
791 %! 1.827245880097156; | |
792 %! -0.763643103133572; | |
793 %! -0.763643068453300]; | |
794 %! | |
7371 | 795 %! obj_opt = 0.0539498477702739; |
7361 | 796 %! |
31493
c05ef94a2bbc
test: Increase tolerance in test for sqp for higher deviations on macOS.
Markus Mützel <markus.muetzel@gmx.de>
parents:
30893
diff
changeset
|
797 %! assert (x, x_opt, 10 * sqrt (eps)); |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
798 %! assert (obj, obj_opt, sqrt (eps)); |
10678
35338deff753
Guarantee equivalent results if sqp called with or wihout bounds
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
799 |
19833
9fc020886ae9
maint: Clean up m-files to follow Octave coding conventions.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
800 ## Test input validation |
28896
90fea9cc9caa
test: Add expected error message <Invalid call> to BIST tests for nargin.
Rik <rik@octave.org>
parents:
28892
diff
changeset
|
801 %!error <Invalid call> sqp () |
90fea9cc9caa
test: Add expected error message <Invalid call> to BIST tests for nargin.
Rik <rik@octave.org>
parents:
28892
diff
changeset
|
802 %!error <Invalid call> sqp (1) |
90fea9cc9caa
test: Add expected error message <Invalid call> to BIST tests for nargin.
Rik <rik@octave.org>
parents:
28892
diff
changeset
|
803 %!error <Invalid call> sqp (1,2,3,4,5) |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
804 %!error sqp (ones (2,2)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
805 %!error sqp (1, cell (4,1)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
806 %!error sqp (1, cell (3,1), cell (3,1)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
807 %!error sqp (1, cell (3,1), cell (2,1), cell (3,1)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
808 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1), ones (2,2),[]) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
809 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),[], ones (2,2)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
810 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),1,-1) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
811 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),[],[], ones (2,2)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
812 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),[],[],-1) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
813 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),[],[],1.5) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
814 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),[],[],[], ones (2,2)) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14327
diff
changeset
|
815 %!error sqp (1, cell (3,1), cell (2,1), cell (2,1),[],[],[],-1) |