octave: libinterp/corefcn/ordqz.cc annotate

annotate libinterp/corefcn/ordqz.cc @ 31605:e88a07dec498 stable

maint: Use macros to begin/end C++ namespaces. * oct-conf-post-public.in.h: Define two macros (OCTAVE_BEGIN_NAMESPACE, OCTAVE_END_NAMESPACE) that can be used to start/end a namespace. * mk-opts.pl, build-env.h, build-env.in.cc, __betainc__.cc, __contourc__.cc, __dsearchn__.cc, __eigs__.cc, __expint__.cc, __ftp__.cc, __gammainc__.cc, __ichol__.cc, __ilu__.cc, __isprimelarge__.cc, __lin_interpn__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc, amd.cc, auto-shlib.cc, auto-shlib.h, balance.cc, base-text-renderer.cc, base-text-renderer.h, besselj.cc, bitfcns.cc, bsxfun.cc, c-file-ptr-stream.cc, c-file-ptr-stream.h, call-stack.cc, call-stack.h, ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, colloc.cc, conv2.cc, daspk.cc, dasrt.cc, dassl.cc, data.cc, data.h, debug.cc, defaults.cc, defaults.h, defun-int.h, defun.cc, det.cc, dirfns.cc, display.cc, display.h, dlmread.cc, dmperm.cc, dot.cc, dynamic-ld.cc, dynamic-ld.h, eig.cc, ellipj.cc, environment.cc, environment.h, error.cc, error.h, errwarn.h, event-manager.cc, event-manager.h, event-queue.cc, event-queue.h, fcn-info.cc, fcn-info.h, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc, ft-text-renderer.cc, ft-text-renderer.h, gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc, gl-render.cc, gl-render.h, gl2ps-print.cc, gl2ps-print.h, graphics-toolkit.cc, graphics-toolkit.h, graphics.cc, graphics.in.h, gsvd.cc, gtk-manager.cc, gtk-manager.h, hash.cc, help.cc, help.h, hess.cc, hex2num.cc, hook-fcn.cc, hook-fcn.h, input.cc, input.h, interpreter-private.cc, interpreter-private.h, interpreter.cc, interpreter.h, inv.cc, jsondecode.cc, jsonencode.cc, kron.cc, latex-text-renderer.cc, latex-text-renderer.h, load-path.cc, load-path.h, load-save.cc, load-save.h, lookup.cc, ls-ascii-helper.cc, ls-ascii-helper.h, ls-oct-text.cc, ls-utils.cc, ls-utils.h, lsode.cc, lu.cc, mappers.cc, matrix_type.cc, max.cc, mex-private.h, mex.cc, mgorth.cc, nproc.cc, oct-fstrm.cc, oct-fstrm.h, oct-hdf5-types.cc, oct-hdf5-types.h, oct-hist.cc, oct-hist.h, oct-iostrm.cc, oct-iostrm.h, oct-opengl.h, oct-prcstrm.cc, oct-prcstrm.h, oct-procbuf.cc, oct-procbuf.h, oct-process.cc, oct-process.h, oct-stdstrm.h, oct-stream.cc, oct-stream.h, oct-strstrm.cc, oct-strstrm.h, oct-tex-lexer.in.ll, oct-tex-parser.yy, ordqz.cc, ordschur.cc, pager.cc, pager.h, pinv.cc, pow2.cc, pr-flt-fmt.cc, pr-output.cc, procstream.cc, procstream.h, psi.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, settings.cc, settings.h, sighandlers.cc, sighandlers.h, sparse-xdiv.cc, sparse-xdiv.h, sparse-xpow.cc, sparse-xpow.h, sparse.cc, spparms.cc, sqrtm.cc, stack-frame.cc, stack-frame.h, stream-euler.cc, strfind.cc, strfns.cc, sub2ind.cc, svd.cc, sylvester.cc, symbfact.cc, syminfo.cc, syminfo.h, symrcm.cc, symrec.cc, symrec.h, symscope.cc, symscope.h, symtab.cc, symtab.h, syscalls.cc, sysdep.cc, sysdep.h, text-engine.cc, text-engine.h, text-renderer.cc, text-renderer.h, time.cc, toplev.cc, tril.cc, tsearch.cc, typecast.cc, url-handle-manager.cc, url-handle-manager.h, urlwrite.cc, utils.cc, utils.h, variables.cc, variables.h, xdiv.cc, xdiv.h, xnorm.cc, xnorm.h, xpow.cc, xpow.h, __delaunayn__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __init_fltk__.cc, __init_gnuplot__.cc, __ode15__.cc, __voronoi__.cc, audiodevinfo.cc, audioread.cc, convhulln.cc, fftw.cc, gzip.cc, mk-build-env-features.sh, mk-builtins.pl, cdef-class.cc, cdef-class.h, cdef-fwd.h, cdef-manager.cc, cdef-manager.h, cdef-method.cc, cdef-method.h, cdef-object.cc, cdef-object.h, cdef-package.cc, cdef-package.h, cdef-property.cc, cdef-property.h, cdef-utils.cc, cdef-utils.h, ov-base.cc, ov-base.h, ov-bool-mat.cc, ov-builtin.h, ov-cell.cc, ov-class.cc, ov-class.h, ov-classdef.cc, ov-classdef.h, ov-complex.cc, ov-fcn-handle.cc, ov-fcn-handle.h, ov-fcn.h, ov-java.cc, ov-java.h, ov-mex-fcn.h, ov-null-mat.cc, ov-oncleanup.cc, ov-struct.cc, ov-typeinfo.cc, ov-typeinfo.h, ov-usr-fcn.cc, ov-usr-fcn.h, ov.cc, ov.h, octave.cc, octave.h, mk-ops.sh, op-b-b.cc, op-b-bm.cc, op-b-sbm.cc, op-bm-b.cc, op-bm-bm.cc, op-bm-sbm.cc, op-cdm-cdm.cc, op-cell.cc, op-chm.cc, op-class.cc, op-cm-cm.cc, op-cm-cs.cc, op-cm-m.cc, op-cm-s.cc, op-cm-scm.cc, op-cm-sm.cc, op-cs-cm.cc, op-cs-cs.cc, op-cs-m.cc, op-cs-s.cc, op-cs-scm.cc, op-cs-sm.cc, op-dm-dm.cc, op-dm-scm.cc, op-dm-sm.cc, op-dm-template.cc, op-dms-template.cc, op-fcdm-fcdm.cc, op-fcm-fcm.cc, op-fcm-fcs.cc, op-fcm-fm.cc, op-fcm-fs.cc, op-fcn.cc, op-fcs-fcm.cc, op-fcs-fcs.cc, op-fcs-fm.cc, op-fcs-fs.cc, op-fdm-fdm.cc, op-fm-fcm.cc, op-fm-fcs.cc, op-fm-fm.cc, op-fm-fs.cc, op-fs-fcm.cc, op-fs-fcs.cc, op-fs-fm.cc, op-fs-fs.cc, op-i16-i16.cc, op-i32-i32.cc, op-i64-i64.cc, op-i8-i8.cc, op-int-concat.cc, op-m-cm.cc, op-m-cs.cc, op-m-m.cc, op-m-s.cc, op-m-scm.cc, op-m-sm.cc, op-mi.cc, op-pm-pm.cc, op-pm-scm.cc, op-pm-sm.cc, op-pm-template.cc, op-range.cc, op-s-cm.cc, op-s-cs.cc, op-s-m.cc, op-s-s.cc, op-s-scm.cc, op-s-sm.cc, op-sbm-b.cc, op-sbm-bm.cc, op-sbm-sbm.cc, op-scm-cm.cc, op-scm-cs.cc, op-scm-m.cc, op-scm-s.cc, op-scm-scm.cc, op-scm-sm.cc, op-sm-cm.cc, op-sm-cs.cc, op-sm-m.cc, op-sm-s.cc, op-sm-scm.cc, op-sm-sm.cc, op-str-m.cc, op-str-s.cc, op-str-str.cc, op-struct.cc, op-ui16-ui16.cc, op-ui32-ui32.cc, op-ui64-ui64.cc, op-ui8-ui8.cc, ops.h, anon-fcn-validator.cc, anon-fcn-validator.h, bp-table.cc, bp-table.h, comment-list.cc, comment-list.h, filepos.h, lex.h, lex.ll, oct-lvalue.cc, oct-lvalue.h, oct-parse.yy, parse.h, profiler.cc, profiler.h, pt-anon-scopes.cc, pt-anon-scopes.h, pt-arg-list.cc, pt-arg-list.h, pt-args-block.cc, pt-args-block.h, pt-array-list.cc, pt-array-list.h, pt-assign.cc, pt-assign.h, pt-binop.cc, pt-binop.h, pt-bp.cc, pt-bp.h, pt-cbinop.cc, pt-cbinop.h, pt-cell.cc, pt-cell.h, pt-check.cc, pt-check.h, pt-classdef.cc, pt-classdef.h, pt-cmd.h, pt-colon.cc, pt-colon.h, pt-const.cc, pt-const.h, pt-decl.cc, pt-decl.h, pt-eval.cc, pt-eval.h, pt-except.cc, pt-except.h, pt-exp.cc, pt-exp.h, pt-fcn-handle.cc, pt-fcn-handle.h, pt-id.cc, pt-id.h, pt-idx.cc, pt-idx.h, pt-jump.h, pt-loop.cc, pt-loop.h, pt-mat.cc, pt-mat.h, pt-misc.cc, pt-misc.h, pt-pr-code.cc, pt-pr-code.h, pt-select.cc, pt-select.h, pt-spmd.cc, pt-spmd.h, pt-stmt.cc, pt-stmt.h, pt-tm-const.cc, pt-tm-const.h, pt-unop.cc, pt-unop.h, pt-vm-eval.cc, pt-walk.cc, pt-walk.h, pt.cc, pt.h, token.cc, token.h, Range.cc, Range.h, idx-vector.cc, idx-vector.h, range-fwd.h, CollocWt.cc, CollocWt.h, aepbalance.cc, aepbalance.h, chol.cc, chol.h, gepbalance.cc, gepbalance.h, gsvd.cc, gsvd.h, hess.cc, hess.h, lo-mappers.cc, lo-mappers.h, lo-specfun.cc, lo-specfun.h, lu.cc, lu.h, oct-convn.cc, oct-convn.h, oct-fftw.cc, oct-fftw.h, oct-norm.cc, oct-norm.h, oct-rand.cc, oct-rand.h, oct-spparms.cc, oct-spparms.h, qr.cc, qr.h, qrp.cc, qrp.h, randgamma.cc, randgamma.h, randmtzig.cc, randmtzig.h, randpoisson.cc, randpoisson.h, schur.cc, schur.h, sparse-chol.cc, sparse-chol.h, sparse-lu.cc, sparse-lu.h, sparse-qr.cc, sparse-qr.h, svd.cc, svd.h, child-list.cc, child-list.h, dir-ops.cc, dir-ops.h, file-ops.cc, file-ops.h, file-stat.cc, file-stat.h, lo-sysdep.cc, lo-sysdep.h, lo-sysinfo.cc, lo-sysinfo.h, mach-info.cc, mach-info.h, oct-env.cc, oct-env.h, oct-group.cc, oct-group.h, oct-password.cc, oct-password.h, oct-syscalls.cc, oct-syscalls.h, oct-time.cc, oct-time.h, oct-uname.cc, oct-uname.h, action-container.cc, action-container.h, base-list.h, cmd-edit.cc, cmd-edit.h, cmd-hist.cc, cmd-hist.h, f77-fcn.h, file-info.cc, file-info.h, lo-array-errwarn.cc, lo-array-errwarn.h, lo-hash.cc, lo-hash.h, lo-ieee.h, lo-regexp.cc, lo-regexp.h, lo-utils.cc, lo-utils.h, oct-base64.cc, oct-base64.h, oct-glob.cc, oct-glob.h, oct-inttypes.h, oct-mutex.cc, oct-mutex.h, oct-refcount.h, oct-shlib.cc, oct-shlib.h, oct-sparse.cc, oct-sparse.h, oct-string.h, octave-preserve-stream-state.h, pathsearch.cc, pathsearch.h, quit.cc, quit.h, unwind-prot.cc, unwind-prot.h, url-transfer.cc, url-transfer.h : Use new macros to begin/end C++ namespaces.

