octave: libinterp/corefcn/balance.cc annotate

annotate libinterp/corefcn/balance.cc @ 21024:6176560b03d9

maint: Remove return stmts after gripe function calls. * gripes.h: Mark gripe functions which call error() wwith GCC_ATTR_NORETURN. * jit-typeinfo.cc (gripe_bad_result): Mark function with GCC_ATTR_NORETURN. * load-save.cc (gripe_file_open): Mark function with GCC_ATTR_NORETURN. * ov-base.cc (gripe_indexed_assignment, gripe_assign_conversion_failed, gripe_no_conversion): Mark functions with GCC_ATTR_NORETURN. * ov-class.cc (gripe_invalid_index1, gripe_invalid_index_for_assignment, gripe_invalid_index_type): Mark functions with GCC_ATTR_NORETURN. * ov-classdef.h (gripe_invalid_object): Mark function with GCC_ATTR_NORETURN. * ov.cc (gripe_cat_op, gripe_cat_op_conv): Mark functions with GCC_ATTR_NORETURN. * balance.cc, det.cc, eig.cc, hess.cc, inv.cc, jit-typeinfo.cc, load-save.cc, qz.cc, schur.cc, sparse-xdiv.cc, sparse-xpow.cc, sylvester.cc, xdiv.cc, xpow.cc, symrcm.cc, ov-base.cc, ov-classdef.h, ov-cs-list.cc, ov-flt-re-mat.cc, ov-re-mat.cc, ov-re-sparse.cc, ov.cc, op-int.h, op-fs-fs.cc, ops.h: Remove statements that occur after a gripe function which does not return. * op-cs-cs.cc, op-dms-template.cc, op-fcs-fcs.cc, op-s-s.cc : Add newline between gripe and return statement for readability.

author	Rik <rik@octave.org>
date	Fri, 01 Jan 2016 13:58:45 -0800
parents	48b2ad5ee801
children	e39e05d90788

