octave: liboctave/numeric/schur.cc annotate

annotate liboctave/numeric/schur.cc @ 31249:de6fc38c78c6

Make Jacobian types offered by dlsode.f accessible by lsode (bug #31626). * liboctave/numeric/LSODE-opts.in: Add options "jacobian type", "lower jacobian subdiagonals", and "upper jacobian subdiagonals". * liboctave/numeric/LSODE.cc (file scope, lsode_j, LSODE::do_integrate (double)): Handle new configurable Jacobian types. * build-aux/mk-opts.pl: Don't implicitly convert to integer in condition.

author	Olaf Till <olaf.till@uni-jena.de>
date	Fri, 12 Nov 2010 08:53:05 +0100
parents	796f54d4ddbf
children	e88a07dec498

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
30564 796f54d4ddbf update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 30073 diff changeset	3 // Copyright (C) 1994-2022 The Octave Project Developers
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	25
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21663 diff changeset	26 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21274 diff changeset	27 # include "config.h"
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	28 #endif
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	29
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	30 #include "Array.h"
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	31 #include "CMatrix.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	32 #include "dMatrix.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	33 #include "fCMatrix.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	34 #include "fMatrix.h"
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	35 #include "lo-error.h"
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22317 diff changeset	36 #include "lo-lapack-proto.h"
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	37 #include "oct-locbuf.h"
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	38 #include "schur.h"
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	39
22317 6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	40 namespace octave
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22135 diff changeset	41 {
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	42 namespace math
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	43 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	44 // For real types.
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	45
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	46 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	47 select_ana (const double& a, const double&)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	48 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	49 return (a < 0.0);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	50 }
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	51
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	52 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	53 select_dig (const double& a, const double& b)
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	54 {
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	55 return (hypot (a, b) < 1.0);
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	56 }
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	57
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	58 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	59 select_ana (const float& a, const float&)
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	60 {
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	61 return (a < 0.0);
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	62 }
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	63
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	64 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	65 select_dig (const float& a, const float& b)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	66 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	67 return (hypot (a, b) < 1.0);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	68 }
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	69
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	70 // For complex types.
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	71
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	72 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	73 select_ana (const F77_DBLE_CMPLX& a_arg)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	74 {
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	75 const Complex a = reinterpret_cast<const Complex&> (a_arg);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	76 return a.real () < 0.0;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	77 }
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	78
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	79 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	80 select_dig (const F77_DBLE_CMPLX& a_arg)
10607 f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	81 {
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	82 const Complex& a = reinterpret_cast<const Complex&> (a_arg);
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	83 return (abs (a) < 1.0);
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	84 }
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	85
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	86 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	87 select_ana (const F77_CMPLX& a_arg)
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	88 {
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	89 const FloatComplex& a = reinterpret_cast<const FloatComplex&> (a_arg);
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	90 return a.real () < 0.0;
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	91 }
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	92
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	93 static F77_INT
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	94 select_dig (const F77_CMPLX& a_arg)
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	95 {
b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	96 const FloatComplex& a = reinterpret_cast<const FloatComplex&> (a_arg);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	97 return (abs (a) < 1.0);
10607 f7501986e42d make schur more Matlab compatible Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	98 }
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	99
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	100 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	101 OCTAVE_API F77_INT
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	102 schur<Matrix>::init (const Matrix& a, const std::string& ord,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	103 bool calc_unitary)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	104 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	105 F77_INT a_nr = to_f77_int (a.rows ());
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	106 F77_INT a_nc = to_f77_int (a.cols ());
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	107
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	108 if (a_nr != a_nc)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	109 (*current_liboctave_error_handler) ("schur: requires square matrix");
3334 5187390bfde6 [project @ 1999-11-03 20:41:16 by jwe] jwe parents: 2847 diff changeset	110
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	111 if (a_nr == 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	112 {
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	113 m_schur_mat.clear ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	114 m_unitary_schur_mat.clear ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	115 return 0;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	116 }
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	117
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	118 // Workspace requirements may need to be fixed if any of the
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	119 // following change.
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	120
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	121 char jobvs;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	122 char sense = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	123 char sort = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	124
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	125 if (calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	126 jobvs = 'V';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	127 else
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	128 jobvs = 'N';
5008 c2bb27ada496 [project @ 2004-09-17 14:45:39 by jwe] jwe parents: 4552 diff changeset	129
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	130 char ord_char = (ord.empty () ? 'U' : ord[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	131
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	132 if (ord_char == 'A' \|\| ord_char == 'D'
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	133 \|\| ord_char == 'a' \|\| ord_char == 'd')
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	134 sort = 'S';
1756 1af643fa00e3 [project @ 1996-01-22 04:55:17 by jwe] jwe parents: 1631 diff changeset	135
23795 980f39c3ab90 Use C++11 nullptr rather than 0 in code (bug #51565). Rik <rik@octave.org> parents: 23696 diff changeset	136 volatile double_selector selector = nullptr;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	137 if (ord_char == 'A' \|\| ord_char == 'a')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	138 selector = select_ana;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	139 else if (ord_char == 'D' \|\| ord_char == 'd')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	140 selector = select_dig;
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	141
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	142 F77_INT n = a_nc;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	143 F77_INT lwork = 8 * n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	144 F77_INT liwork = 1;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	145 F77_INT info;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	146 F77_INT sdim;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	147 double rconde;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	148 double rcondv;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	149
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	150 m_schur_mat = a;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	151
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	152 if (calc_unitary)
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	153 m_unitary_schur_mat.clear (n, n);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	154
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	155 double *s = m_schur_mat.fortran_vec ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	156 double *q = m_unitary_schur_mat.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	157
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	158 Array<double> wr (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	159 double *pwr = wr.fortran_vec ();
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	160
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	161 Array<double> wi (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	162 double *pwi = wi.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	163
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	164 Array<double> work (dim_vector (lwork, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	165 double *pwork = work.fortran_vec ();
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	166
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	167 // BWORK is not referenced for the non-ordered Schur routine.
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	168 F77_INT ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	169 Array<F77_INT> bwork (dim_vector (ntmp, 1));
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	170 F77_INT *pbwork = bwork.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	171
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	172 Array<F77_INT> iwork (dim_vector (liwork, 1));
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	173 F77_INT *piwork = iwork.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	174
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	175 F77_XFCN (dgeesx, DGEESX, (F77_CONST_CHAR_ARG2 (&jobvs, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	176 F77_CONST_CHAR_ARG2 (&sort, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	177 selector,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	178 F77_CONST_CHAR_ARG2 (&sense, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	179 n, s, n, sdim, pwr, pwi, q, n, rconde, rcondv,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	180 pwork, lwork, piwork, liwork, pbwork, info
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	181 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	182 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	183 F77_CHAR_ARG_LEN (1)));
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	184
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	185 return info;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	186 }
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	187
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	188 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	189 OCTAVE_API F77_INT
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	190 schur<FloatMatrix>::init (const FloatMatrix& a, const std::string& ord,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	191 bool calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	192 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	193 F77_INT a_nr = to_f77_int (a.rows ());
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	194 F77_INT a_nc = to_f77_int (a.cols ());
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	195
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	196 if (a_nr != a_nc)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	197 (*current_liboctave_error_handler) ("SCHUR requires square matrix");
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	198
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	199 if (a_nr == 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	200 {
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	201 m_schur_mat.clear ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	202 m_unitary_schur_mat.clear ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	203 return 0;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	204 }
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	205
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	206 // Workspace requirements may need to be fixed if any of the
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	207 // following change.
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	208
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	209 char jobvs;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	210 char sense = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	211 char sort = 'N';
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	212
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	213 if (calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	214 jobvs = 'V';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	215 else
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	216 jobvs = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	217
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	218 char ord_char = (ord.empty () ? 'U' : ord[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	219
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	220 if (ord_char == 'A' \|\| ord_char == 'D'
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	221 \|\| ord_char == 'a' \|\| ord_char == 'd')
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	222 sort = 'S';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	223
23795 980f39c3ab90 Use C++11 nullptr rather than 0 in code (bug #51565). Rik <rik@octave.org> parents: 23696 diff changeset	224 volatile float_selector selector = nullptr;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	225 if (ord_char == 'A' \|\| ord_char == 'a')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	226 selector = select_ana;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	227 else if (ord_char == 'D' \|\| ord_char == 'd')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	228 selector = select_dig;
1929 908f5b6676d7 [project @ 1996-02-11 22:05:08 by jwe] jwe parents: 1882 diff changeset	229
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	230 F77_INT n = a_nc;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	231 F77_INT lwork = 8 * n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	232 F77_INT liwork = 1;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	233 F77_INT info;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	234 F77_INT sdim;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	235 float rconde;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	236 float rcondv;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	237
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	238 m_schur_mat = a;
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	239
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	240 if (calc_unitary)
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	241 m_unitary_schur_mat.clear (n, n);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	242
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	243 float *s = m_schur_mat.fortran_vec ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	244 float *q = m_unitary_schur_mat.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	245
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	246 Array<float> wr (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	247 float *pwr = wr.fortran_vec ();
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	248
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	249 Array<float> wi (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	250 float *pwi = wi.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	251
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	252 Array<float> work (dim_vector (lwork, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	253 float *pwork = work.fortran_vec ();
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	254
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	255 // BWORK is not referenced for the non-ordered Schur routine.
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	256 F77_INT ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	257 Array<F77_INT> bwork (dim_vector (ntmp, 1));
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	258 F77_INT *pbwork = bwork.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	259
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	260 Array<F77_INT> iwork (dim_vector (liwork, 1));
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	261 F77_INT *piwork = iwork.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	262
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	263 F77_XFCN (sgeesx, SGEESX, (F77_CONST_CHAR_ARG2 (&jobvs, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	264 F77_CONST_CHAR_ARG2 (&sort, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	265 selector,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	266 F77_CONST_CHAR_ARG2 (&sense, 1),
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	267 n, s, n, sdim, pwr, pwi, q, n, rconde, rcondv,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	268 pwork, lwork, piwork, liwork, pbwork, info
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	269 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	270 F77_CHAR_ARG_LEN (1)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	271 F77_CHAR_ARG_LEN (1)));
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	272
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	273 return info;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	274 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	275
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	276 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	277 OCTAVE_API F77_INT
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	278 schur<ComplexMatrix>::init (const ComplexMatrix& a, const std::string& ord,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	279 bool calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	280 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	281 F77_INT a_nr = to_f77_int (a.rows ());
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	282 F77_INT a_nc = to_f77_int (a.cols ());
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	283
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	284 if (a_nr != a_nc)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	285 (*current_liboctave_error_handler) ("SCHUR requires square matrix");
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	286
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	287 if (a_nr == 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	288 {
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	289 m_schur_mat.clear ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	290 m_unitary_schur_mat.clear ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	291 return 0;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	292 }
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	293
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	294 // Workspace requirements may need to be fixed if any of the
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	295 // following change.
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	296
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	297 char jobvs;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	298 char sense = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	299 char sort = 'N';
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	300
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	301 if (calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	302 jobvs = 'V';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	303 else
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	304 jobvs = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	305
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	306 char ord_char = (ord.empty () ? 'U' : ord[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	307
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	308 if (ord_char == 'A' \|\| ord_char == 'D'
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	309 \|\| ord_char == 'a' \|\| ord_char == 'd')
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	310 sort = 'S';
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	311
23795 980f39c3ab90 Use C++11 nullptr rather than 0 in code (bug #51565). Rik <rik@octave.org> parents: 23696 diff changeset	312 volatile complex_selector selector = nullptr;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	313 if (ord_char == 'A' \|\| ord_char == 'a')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	314 selector = select_ana;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	315 else if (ord_char == 'D' \|\| ord_char == 'd')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	316 selector = select_dig;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	317
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	318 F77_INT n = a_nc;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	319 F77_INT lwork = 8 * n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	320 F77_INT info;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	321 F77_INT sdim;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	322 double rconde;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	323 double rcondv;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	324
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	325 m_schur_mat = a;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	326 if (calc_unitary)
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	327 m_unitary_schur_mat.clear (n, n);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	328
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	329 Complex *s = m_schur_mat.fortran_vec ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	330 Complex *q = m_unitary_schur_mat.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	331
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	332 Array<double> rwork (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	333 double *prwork = rwork.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	334
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	335 Array<Complex> w (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	336 Complex *pw = w.fortran_vec ();
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	337
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	338 Array<Complex> work (dim_vector (lwork, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	339 Complex *pwork = work.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	340
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	341 // BWORK is not referenced for non-ordered Schur.
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	342 F77_INT ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	343 Array<F77_INT> bwork (dim_vector (ntmp, 1));
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	344 F77_INT *pbwork = bwork.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	345
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	346 F77_XFCN (zgeesx, ZGEESX,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	347 (F77_CONST_CHAR_ARG2 (&jobvs, 1),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	348 F77_CONST_CHAR_ARG2 (&sort, 1),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	349 selector,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	350 F77_CONST_CHAR_ARG2 (&sense, 1),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	351 n, F77_DBLE_CMPLX_ARG (s), n, sdim, F77_DBLE_CMPLX_ARG (pw),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	352 F77_DBLE_CMPLX_ARG (q), n, rconde, rcondv,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	353 F77_DBLE_CMPLX_ARG (pwork), lwork, prwork, pbwork, info
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	354 F77_CHAR_ARG_LEN (1)
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	355 F77_CHAR_ARG_LEN (1)
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	356 F77_CHAR_ARG_LEN (1)));
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	357
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	358 return info;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	359 }
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	360
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	361 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	362 OCTAVE_API schur<ComplexMatrix>
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	363 rsf2csf<ComplexMatrix, Matrix> (const Matrix& s_arg, const Matrix& u_arg)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	364 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	365 ComplexMatrix s (s_arg);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	366 ComplexMatrix u (u_arg);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	367
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	368 F77_INT n = to_f77_int (s.rows ());
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	369
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	370 if (s.columns () != n \|\| u.rows () != n \|\| u.columns () != n)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	371 (*current_liboctave_error_handler)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	372 ("rsf2csf: inconsistent matrix dimensions");
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	373
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	374 if (n > 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	375 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	376 OCTAVE_LOCAL_BUFFER (double, c, n-1);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	377 OCTAVE_LOCAL_BUFFER (double, sx, n-1);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	378
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	379 F77_XFCN (zrsf2csf, ZRSF2CSF,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	380 (n, F77_DBLE_CMPLX_ARG (s.fortran_vec ()),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	381 F77_DBLE_CMPLX_ARG (u.fortran_vec ()), c, sx));
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	382 }
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	383
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	384 return schur<ComplexMatrix> (s, u);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	385 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	386
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	387 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	388 OCTAVE_API F77_INT
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	389 schur<FloatComplexMatrix>::init (const FloatComplexMatrix& a,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	390 const std::string& ord, bool calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	391 {
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	392 F77_INT a_nr = to_f77_int (a.rows ());
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	393 F77_INT a_nc = to_f77_int (a.cols ());
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	394
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	395 if (a_nr != a_nc)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	396 (*current_liboctave_error_handler) ("SCHUR requires square matrix");
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	397
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	398 if (a_nr == 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	399 {
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	400 m_schur_mat.clear ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	401 m_unitary_schur_mat.clear ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	402 return 0;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	403 }
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	404
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	405 // Workspace requirements may need to be fixed if any of the
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	406 // following change.
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	407
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	408 char jobvs;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	409 char sense = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	410 char sort = 'N';
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	411
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	412 if (calc_unitary)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	413 jobvs = 'V';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	414 else
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	415 jobvs = 'N';
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	416
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23220 diff changeset	417 char ord_char = (ord.empty () ? 'U' : ord[0]);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	418
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	419 if (ord_char == 'A' \|\| ord_char == 'D'
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	420 \|\| ord_char == 'a' \|\| ord_char == 'd')
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	421 sort = 'S';
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	422
23795 980f39c3ab90 Use C++11 nullptr rather than 0 in code (bug #51565). Rik <rik@octave.org> parents: 23696 diff changeset	423 volatile float_complex_selector selector = nullptr;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	424 if (ord_char == 'A' \|\| ord_char == 'a')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	425 selector = select_ana;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	426 else if (ord_char == 'D' \|\| ord_char == 'd')
22839 b27125252bb0 use Fortran data types for Fortran function pointer typedefs John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	427 selector = select_dig;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	428
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	429 F77_INT n = a_nc;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	430 F77_INT lwork = 8 * n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	431 F77_INT info;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	432 F77_INT sdim;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	433 float rconde;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	434 float rcondv;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	435
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	436 m_schur_mat = a;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	437 if (calc_unitary)
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	438 m_unitary_schur_mat.clear (n, n);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	439
30060 ae47678c76e8 maint: use "m_" prefix for member variables in class schur. Rik <rik@octave.org> parents: 29359 diff changeset	440 FloatComplex *s = m_schur_mat.fortran_vec ();
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	441 FloatComplex *q = m_unitary_schur_mat.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	442
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	443 Array<float> rwork (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	444 float *prwork = rwork.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	445
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	446 Array<FloatComplex> w (dim_vector (n, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	447 FloatComplex *pw = w.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	448
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	449 Array<FloatComplex> work (dim_vector (lwork, 1));
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	450 FloatComplex *pwork = work.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	451
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	452 // BWORK is not referenced for non-ordered Schur.
22844 7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	453 F77_INT ntmp = (ord_char == 'N' \|\| ord_char == 'n') ? 0 : n;
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	454 Array<F77_INT> bwork (dim_vector (ntmp, 1));
7bb977866f55 use F77_INT instead of octave_idx_type for liboctave schur class John W. Eaton <jwe@octave.org> parents: 22839 diff changeset	455 F77_INT *pbwork = bwork.fortran_vec ();
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	456
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	457 F77_XFCN (cgeesx, CGEESX,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	458 (F77_CONST_CHAR_ARG2 (&jobvs, 1),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	459 F77_CONST_CHAR_ARG2 (&sort, 1),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	460 selector,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	461 F77_CONST_CHAR_ARG2 (&sense, 1),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	462 n, F77_CMPLX_ARG (s), n, sdim, F77_CMPLX_ARG (pw),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	463 F77_CMPLX_ARG (q), n,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	464 rconde, rcondv,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	465 F77_CMPLX_ARG (pwork), lwork, prwork, pbwork, info
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	466 F77_CHAR_ARG_LEN (1)
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	467 F77_CHAR_ARG_LEN (1)
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	468 F77_CHAR_ARG_LEN (1)));
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	469
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	470 return info;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	471 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	472
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	473 template <>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	474 OCTAVE_API schur<FloatComplexMatrix>
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22329 diff changeset	475 rsf2csf<FloatComplexMatrix, FloatMatrix> (const FloatMatrix& s_arg,
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22329 diff changeset	476 const FloatMatrix& u_arg)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	477 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	478 FloatComplexMatrix s (s_arg);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	479 FloatComplexMatrix u (u_arg);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	480
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	481 F77_INT n = to_f77_int (s.rows ());
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	482
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	483 if (s.columns () != n \|\| u.rows () != n \|\| u.columns () != n)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	484 (*current_liboctave_error_handler)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	485 ("rsf2csf: inconsistent matrix dimensions");
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	486
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	487 if (n > 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	488 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	489 OCTAVE_LOCAL_BUFFER (float, c, n-1);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	490 OCTAVE_LOCAL_BUFFER (float, sx, n-1);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	491
30073 1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	492 F77_XFCN (crsf2csf, CRSF2CSF,
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	493 (n, F77_CMPLX_ARG (s.fortran_vec ()),
1e277c6b6626 maint: rename schur member function from "unitary_matrix" to "unitary_schur_matrix" Rik <rik@octave.org> parents: 30060 diff changeset	494 F77_CMPLX_ARG (u.fortran_vec ()), c, sx));
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	495 }
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	496
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	497 return schur<FloatComplexMatrix> (s, u);
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	498 }
e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	499
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	500 // Instantiations we need.
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	501
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	502 template class schur<ComplexMatrix>;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	503
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	504 template class schur<FloatComplexMatrix>;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	505
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	506 template class schur<FloatMatrix>;
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	507
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	508 template class schur<Matrix>;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	509 }
21266 e69eaee28737 make better use of templates for Schur decomposition John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	510 }

Mercurial > octave

annotate liboctave/numeric/schur.cc @ 31249:de6fc38c78c6