Mercurial > octave-nkf
annotate libinterp/corefcn/hess.cc @ 19074:0850b5212619
doc: Add @nospell macro around proper names in documentation.
* aspell-octave.en.pws: Remove proper names from Octave's custom dictionary.
* doc/interpreter/audio.txi, doc/interpreter/diagperm.txi,
doc/interpreter/diffeq.txi, doc/interpreter/emacs.txi,
doc/interpreter/func.txi, doc/interpreter/geometry.txi,
doc/interpreter/intro.txi, doc/interpreter/linalg.txi,
doc/interpreter/numbers.txi, doc/interpreter/preface.txi,
doc/interpreter/quad.txi, doc/interpreter/sparse.txi,
libinterp/corefcn/colloc.cc, libinterp/corefcn/ellipj.cc,
libinterp/corefcn/gammainc.cc, libinterp/corefcn/hess.cc,
libinterp/corefcn/quadcc.cc, libinterp/corefcn/rand.cc,
libinterp/corefcn/sqrtm.cc, libinterp/dldfcn/ccolamd.cc,
libinterp/dldfcn/colamd.cc, libinterp/dldfcn/dmperm.cc,
libinterp/dldfcn/symrcm.cc, scripts/general/quadgk.m, scripts/general/quadl.m,
scripts/linear-algebra/commutation_matrix.m, scripts/linear-algebra/condest.m,
scripts/linear-algebra/duplication_matrix.m, scripts/linear-algebra/expm.m,
scripts/linear-algebra/housh.m, scripts/linear-algebra/krylov.m,
scripts/linear-algebra/onenormest.m, scripts/linear-algebra/qzhess.m,
scripts/optimization/fminsearch.m, scripts/optimization/glpk.m,
scripts/signal/arch_fit.m, scripts/signal/bartlett.m,
scripts/signal/blackman.m, scripts/signal/diffpara.m,
scripts/signal/durbinlevinson.m, scripts/signal/hamming.m,
scripts/signal/hanning.m, scripts/sparse/eigs.m, scripts/sparse/pcg.m,
scripts/sparse/pcr.m, scripts/specfun/ellipke.m,
scripts/special-matrix/gallery.m, scripts/special-matrix/invhilb.m,
scripts/statistics/base/quantile.m,
scripts/statistics/tests/kruskal_wallis_test.m,
scripts/statistics/tests/manova.m, scripts/statistics/tests/mcnemar_test.m:
Add @nospell macro around proper names in documentation.
author | Rik <rik@octave.org> |
---|---|
date | Thu, 21 Aug 2014 15:57:18 -0700 |
parents | 175b392e91fe |
children | 4197fc428c7d |
rev | line source |
---|---|
2928 | 1 /* |
2 | |
17744
d63878346099
maint: Update copyright notices for release.
John W. Eaton <jwe@octave.org>
parents:
16920
diff
changeset
|
3 Copyright (C) 1996-2013 John W. Eaton |
2928 | 4 |
5 This file is part of Octave. | |
6 | |
7 Octave is free software; you can redistribute it and/or modify it | |
8 under the terms of the GNU General Public License as published by the | |
7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
10 option) any later version. | |
2928 | 11 |
12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
15 for more details. | |
16 | |
17 You should have received a copy of the GNU General Public License | |
7016 | 18 along with Octave; see the file COPYING. If not, see |
19 <http://www.gnu.org/licenses/>. | |
2928 | 20 |
21 */ | |
22 | |
23 #ifdef HAVE_CONFIG_H | |
24 #include <config.h> | |
25 #endif | |
26 | |
27 #include "CmplxHESS.h" | |
28 #include "dbleHESS.h" | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
29 #include "fCmplxHESS.h" |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
30 #include "floatHESS.h" |
2928 | 31 |
15039
e753177cde93
maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents:
14501
diff
changeset
|
32 #include "defun.h" |
2928 | 33 #include "error.h" |
34 #include "gripes.h" | |
35 #include "oct-obj.h" | |
36 #include "utils.h" | |
37 | |
15039
e753177cde93
maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents:
14501
diff
changeset
|
38 DEFUN (hess, args, nargout, |
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
39 "-*- texinfo -*-\n\ |
15039
e753177cde93
maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents:
14501
diff
changeset
|
40 @deftypefn {Built-in Function} {@var{H} =} hess (@var{A})\n\ |
e753177cde93
maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents:
14501
diff
changeset
|
41 @deftypefnx {Built-in Function} {[@var{P}, @var{H}] =} hess (@var{A})\n\ |
3372 | 42 @cindex Hessenberg decomposition\n\ |
11553
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
43 Compute the Hessenberg decomposition of the matrix @var{A}.\n\ |
3372 | 44 \n\ |
11553
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
45 The Hessenberg decomposition is\n\ |
3372 | 46 @tex\n\ |
47 $$\n\ | |
48 A = PHP^T\n\ | |
49 $$\n\ | |
11553
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
50 where $P$ is a square unitary matrix ($P^TP = I$), and $H$\n\ |
3372 | 51 is upper Hessenberg ($H_{i,j} = 0, \\forall i \\ge j+1$).\n\ |
52 @end tex\n\ | |
8517
81d6ab3ac93c
Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents:
7814
diff
changeset
|
53 @ifnottex\n\ |
12584
7ef7e20057fa
Improve documentation strings in Linear Algebra chapter.
Rik <octave@nomad.inbox5.com>
parents:
11553
diff
changeset
|
54 @code{@var{P} * @var{H} * @var{P}' = @var{A}} where @var{P} is a square\n\ |
7ef7e20057fa
Improve documentation strings in Linear Algebra chapter.
Rik <octave@nomad.inbox5.com>
parents:
11553
diff
changeset
|
55 unitary matrix (@code{@var{P}' * @var{P} = I}, using complex-conjugate\n\ |
11553
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
56 transposition) and @var{H} is upper Hessenberg\n\ |
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
57 (@code{@var{H}(i, j) = 0 forall i >= j+1)}.\n\ |
8517
81d6ab3ac93c
Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents:
7814
diff
changeset
|
58 @end ifnottex\n\ |
11553
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
59 \n\ |
01f703952eff
Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents:
11523
diff
changeset
|
60 The Hessenberg decomposition is usually used as the first step in an\n\ |
19074
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17787
diff
changeset
|
61 eigenvalue computation, but has other applications as well\n\ |
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17787
diff
changeset
|
62 (see @nospell{Golub, Nash, and Van Loan},\n\ |
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17787
diff
changeset
|
63 IEEE Transactions on Automatic Control, 1979).\n\ |
16920
53eaa83e4181
doc: Add seealso links between various factorization forms.
Rik <rik@octave.org>
parents:
15195
diff
changeset
|
64 @seealso{eig, chol, lu, qr, qz, schur, svd}\n\ |
3372 | 65 @end deftypefn") |
2928 | 66 { |
67 octave_value_list retval; | |
68 | |
69 int nargin = args.length (); | |
70 | |
71 if (nargin != 1 || nargout > 2) | |
72 { | |
5823 | 73 print_usage (); |
2928 | 74 return retval; |
75 } | |
76 | |
77 octave_value arg = args(0); | |
78 | |
5275 | 79 octave_idx_type nr = arg.rows (); |
80 octave_idx_type nc = arg.columns (); | |
2928 | 81 |
82 int arg_is_empty = empty_arg ("hess", nr, nc); | |
83 | |
84 if (arg_is_empty < 0) | |
85 return retval; | |
86 else if (arg_is_empty > 0) | |
87 return octave_value_list (2, Matrix ()); | |
88 | |
89 if (nr != nc) | |
90 { | |
91 gripe_square_matrix_required ("hess"); | |
92 return retval; | |
93 } | |
94 | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
95 if (arg.is_single_type ()) |
2928 | 96 { |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
97 if (arg.is_real_type ()) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
98 { |
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
99 FloatMatrix tmp = arg.float_matrix_value (); |
2928 | 100 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
101 if (! error_state) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
102 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
103 FloatHESS result (tmp); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
104 |
10539
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
105 if (nargout <= 1) |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
106 retval(0) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
107 else |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
108 { |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
109 retval(1) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
110 retval(0) = result.unitary_hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
111 } |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
112 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
113 } |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
114 else if (arg.is_complex_type ()) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
115 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
116 FloatComplexMatrix ctmp = arg.float_complex_matrix_value (); |
2928 | 117 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
118 if (! error_state) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
119 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
120 FloatComplexHESS result (ctmp); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
121 |
10539
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
122 if (nargout <= 1) |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
123 retval(0) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
124 else |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
125 { |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
126 retval(1) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
127 retval(0) = result.unitary_hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
128 } |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
129 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
130 } |
2928 | 131 } |
132 else | |
133 { | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
134 if (arg.is_real_type ()) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
135 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
136 Matrix tmp = arg.matrix_value (); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
137 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
138 if (! error_state) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
139 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
140 HESS result (tmp); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
141 |
10539
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
142 if (nargout <= 1) |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
143 retval(0) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
144 else |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
145 { |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
146 retval(1) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
147 retval(0) = result.unitary_hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
148 } |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
149 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
150 } |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
151 else if (arg.is_complex_type ()) |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
152 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
153 ComplexMatrix ctmp = arg.complex_matrix_value (); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
154 |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
155 if (! error_state) |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
156 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
157 ComplexHESS result (ctmp); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
158 |
10539
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
159 if (nargout <= 1) |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
160 retval(0) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
161 else |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
162 { |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
163 retval(1) = result.hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
164 retval(0) = result.unitary_hess_matrix (); |
8db807a0eb78
Return correct Hessenberg matrix when called with 0 or 1 output arguments
Rik <code@nomad.inbox5.com>
parents:
10155
diff
changeset
|
165 } |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
166 } |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
167 } |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
168 else |
10154
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
169 { |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
170 gripe_wrong_type_arg ("hess", arg); |
40dfc0c99116
DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents:
9758
diff
changeset
|
171 } |
2928 | 172 } |
173 | |
174 return retval; | |
175 } | |
176 | |
177 /* | |
7814
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
178 %!test |
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
179 %! a = [1, 2, 3; 5, 4, 6; 8, 7, 9]; |
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
180 %! [p, h] = hess (a); |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
181 %! assert (p * h * p', a, sqrt (eps)); |
7814
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
182 |
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
183 %!test |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
184 %! a = single ([1, 2, 3; 5, 4, 6; 8, 7, 9]); |
7814
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
185 %! [p, h] = hess (a); |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
186 %! assert (p * h * p', a, sqrt (eps ("single"))); |
7814
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
187 |
14501
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
188 %!error hess () |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
189 %!error hess ([1, 2; 3, 4], 2) |
60e5cf354d80
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
190 %!error <argument must be a square matrix> hess ([1, 2; 3, 4; 5, 6]) |
7814
87865ed7405f
Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents:
7789
diff
changeset
|
191 */ |