octave: liboctave/numeric/aepbalance.cc annotate

annotate liboctave/numeric/aepbalance.cc @ 27933:863ae57eee69

maint: Use Octave coding conventions in liboctave/ * Array-C.cc, Array-d.cc, Array-f.cc, MSparse.cc, Sparse.cc, dim-vector.h, xerbla.cc, aepbalance.cc, eigs-base.cc, gepbalance.cc, oct-fftw.cc, randmtzig.cc, mx-inlines.cc, lo-sysdep.cc, base-list.h, cmd-edit.h, lo-regexp.cc, oct-atomic.h, oct-binmap.h, oct-inttypes.cc, oct-inttypes.h, quit.h, url-transfer.cc: Use Octave coding conventions in liboctave.

author	Rik <rik@octave.org>
date	Sat, 11 Jan 2020 12:53:20 -0800
parents	bd51beb6205e
children	5c14f81e0937 0a5b15007766

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

Mercurial > octave

annotate liboctave/numeric/aepbalance.cc @ 27933:863ae57eee69