octave: liboctave/numeric/gsvd.cc annotate

author	Barbara Locsi <locsi.barbara@gmail.com>
date	Tue, 09 Aug 2016 18:02:11 +0200
parents
children	da201af35c97

rev	line source
22236 065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	1 // Copyright (C) 1996, 1997 John W. Eaton
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	2 // Copyright (C) 2006 Pascal Dupuis <Pascal.Dupuis@uclouvain.be>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	3 // Copyright (C) 2016 Barbara Lócsi
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	4 //
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	5 // This program is free software; you can redistribute it and/or modify it under
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	6 // the terms of the GNU General Public License as published by the Free Software
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	7 // Foundation; either version 3 of the License, or (at your option) any later
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	8 // version.
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	9 //
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	10 // This program is distributed in the hope that it will be useful, but WITHOUT
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	11 // ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	12 // FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	13 // details.
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	14 //
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	15 // You should have received a copy of the GNU General Public License along with
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	16 // this program; if not, see <http://www.gnu.org/licenses/>.
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	17
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	18 #ifdef HAVE_CONFIG_H
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	19 # include <config.h>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	20 #endif
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	21
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	22 #include "gsvd.h"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	23 #include "f77-fcn.h"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	24 #include "lo-error.h"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	25 #include "CMatrix.h"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	26 #include "dDiagMatrix.h"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	27 #include "dMatrix.h"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	28
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	29 extern "C"
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	30 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	31 F77_RET_T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	32 F77_FUNC (dggsvd, DGGSVD)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	33 (
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	34 F77_CONST_CHAR_ARG_DECL, // JOBU (input) CHARACTER*1
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	35 F77_CONST_CHAR_ARG_DECL, // JOBV (input) CHARACTER*1
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	36 F77_CONST_CHAR_ARG_DECL, // JOBQ (input) CHARACTER*1
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	37 const F77_INT&, // M (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	38 const F77_INT&, // N (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	39 const F77_INT&, // P (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	40 F77_INT &, // K (output) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	41 F77_INT &, // L (output) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	42 F77_DBLE*, // A (input/output) DOUBLE PRECISION array, dimension (LDA,N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	43 const F77_INT&, // LDA (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	44 F77_DBLE*, // B (input/output) DOUBLE PRECISION array, dimension (LDB,N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	45 const F77_INT&, // LDB (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	46 F77_DBLE*, // ALPHA (output) DOUBLE PRECISION array, dimension (N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	47 F77_DBLE*, // BETA (output) DOUBLE PRECISION array, dimension (N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	48 F77_DBLE*, // U (output) DOUBLE PRECISION array, dimension (LDU,M)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	49 const F77_INT&, // LDU (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	50 F77_DBLE*, // V (output) DOUBLE PRECISION array, dimension (LDV,P)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	51 const F77_INT&, // LDV (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	52 F77_DBLE*, // Q (output) DOUBLE PRECISION array, dimension (LDQ,N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	53 const F77_INT&, // LDQ (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	54 F77_DBLE*, // WORK (workspace) DOUBLE PRECISION array
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	55 int*, // IWORK (workspace/output) INTEGER array, dimension (N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	56 F77_INT& // INFO (output)INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	57 F77_CHAR_ARG_LEN_DECL
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	58 F77_CHAR_ARG_LEN_DECL
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	59 F77_CHAR_ARG_LEN_DECL
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	60 );
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	61
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	62 F77_RET_T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	63 F77_FUNC (zggsvd, ZGGSVD)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	64 (
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	65 F77_CONST_CHAR_ARG_DECL, // JOBU (input) CHARACTER*1
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	66 F77_CONST_CHAR_ARG_DECL, // JOBV (input) CHARACTER*1
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	67 F77_CONST_CHAR_ARG_DECL, // JOBQ (input) CHARACTER*1
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	68 const F77_INT&, // M (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	69 const F77_INT&, // N (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	70 const F77_INT&, // P (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	71 F77_INT &, // K (output) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	72 F77_INT &, // L (output) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	73 F77_DBLE_CMPLX, // A (input/output) COMPLEX16 array, dimension (LDA,N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	74 const F77_INT&, // LDA (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	75 F77_DBLE_CMPLX, // B (input/output) COMPLEX16 array, dimension (LDB,N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	76 const F77_INT&, // LDB (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	77 F77_DBLE*, // ALPHA (output) DOUBLE PRECISION array, dimension (N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	78 F77_DBLE*, // BETA (output) DOUBLE PRECISION array, dimension (N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	79 F77_DBLE_CMPLX, // U (output) COMPLEX16 array, dimension (LDU,M)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	80 const F77_INT&, // LDU (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	81 F77_DBLE_CMPLX, // V (output) COMPLEX16 array, dimension (LDV,P)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	82 const F77_INT&, // LDV (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	83 F77_DBLE_CMPLX, // Q (output) COMPLEX16 array, dimension (LDQ,N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	84 const F77_INT&, // LDQ (input) INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	85 F77_DBLE_CMPLX, // WORK (workspace) COMPLEX16 array
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	86 F77_DBLE*, // RWORK (workspace) DOUBLE PRECISION array
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	87 int*, // IWORK (workspace/output) INTEGER array, dimension (N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	88 F77_INT& // INFO (output)INTEGER
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	89 F77_CHAR_ARG_LEN_DECL
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	90 F77_CHAR_ARG_LEN_DECL
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	91 F77_CHAR_ARG_LEN_DECL
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	92 );
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	93 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	94
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	95 template <>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	96 void
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	97 gsvd<Matrix>::ggsvd (char& jobu, char& jobv, char& jobq, octave_idx_type m,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	98 octave_idx_type n, octave_idx_type p, octave_idx_type& k,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	99 octave_idx_type& l, double *tmp_dataA, octave_idx_type m1,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	100 double *tmp_dataB, octave_idx_type p1, Matrix& alpha,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	101 Matrix& beta, double u, octave_idx_type nrow_u, double v,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	102 octave_idx_type nrow_v, double *q, octave_idx_type nrow_q,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	103 Matrix& work, octave_idx_type* iwork, octave_idx_type& info)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	104 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	105 F77_XFCN (dggsvd, DGGSVD, (F77_CONST_CHAR_ARG2 (&jobu, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	106 F77_CONST_CHAR_ARG2 (&jobv, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	107 F77_CONST_CHAR_ARG2 (&jobq, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	108 m, n, p, k, l, tmp_dataA, m1,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	109 tmp_dataB, p1, alpha.fortran_vec (),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	110 beta.fortran_vec (), u, nrow_u,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	111 v, nrow_v, q, nrow_q, work.fortran_vec (),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	112 iwork, info
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	113 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	114 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	115 F77_CHAR_ARG_LEN (1)));
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	116 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	117
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	118
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	119 template <>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	120 void
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	121 gsvd<ComplexMatrix>::ggsvd (char& jobu, char& jobv, char& jobq,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	122 octave_idx_type m, octave_idx_type n,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	123 octave_idx_type p, octave_idx_type& k,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	124 octave_idx_type& l, Complex *tmp_dataA,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	125 octave_idx_type m1, Complex *tmp_dataB,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	126 octave_idx_type p1, Matrix& alpha,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	127 Matrix& beta, Complex *u, octave_idx_type nrow_u,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	128 Complex v, octave_idx_type nrow_v, Complex q,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	129 octave_idx_type nrow_q, ComplexMatrix& work,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	130 octave_idx_type* iwork, octave_idx_type& info)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	131 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	132 Matrix rwork(2*n, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	133 F77_XFCN (zggsvd, ZGGSVD, (F77_CONST_CHAR_ARG2 (&jobu, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	134 F77_CONST_CHAR_ARG2 (&jobv, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	135 F77_CONST_CHAR_ARG2 (&jobq, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	136 m, n, p, k, l, F77_DBLE_CMPLX_ARG (tmp_dataA), m1,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	137 F77_DBLE_CMPLX_ARG (tmp_dataB), p1,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	138 alpha.fortran_vec (), beta.fortran_vec (),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	139 F77_DBLE_CMPLX_ARG (u), nrow_u,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	140 F77_DBLE_CMPLX_ARG (v), nrow_v,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	141 F77_DBLE_CMPLX_ARG (q), nrow_q,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	142 F77_DBLE_CMPLX_ARG (work.fortran_vec ()),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	143 rwork.fortran_vec (), iwork, info
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	144 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	145 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	146 F77_CHAR_ARG_LEN (1)));
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	147 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	148
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	149 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	150 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	151 gsvd<T>::left_singular_matrix_A (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	152 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	153 if (type == gsvd::Type::sigma_only)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	154 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	155 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	156 ("gsvd: U not computed because type == gsvd::sigma_only");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	157 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	158 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	159 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	160 return left_smA;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	161 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	162
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	163 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	164 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	165 gsvd<T>::left_singular_matrix_B (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	166 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	167 if (type == gsvd::Type::sigma_only)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	168 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	169 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	170 ("gsvd: V not computed because type == gsvd::sigma_only");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	171 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	172 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	173 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	174 return left_smB;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	175 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	176
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	177 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	178 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	179 gsvd<T>::right_singular_matrix (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	180 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	181 if (type == gsvd::Type::sigma_only)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	182 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	183 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	184 ("gsvd: X not computed because type == gsvd::sigma_only");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	185 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	186 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	187 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	188 return right_sm;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	189 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	190
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	191 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	192 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	193 gsvd<T>::R_matrix (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	194 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	195 if (type != gsvd::Type::std)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	196 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	197 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	198 ("gsvd: R not computed because type != gsvd::std");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	199 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	200 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	201 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	202 return R;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	203 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	204
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	205 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	206 gsvd<T>::gsvd (const T& a, const T& b, gsvd::Type gsvd_type)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	207 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	208 octave_idx_type info;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	209
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	210 octave_idx_type m = a.rows ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	211 octave_idx_type n = a.cols ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	212 octave_idx_type p = b.rows ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	213
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	214 T atmp = a;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	215 P *tmp_dataA = atmp.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	216
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	217 T btmp = b;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	218 P *tmp_dataB = btmp.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	219
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	220 char jobu = 'U';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	221 char jobv = 'V';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	222 char jobq = 'Q';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	223
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	224 octave_idx_type nrow_u = m;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	225 octave_idx_type nrow_v = p;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	226 octave_idx_type nrow_q = n;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	227
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	228 octave_idx_type k, l;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	229
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	230 switch (gsvd_type)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	231 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	232 case gsvd<T>::Type::sigma_only:
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	233
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	234 // Note: for this case, both jobu and jobv should be 'N', but
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	235 // there seems to be a bug in dgesvd from Lapack V2.0. To
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	236 // demonstrate the bug, set both jobu and jobv to 'N' and find
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	237 // the singular values of [eye(3), eye(3)]. The result is
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	238 // [-sqrt(2), -sqrt(2), -sqrt(2)].
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	239 //
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	240 // For Lapack 3.0, this problem seems to be fixed.
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	241
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	242 jobu = 'N';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	243 jobv = 'N';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	244 jobq = 'N';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	245 nrow_u = nrow_v = nrow_q = 1;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	246 break;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	247
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	248 default:
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	249 break;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	250 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	251
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	252 type = gsvd_type;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	253
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	254 if (! (jobu == 'N' \|\| jobu == 'O'))
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	255 left_smA.resize (nrow_u, m);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	256
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	257 P *u = left_smA.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	258
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	259 if (! (jobv == 'N' \|\| jobv == 'O'))
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	260 left_smB.resize (nrow_v, p);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	261
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	262 P *v = left_smB.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	263
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	264 if (! (jobq == 'N' \|\| jobq == 'O'))
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	265 right_sm.resize (nrow_q, n);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	266
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	267 P *q = right_sm.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	268
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	269 octave_idx_type lwork = 3*n;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	270 lwork = lwork > m ? lwork : m;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	271 lwork = (lwork > p ? lwork : p) + n;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	272
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	273 T work (lwork, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	274 Matrix alpha (n, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	275 Matrix beta (n, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	276
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	277 std::vector<octave_idx_type> iwork (n);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	278
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	279 gsvd<T>::ggsvd (jobu, jobv, jobq, m, n, p, k, l,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	280 tmp_dataA, m, tmp_dataB, p, alpha, beta, u,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	281 nrow_u, v, nrow_v, q, nrow_q, work, iwork.data (), info);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	282
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	283 if (f77_exception_encountered)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	284 (current_liboctave_error_handler) ("unrecoverable error in ggsvd");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	285
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	286 if (info < 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	287 (current_liboctave_error_handler) ("ggsvd.f: argument %d illegal", -info);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	288 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	289 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	290 if (info > 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	291 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	292 ("*ggsvd.f: Jacobi-type procedure failed to converge.");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	293 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	294 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	295 octave_idx_type i, j;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	296
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	297 if (gsvd<T>::Type::std == gsvd_type)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	298 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	299 R.resize(k+l, k+l);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	300 int astart = n-k-l;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	301 if (m - k - l >= 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	302 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	303 astart = n-k-l;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	304 // R is stored in A(1:K+L,N-K-L+1:N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	305 for (i = 0; i < k+l; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	306 for (j = 0; j < k+l; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	307 R.xelem (i, j) = atmp.xelem (i, astart + j);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	308 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	309 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	310 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	311 // (R11 R12 R13 ) is stored in A(1:M, N-K-L+1:N),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	312 // ( 0 R22 R23 )
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	313
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	314 for (i = 0; i < m; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	315 for (j = 0; j < k+l; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	316 R.xelem (i, j) = atmp.xelem (i, astart + j);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	317 // and R33 is stored in B(M-K+1:L,N+M-K-L+1:N)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	318 for (i = k+l-1; i >=m; i--) {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	319 for (j = 0; j < m; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	320 R.xelem(i, j) = 0.0;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	321 for (j = m; j < k+l; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	322 R.xelem (i, j) = btmp.xelem (i - k, astart + j);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	323 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	324 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	325 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	326
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	327 if (m-k-l >= 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	328 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	329 // Fills in C and S
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	330 sigmaA.resize (l, l);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	331 sigmaB.resize (l, l);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	332 for (i = 0; i < l; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	333 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	334 sigmaA.dgxelem(i) = alpha.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	335 sigmaB.dgxelem(i) = beta.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	336 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	337 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	338 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	339 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	340 // Fills in C and S
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	341 sigmaA.resize (m-k, m-k);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	342 sigmaB.resize (m-k, m-k);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	343 for (i = 0; i < m-k; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	344 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	345 sigmaA.dgxelem(i) = alpha.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	346 sigmaB.dgxelem(i) = beta.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	347 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	348 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	349 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	350 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	351 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	352
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	353 // Instantiations we need.
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	354
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	355 template class gsvd<Matrix>;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	356
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	357 template class gsvd<ComplexMatrix>;

22236

065a44375723 gsvd: reduce code duplication with templates.