Mercurial > octave

annotate liboctave/numeric/aepbalance.cc @ 27919:1891570abac8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2020.
author John W. Eaton <jwe@octave.org>
date Mon, 06 Jan 2020 22:29:51 -0500
parents b442ec6dda5c
children bd51beb6205e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
1 /*
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
2
27919
1891570abac8 update Octave Project Developers copyright for the new year
John W. Eaton <jwe@octave.org>
parents: 27918
diff changeset
3 Copyright (C) 1994-2020 The Octave Project Developers
27918
b442ec6dda5c use centralized file for copyright info for individual contributors
John W. Eaton <jwe@octave.org>
parents: 27101
diff changeset
4
b442ec6dda5c use centralized file for copyright info for individual contributors
John W. Eaton <jwe@octave.org>
parents: 27101
diff changeset
5 See the file COPYRIGHT.md in the top-level directory of this distribution
b442ec6dda5c use centralized file for copyright info for individual contributors
John W. Eaton <jwe@octave.org>
parents: 27101
diff changeset
6 or <https://octave.org/COPYRIGHT.html/>.
b442ec6dda5c use centralized file for copyright info for individual contributors
John W. Eaton <jwe@octave.org>
parents: 27101
diff changeset
7
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
8
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
9 This file is part of Octave.
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
10
24534
194eb4bd202b maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents: 23696
diff changeset
11 Octave is free software: you can redistribute it and/or modify it
22755
3a2b891d0b33 maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents: 22402
diff changeset
12 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: 23696
diff changeset
13 the Free Software Foundation, either version 3 of the License, or
22755
3a2b891d0b33 maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents: 22402
diff changeset
14 (at your option) any later version.
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
15
22755
3a2b891d0b33 maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents: 22402
diff changeset
16 Octave is distributed in the hope that it will be useful, but
3a2b891d0b33 maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents: 22402
diff changeset
17 WITHOUT ANY WARRANTY; without even the implied warranty of
3a2b891d0b33 maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents: 22402
diff changeset
18 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
3a2b891d0b33 maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents: 22402
diff changeset
19 GNU General Public License for more details.
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
20
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
21 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
22 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: 23696
diff changeset
23 <https://www.gnu.org/licenses/>.
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
24
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
25 */
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
26
21724
aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents: 21301
diff changeset
27 #if defined (HAVE_CONFIG_H)
21301
40de9f8f23a6 Use '#include "config.h"' rather than <config.h>.
Rik <rik@octave.org>
parents: 21268
diff changeset
28 # include "config.h"
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
29 #endif
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
30
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
31 #include "CMatrix.h"
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
32 #include "aepbalance.h"
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
33 #include "dColVector.h"
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
34 #include "dMatrix.h"
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
35 #include "fCMatrix.h"
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
36 #include "fColVector.h"
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
37 #include "fMatrix.h"
23475
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
38 #include "lo-error.h"
22322
93b3cdd36854 move most f77 function decls to separate header files
John W. Eaton <jwe@octave.org>
parents: 22317
diff changeset
39 #include "lo-lapack-proto.h"
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
40
22317
6ca3acf5fad8 move some new numeric classes to namespace octave::math
John W. Eaton <jwe@octave.org>
parents: 22135
diff changeset
41 namespace octave
6ca3acf5fad8 move some new numeric classes to namespace octave::math
John W. Eaton <jwe@octave.org>
parents: 22135
diff changeset
42 {
27101
097774bed4ed move some static functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
43 static inline char
097774bed4ed move some static functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
44 get_job (bool noperm, bool noscal)
097774bed4ed move some static functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
45 {
097774bed4ed move some static functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
46 return noperm ? (noscal ? 'N' : 'S') : (noscal ? 'P' : 'B');
097774bed4ed move some static functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
47 }
097774bed4ed move some static functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
48
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
49 namespace math
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
50 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
51 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
52 aepbalance<Matrix>::aepbalance (const Matrix& a, bool noperm, bool noscal)
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
53 : balanced_mat (a), scale (), ilo (), ihi (),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
54 job (get_job (noperm, noscal))
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
55 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
56 F77_INT n = to_f77_int (a.cols ());
22317
6ca3acf5fad8 move some new numeric classes to namespace octave::math
John W. Eaton <jwe@octave.org>
parents: 22135
diff changeset
57
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
58 if (a.rows () != n)
23475
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
59 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
60 ("aepbalance: requires square matrix");
1730
a744f4d0ba59 [project @ 1996-01-09 12:16:29 by jwe]
jwe
parents: 1360
diff changeset
61
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
62 scale = ColumnVector (n);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
63
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
64 F77_INT info, t_ilo, t_ihi;
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
65
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
66 F77_XFCN (dgebal, DGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1), n,
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
67 balanced_mat.fortran_vec (), n,
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
68 t_ilo, t_ihi, scale.fortran_vec (), info
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
69 F77_CHAR_ARG_LEN (1)));
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
70
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
71 ilo = t_ilo;
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
72 ihi = t_ihi;
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
73 }
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
74
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
75 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
76 Matrix
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
77 aepbalance<Matrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
78 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
79 F77_INT n = to_f77_int (balanced_mat.rows ());
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
80
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
81 Matrix balancing_mat (n, n, 0.0);
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
82 for (F77_INT i = 0; i < n; i++)
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
83 balancing_mat.elem (i ,i) = 1.0;
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
84
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
85 F77_INT info;
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
86 F77_INT t_ilo = to_f77_int (ilo);
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
87 F77_INT t_ihi = to_f77_int (ihi);
1933
333923894848 [project @ 1996-02-12 03:07:39 by jwe]
jwe
parents: 1882
diff changeset
88
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
89 char side = 'R';
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
90
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
91 F77_XFCN (dgebak, DGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
92 F77_CONST_CHAR_ARG2 (&side, 1),
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
93 n, t_ilo, t_ihi, scale.data (), n,
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
94 balancing_mat.fortran_vec (), n, info
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
95 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
96 F77_CHAR_ARG_LEN (1)));
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
97
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
98 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
99 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
100
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
101 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
102 aepbalance<FloatMatrix>::aepbalance (const FloatMatrix& a, bool noperm,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
103 bool noscal)
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
104 : balanced_mat (a), scale (), ilo (), ihi (),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
105 job (get_job (noperm, noscal))
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
106 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
107 F77_INT n = to_f77_int (a.cols ());
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
108
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
109 if (a.rows () != n)
23475
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
110 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
111 ("aepbalance: requires square matrix");
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
112
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
113 scale = FloatColumnVector (n);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
114
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
115 F77_INT info, t_ilo, t_ihi;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
116
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
117 F77_XFCN (sgebal, SGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1), n,
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
118 balanced_mat.fortran_vec (), n, t_ilo,
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
119 t_ihi, scale.fortran_vec (), info
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
120 F77_CHAR_ARG_LEN (1)));
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
121
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
122 ilo = t_ilo;
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
123 ihi = t_ihi;
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
124 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
125
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
126 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
127 FloatMatrix
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
128 aepbalance<FloatMatrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
129 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
130 F77_INT n = to_f77_int (balanced_mat.rows ());
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
131
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
132 FloatMatrix balancing_mat (n, n, 0.0);
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
133 for (F77_INT i = 0; i < n; i++)
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
134 balancing_mat.elem (i ,i) = 1.0;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
135
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
136 F77_INT info;
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
137 F77_INT t_ilo = to_f77_int (ilo);
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
138 F77_INT t_ihi = to_f77_int (ihi);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
139
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
140 char side = 'R';
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
141
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
142 F77_XFCN (sgebak, SGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
143 F77_CONST_CHAR_ARG2 (&side, 1),
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
144 n, t_ilo, t_ihi, scale.data (), n,
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
145 balancing_mat.fortran_vec (), n, info
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
146 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
147 F77_CHAR_ARG_LEN (1)));
457
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe]
jwe
parents:
diff changeset
148
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
149 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
150 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
151
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
152 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
153 aepbalance<ComplexMatrix>::aepbalance (const ComplexMatrix& a, bool noperm,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
154 bool noscal)
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
155 : balanced_mat (a), scale (), ilo (), ihi (),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
156 job (get_job (noperm, noscal))
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
157 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
158 F77_INT n = to_f77_int (a.cols ());
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
159
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
160 if (a.rows () != n)
23475
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
161 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
162 ("aepbalance: requires square matrix");
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
163
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
164 scale = ColumnVector (n);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
165
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
166 F77_INT info, t_ilo, t_ihi;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
167
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
168 F77_XFCN (zgebal, ZGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1), n,
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
169 F77_DBLE_CMPLX_ARG (balanced_mat.fortran_vec ()),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
170 n, t_ilo, t_ihi, scale.fortran_vec (), info
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
171 F77_CHAR_ARG_LEN (1)));
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
172
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
173 ilo = t_ilo;
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
174 ihi = t_ihi;
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
175 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
176
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
177 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
178 ComplexMatrix
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
179 aepbalance<ComplexMatrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
180 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
181 F77_INT n = to_f77_int (balanced_mat.rows ());
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
182
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
183 ComplexMatrix balancing_mat (n, n, 0.0);
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
184 for (F77_INT i = 0; i < n; i++)
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
185 balancing_mat.elem (i, i) = 1.0;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
186
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
187 F77_INT info;
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
188 F77_INT t_ilo = to_f77_int (ilo);
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
189 F77_INT t_ihi = to_f77_int (ihi);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
190
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
191 char side = 'R';
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
192
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
193 F77_XFCN (zgebak, ZGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
194 F77_CONST_CHAR_ARG2 (&side, 1),
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
195 n, t_ilo, t_ihi, scale.data (), n,
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
196 F77_DBLE_CMPLX_ARG (balancing_mat.fortran_vec ()),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
197 n, info
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
198 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
199 F77_CHAR_ARG_LEN (1)));
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
200
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
201 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
202 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
203
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
204 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
205 aepbalance<FloatComplexMatrix>::aepbalance (const FloatComplexMatrix& a,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
206 bool noperm, bool noscal)
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
207 : balanced_mat (a), scale (), ilo (), ihi (),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
208 job (get_job (noperm, noscal))
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
209 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
210 F77_INT n = to_f77_int (a.cols ());
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
211
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
212 if (a.rows () != n)
23475
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
213 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents: 23220
diff changeset
214 ("aepbalance: requires square matrix");
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
215
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
216 scale = FloatColumnVector (n);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
217
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
218 F77_INT info, t_ilo, t_ihi;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
219
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
220 F77_XFCN (cgebal, CGEBAL, (F77_CONST_CHAR_ARG2 (&job, 1), n,
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
221 F77_CMPLX_ARG (balanced_mat.fortran_vec ()),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
222 n, t_ilo, t_ihi, scale.fortran_vec (), info
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
223 F77_CHAR_ARG_LEN (1)));
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
224
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
225 ilo = t_ilo;
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
226 ihi = t_ihi;
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
227 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
228
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
229 template <>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
230 FloatComplexMatrix
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
231 aepbalance<FloatComplexMatrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
232 {
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
233 F77_INT n = to_f77_int (balanced_mat.rows ());
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
234
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
235 FloatComplexMatrix balancing_mat (n, n, 0.0);
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
236 for (F77_INT i = 0; i < n; i++)
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
237 balancing_mat.elem (i, i) = 1.0;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
238
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
239 F77_INT info;
23696
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
240 F77_INT t_ilo = to_f77_int (ilo);
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 23475
diff changeset
241 F77_INT t_ihi = to_f77_int (ihi);
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
242
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
243 char side = 'R';
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
244
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
245 F77_XFCN (cgebak, CGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
246 F77_CONST_CHAR_ARG2 (&side, 1),
22947
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
247 n, t_ilo, t_ihi, scale.data (), n,
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
248 F77_CMPLX_ARG (balancing_mat.fortran_vec ()),
0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class
John W. Eaton <jwe@octave.org>
parents: 22755
diff changeset
249 n, info
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
250 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
251 F77_CHAR_ARG_LEN (1)));
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
252
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
253 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
254 }
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
255
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
256 // Instantiations we need.
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
257
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
258 template class aepbalance<Matrix>;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
259
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
260 template class aepbalance<FloatMatrix>;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
261
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
262 template class aepbalance<ComplexMatrix>;
21268
f08ae27289e4 better use of templates for balance classes
John W. Eaton <jwe@octave.org>
parents: 21202
diff changeset
263
22329
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
264 template class aepbalance<FloatComplexMatrix>;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents: 22323
diff changeset
265 }
22317
6ca3acf5fad8 move some new numeric classes to namespace octave::math
John W. Eaton <jwe@octave.org>
parents: 22135
diff changeset
266 }