octave: liboctave/numeric/gsvd.cc annotate

annotate liboctave/numeric/gsvd.cc @ 22322:93b3cdd36854

move most f77 function decls to separate header files * liboctave/numeric/lo-amos-proto.h, liboctave/numeric/lo-arpack-proto.h, liboctave/numeric/lo-blas-proto.h, liboctave/numeric/lo-fftpack-proto.h, liboctave/numeric/lo-lapack-proto.h, liboctave/numeric/lo-qrupdate-proto.h, liboctave/numeric/lo-ranlib-proto.h, liboctave/numeric/lo-slatec-proto.h: New files. * liboctave/numeric/module.mk: Update. * __pchip_deriv__.cc, dot.cc, interpreter.cc, ordschur.cc, qz.cc, CColVector.cc, CMatrix.cc, CNDArray.cc, CRowVector.cc, CSparse.cc, dColVector.cc, dMatrix.cc, dNDArray.cc, dRowVector.cc, dSparse.cc, fCColVector.cc, fCMatrix.cc, fCNDArray.cc, fCRowVector.cc, fColVector.cc, fMatrix.cc, fNDArray.cc, fRowVector.cc, EIG.cc, aepbalance.cc, chol.cc, eigs-base.cc, fEIG.cc, gepbalance.cc, gsvd.cc, hess.cc, lo-specfun.cc, lu.cc, oct-rand.cc, qr.cc, qrp.cc, randpoisson.cc, schur.cc, sparse-qr.cc, svd.cc: Use new header files.