author	Rik <rik@octave.org>
date	Thu, 01 Dec 2022 14:23:45 -0800
parents	796f54d4ddbf
children	aac27ad79be6

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 //
30564 796f54d4ddbf update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 30390 diff changeset	3 // Copyright (C) 2020-2022 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 (),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	166 [] (char & c) { c = std::tolower (c); });
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;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	332 default:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	333 // default: case just here to suppress compiler warning.
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	334 panic_impossible ();
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	335 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	336 }
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	337
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	338 F77_LOGICAL wantq, wantz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	339 wantq = 1, wantz = 1;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	340 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	341 ijob = 0, lrwork3 = 1, liwork = 1;
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	342 F77_DBLE pl, pr;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	343 ComplexRowVector work3 (lrwork3);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	344 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	345
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	346 F77_XFCN (ztgsen, ZTGSEN,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	347 (ijob, wantq, wantz,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	348 select.fortran_vec (), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	349 F77_DBLE_CMPLX_ARG (caa.fortran_vec ()), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	350 F77_DBLE_CMPLX_ARG (cbb.fortran_vec ()), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	351 F77_DBLE_CMPLX_ARG (alpha.fortran_vec ()),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	352 F77_DBLE_CMPLX_ARG (beta.fortran_vec ()),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	353 F77_DBLE_CMPLX_ARG (cqq.fortran_vec ()), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	354 F77_DBLE_CMPLX_ARG (czz.fortran_vec ()), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	355 mm,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	356 pl, pr,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	357 nullptr,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	358 F77_DBLE_CMPLX_ARG (work3.fortran_vec ()), lrwork3,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	359 iwork.fortran_vec (), liwork,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	360 info));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	361 if (info != 0)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	362 error_with_id ("Octave:ordqz:ztgsen_failed",
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	363 "ordqz: failed to reorder eigenvalues");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	364 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	365 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	366 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	367 // Extract eigenvalues
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	368 RowVector alphar (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	369 RowVector alphai (dim_vector (nn, 1));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	370 RowVector beta (dim_vector (nn, 1));
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 octave_idx_type k;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	373
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	374 k = 0;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	375 while (k < nn)
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 #ifdef DEBUG
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	378 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	379 #endif
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	380 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	381 {
30390 a61e1a0f6024 maint: style check C++ files in libinterp/ ahead of 7.1 release. Rik <rik@octave.org> parents: 30371 diff changeset	382 alphar(k) = aa(k, k);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	383 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	384 beta(k) = bb(k, k);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	385 k++;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	386 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	387 else
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	388 {
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	389 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	390 char qz_job = 'E';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	391 char comp_q = 'N';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	392 char comp_z = 'N';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	393 F77_INT nl = 2;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	394 F77_INT ilo = 1;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	395 F77_INT ihi = 2;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	396 F77_INT lwork = 4;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	397 F77_INT info = 0;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	398 double * aa_vec = aa.fortran_vec ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	399 double * bb_vec = bb.fortran_vec ();
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	400
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	401 F77_XFCN (dhgeqz, DHGEQZ,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	402 (F77_CONST_CHAR_ARG2 (&qz_job, 1),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	403 F77_CONST_CHAR_ARG2 (&comp_q, 1),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	404 F77_CONST_CHAR_ARG2 (&comp_z, 1),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	405 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	406 &aa_vec[k+k*nn], nn,
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	407 &bb_vec[k+k*nn], nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	408 ar, ai, b,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	409 nullptr, nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	410 nullptr, nn, work, lwork, info
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	411 F77_CHAR_ARG_LEN (1)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	412 F77_CHAR_ARG_LEN (1)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	413 F77_CHAR_ARG_LEN (1)));
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	414 if (info != 0)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	415 error("ordqz: failed to extract eigenvalues");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	416
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	417 alphar(k) = ar[0];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	418 alphar(k+1) = ar[1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	419 alphai(k) = ai[0];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	420 alphai(k+1) = ai[1];
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	421 beta(k) = b[0];
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	422 beta(k+1) = b[1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	423
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	424 k += 2;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	425 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	426
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	427 }
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 for (k = 0; k < nn; k++)
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 switch (select_mode)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	432 {
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	433 case LHP:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	434 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	435 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	436 case RHP:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	437 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	438 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	439 case UDI:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	440 select(k) = alphar(k)*alphar(k)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	441 + 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	442 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	443 case UDO:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	444 select(k) = alphar(k)*alphar(k)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	445 + 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	446 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	447 case VEC:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	448 if (select_array(k) != 0.0)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	449 select(k) = true;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	450 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	451 select(k) = false;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	452 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	453 default:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	454 // default: case just here to suppress compiler warning.
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	455 panic_impossible();
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	456 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	457 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	458
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	459 F77_LOGICAL wantq, wantz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	460 wantq = 1, wantz = 1;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	461 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	462 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	463 F77_DBLE pl, pr;
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	464 RowVector rwork3 (lrwork3);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	465 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	466
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	467 F77_XFCN (dtgsen, DTGSEN,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	468 (ijob, wantq, wantz,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	469 select.fortran_vec (), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	470 aa.fortran_vec (), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	471 bb.fortran_vec (), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	472 alphar.fortran_vec (),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	473 alphai.fortran_vec (),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	474 beta.fortran_vec (),
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	475 qq.fortran_vec (), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	476 zz.fortran_vec (), nn,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	477 mm,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	478 pl, pr,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	479 nullptr,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	480 rwork3.fortran_vec (), lrwork3,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	481 iwork.fortran_vec (), liwork,
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	482 info));
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	483 if (info != 0)
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	484 error("ordqz: failed to reorder eigenvalues");
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	485 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	486
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	487 octave_value_list retval (nargout);
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	488 switch (nargout)
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	489 {
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	490 case 4:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	491 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	492 retval(3) = czz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	493 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	494 retval(3) = zz;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	495 OCTAVE_FALLTHROUGH;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	496 case 3:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	497 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	498 retval(2) = cqq.hermitian();
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	499 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	500 retval(2) = qq.transpose();
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	501 OCTAVE_FALLTHROUGH;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	502 case 2:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	503 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	504 retval(1) = cbb;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	505 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	506 retval(1) = bb;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	507 OCTAVE_FALLTHROUGH;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	508 case 1:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	509 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	510 retval(0) = caa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	511 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	512 retval(0) = aa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	513 break;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	514 case 0:
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	515 if (complex_case)
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	516 retval(0) = caa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	517 else
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	518 retval(0) = aa;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	519 break;
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	520 }
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	521
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	522 return retval;
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	523 }
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	524
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	525
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	526 /*
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	527 %!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	528 %! 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	529 %! 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	530 %! -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	531 %! -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	532 %! 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	533 %! 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	534 %! 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	535 %! 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	536 %! 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	537 %! 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	538 %! 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	539 %! 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	540 %! 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	541 %! 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	542 %! 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	543 %! 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	544 %! [AA, BB, QQ, ZZ] = qz (A, B);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	545 %! [AAC, BBC, QQC, ZZC] = qz (AC, BC);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	546 %! select = [0 0 1 1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	547 %! selectc = [0 0 0 1];
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	548
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	549 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	550 %! [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	551 %! 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	552 %! 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	553 %! 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	554 %! 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	555
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	556 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	557 %! [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	558 %! 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	559 %! 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	560
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	561 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	562 %! [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	563 %! 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	564 %! 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	565 %! 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	566 %! 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	567
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	568 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	569 %! [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	570 %! 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	571 %! 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	572
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	573 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	574 %! [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	575 %! 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	576 %! 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	577 %! 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	578 %! 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	579
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	580 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	581 %! [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	582 %! 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	583 %! 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	584
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	585 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	586 %! [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	587 %! 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	588 %! 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	589 %! 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	590 %! 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	591
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	592 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	593 %! [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	594 %! 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	595 %! 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	596
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	597 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	598 %! [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	599 %! 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	600 %! 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	601 %! 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	602
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	603 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	604 %! [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	605 %! 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	606 %! 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	607 %! 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	608 %! 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	609
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	610 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	611 %! [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	612 %! 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	613 %! 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	614 %! 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	615 %! 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	616
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	617 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	618 %! [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	619 %! 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	620 %! 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	621 %! 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	622 %! 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	623
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	624 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	625 %! [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	626 %! 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	627 %! 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	628 %! 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	629 %! 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	630
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	631 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	632 %! [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	633 %! 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	634 %! assert(ev > 0.6 && ev < 0.7);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	635 %! 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	636 %! 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	637
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	638 %!test
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	639 %! A = toeplitz ([1,2,3,4]);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	640 %! [B, A] = qr (A);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	641 %! B = B';
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	642 %! [AA, BB, Q, Z] = qz (A, B);
69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	643 %! [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	644 %! E2 = ordeig (AAS, BBS);
28995 0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	645 %! ECOMP = [-3.414213562373092; -1.099019513592784;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	646 %! -0.5857864376269046; 9.099019513592784];
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	647 %! 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	648
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	649 ## Test input validation
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	650 %!error <Invalid call> ordqz ()
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	651 %!error <Invalid call> ordqz (eye (2))
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	652 %!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	653 %!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	654 %!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	655 %!error id=Octave:ordqz:unknown-keyword
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	656 %! 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	657 %!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	658 %! 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	659 %!error id=Octave:ordqz:unknown-arg
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	660 %! 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	661 %!error id=Octave:ordqz:nargout
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	662 %! [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	663 %!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	664 %!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	665 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	666 %! 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	667 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	668 %! 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	669 %!error <nonconformant matrices>
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 (3), eye (2), "udi");
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	671 %!error <nonconformant matrices>
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	672 %! 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	673 %!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	674 %! 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	675 %!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	676 %! AA = zeros (2, 2);
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	677 %! AA(2,1) = i;
0e86e50bebeb ordqz.cc: Review and update for Octave coding conventions. Rik <rik@octave.org> parents: 28986 diff changeset	678 %! 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	679
28986 69ec8d9e769b Add function ordqz (patch #9897). Martin Köhler <koehlerm@mpi-magdeburg.mpg.de> parents: diff changeset	680 */
29958 32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	681
31605 e88a07dec498 maint: Use macros to begin/end C++ namespaces. Rik <rik@octave.org> parents: 30564 diff changeset	682 OCTAVE_END_NAMESPACE(octave)

Mercurial > octave

annotate libinterp/corefcn/ordqz.cc @ 31605:e88a07dec498 stable