octave: libinterp/corefcn/ordqz.cc annotate

annotate libinterp/corefcn/ordqz.cc @ 33617:ec2635a02328 bytecode-interpreter tip

maint: Merge default to bytecode-interpreter.

author	Markus Mützel <markus.muetzel@gmx.de>
date	Tue, 21 May 2024 18:29:03 +0200
parents	3606817b9994
children

rev	line source
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	1 ////////////////////////////////////////////////////////////////////////
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	2 //
32632 2e484f9f1f18 maint: update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 31706 diff changeset	3 // Copyright (C) 2020-2024 The Octave Project Developers
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	4 //
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	6 // distribution or <https://octave.org/copyright/>.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	7 //
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	8 // This file is part of Octave.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	9 //
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	11 // under the terms of the GNU General Public License as published by
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	13 // (at your option) any later version.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	14 //
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	15 // Octave is distributed in the hope that it will be useful, but
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	18 // GNU General Public License for more details.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	19 //
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	20 // You should have received a copy of the GNU General Public License
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	21 // along with Octave; see the file COPYING. If not, see
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	22 // <https://www.gnu.org/licenses/>.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	23 //
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	24 ////////////////////////////////////////////////////////////////////////
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	25
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	26 // Generalized eigenvalue reordering via LAPACK
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	27
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	28 // Originally written by M. Koehler <koehlerm(AT)mpi-magdeburg.mpg.de>
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	29
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	30 #undef DEBUG
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	31
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	32 #if defined (HAVE_CONFIG_H)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	33 # include "config.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	34 #endif
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	35
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	36 #include <cctype>
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	37 #include <cmath>
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	38
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	39 #include "f77-fcn.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	40 #include "lo-lapack-proto.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	41 #include "qr.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	42 #include "quit.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	43
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	44 #include "defun.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	45 #include "error.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	46 #include "errwarn.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	47 #include "ovl.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	48
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	49
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	50 #if defined (DEBUG)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	51 # include "pager.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	52 # include "pr-output.h"
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	53 #endif
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	54
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	55
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 30564 diff changeset	56 OCTAVE_BEGIN_NAMESPACE(octave)
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	57
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	58 DEFUN (ordqz, args, nargout,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	59 doc: /* -- texinfo --
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	60 @deftypefn {} {[@var{AR}, @var{BR}, @var{QR}, @var{ZR}] =} ordqz (@var{AA}, @var{BB}, @var{Q}, @var{Z}, @var{keyword})
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	61 @deftypefnx {} {[@var{AR}, @var{BR}, @var{QR}, @var{ZR}] =} ordqz (@var{AA}, @var{BB}, @var{Q}, @var{Z}, @var{select})
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	62 Reorder the QZ@tie{}decomposition of a generalized eigenvalue problem.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	63
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	64 The generalized eigenvalue problem is defined as
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	65
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	66 @tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	67 $$A x = \lambda B x$$
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	68 @end tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	69 @ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	70
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	71 @math{A x = @var{lambda} B x}
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	72
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	73 @end ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	74
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	75 Its generalized Schur decomposition is computed using the @code{qz} algorithm:
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	76
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	77 @code{[@var{AA}, @var{BB}, @var{Q}, @var{Z}] = qz (@var{A}, @var{B})}
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	78
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	79 where @var{AA}, @var{BB}, @var{Q}, and @var{Z} fulfill
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	80 @tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	81 $$ AA = Q \cdot A \cdot Z, BB = Q \cdot B \cdot Z $$
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	82 @end tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	83 @ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	84
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	85 @example
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	86 @group
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	87
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	88 @var{AA} = @var{Q} * @var{A} * @var{Z}, @var{BB} = @var{Q} * @var{B} * @var{Z}
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	89
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	90 @end group
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	91 @end example
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	92
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	93 @end ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	94
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	95 The @code{ordqz} function computes a unitary transformation @var{QR} and
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	96 @var{ZR} such that the order of the eigenvalue on the diagonal of @var{AA} and
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	97 @var{BB} is changed. The resulting reordered matrices @var{AR} and @var{BR}
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	98 fulfill:
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	99
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	100 @tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	101 $$ A_R = Q_R \cdot A \cdot Z_R, B_R = Q_R \cdot B \cdot Z_R $$
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	102 @end tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	103 @ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	104
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	105 @example
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	106 @group
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	107
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	108 @var{AR} = @var{QR} * @var{A} * @var{ZR}, @var{BR} = @var{QR} * @var{B} * @var{ZR}
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	109
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	110 @end group
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	111 @end example
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	112
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	113 @end ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	114
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	115 The function can either be called with the @var{keyword} argument which
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	116 selects the eigenvalues in the top left block of @var{AR} and @var{BR} in the
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	117 following way:
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	118
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	119 @table @asis
30371 1c114a55f4fb doc: Spellcheck documentation ahead of 7.1 release. Rik <rik@octave.org> parents: 29961 diff changeset	120 @item @qcode{"S"}, @nospell{@qcode{"udi"}}
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	121 small: leading block has all
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	122 @tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	123 $\|\lambda\| < 1$
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	124 @end tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	125 @ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	126 \|@var{lambda}\| < 1
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	127 @end ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	128
30371 1c114a55f4fb doc: Spellcheck documentation ahead of 7.1 release. Rik <rik@octave.org> parents: 29961 diff changeset	129 @item @qcode{"B"}, @nospell{@qcode{"udo"}}
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	130 big: leading block has all
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	131 @tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	132 $\|\lambda\| \geq 1$
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	133 @end tex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	134 @ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	135 \|@var{lambda}\| @geq{} 1
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	136 @end ifnottex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	137
30371 1c114a55f4fb doc: Spellcheck documentation ahead of 7.1 release. Rik <rik@octave.org> parents: 29961 diff changeset	138 @item @qcode{"-"}, @nospell{@qcode{"lhp"}}
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	139 negative real part: leading block has all eigenvalues in the open left
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	140 half-plane
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	141
30371 1c114a55f4fb doc: Spellcheck documentation ahead of 7.1 release. Rik <rik@octave.org> parents: 29961 diff changeset	142 @item @qcode{"+"}, @nospell{@qcode{"rhp"}}
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	143 non-negative real part: leading block has all eigenvalues in the closed right
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	144 half-plane
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	145 @end table
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	146
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	147 If a logical vector @var{select} is given instead of a keyword the @code{ordqz}
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	148 function reorders all eigenvalues @code{k} to the left block for which
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	149 @code{select(k)} is true.
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	150
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	151 Note: The keywords are compatible with the ones from @code{qr}.
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	152
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	153 @seealso{eig, ordeig, qz, schur, ordschur}
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	154 @end deftypefn */)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	155 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	156 enum { LHP, RHP, UDI, UDO, VEC, NONE } select_mode = NONE;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	157
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	158 if (args.length () != 5)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	159 print_usage ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	160
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	161 // Check select argument
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	162 if (args(4).is_string())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	163 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	164 std::string opts = args(4).string_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	165 std::for_each (opts.begin (), opts.end (),
31607 aac27ad79be6 maint: Re-indent code after switch to using namespace macros. Rik <rik@octave.org> parents: 31605 diff changeset	166 [] (char& c) { c = std::tolower (c); });
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	167 if (opts == "lhp" \|\| opts == "-")
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	168 select_mode = LHP;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	169 else if (opts == "rhp" \|\| opts == "+")
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	170 select_mode = RHP;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	171 else if (opts == "udi" \|\| opts == "s")
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	172 select_mode = UDI;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	173 else if (opts == "udo" \|\| opts == "b")
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	174 select_mode = UDO;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	175 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	176 error_with_id ("Octave:ordqz:unknown-keyword",
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	177 "ordqz: unknown KEYWORD, possible values: "
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	178 "lhp, rhp, udi, udo");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	179 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	180 else if (args(4).isreal () \|\| args(4).isinteger () \|\| args(4).islogical ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	181 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	182 if (args(4).rows () > 1 && args(4).columns () > 1)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	183 error_with_id ("Octave:ordqz:select-not-vector",
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	184 "ordqz: SELECT argument must be a vector");
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	185 select_mode = VEC;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	186 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	187 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	188 error_with_id ("Octave:ordqz:unknown-arg",
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	189 "ordqz: OPT must be string or a logical vector");
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	190
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	191 if (nargout > 4)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	192 error_with_id ("Octave:ordqz:nargout",
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	193 "ordqz: at most four output arguments possible");
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	194
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	195 // Matrix A: check dimensions.
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	196 F77_INT nn = to_f77_int (args(0).rows ());
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	197 F77_INT nc = to_f77_int (args(0).columns ());
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	198
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	199 if (args(0).isempty ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	200 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	201 warn_empty_arg ("qz: A");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	202 return octave_value_list (2, Matrix ());
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	203 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	204 else if (nc != nn)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	205 err_square_matrix_required ("qz", "A");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	206
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	207 // Matrix A: get value.
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	208 Matrix aa;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	209 ComplexMatrix caa;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	210
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	211 if (args(0).iscomplex ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	212 caa = args(0).complex_matrix_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	213 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	214 aa = args(0).matrix_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	215
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	216 // Extract argument 2 (bb, or cbb if complex).
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	217 F77_INT b_nr = to_f77_int (args(1).rows ());
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	218 F77_INT b_nc = to_f77_int (args(1).columns ());
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	219
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	220 if (nn != b_nc \|\| nn != b_nr)
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	221 ::err_nonconformant ();
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	222
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	223 Matrix bb;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	224 ComplexMatrix cbb;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	225
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	226 if (args(1).iscomplex ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	227 cbb = args(1).complex_matrix_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	228 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	229 bb = args(1).matrix_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	230
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	231 // Extract argument 3 (qq, or cqq if complex).
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	232 F77_INT q_nr = to_f77_int (args(2).rows ());
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	233 F77_INT q_nc = to_f77_int (args(2).columns ());
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	234
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	235 if (nn != q_nc \|\| nn != q_nr)
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	236 ::err_nonconformant ();
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	237
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	238 Matrix qq;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	239 ComplexMatrix cqq;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	240
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	241 if (args(2).iscomplex ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	242 cqq = args(2).complex_matrix_value ().hermitian ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	243 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	244 qq = args(2).matrix_value ().transpose ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	245
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	246 // Extract argument 4 (zz, or czz if complex).
29961 7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	247 F77_INT z_nr = to_f77_int (args(3).rows ());
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more John W. Eaton <jwe@octave.org> parents: 29958 diff changeset	248 F77_INT z_nc = to_f77_int (args(3).columns ());
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	249
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	250 if (nn != z_nc \|\| nn != z_nr)
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	251 ::err_nonconformant ();
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	252
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	253 Matrix zz;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	254 ComplexMatrix czz;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	255
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	256 if (args(3).iscomplex ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	257 czz = args(3).complex_matrix_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	258 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	259 zz = args(3).matrix_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	260
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	261 bool complex_case = (args(0).iscomplex () \|\| args(1).iscomplex ()
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	262 \|\| args(2).iscomplex () \|\| args(3).iscomplex ());
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	263
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	264 if (select_mode == VEC && args(4).rows () != nn && args(4).columns () != nn)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	265 error_with_id ("Octave:ordqz:numel_select",
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	266 "ordqz: SELECT vector has the wrong number of elements");
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	267
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	268 Array<double> select_array (dim_vector (nn, 1));
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	269 if (select_mode == VEC)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	270 select_array = args(4).vector_value ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	271
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	272 Array<F77_LOGICAL> select (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	273
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	274 if (complex_case)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	275 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	276 // Complex
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	277 if (args(0).isreal ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	278 caa = ComplexMatrix (aa);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	279 if (args(1).isreal ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	280 cbb = ComplexMatrix (bb);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	281 if (args(2).isreal ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	282 cqq = ComplexMatrix (qq);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	283 if (args(3).isreal ())
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	284 czz = ComplexMatrix (zz);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	285
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	286 ComplexRowVector alpha (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	287 ComplexRowVector beta (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	288 octave_idx_type k;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	289
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	290 for (k = 0; k < nn-1; k++)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	291 {
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	292 if (caa(k+1, k) != 0.0)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	293 error_with_id ("Octave:ordqz:unsupported_AA",
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	294 "ordqz: quasi upper triangular matrices are not "
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	295 "allowed with complex data");
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	296 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	297
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	298 for (k = 0; k < nn; k++)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	299 {
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	300 alpha(k) = caa(k, k);
a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	301 beta(k) = cbb(k, k);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	302 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	303
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	304 for (k = 0; k < nn; k++)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	305 {
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	306 switch (select_mode)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	307 {
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	308 case LHP:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	309 select(k) = real (alpha(k) * beta(k)) < 0;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	310 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	311 case RHP:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	312 select(k) = real (alpha(k) * beta(k)) > 0;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	313 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	314 case UDI:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	315 if (beta(k) != 0.0)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	316 select(k) = abs (alpha(k)/beta(k)) < 1.0;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	317 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	318 select(k) = false;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	319 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	320 case UDO:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	321 if (beta(k) != 0.0)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	322 select(k) = abs (alpha(k)/beta(k)) > 1.0;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	323 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	324 select(k) = true;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	325 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	326 case VEC:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	327 if (select_array(k) != 0.0)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	328 select(k) = true;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	329 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	330 select(k) = false;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	331 break;
33344 3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	332 case NONE:
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	333 error ("ordqz: invalid select mode NONE");
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	334 break;
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	335
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	336 // We should have handled all possible enum values above.
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	337 // Rely on compiler diagnostics to warn if we haven't.
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	338 // For example, GCC's -Wswitch option, enabled by -Wall,
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	339 // will provide a warning.
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	340 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	341 }
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	342
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	343 F77_LOGICAL wantq, wantz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	344 wantq = 1, wantz = 1;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	345 F77_INT ijob, mm, lrwork3, liwork, info;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	346 ijob = 0, lrwork3 = 1, liwork = 1;
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	347 F77_DBLE pl, pr;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	348 ComplexRowVector work3 (lrwork3);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	349 Array<F77_INT> iwork (dim_vector (liwork, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	350
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	351 F77_XFCN (ztgsen, ZTGSEN,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	352 (ijob, wantq, wantz,
32660 f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	353 select.rwdata (), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	354 F77_DBLE_CMPLX_ARG (caa.rwdata ()), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	355 F77_DBLE_CMPLX_ARG (cbb.rwdata ()), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	356 F77_DBLE_CMPLX_ARG (alpha.rwdata ()),
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	357 F77_DBLE_CMPLX_ARG (beta.rwdata ()),
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	358 F77_DBLE_CMPLX_ARG (cqq.rwdata ()), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	359 F77_DBLE_CMPLX_ARG (czz.rwdata ()), nn,
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	360 mm,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	361 pl, pr,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	362 nullptr,
32660 f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	363 F77_DBLE_CMPLX_ARG (work3.rwdata ()), lrwork3,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	364 iwork.rwdata (), liwork,
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	365 info));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	366 if (info != 0)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	367 error_with_id ("Octave:ordqz:ztgsen_failed",
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	368 "ordqz: failed to reorder eigenvalues");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	369 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	370 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	371 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	372 // Extract eigenvalues
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	373 RowVector alphar (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	374 RowVector alphai (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	375 RowVector beta (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	376
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	377 octave_idx_type k;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	378
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	379 k = 0;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	380 while (k < nn)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	381 {
33167 8f8bb94a1b83 for consistent coding style, use #if defined, not #ifdef or #ifndef John W. Eaton <jwe@octave.org> parents: 32660 diff changeset	382 #if defined (DEBUG)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	383 octave_stdout << "ordqz: k = " << k << " nn = " << nn << " \n";
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	384 #endif
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	385 if ((k < nn-1 && aa(k+1, k) == 0.0) \|\| k == nn-1)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	386 {
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	387 alphar(k) = aa(k, k);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	388 alphai(k) = 0.0;
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	389 beta(k) = bb(k, k);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	390 k++;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	391 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	392 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	393 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	394 double ar[2], ai[2], b[2], work[4];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	395 char qz_job = 'E';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	396 char comp_q = 'N';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	397 char comp_z = 'N';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	398 F77_INT nl = 2;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	399 F77_INT ilo = 1;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	400 F77_INT ihi = 2;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	401 F77_INT lwork = 4;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	402 F77_INT info = 0;
32660 f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	403 double *aa_vec = aa.rwdata ();
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	404 double *bb_vec = bb.rwdata ();
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	405
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	406 F77_XFCN (dhgeqz, DHGEQZ,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	407 (F77_CONST_CHAR_ARG2 (&qz_job, 1),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	408 F77_CONST_CHAR_ARG2 (&comp_q, 1),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	409 F77_CONST_CHAR_ARG2 (&comp_z, 1),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	410 nl, ilo, ihi,
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	411 &aa_vec[k+k*nn], nn,
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	412 &bb_vec[k+k*nn], nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	413 ar, ai, b,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	414 nullptr, nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	415 nullptr, nn, work, lwork, info
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	416 F77_CHAR_ARG_LEN (1)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	417 F77_CHAR_ARG_LEN (1)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	418 F77_CHAR_ARG_LEN (1)));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	419 if (info != 0)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	420 error("ordqz: failed to extract eigenvalues");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	421
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	422 alphar(k) = ar[0];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	423 alphar(k+1) = ar[1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	424 alphai(k) = ai[0];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	425 alphai(k+1) = ai[1];
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	426 beta(k) = b[0];
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	427 beta(k+1) = b[1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	428
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	429 k += 2;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	430 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	431
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	432 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	433
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	434 for (k = 0; k < nn; k++)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	435 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	436 switch (select_mode)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	437 {
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	438 case LHP:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	439 select(k) = alphar(k) * beta(k) < 0;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	440 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	441 case RHP:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	442 select(k) = alphar(k) * beta(k) > 0;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	443 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	444 case UDI:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	445 select(k) = alphar(k)*alphar(k)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	446 + alphai(k)alphai(k) < beta(k)beta(k);
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	447 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	448 case UDO:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	449 select(k) = alphar(k)*alphar(k)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	450 + alphai(k)alphai(k) > beta(k)beta(k);
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	451 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	452 case VEC:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	453 if (select_array(k) != 0.0)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	454 select(k) = true;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	455 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	456 select(k) = false;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	457 break;
33344 3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	458 case NONE:
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	459 error ("ordqz: invalid select mode NONE");
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	460 break;
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	461
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	462 // We should have handled all possible enum values above.
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	463 // Rely on compiler diagnostics to warn if we haven't.
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	464 // For example, GCC's -Wswitch option, enabled by -Wall,
3606817b9994 take adavantage of compiler warnings for unhandled enum switch cases John W. Eaton <jwe@octave.org> parents: 33167 diff changeset	465 // will provide a warning.
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	466 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	467 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	468
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	469 F77_LOGICAL wantq, wantz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	470 wantq = 1, wantz = 1;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	471 F77_INT ijob, mm, lrwork3, liwork, info;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	472 ijob = 0, lrwork3 = 4*nn+16, liwork = nn;
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	473 F77_DBLE pl, pr;
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	474 RowVector rwork3 (lrwork3);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	475 Array<F77_INT> iwork (dim_vector (liwork, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	476
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	477 F77_XFCN (dtgsen, DTGSEN,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	478 (ijob, wantq, wantz,
32660 f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	479 select.rwdata (), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	480 aa.rwdata (), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	481 bb.rwdata (), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	482 alphar.rwdata (),
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	483 alphai.rwdata (),
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	484 beta.rwdata (),
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	485 qq.rwdata (), nn,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	486 zz.rwdata (), nn,
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	487 mm,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	488 pl, pr,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	489 nullptr,
32660 f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	490 rwork3.rwdata (), lrwork3,
f53ac65ffba6 maint: New method rwdata() as clearer alternative to fortran_vec(). Rik <rik@octave.org> parents: 32632 diff changeset	491 iwork.rwdata (), liwork,
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	492 info));
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	493 if (info != 0)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	494 error("ordqz: failed to reorder eigenvalues");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	495 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	496
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	497 octave_value_list retval (nargout);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	498 switch (nargout)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	499 {
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	500 case 4:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	501 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	502 retval(3) = czz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	503 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	504 retval(3) = zz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	505 OCTAVE_FALLTHROUGH;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	506 case 3:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	507 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	508 retval(2) = cqq.hermitian();
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	509 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	510 retval(2) = qq.transpose();
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	511 OCTAVE_FALLTHROUGH;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	512 case 2:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	513 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	514 retval(1) = cbb;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	515 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	516 retval(1) = bb;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	517 OCTAVE_FALLTHROUGH;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	518 case 1:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	519 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	520 retval(0) = caa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	521 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	522 retval(0) = aa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	523 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	524 case 0:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	525 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	526 retval(0) = caa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	527 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	528 retval(0) = aa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	529 break;
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	530 }
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	531
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	532 return retval;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	533 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	534
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	535
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	536 /*
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	537 %!shared A, B, AA, BB, QQ, ZZ, AC, BC, AAC, BBC, QQC, ZZC, select, selectc
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	538 %! A = [ -1.03428 0.24929 0.43205 -0.12860;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	539 %! 1.16228 0.27870 2.12954 0.69250;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	540 %! -0.51524 -0.34939 -0.77820 2.13721;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	541 %! -1.32941 2.11870 0.72005 1.00835 ];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	542 %! B = [ 1.407302 -0.632956 -0.360628 0.068534;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	543 %! 0.149898 0.298248 0.991777 0.023652;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	544 %! 0.169281 -0.405205 -1.775834 1.511730;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	545 %! 0.717770 1.291390 -1.766607 -0.531352 ];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	546 %! AC = [ 0.4577 + 0.7199i 0.1476 + 0.6946i 0.6202 + 0.2092i 0.7559 + 0.2759i;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	547 %! 0.5868 + 0.7275i 0.9174 + 0.8781i 0.6741 + 0.1985i 0.4320 + 0.7023i;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	548 %! 0.2408 + 0.6359i 0.2959 + 0.8501i 0.3904 + 0.5613i 0.5000 + 0.1428i;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	549 %! 0.8177 + 0.8581i 0.2583 + 0.8970i 0.7706 + 0.5451i 0.1068 + 0.1650i];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	550 %! BC = [ 0.089898 + 0.209257i 0.157769 + 0.311387i 0.018926 + 0.622517i 0.058825 + 0.374647i;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	551 %! 0.009367 + 0.098211i 0.736087 + 0.095797i 0.973192 + 0.583765i 0.434018 + 0.461909i;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	552 %! 0.880784 + 0.868215i 0.032839 + 0.569461i 0.873437 + 0.266081i 0.739426 + 0.362017i;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	553 %! 0.121649 + 0.115111i 0.426695 + 0.492222i 0.247670 + 0.034414i 0.771629 + 0.078153i];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	554 %! [AA, BB, QQ, ZZ] = qz (A, B);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	555 %! [AAC, BBC, QQC, ZZC] = qz (AC, BC);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	556 %! select = [0 0 1 1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	557 %! selectc = [0 0 0 1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	558
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	559 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	560 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "rhp");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	561 %! assert (all (real (eig (AAX(1:3,1:3), BBX(1:3,1:3))) >= 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	562 %! assert (all (real (eig (AAX(4:4,4:4), BBX(4:4,4:4))) < 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	563 %! assert (norm (QQX'AAXZZX' - A, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	564 %! assert (norm (QQX'BBXZZX' - B, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	565
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	566 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	567 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "+");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	568 %! assert (all (real (eig (AAX(1:3,1:3), BBX(1:3,1:3))) >= 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	569 %! assert (all (real (eig (AAX(4:4,4:4), BBX(4:4,4:4))) < 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	570
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	571 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	572 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "lhp");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	573 %! assert (all (real (eig (AAX(2:4,2:4), BBX(2:4,2:4))) >= 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	574 %! assert (all (real (eig (AAX(1:1,1:1), BBX(1:1,1:1))) < 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	575 %! assert (norm (QQX'AAXZZX' - A, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	576 %! assert (norm (QQX'BBXZZX' - B, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	577
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	578 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	579 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "-");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	580 %! assert (all (real (eig (AAX(2:4,2:4), BBX(2:4,2:4))) >= 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	581 %! assert (all (real (eig (AAX(1:1,1:1), BBX(1:1,1:1))) < 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	582
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	583 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	584 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "udi");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	585 %! assert (all (abs (eig (AAX(1:1,1:1), BBX(1:1,1:1))) < 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	586 %! assert (all (abs (eig (AAX(2:4,2:4), BBX(2:4,2:4))) > 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	587 %! assert (norm (QQX'AAXZZX' - A, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	588 %! assert (norm (QQX'BBXZZX' - B, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	589
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	590 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	591 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "S");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	592 %! assert (all (abs (eig (AAX(1:1,1:1), BBX(1:1,1:1))) < 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	593 %! assert (all (abs (eig (AAX(2:4,2:4), BBX(2:4,2:4))) > 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	594
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	595 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	596 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "udo");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	597 %! assert (all (abs (eig (AAX(1:3,1:3), BBX(1:3,1:3))) >= 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	598 %! assert (all (abs (eig (AAX(4:4,4:4), BBX(4:4,4:4))) < 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	599 %! assert (norm (QQX'AAXZZX' - A, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	600 %! assert (norm (QQX'BBXZZX' - B, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	601
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	602 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	603 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, "B");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	604 %! assert (all (abs (eig (AAX(1:3,1:3), BBX(1:3,1:3))) >= 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	605 %! assert (all (abs (eig (AAX(4:4,4:4), BBX(4:4,4:4))) < 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	606
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	607 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	608 %! [AAX, BBX, QQX, ZZX] = ordqz (AA, BB, QQ, ZZ, select);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	609 %! assert (all (iscomplex (eig (AAX(1:2,1:2), BBX(1:2,1:2)))));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	610 %! assert (norm (QQX'AAXZZX' - A, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	611 %! assert (norm (QQX'BBXZZX' - B, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	612
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	613 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	614 %! [AACX, BBCX, QQCX, ZZCX] = ordqz (AAC, BBC, QQC, ZZC, "rhp");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	615 %! assert (all (real (eig (AACX(1:2,1:2), BBCX(1:2,1:2))) >= 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	616 %! assert (all (real (eig (AACX(3:4,3:4), BBCX(3:4,3:4))) < 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	617 %! assert (norm (QQCX'AACXZZCX' - AC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	618 %! assert (norm (QQCX'BBCXZZCX' - BC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	619
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	620 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	621 %! [AACX, BBCX, QQCX, ZZCX] = ordqz (AAC, BBC, QQC, ZZC, "lhp");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	622 %! assert (all (real (eig (AACX(1:2,1:2), BBCX(1:2,1:2))) < 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	623 %! assert (all (real (eig (AACX(3:4,3:4), BBCX(3:4,3:4))) >= 0));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	624 %! assert (norm (QQCX'AACXZZCX' - AC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	625 %! assert (norm (QQCX'BBCXZZCX' - BC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	626
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	627 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	628 %! [AACX, BBCX, QQCX, ZZCX] = ordqz (AAC, BBC, QQC, ZZC, "udi");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	629 %! assert (all (abs (eig (AACX(1:2,1:2), BBCX(1:2,1:2))) < 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	630 %! assert (all (abs (eig (AACX(3:4,3:4), BBCX(3:4,3:4))) >= 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	631 %! assert (norm (QQCX'AACXZZCX' - AC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	632 %! assert (norm (QQCX'BBCXZZCX' - BC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	633
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	634 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	635 %! [AACX, BBCX, QQCX, ZZCX] = ordqz (AAC, BBC, QQC, ZZC, "udo");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	636 %! assert (all (abs (eig (AACX(1:2,1:2), BBCX(1:2,1:2))) >= 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	637 %! assert (all (abs (eig (AACX(3:4,3:4), BBCX(3:4,3:4))) < 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	638 %! assert (norm (QQCX'AACXZZCX' - AC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	639 %! assert (norm (QQCX'BBCXZZCX' - BC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	640
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	641 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	642 %! [AACX, BBCX, QQCX, ZZCX] = ordqz (AAC, BBC, QQC, ZZC, selectc);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	643 %! ev = abs (eig (AACX(1:1,1:1), BBCX(1:1,1:1)));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	644 %! assert(ev > 0.6 && ev < 0.7);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	645 %! assert (norm (QQCX'AACXZZCX' - AC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	646 %! assert (norm (QQCX'BBCXZZCX' - BC, "fro"), 0, 1e-12);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	647
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	648 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	649 %! A = toeplitz ([1,2,3,4]);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	650 %! [B, A] = qr (A);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	651 %! B = B';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	652 %! [AA, BB, Q, Z] = qz (A, B);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	653 %! [AAS, BBS, QS, ZS] = ordqz (AA, BB, Q, Z, "lhp");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	654 %! E2 = ordeig (AAS, BBS);
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	655 %! ECOMP = [-3.414213562373092; -1.099019513592784;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	656 %! -0.5857864376269046; 9.099019513592784];
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	657 %! assert (norm (ECOMP - E2, "Inf"), 0, 1e-8);
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	658
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	659 ## Test input validation
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	660 %!error <Invalid call> ordqz ()
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	661 %!error <Invalid call> ordqz (eye (2))
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	662 %!error <Invalid call> ordqz (eye (2), eye (2))
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	663 %!error <Invalid call> ordqz (eye (2), eye (2), eye (2))
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	664 %!error <Invalid call> ordqz (eye (2), eye (2), eye (2), eye (2))
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	665 %!error id=Octave:ordqz:unknown-keyword
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	666 %! ordqz (eye (2), eye (2), eye (2), eye (2), "foobar");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	667 %!error id=Octave:ordqz:select-not-vector
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	668 %! ordqz (eye (2), eye (2), eye (2), eye (2), eye (2));
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	669 %!error id=Octave:ordqz:unknown-arg
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	670 %! ordqz (eye (2), eye (2), eye (2), eye (2), {"foobar"});
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	671 %!error id=Octave:ordqz:nargout
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	672 %! [a,b,c,d,e] = ordqz (eye (2), eye (2), eye (2), eye (2), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	673 %!warning <A: argument is empty matrix> ordqz ([], [], [], [], "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	674 %!error <A must be a square matrix> ordqz (ones (1,2), [], [], [], "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	675 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	676 %! ordqz (eye (3), eye (2), eye (2), eye (2), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	677 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	678 %! ordqz (eye (2), eye (3), eye (2), eye (2), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	679 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	680 %! ordqz (eye (2), eye (2), eye (3), eye (2), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	681 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	682 %! ordqz (eye (2), eye (2), eye (2), eye (3), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	683 %!error <SELECT vector .* wrong number of elements>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	684 %! ordqz (eye (2), eye (2), eye (2), eye (2), ones (1,5));
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	685 %!error <quasi upper triangular matrices are not allowed with complex data>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	686 %! AA = zeros (2, 2);
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	687 %! AA(2,1) = i;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	688 %! ordqz (AA, eye (2), eye (2), eye (2), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	689
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	690 */
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	691
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 30564 diff changeset	692 OCTAVE_END_NAMESPACE(octave)

Mercurial > octave

annotate libinterp/corefcn/ordqz.cc @ 33617:ec2635a02328 bytecode-interpreter tip