author	John W. Eaton <jwe@octave.org>
date	Wed, 17 Aug 2016 00:18:08 -0400
parents	bc2a5db96754
children	bac0d6f07a3e

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 "lo-error.h"
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22297 diff changeset	24 #include "lo-lapack-proto.h"
22236 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
22249 da201af35c97 maint: properly include <vector> in svd classes. Kai T. Ohlhus <k.ohlhus@gmail.com> parents: 22236 diff changeset	29 #include <vector>
da201af35c97 maint: properly include <vector> in svd classes. Kai T. Ohlhus <k.ohlhus@gmail.com> parents: 22236 diff changeset	30
22236 065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	31 template <>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	32 void
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	33 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	34 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	35 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	36 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	37 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	38 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	39 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	40 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	41 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	42 F77_CONST_CHAR_ARG2 (&jobv, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	43 F77_CONST_CHAR_ARG2 (&jobq, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	44 m, n, p, k, l, tmp_dataA, m1,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	45 tmp_dataB, p1, alpha.fortran_vec (),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	46 beta.fortran_vec (), u, nrow_u,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	47 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	48 iwork, info
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	49 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	50 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	51 F77_CHAR_ARG_LEN (1)));
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	52 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	53
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	54
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	55 template <>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	56 void
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	57 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	58 octave_idx_type m, octave_idx_type n,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	59 octave_idx_type p, octave_idx_type& k,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	60 octave_idx_type& l, Complex *tmp_dataA,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	61 octave_idx_type m1, Complex *tmp_dataB,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	62 octave_idx_type p1, Matrix& alpha,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	63 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	64 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	65 octave_idx_type nrow_q, ComplexMatrix& work,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	66 octave_idx_type* iwork, octave_idx_type& info)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	67 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	68 Matrix rwork(2*n, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	69 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	70 F77_CONST_CHAR_ARG2 (&jobv, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	71 F77_CONST_CHAR_ARG2 (&jobq, 1),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	72 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	73 F77_DBLE_CMPLX_ARG (tmp_dataB), p1,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	74 alpha.fortran_vec (), beta.fortran_vec (),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	75 F77_DBLE_CMPLX_ARG (u), nrow_u,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	76 F77_DBLE_CMPLX_ARG (v), nrow_v,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	77 F77_DBLE_CMPLX_ARG (q), nrow_q,
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	78 F77_DBLE_CMPLX_ARG (work.fortran_vec ()),
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	79 rwork.fortran_vec (), iwork, info
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	80 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	81 F77_CHAR_ARG_LEN (1)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	82 F77_CHAR_ARG_LEN (1)));
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	83 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	84
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	85 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	86 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	87 gsvd<T>::left_singular_matrix_A (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	88 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	89 if (type == gsvd::Type::sigma_only)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	90 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	91 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	92 ("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	93 return T ();
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 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	96 return left_smA;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	97 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	98
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	99 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	100 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	101 gsvd<T>::left_singular_matrix_B (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	102 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	103 if (type == gsvd::Type::sigma_only)
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 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	106 ("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	107 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	108 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	109 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	110 return left_smB;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	111 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	112
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	113 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	114 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	115 gsvd<T>::right_singular_matrix (void) const
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 if (type == gsvd::Type::sigma_only)
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 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	120 ("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	121 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	122 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	123 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	124 return right_sm;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	125 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	126
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	127 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	128 T
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	129 gsvd<T>::R_matrix (void) const
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	130 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	131 if (type != gsvd::Type::std)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	132 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	133 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	134 ("gsvd: R not computed because type != gsvd::std");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	135 return T ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	136 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	137 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	138 return R;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	139 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	140
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	141 template <typename T>
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	142 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	143 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	144 octave_idx_type info;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	145
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	146 octave_idx_type m = a.rows ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	147 octave_idx_type n = a.cols ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	148 octave_idx_type p = b.rows ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	149
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	150 T atmp = a;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	151 P *tmp_dataA = atmp.fortran_vec ();
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 T btmp = b;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	154 P *tmp_dataB = btmp.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	155
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	156 char jobu = 'U';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	157 char jobv = 'V';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	158 char jobq = 'Q';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	159
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	160 octave_idx_type nrow_u = m;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	161 octave_idx_type nrow_v = p;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	162 octave_idx_type nrow_q = n;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	163
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	164 octave_idx_type k, l;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	165
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	166 switch (gsvd_type)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	167 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	168 case gsvd<T>::Type::sigma_only:
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	169
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	170 // 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	171 // 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	172 // 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	173 // 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	174 // [-sqrt(2), -sqrt(2), -sqrt(2)].
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 // 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	177
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	178 jobu = 'N';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	179 jobv = 'N';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	180 jobq = 'N';
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	181 nrow_u = nrow_v = nrow_q = 1;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	182 break;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	183
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	184 default:
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	185 break;
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
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	188 type = gsvd_type;
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 if (! (jobu == 'N' \|\| jobu == 'O'))
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	191 left_smA.resize (nrow_u, m);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	192
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	193 P *u = left_smA.fortran_vec ();
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 (! (jobv == 'N' \|\| jobv == 'O'))
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	196 left_smB.resize (nrow_v, p);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	197
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	198 P *v = left_smB.fortran_vec ();
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	199
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	200 if (! (jobq == 'N' \|\| jobq == 'O'))
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	201 right_sm.resize (nrow_q, n);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	202
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	203 P *q = right_sm.fortran_vec ();
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 octave_idx_type lwork = 3*n;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	206 lwork = lwork > m ? lwork : m;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	207 lwork = (lwork > p ? lwork : p) + n;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	208
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	209 T work (lwork, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	210 Matrix alpha (n, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	211 Matrix beta (n, 1);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	212
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	213 std::vector<octave_idx_type> iwork (n);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	214
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	215 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	216 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	217 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	218
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	219 if (f77_exception_encountered)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	220 (current_liboctave_error_handler) ("unrecoverable error in ggsvd");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	221
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	222 if (info < 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	223 (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	224 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	225 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	226 if (info > 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	227 (*current_liboctave_error_handler)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	228 ("*ggsvd.f: Jacobi-type procedure failed to converge.");
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	229 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	230 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	231 octave_idx_type i, j;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	232
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	233 if (gsvd<T>::Type::std == gsvd_type)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	234 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	235 R.resize(k+l, k+l);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	236 int astart = n-k-l;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	237 if (m - k - l >= 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	238 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	239 astart = n-k-l;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	240 // 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	241 for (i = 0; i < k+l; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	242 for (j = 0; j < k+l; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	243 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	244 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	245 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	246 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	247 // (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	248 // ( 0 R22 R23 )
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	249
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	250 for (i = 0; i < m; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	251 for (j = 0; j < k+l; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	252 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	253 // 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	254 for (i = k+l-1; i >=m; i--) {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	255 for (j = 0; j < m; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	256 R.xelem(i, j) = 0.0;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	257 for (j = m; j < k+l; j++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	258 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	259 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	260 }
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
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	263 if (m-k-l >= 0)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	264 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	265 // Fills in C and S
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	266 sigmaA.resize (l, l);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	267 sigmaB.resize (l, l);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	268 for (i = 0; i < l; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	269 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	270 sigmaA.dgxelem(i) = alpha.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	271 sigmaB.dgxelem(i) = beta.elem(k+i);
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 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	274 else
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	275 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	276 // Fills in C and S
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	277 sigmaA.resize (m-k, m-k);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	278 sigmaB.resize (m-k, m-k);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	279 for (i = 0; i < m-k; i++)
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	280 {
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	281 sigmaA.dgxelem(i) = alpha.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	282 sigmaB.dgxelem(i) = beta.elem(k+i);
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	283 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	284 }
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 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	287 }
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	288
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	289 // Instantiations we need.
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	290
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	291 template class gsvd<Matrix>;
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	292
065a44375723 gsvd: reduce code duplication with templates. Barbara Locsi <locsi.barbara@gmail.com> parents: diff changeset	293 template class gsvd<ComplexMatrix>;

Mercurial > octave

annotate liboctave/numeric/gsvd.cc @ 22322:93b3cdd36854