octave: liboctave/numeric/aepbalance.cc annotate

annotate liboctave/numeric/aepbalance.cc @ 30564:796f54d4ddbf stable

update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2021. In all .txi and .texi files except gpl.txi and gpl.texi in the doc/liboctave and doc/interpreter directories, change the copyright to "Octave Project Developers", the same as used for other source files. Update copyright notices for 2022 (not done since 2019). For gpl.txi and gpl.texi, change the copyright notice to be "Free Software Foundation, Inc." and leave the date at 2007 only because this file only contains the text of the GPL, not anything created by the Octave Project Developers. Add Paul Thomas to contributors.in.

author	John W. Eaton <jwe@octave.org>
date	Tue, 28 Dec 2021 18:22:40 -0500
parents	f3f3e3793fb5
children	e88a07dec498

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
30564 796f54d4ddbf update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 30394 diff changeset	3 // Copyright (C) 1994-2022 The Octave Project Developers
27923 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 <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	51 OCTAVE_API
22329 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)
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	53 : m_balanced_mat (a), m_scale (), m_ilo (), m_ihi (),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	54 m_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
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	62 m_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
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	66 F77_XFCN (dgebal, DGEBAL, (F77_CONST_CHAR_ARG2 (&m_job, 1), n,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	67 m_balanced_mat.fortran_vec (), n,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	68 t_ilo, t_ihi, m_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
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	71 m_ilo = t_ilo;
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	72 m_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 <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	76 OCTAVE_API Matrix
22329 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 {
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	79 F77_INT n = to_f77_int (m_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++)
30394 f3f3e3793fb5 maint: style check C++ files in liboctave/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30078 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;
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	86 F77_INT t_ilo = to_f77_int (m_ilo);
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	87 F77_INT t_ihi = to_f77_int (m_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
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	91 F77_XFCN (dgebak, DGEBAK, (F77_CONST_CHAR_ARG2 (&m_job, 1),
22329 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),
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	93 n, t_ilo, t_ihi, m_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 <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	102 OCTAVE_API
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	103 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	104 bool noscal)
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	105 : m_balanced_mat (a), m_scale (), m_ilo (), m_ihi (),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	106 m_job (get_job (noperm, noscal))
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	107 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	108 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	109
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	110 if (a.rows () != n)
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23220 diff changeset	111 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23220 diff changeset	112 ("aepbalance: requires square matrix");
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	113
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	114 m_scale = FloatColumnVector (n);
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	115
22947 0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	116 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	117
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	118 F77_XFCN (sgebal, SGEBAL, (F77_CONST_CHAR_ARG2 (&m_job, 1), n,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	119 m_balanced_mat.fortran_vec (), n, t_ilo,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	120 t_ihi, m_scale.fortran_vec (), info
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	121 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	122
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	123 m_ilo = t_ilo;
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	124 m_ihi = t_ihi;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	125 }
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	126
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	127 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	128 OCTAVE_API FloatMatrix
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	129 aepbalance<FloatMatrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	130 {
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	131 F77_INT n = to_f77_int (m_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	132
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	133 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	134 for (F77_INT i = 0; i < n; i++)
30394 f3f3e3793fb5 maint: style check C++ files in liboctave/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30078 diff changeset	135 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	136
22947 0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	137 F77_INT info;
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	138 F77_INT t_ilo = to_f77_int (m_ilo);
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	139 F77_INT t_ihi = to_f77_int (m_ihi);
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 char side = 'R';
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	142
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	143 F77_XFCN (sgebak, SGEBAK, (F77_CONST_CHAR_ARG2 (&m_job, 1),
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	144 F77_CONST_CHAR_ARG2 (&side, 1),
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	145 n, t_ilo, t_ihi, m_scale.data (), n,
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	146 balancing_mat.fortran_vec (), n, info
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	147 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	148 F77_CHAR_ARG_LEN (1)));
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	149
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	150 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	151 }
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	152
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	153 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	154 OCTAVE_API
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	155 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	156 bool noscal)
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	157 : m_balanced_mat (a), m_scale (), m_ilo (), m_ihi (),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	158 m_job (get_job (noperm, noscal))
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	159 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	160 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	161
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	162 if (a.rows () != n)
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23220 diff changeset	163 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23220 diff changeset	164 ("aepbalance: requires square matrix");
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	165
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	166 m_scale = ColumnVector (n);
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	167
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_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	169
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	170 F77_XFCN (zgebal, ZGEBAL,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	171 (F77_CONST_CHAR_ARG2 (&m_job, 1), n,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	172 F77_DBLE_CMPLX_ARG (m_balanced_mat.fortran_vec ()),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	173 n, t_ilo, t_ihi, m_scale.fortran_vec (), info
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	174 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	175
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	176 m_ilo = t_ilo;
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	177 m_ihi = t_ihi;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	178 }
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	179
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	180 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	181 OCTAVE_API ComplexMatrix
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	182 aepbalance<ComplexMatrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	183 {
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	184 F77_INT n = to_f77_int (m_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	185
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	186 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	187 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	188 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	189
22947 0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	190 F77_INT info;
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	191 F77_INT t_ilo = to_f77_int (m_ilo);
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	192 F77_INT t_ihi = to_f77_int (m_ihi);
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	193
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	194 char side = 'R';
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	195
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	196 F77_XFCN (zgebak, ZGEBAK,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	197 (F77_CONST_CHAR_ARG2 (&m_job, 1),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	198 F77_CONST_CHAR_ARG2 (&side, 1),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	199 n, t_ilo, t_ihi, m_scale.data (), n,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	200 F77_DBLE_CMPLX_ARG (balancing_mat.fortran_vec ()),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	201 n, info
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	202 F77_CHAR_ARG_LEN (1)
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	203 F77_CHAR_ARG_LEN (1)));
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	204
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	205 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	206 }
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	207
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	208 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	209 OCTAVE_API
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	210 aepbalance<FloatComplexMatrix>::aepbalance (const FloatComplexMatrix& a,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	211 bool noperm, bool noscal)
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	212 : m_balanced_mat (a), m_scale (), m_ilo (), m_ihi (),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	213 m_job (get_job (noperm, noscal))
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	214 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	215 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	216
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	217 if (a.rows () != n)
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23220 diff changeset	218 (*current_liboctave_error_handler)
d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23220 diff changeset	219 ("aepbalance: requires square matrix");
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	220
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	221 m_scale = FloatColumnVector (n);
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	222
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 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	224
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	225 F77_XFCN (cgebal, CGEBAL, (F77_CONST_CHAR_ARG2 (&m_job, 1), n,
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	226 F77_CMPLX_ARG (m_balanced_mat.fortran_vec ()),
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	227 n, t_ilo, t_ihi, m_scale.fortran_vec (), info
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	228 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	229
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	230 m_ilo = t_ilo;
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	231 m_ihi = t_ihi;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	232 }
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	233
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	234 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27933 diff changeset	235 OCTAVE_API FloatComplexMatrix
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	236 aepbalance<FloatComplexMatrix>::balancing_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	237 {
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	238 F77_INT n = to_f77_int (m_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	239
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	240 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	241 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	242 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	243
22947 0b2ef933b738 use F77_INT instead of octave_idx_type for liboctave aepbalance class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	244 F77_INT info;
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	245 F77_INT t_ilo = to_f77_int (m_ilo);
443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	246 F77_INT t_ihi = to_f77_int (m_ihi);
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	247
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	248 char side = 'R';
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	249
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	250 F77_XFCN (cgebak, CGEBAK, (F77_CONST_CHAR_ARG2 (&m_job, 1),
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	251 F77_CONST_CHAR_ARG2 (&side, 1),
30078 443571d367fe maint: use "m_" prefix for member variables in class aepbalance. Rik <rik@octave.org> parents: 29359 diff changeset	252 n, t_ilo, t_ihi, m_scale.data (), 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	253 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	254 n, info
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	255 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	256 F77_CHAR_ARG_LEN (1)));
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 return balancing_mat;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	259 }
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 // Instantiations we need.
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<Matrix>;
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	264
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	265 template class aepbalance<FloatMatrix>;
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	266
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	267 template class aepbalance<ComplexMatrix>;
21268 f08ae27289e4 better use of templates for balance classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	268
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	269 template class aepbalance<FloatComplexMatrix>;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	270 }
22317 6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	271 }

Mercurial > octave

annotate liboctave/numeric/aepbalance.cc @ 30564:796f54d4ddbf stable