rev	line source
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	1 /*
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	2
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 18812 diff changeset	3 Copyright (C) 1996-2015 John W. Eaton
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 10846 diff changeset	4 Copyright (C) 2008-2009 Jaroslav Hajek
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	5
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	7
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5823 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5823 diff changeset	11 option) any later version.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	12
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	16 for more details.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	17
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	18 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: 5823 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5823 diff changeset	20 <http://www.gnu.org/licenses/>.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	21
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	22 */
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	23
3911 8389e78e67d4 [project @ 2002-04-28 02:15:38 by jwe] jwe parents: 3887 diff changeset	24 // Author: A. S. Hodel <scotte@eng.auburn.edu>
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	25
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	26 #ifdef HAVE_CONFIG_H
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	27 #include <config.h>
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	28 #endif
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	29
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	30 #include <string>
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	31
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	32 #include "CmplxAEPBAL.h"
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	33 #include "fCmplxAEPBAL.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	34 #include "dbleAEPBAL.h"
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	35 #include "floatAEPBAL.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	36 #include "CmplxGEPBAL.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	37 #include "fCmplxGEPBAL.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	38 #include "dbleGEPBAL.h"
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	39 #include "floatGEPBAL.h"
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3911 diff changeset	40 #include "quit.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	41
15039 e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory Rik <rik@octave.org> parents: 14846 diff changeset	42 #include "defun.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	43 #include "error.h"
3181 3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	44 #include "f77-fcn.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	45 #include "gripes.h"
20940 48b2ad5ee801 maint: Rename oct-obj.[cc\|h] to ovl.[cc\|h] for clarity. Rik <rik@octave.org> parents: 20918 diff changeset	46 #include "ovl.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	47 #include "utils.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	48
15039 e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory Rik <rik@octave.org> parents: 14846 diff changeset	49 DEFUN (balance, args, nargout,
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	50 "-- texinfo --\n\
20853 1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20802 diff changeset	51 @deftypefn {} {@var{AA} =} balance (@var{A})\n\
1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20802 diff changeset	52 @deftypefnx {} {@var{AA} =} balance (@var{A}, @var{opt})\n\
1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20802 diff changeset	53 @deftypefnx {} {[@var{DD}, @var{AA}] =} balance (@var{A}, @var{opt})\n\
1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20802 diff changeset	54 @deftypefnx {} {[@var{D}, @var{P}, @var{AA}] =} balance (@var{A}, @var{opt})\n\
1142cf6abc0d 2015 Code Sprint: remove class of function from docstring for all C++ files. Rik <rik@octave.org> parents: 20802 diff changeset	55 @deftypefnx {} {[@var{CC}, @var{DD}, @var{AA}, @var{BB}] =} balance (@var{A}, @var{B}, @var{opt})\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	56 \n\
20172 4f45eaf83908 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 19697 diff changeset	57 Balance the matrix @var{A} to reduce numerical errors in future\n\
4f45eaf83908 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 19697 diff changeset	58 calculations.\n\
4f45eaf83908 doc: Update more docstrings to have one sentence summary as first line. Rik <rik@octave.org> parents: 19697 diff changeset	59 \n\
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	60 Compute @code{@var{AA} = @var{DD} \\ @var{A} * @var{DD}} in which @var{AA}\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	61 is a matrix whose row and column norms are roughly equal in magnitude, and\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	62 @code{@var{DD} = @var{P} * @var{D}}, in which @var{P} is a permutation\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	63 matrix and @var{D} is a diagonal matrix of powers of two. This allows the\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	64 equilibration to be computed without round-off. Results of eigenvalue\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	65 calculation are typically improved by balancing first.\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	66 \n\
12642 f96b9b9f141b doc: Periodic grammarcheck and spellcheck of documentation. Rik <octave@nomad.inbox5.com> parents: 12584 diff changeset	67 If two output values are requested, @code{balance} returns\n\
12584 7ef7e20057fa Improve documentation strings in Linear Algebra chapter. Rik <octave@nomad.inbox5.com> parents: 11586 diff changeset	68 the diagonal @var{D} and the permutation @var{P} separately as vectors.\n\
7ef7e20057fa Improve documentation strings in Linear Algebra chapter. Rik <octave@nomad.inbox5.com> parents: 11586 diff changeset	69 In this case, @code{@var{DD} = eye(n)(:,@var{P}) * diag (@var{D})}, where\n\
12642 f96b9b9f141b doc: Periodic grammarcheck and spellcheck of documentation. Rik <octave@nomad.inbox5.com> parents: 12584 diff changeset	70 @math{n} is the matrix size.\n\
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	71 \n\
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	72 If four output values are requested, compute @code{@var{AA} =\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	73 @var{CC}@var{A}@var{DD}} and @code{@var{BB} = @var{CC}@var{B}@var{DD}},\n\
18812 9ac2357f19bc doc: Replace "non-zero" with "nonzero" to match existing usage. Rik <rik@octave.org> parents: 18100 diff changeset	74 in which @var{AA} and @var{BB} have nonzero elements of approximately the\n\
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	75 same magnitude and @var{CC} and @var{DD} are permuted diagonal matrices as\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	76 in @var{DD} for the algebraic eigenvalue problem.\n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	77 \n\
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	78 The eigenvalue balancing option @var{opt} may be one of:\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	79 \n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	80 @table @asis\n\
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 15195 diff changeset	81 @item @qcode{\"noperm\"}, @qcode{\"S\"}\n\
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	82 Scale only; do not permute.\n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	83 \n\
17281 bc924baa2c4e doc: Add new @qcode macro for code samples which are quoted. Rik <rik@octave.org> parents: 15195 diff changeset	84 @item @qcode{\"noscal\"}, @qcode{\"P\"}\n\
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	85 Permute only; do not scale.\n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	86 @end table\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	87 \n\
9209 923c7cb7f13f Simplify TeXinfo files by eliminating redundant @iftex followed by @tex construction. Rik <rdrider0-list@yahoo.com> parents: 9064 diff changeset	88 Algebraic eigenvalue balancing uses standard @sc{lapack} routines.\n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	89 \n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	90 Generalized eigenvalue problem balancing uses Ward's algorithm\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	91 (SIAM Journal on Scientific and Statistical Computing, 1981).\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3221 diff changeset	92 @end deftypefn")
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	93 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	94 int nargin = args.length ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	95
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	96 if (nargin < 1 \|\| nargin > 3 \|\| nargout < 0)
20802 8bb38ba1bad6 eliminate return statements after calls to print_usage John W. Eaton <jwe@octave.org> parents: 20700 diff changeset	97 print_usage ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	98
20918 6f0bd96f93c0 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20909 diff changeset	99 octave_value_list retval;
6f0bd96f93c0 maint: Use new C++ archetype in more files. Rik <rik@octave.org> parents: 20909 diff changeset	100
3181 3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	101 // determine if it's AEP or GEP
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	102 bool AEPcase = nargin == 1 \|\| args(1).is_string ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	103
3181 3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	104 // problem dimension
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4628 diff changeset	105 octave_idx_type nn = args(0).rows ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	106
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	107 if (nn != args(0).columns ())
21024 6176560b03d9 maint: Remove return stmts after gripe function calls. Rik <rik@octave.org> parents: 20940 diff changeset	108 gripe_square_matrix_required ("balance");
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	109
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	110 bool isfloat = args(0).is_single_type ()
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	111 \|\| (! AEPcase && args(1).is_single_type ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	112
17787 175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	113 bool complex_case = args(0).is_complex_type ()
175b392e91fe Use GNU style coding conventions for code in libinterp/ Rik <rik@octave.org> parents: 17744 diff changeset	114 \|\| (! AEPcase && args(1).is_complex_type ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	115
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	116 // Extract argument 1 parameter for both AEP and GEP.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	117 Matrix aa;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	118 ComplexMatrix caa;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	119 FloatMatrix faa;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	120 FloatComplexMatrix fcaa;
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	121
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	122 if (isfloat)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	123 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	124 if (complex_case)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	125 fcaa = args(0).float_complex_matrix_value ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	126 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	127 faa = args(0).float_matrix_value ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	128 }
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	129 else
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	130 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	131 if (complex_case)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	132 caa = args(0).complex_matrix_value ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	133 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	134 aa = args(0).matrix_value ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	135 }
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	136
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	137 // Treat AEP/GEP cases.
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	138 if (AEPcase)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11553 diff changeset	139 {
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	140 // Algebraic eigenvalue problem.
18100 6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	141 bool noperm = false;
6a71e5030df5 Follow coding convention of defining and initializing only 1 variable per line in liboctinterp. Rik <rik@octave.org> parents: 17787 diff changeset	142 bool noscal = false;
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	143 if (nargin > 1)
a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	144 {
a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	145 std::string a1s = args(1).string_value ();
a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	146 noperm = a1s == "noperm" \|\| a1s == "S";
a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	147 noscal = a1s == "noscal" \|\| a1s == "P";
a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	148 }
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	149
9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	150 // balance the AEP
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	151 if (isfloat)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	152 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	153 if (complex_case)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	154 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	155 FloatComplexAEPBALANCE result (fcaa, noperm, noscal);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	156
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	157 if (nargout == 0 \|\| nargout == 1)
20884 f1b2a2dbc0e1 2015 Code Sprint: use ovl () in C++ files. José Luis García Pallero <jgpallero@gmail.com> parents: 20853 diff changeset	158 retval = ovl (result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	159 else if (nargout == 2)
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	160 retval = ovl (result.balancing_matrix (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	161 result.balanced_matrix ());
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	162 else
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	163 retval = ovl (result.scaling_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	164 result.permuting_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	165 result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	166 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	167 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	168 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	169 FloatAEPBALANCE result (faa, noperm, noscal);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	170
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	171 if (nargout == 0 \|\| nargout == 1)
20884 f1b2a2dbc0e1 2015 Code Sprint: use ovl () in C++ files. José Luis García Pallero <jgpallero@gmail.com> parents: 20853 diff changeset	172 retval = ovl (result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	173 else if (nargout == 2)
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	174 retval = ovl (result.balancing_matrix (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	175 result.balanced_matrix ());
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	176 else
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	177 retval = ovl (result.scaling_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	178 result.permuting_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	179 result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	180 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	181 }
3181 3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	182 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	183 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	184 if (complex_case)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	185 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	186 ComplexAEPBALANCE result (caa, noperm, noscal);
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	187
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	188 if (nargout == 0 \|\| nargout == 1)
20884 f1b2a2dbc0e1 2015 Code Sprint: use ovl () in C++ files. José Luis García Pallero <jgpallero@gmail.com> parents: 20853 diff changeset	189 retval = ovl (result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	190 else if (nargout == 2)
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	191 retval = ovl (result.balancing_matrix (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	192 result.balanced_matrix ());
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	193 else
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	194 retval = ovl (result.scaling_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	195 result.permuting_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	196 result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	197 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	198 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	199 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	200 AEPBALANCE result (aa, noperm, noscal);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	201
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	202 if (nargout == 0 \|\| nargout == 1)
20884 f1b2a2dbc0e1 2015 Code Sprint: use ovl () in C++ files. José Luis García Pallero <jgpallero@gmail.com> parents: 20853 diff changeset	203 retval = ovl (result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	204 else if (nargout == 2)
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	205 retval = ovl (result.balancing_matrix (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	206 result.balanced_matrix ());
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	207 else
20909 03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	208 retval = ovl (result.scaling_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	209 result.permuting_vector (),
03e4ddd49396 omit unnecessary nargout checks John W. Eaton <jwe@octave.org> parents: 20884 diff changeset	210 result.balanced_matrix ());
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	211 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	212 }
3181 3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	213 }
3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	214 else
3f6a813eb09e [project @ 1998-09-25 02:50:29 by jwe] jwe parents: 3179 diff changeset	215 {
8386 a5e080076778 make balance more Matlab-compatible Jaroslav Hajek <highegg@gmail.com> parents: 7789 diff changeset	216 std::string bal_job;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	217 if (nargout == 1)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	218 warning ("balance: used GEP, should have two output arguments");
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	219
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	220 // Generalized eigenvalue problem.
9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	221 if (nargin == 2)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	222 bal_job = "B";
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	223 else
20700 68e3a747ca02 rename octave_value value extractors that accept error message args John W. Eaton <jwe@octave.org> parents: 20582 diff changeset	224 bal_job = args(2).xstring_value ("balance: OPT argument must be a string");
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	225
9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	226 if ((nn != args(1).columns ()) \|\| (nn != args(1).rows ()))
21024 6176560b03d9 maint: Remove return stmts after gripe function calls. Rik <rik@octave.org> parents: 20940 diff changeset	227 gripe_nonconformant ();
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	228
9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	229 Matrix bb;
9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	230 ComplexMatrix cbb;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	231 FloatMatrix fbb;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	232 FloatComplexMatrix fcbb;
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	233
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	234 if (isfloat)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	235 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	236 if (complex_case)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	237 fcbb = args(1).float_complex_matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	238 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	239 fbb = args(1).float_matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	240 }
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	241 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	242 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	243 if (complex_case)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	244 cbb = args(1).complex_matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	245 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	246 bb = args(1).matrix_value ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	247 }
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3911 diff changeset	248
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	249 // balance the GEP
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	250 if (isfloat)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	251 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	252 if (complex_case)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	253 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	254 FloatComplexGEPBALANCE result (fcaa, fcbb, bal_job);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	255
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	256 switch (nargout)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	257 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	258 case 4:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	259 retval(3) = result.balanced_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	260 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	261 case 3:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	262 retval(2) = result.balanced_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	263 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	264 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	265 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	266 case 2:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	267 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	268 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	269 case 1:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	270 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	271 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	272 default:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	273 error ("balance: invalid number of output arguments");
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	274 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	275 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	276 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	277 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	278 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	279 FloatGEPBALANCE result (faa, fbb, bal_job);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	280
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	281 switch (nargout)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	282 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	283 case 4:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	284 retval(3) = result.balanced_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	285 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	286 case 3:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	287 retval(2) = result.balanced_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	288 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	289 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	290 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	291 case 2:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	292 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	293 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	294 case 1:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	295 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	296 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	297 default:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	298 error ("balance: invalid number of output arguments");
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	299 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	300 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	301 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	302 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	303 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	304 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	305 if (complex_case)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	306 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	307 ComplexGEPBALANCE result (caa, cbb, bal_job);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	308
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	309 switch (nargout)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	310 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	311 case 4:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	312 retval(3) = result.balanced_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	313 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	314 case 3:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	315 retval(2) = result.balanced_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	316 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	317 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	318 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	319 case 2:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	320 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	321 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	322 case 1:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	323 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	324 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	325 default:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	326 error ("balance: invalid number of output arguments");
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	327 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	328 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	329 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	330 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	331 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	332 GEPBALANCE result (aa, bb, bal_job);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7482 diff changeset	333
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	334 switch (nargout)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	335 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	336 case 4:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	337 retval(3) = result.balanced_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	338 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	339 case 3:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	340 retval(2) = result.balanced_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	341 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	342 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	343 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	344 case 2:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	345 retval(1) = result.balancing_matrix2 ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	346 // fall through
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	347 case 1:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	348 retval(0) = result.balancing_matrix ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	349 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	350 default:
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	351 error ("balance: invalid number of output arguments");
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	352 break;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	353 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	354 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9821 diff changeset	355 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	356 }
3185 9580887dd160 [project @ 1998-09-26 02:45:55 by jwe] jwe parents: 3181 diff changeset	357
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	358 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	359 }

Mercurial > octave

annotate libinterp/corefcn/balance.cc @ 21024:6176560b03d9