octave: liboctave/numeric/schur.cc annotate

annotate liboctave/numeric/schur.cc @ 22317:6ca3acf5fad8

move some new numeric classes to namespace octave::math * aepbalance.cc, aepbalance.h, chol.cc, chol.h, gepbalance.cc, gepbalance.h, hess.cc, hess.h, lu.cc, lu.h, qr.cc, qr.h, qrp.cc, qrp.h, schur.cc, schur.h, sparse-chol.cc, sparse-chol.h, sparse-dmsolve.cc,sparse-lu.cc, sparse-lu.h, sparse-qr.cc, sparse-qr.h, svd.cc, svd.h: Move classes to namespace octave::math. * __luinc__.cc, __qp__.cc, balance.cc, hess.cc, lu.cc, qz.cc, schur.cc, sqrtm.cc, svd.cc, chol.cc, dmperm.cc, qr.cc, lex.h, CMatrix.cc, CSparse.cc, dMatrix.cc, dSparse.cc, fCMatrix.cc, fMatrix.cc, eigs-base.cc, oct-norm.cc: Update for new namespaces.

author	John W. Eaton <jwe@octave.org>
date	Tue, 16 Aug 2016 15:48:22 -0400
parents	407c66ae1e20
children	93b3cdd36854

rev	line source
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	1 /*
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	2
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 17769 diff changeset	3 Copyright (C) 1994-2015 John W. Eaton
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	4
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	5 This file is part of Octave.
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	6
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	7 Octave is free software; you can redistribute it and/or modify it
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	8 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	10 option) any later version.
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	11
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	12 Octave is distributed in the hope that it will be useful, but WITHOUT
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	15 for more details.
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	16
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	17 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	19 <http://www.gnu.org/licenses/>.
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	20
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	21 */
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	22
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21663 diff changeset	23 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21274 diff changeset	24 # include "config.h"
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	25 #endif
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	26
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	27 #include "CMatrix.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	28 #include "dMatrix.h"
1847 2ffe49eb95a5 [project @ 1996-02-03 12:47:55 by jwe] jwe parents: 1756 diff changeset	29 #include "f77-fcn.h"
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	30 #include "fCMatrix.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	31 #include "fMatrix.h"
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	32 #include "lo-error.h"
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	33 #include "schur.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	34
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	35 typedef octave_idx_type (*double_selector) (const double&, const double&);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	36 typedef octave_idx_type (*float_selector) (const float&, const float&);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	37 typedef octave_idx_type (*complex_selector) (const Complex&);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	38 typedef octave_idx_type (*float_complex_selector) (const FloatComplex&);
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	39
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	40 extern "C"
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	41 {
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4192 diff changeset	42 F77_RET_T
6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4192 diff changeset	43 F77_FUNC (dgeesx, DGEESX) (F77_CONST_CHAR_ARG_DECL,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	44 F77_CONST_CHAR_ARG_DECL,
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	45 double_selector,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	46 F77_CONST_CHAR_ARG_DECL,
22133 59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	47 const F77_INT&, F77_DBLE*,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	48 const F77_INT&, F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	49 F77_DBLE, F77_DBLE, F77_DBLE*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	50 F77_DBLE&, F77_DBLE&, F77_DBLE*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	51 F77_INT*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	52 F77_INT*, F77_INT&
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	53 F77_CHAR_ARG_LEN_DECL
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	54 F77_CHAR_ARG_LEN_DECL
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	55 F77_CHAR_ARG_LEN_DECL);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	56
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	57 F77_RET_T
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	58 F77_FUNC (sgeesx, SGEESX) (F77_CONST_CHAR_ARG_DECL,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	59 F77_CONST_CHAR_ARG_DECL,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	60 float_selector,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	61 F77_CONST_CHAR_ARG_DECL,
22133 59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	62 const F77_INT&, F77_REAL*,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	63 const F77_INT&, F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	64 F77_REAL, F77_REAL, F77_REAL*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	65 F77_REAL&, F77_REAL&, F77_REAL*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	66 F77_INT*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	67 F77_INT*, F77_INT&
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	68 F77_CHAR_ARG_LEN_DECL
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	69 F77_CHAR_ARG_LEN_DECL
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	70 F77_CHAR_ARG_LEN_DECL);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	71
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	72 F77_RET_T
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	73 F77_FUNC (zgeesx, ZGEESX) (F77_CONST_CHAR_ARG_DECL,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	74 F77_CONST_CHAR_ARG_DECL,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	75 complex_selector,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	76 F77_CONST_CHAR_ARG_DECL,
22133 59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	77 const F77_INT&, F77_DBLE_CMPLX*,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	78 const F77_INT&, F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	79 F77_DBLE_CMPLX, F77_DBLE_CMPLX, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	80 F77_DBLE&, F77_DBLE&, F77_DBLE_CMPLX*,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	81 const F77_INT&, F77_DBLE*,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	82 F77_INT*, F77_INT&
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	83 F77_CHAR_ARG_LEN_DECL
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	84 F77_CHAR_ARG_LEN_DECL
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	85 F77_CHAR_ARG_LEN_DECL);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	86
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	87 F77_RET_T
22133 59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	88 F77_FUNC (zrsf2csf, ZRSF2CSF) (const F77_INT&, F77_DBLE_CMPLX *,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	89 F77_DBLE_CMPLX , F77_DBLE , F77_DBLE *);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	90 F77_RET_T
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	91 F77_FUNC (cgeesx, CGEESX) (F77_CONST_CHAR_ARG_DECL,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	92 F77_CONST_CHAR_ARG_DECL,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	93 float_complex_selector,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	94 F77_CONST_CHAR_ARG_DECL,
22133 59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	95 const F77_INT&, F77_CMPLX*,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	96 const F77_INT&, F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	97 F77_CMPLX, F77_CMPLX,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	98 const F77_INT&, F77_REAL&, F77_REAL&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	99 F77_CMPLX*, const F77_INT&,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	100 F77_REAL, F77_INT, F77_INT&
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	101 F77_CHAR_ARG_LEN_DECL
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	102 F77_CHAR_ARG_LEN_DECL
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	103 F77_CHAR_ARG_LEN_DECL);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	104
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	105 F77_RET_T
22133 59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	106 F77_FUNC (crsf2csf, CRSF2CSF) (const F77_INT&, F77_CMPLX *,
59cadee1c74b new macros for F77 data types John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	107 F77_CMPLX , F77_REAL , F77_REAL *);
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	108 }
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	109
22317 6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	110 namespace octave
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	111 {
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	112 namespace math
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	113 {
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	114
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	115 // For real types.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	116
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	117 template <typename T>
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	118 static octave_idx_type
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	119 select_ana (const T& a, const T&)
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	120 {
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	121 return (a < 0.0);
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	122 }
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	123
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	124 template <typename T>
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	125 static octave_idx_type
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	126 select_dig (const T& a, const T& b)
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	127 {
1251 97eac19837dc [project @ 1995-04-11 15:58:32 by jwe] jwe parents: 1192 diff changeset	128 return (hypot (a, b) < 1.0);
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	129 }
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	130
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	131 // For complex types.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	132
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	133 template <typename T>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	134 static octave_idx_type
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	135 select_ana (const T& a)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	136 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	137 return a.real () < 0.0;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	138 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	139
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	140 template <typename T>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	141 static octave_idx_type
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	142 select_dig (const T& a)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	143 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	144 return (abs (a) < 1.0);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	145 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	146
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	147 template <>
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	148 octave_idx_type
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	149 schur<Matrix>::init (const Matrix& a, const std::string& ord, bool calc_unitary)
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	150 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	151 octave_idx_type a_nr = a.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	152 octave_idx_type a_nc = a.cols ();
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	153
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	154 if (a_nr != a_nc)
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	155 (*current_liboctave_error_handler) ("schur: requires square matrix");
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 19697 diff changeset	156
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 19697 diff changeset	157 if (a_nr == 0)
10607 f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	158 {
f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	159 schur_mat.clear ();
f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	160 unitary_mat.clear ();
f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	161 return 0;
f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	162 }
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	163
3334 5187390bfde6 [project @ 1999-11-03 20:41:16 by jwe] jwe parents: 2847 diff changeset	164 // Workspace requirements may need to be fixed if any of the
5187390bfde6 [project @ 1999-11-03 20:41:16 by jwe] jwe parents: 2847 diff changeset	165 // following change.
5187390bfde6 [project @ 1999-11-03 20:41:16 by jwe] jwe parents: 2847 diff changeset	166
5008 c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	167 char jobvs;
1930 d20ab06301e8 [project @ 1996-02-11 22:30:18 by jwe] jwe parents: 1929 diff changeset	168 char sense = 'N';
d20ab06301e8 [project @ 1996-02-11 22:30:18 by jwe] jwe parents: 1929 diff changeset	169 char sort = 'N';
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	170
5008 c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	171 if (calc_unitary)
c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	172 jobvs = 'V';
c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	173 else
c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	174 jobvs = 'N';
c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	175
1756 1af643fa00e3 [project @ 1996-01-22 04:55:17 by jwe] jwe parents: 1631 diff changeset	176 char ord_char = ord.empty () ? 'U' : ord[0];
1af643fa00e3 [project @ 1996-01-22 04:55:17 by jwe] jwe parents: 1631 diff changeset	177
1af643fa00e3 [project @ 1996-01-22 04:55:17 by jwe] jwe parents: 1631 diff changeset	178 if (ord_char == 'A' \|\| ord_char == 'D' \|\| ord_char == 'a' \|\| ord_char == 'd')
1930 d20ab06301e8 [project @ 1996-02-11 22:30:18 by jwe] jwe parents: 1929 diff changeset	179 sort = 'S';
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	180
21663 79414f8e6f24 schur.cc: Declare selector variable as volatile to suppress gcc warning. Rik <rik@octave.org> parents: 21301 diff changeset	181 volatile double_selector selector = 0;
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	182 if (ord_char == 'A' \|\| ord_char == 'a')
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	183 selector = select_ana<double>;
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	184 else if (ord_char == 'D' \|\| ord_char == 'd')
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	185 selector = select_dig<double>;
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	186
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	187 octave_idx_type n = a_nc;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	188 octave_idx_type lwork = 8 * n;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	189 octave_idx_type liwork = 1;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	190 octave_idx_type info;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	191 octave_idx_type sdim;
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	192 double rconde;
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	193 double rcondv;
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	194
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	195 schur_mat = a;
5008 c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	196
c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	197 if (calc_unitary)
10607 f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	198 unitary_mat.clear (n, n);
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	199
908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	200 double *s = schur_mat.fortran_vec ();
908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	201 double *q = unitary_mat.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	202
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	203 Array<double> wr (dim_vector (n, 1));
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	204 double *pwr = wr.fortran_vec ();
908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	205
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	206 Array<double> wi (dim_vector (n, 1));
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	207 double *pwi = wi.fortran_vec ();
908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	208
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	209 Array<double> work (dim_vector (lwork, 1));
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	210 double *pwork = work.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	211
3334 5187390bfde6 [project @ 1999-11-03 20:41:16 by jwe] jwe parents: 2847 diff changeset	212 // BWORK is not referenced for the non-ordered Schur routine.
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	213 octave_idx_type ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	214 Array<octave_idx_type> bwork (dim_vector (ntmp, 1));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	215 octave_idx_type *pbwork = bwork.fortran_vec ();
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	216
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	217 Array<octave_idx_type> iwork (dim_vector (liwork, 1));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	218 octave_idx_type *piwork = iwork.fortran_vec ();
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	219
4552 6f3382e08a52 [project @ 2003-10-27 20:38:02 by jwe] jwe parents: 4192 diff changeset	220 F77_XFCN (dgeesx, DGEESX, (F77_CONST_CHAR_ARG2 (&jobvs, 1),
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	221 F77_CONST_CHAR_ARG2 (&sort, 1),
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	222 selector,
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	223 F77_CONST_CHAR_ARG2 (&sense, 1),
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	224 n, s, n, sdim, pwr, pwi, q, n, rconde, rcondv,
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	225 pwork, lwork, piwork, liwork, pbwork, info
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	226 F77_CHAR_ARG_LEN (1)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	227 F77_CHAR_ARG_LEN (1)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	228 F77_CHAR_ARG_LEN (1)));
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	229
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	230 return info;
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	231 }
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	232
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	233 template <>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	234 octave_idx_type
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	235 schur<FloatMatrix>::init (const FloatMatrix& a, const std::string& ord,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	236 bool calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	237 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	238 octave_idx_type a_nr = a.rows ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	239 octave_idx_type a_nc = a.cols ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	240
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	241 if (a_nr != a_nc)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	242 (*current_liboctave_error_handler) ("SCHUR requires square matrix");
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	243
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	244 if (a_nr == 0)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	245 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	246 schur_mat.clear ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	247 unitary_mat.clear ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	248 return 0;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	249 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	250
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	251 // Workspace requirements may need to be fixed if any of the
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	252 // following change.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	253
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	254 char jobvs;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	255 char sense = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	256 char sort = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	257
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	258 if (calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	259 jobvs = 'V';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	260 else
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	261 jobvs = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	262
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	263 char ord_char = ord.empty () ? 'U' : ord[0];
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	264
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	265 if (ord_char == 'A' \|\| ord_char == 'D' \|\| ord_char == 'a' \|\| ord_char == 'd')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	266 sort = 'S';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	267
21663 79414f8e6f24 schur.cc: Declare selector variable as volatile to suppress gcc warning. Rik <rik@octave.org> parents: 21301 diff changeset	268 volatile float_selector selector = 0;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	269 if (ord_char == 'A' \|\| ord_char == 'a')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	270 selector = select_ana<float>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	271 else if (ord_char == 'D' \|\| ord_char == 'd')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	272 selector = select_dig<float>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	273
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	274 octave_idx_type n = a_nc;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	275 octave_idx_type lwork = 8 * n;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	276 octave_idx_type liwork = 1;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	277 octave_idx_type info;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	278 octave_idx_type sdim;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	279 float rconde;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	280 float rcondv;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	281
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	282 schur_mat = a;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	283
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	284 if (calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	285 unitary_mat.clear (n, n);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	286
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	287 float *s = schur_mat.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	288 float *q = unitary_mat.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	289
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	290 Array<float> wr (dim_vector (n, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	291 float *pwr = wr.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	292
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	293 Array<float> wi (dim_vector (n, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	294 float *pwi = wi.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	295
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	296 Array<float> work (dim_vector (lwork, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	297 float *pwork = work.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	298
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	299 // BWORK is not referenced for the non-ordered Schur routine.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	300 octave_idx_type ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	301 Array<octave_idx_type> bwork (dim_vector (ntmp, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	302 octave_idx_type *pbwork = bwork.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	303
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	304 Array<octave_idx_type> iwork (dim_vector (liwork, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	305 octave_idx_type *piwork = iwork.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	306
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	307 F77_XFCN (sgeesx, SGEESX, (F77_CONST_CHAR_ARG2 (&jobvs, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	308 F77_CONST_CHAR_ARG2 (&sort, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	309 selector,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	310 F77_CONST_CHAR_ARG2 (&sense, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	311 n, s, n, sdim, pwr, pwi, q, n, rconde, rcondv,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	312 pwork, lwork, piwork, liwork, pbwork, info
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	313 F77_CHAR_ARG_LEN (1)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	314 F77_CHAR_ARG_LEN (1)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	315 F77_CHAR_ARG_LEN (1)));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	316
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	317 return info;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	318 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	319
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	320 template <>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	321 octave_idx_type
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	322 schur<ComplexMatrix>::init (const ComplexMatrix& a, const std::string& ord,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	323 bool calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	324 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	325 octave_idx_type a_nr = a.rows ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	326 octave_idx_type a_nc = a.cols ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	327
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	328 if (a_nr != a_nc)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	329 (*current_liboctave_error_handler) ("SCHUR requires square matrix");
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	330
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	331 if (a_nr == 0)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	332 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	333 schur_mat.clear ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	334 unitary_mat.clear ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	335 return 0;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	336 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	337
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	338 // Workspace requirements may need to be fixed if any of the
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	339 // following change.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	340
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	341 char jobvs;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	342 char sense = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	343 char sort = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	344
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	345 if (calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	346 jobvs = 'V';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	347 else
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	348 jobvs = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	349
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	350 char ord_char = ord.empty () ? 'U' : ord[0];
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	351
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	352 if (ord_char == 'A' \|\| ord_char == 'D' \|\| ord_char == 'a' \|\| ord_char == 'd')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	353 sort = 'S';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	354
21663 79414f8e6f24 schur.cc: Declare selector variable as volatile to suppress gcc warning. Rik <rik@octave.org> parents: 21301 diff changeset	355 volatile complex_selector selector = 0;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	356 if (ord_char == 'A' \|\| ord_char == 'a')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	357 selector = select_ana<Complex>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	358 else if (ord_char == 'D' \|\| ord_char == 'd')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	359 selector = select_dig<Complex>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	360
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	361 octave_idx_type n = a_nc;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	362 octave_idx_type lwork = 8 * n;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	363 octave_idx_type info;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	364 octave_idx_type sdim;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	365 double rconde;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	366 double rcondv;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	367
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	368 schur_mat = a;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	369 if (calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	370 unitary_mat.clear (n, n);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	371
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	372 Complex *s = schur_mat.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	373 Complex *q = unitary_mat.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	374
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	375 Array<double> rwork (dim_vector (n, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	376 double *prwork = rwork.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	377
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	378 Array<Complex> w (dim_vector (n, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	379 Complex *pw = w.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	380
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	381 Array<Complex> work (dim_vector (lwork, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	382 Complex *pwork = work.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	383
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	384 // BWORK is not referenced for non-ordered Schur.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	385 octave_idx_type ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	386 Array<octave_idx_type> bwork (dim_vector (ntmp, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	387 octave_idx_type *pbwork = bwork.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	388
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	389 F77_XFCN (zgeesx, ZGEESX, (F77_CONST_CHAR_ARG2 (&jobvs, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	390 F77_CONST_CHAR_ARG2 (&sort, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	391 selector,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	392 F77_CONST_CHAR_ARG2 (&sense, 1),
22135 407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	393 n, F77_DBLE_CMPLX_ARG (s), n, sdim, F77_DBLE_CMPLX_ARG (pw), F77_DBLE_CMPLX_ARG (q), n, rconde, rcondv,
407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	394 F77_DBLE_CMPLX_ARG (pwork), lwork, prwork, pbwork, info
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	395 F77_CHAR_ARG_LEN (1)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	396 F77_CHAR_ARG_LEN (1)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	397 F77_CHAR_ARG_LEN (1)));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	398
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	399 return info;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	400 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	401
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	402 template <>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	403 schur<ComplexMatrix>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	404 rsf2csf<ComplexMatrix, Matrix> (const Matrix& s_arg, const Matrix& u_arg)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	405 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	406 ComplexMatrix s (s_arg);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	407 ComplexMatrix u (u_arg);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	408
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	409 octave_idx_type n = s.rows ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	410
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	411 if (s.columns () != n \|\| u.rows () != n \|\| u.columns () != n)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	412 (*current_liboctave_error_handler)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	413 ("rsf2csf: inconsistent matrix dimensions");
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	414
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	415 if (n > 0)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	416 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	417 OCTAVE_LOCAL_BUFFER (double, c, n-1);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	418 OCTAVE_LOCAL_BUFFER (double, sx, n-1);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	419
22135 407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	420 F77_XFCN (zrsf2csf, ZRSF2CSF, (n, F77_DBLE_CMPLX_ARG (s.fortran_vec ()),
407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	421 F77_DBLE_CMPLX_ARG (u.fortran_vec ()), c, sx));
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	422 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	423
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	424 return schur<ComplexMatrix> (s, u);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	425 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	426
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	427 template <>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	428 octave_idx_type
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	429 schur<FloatComplexMatrix>::init (const FloatComplexMatrix& a,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	430 const std::string& ord, bool calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	431 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	432 octave_idx_type a_nr = a.rows ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	433 octave_idx_type a_nc = a.cols ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	434
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	435 if (a_nr != a_nc)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	436 (*current_liboctave_error_handler) ("SCHUR requires square matrix");
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	437
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	438 if (a_nr == 0)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	439 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	440 schur_mat.clear ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	441 unitary_mat.clear ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	442 return 0;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	443 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	444
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	445 // Workspace requirements may need to be fixed if any of the
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	446 // following change.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	447
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	448 char jobvs;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	449 char sense = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	450 char sort = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	451
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	452 if (calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	453 jobvs = 'V';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	454 else
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	455 jobvs = 'N';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	456
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	457 char ord_char = ord.empty () ? 'U' : ord[0];
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	458
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	459 if (ord_char == 'A' \|\| ord_char == 'D' \|\| ord_char == 'a' \|\| ord_char == 'd')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	460 sort = 'S';
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	461
21663 79414f8e6f24 schur.cc: Declare selector variable as volatile to suppress gcc warning. Rik <rik@octave.org> parents: 21301 diff changeset	462 volatile float_complex_selector selector = 0;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	463 if (ord_char == 'A' \|\| ord_char == 'a')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	464 selector = select_ana<FloatComplex>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	465 else if (ord_char == 'D' \|\| ord_char == 'd')
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	466 selector = select_dig<FloatComplex>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	467
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	468 octave_idx_type n = a_nc;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	469 octave_idx_type lwork = 8 * n;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	470 octave_idx_type info;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	471 octave_idx_type sdim;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	472 float rconde;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	473 float rcondv;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	474
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	475 schur_mat = a;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	476 if (calc_unitary)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	477 unitary_mat.clear (n, n);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	478
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	479 FloatComplex *s = schur_mat.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	480 FloatComplex *q = unitary_mat.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	481
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	482 Array<float> rwork (dim_vector (n, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	483 float *prwork = rwork.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	484
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	485 Array<FloatComplex> w (dim_vector (n, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	486 FloatComplex *pw = w.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	487
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	488 Array<FloatComplex> work (dim_vector (lwork, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	489 FloatComplex *pwork = work.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	490
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	491 // BWORK is not referenced for non-ordered Schur.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	492 octave_idx_type ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	493 Array<octave_idx_type> bwork (dim_vector (ntmp, 1));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	494 octave_idx_type *pbwork = bwork.fortran_vec ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	495
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	496 F77_XFCN (cgeesx, CGEESX, (F77_CONST_CHAR_ARG2 (&jobvs, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	497 F77_CONST_CHAR_ARG2 (&sort, 1),
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	498 selector,
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	499 F77_CONST_CHAR_ARG2 (&sense, 1),
22135 407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	500 n, F77_CMPLX_ARG (s), n, sdim, F77_CMPLX_ARG (pw), F77_CMPLX_ARG (q), n, rconde, rcondv,
407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	501 F77_CMPLX_ARG (pwork), lwork, prwork, pbwork, info
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	502 F77_CHAR_ARG_LEN (1)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	503 F77_CHAR_ARG_LEN (1)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	504 F77_CHAR_ARG_LEN (1)));
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	505
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	506 return info;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	507 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	508
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	509 template <>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	510 schur<FloatComplexMatrix>
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	511 rsf2csf<FloatComplexMatrix, FloatMatrix> (const FloatMatrix& s_arg, const FloatMatrix& u_arg)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	512 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	513 FloatComplexMatrix s (s_arg);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	514 FloatComplexMatrix u (u_arg);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	515
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	516 octave_idx_type n = s.rows ();
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	517
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	518 if (s.columns () != n \|\| u.rows () != n \|\| u.columns () != n)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	519 (*current_liboctave_error_handler)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	520 ("rsf2csf: inconsistent matrix dimensions");
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	521
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	522 if (n > 0)
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	523 {
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	524 OCTAVE_LOCAL_BUFFER (float, c, n-1);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	525 OCTAVE_LOCAL_BUFFER (float, sx, n-1);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	526
22135 407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	527 F77_XFCN (crsf2csf, CRSF2CSF, (n, F77_CMPLX_ARG (s.fortran_vec ()),
407c66ae1e20 reduce warnings from GCC's link-time optimization feature (bug #48531) John W. Eaton <jwe@octave.org> parents: 22133 diff changeset	528 F77_CMPLX_ARG (u.fortran_vec ()), c, sx));
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	529 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	530
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	531 return schur<FloatComplexMatrix> (s, u);
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	532 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	533
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	534 // Instantiations we need.
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	535
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	536 template class schur<ComplexMatrix>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	537
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	538 template class schur<FloatComplexMatrix>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	539
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	540 template class schur<FloatMatrix>;
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	541
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	542 template class schur<Matrix>;
22317 6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	543
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	544 }
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	545 }

Mercurial > octave

annotate liboctave/numeric/schur.cc @ 22317:6ca3acf5fad8