octave: liboctave/numeric/svd.cc annotate

annotate liboctave/numeric/svd.cc @ 29416:3f45f5b6fa7f

svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). * liboctave/numeric/svd.cc (gejsv_lwork<T>::optimal, svd<T>::svd): Use std::max<F77_INT> for Fortran integers.

author	Markus Mützel <markus.muetzel@gmx.de>
date	Tue, 09 Mar 2021 20:15:27 +0100
parents	9b6bf68ea663
children	9245ae55b6bd

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 //
29358 0a5b15007766 update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 27923 diff changeset	3 // Copyright (C) 1994-2021 The Octave Project Developers
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	4 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	7 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	9 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	14 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	19 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 //
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ////////////////////////////////////////////////////////////////////////
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	25
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21301 diff changeset	26 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21273 diff changeset	27 # include "config.h"
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	28 #endif
3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	29
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	30 #include <cassert>
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	31
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	32 #include <algorithm>
1631 b0b22b6ce22f [project @ 1995-11-14 20:24:15 by jwe] jwe parents: 1545 diff changeset	33
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	34 #include "CMatrix.h"
cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	35 #include "dDiagMatrix.h"
cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	36 #include "dMatrix.h"
cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	37 #include "fCMatrix.h"
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	38 #include "fDiagMatrix.h"
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	39 #include "fMatrix.h"
cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	40 #include "lo-error.h"
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22317 diff changeset	41 #include "lo-lapack-proto.h"
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	42 #include "svd.h"
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	43
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	44 // class to compute optimal work space size (lwork) for DGEJSV and SGEJSV
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	45 template<typename T>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	46 class
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	47 gejsv_lwork
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	48 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	49 public:
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	50 gejsv_lwork () = delete;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	51
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	52 // Unfortunately, dgejsv and sgejsv do not provide estimation of 'lwork'.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	53 // Thus, we have to estimate it according to corresponding LAPACK
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	54 // documentation and related source codes (e.g. cgejsv).
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	55 // In LAPACKE (C interface to LAPACK), the memory handling code in
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	56 // LAPACKE_dgejsv() (lapacke_dgejsv.c, last visit 2019-02-17) uses
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	57 // the minimum required working space. In contrast, here the optimal
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	58 // working space size is computed, at the cost of much longer code.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	59
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	60 static F77_INT optimal (char& joba, char& jobu, char& jobv,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	61 F77_INT m, F77_INT n);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	62
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	63 private:
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	64 typedef typename T::element_type P;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	65
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	66 // functions could be called from GEJSV
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	67 static F77_INT geqp3_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	68 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	69 F77_INT* jpvt, P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	70 F77_INT lwork, F77_INT& info);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	71
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	72 static F77_INT geqrf_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	73 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	74 P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	75 F77_INT lwork, F77_INT& info);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	76
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	77 static F77_INT gelqf_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	78 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	79 P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	80 F77_INT lwork, F77_INT& info);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	81
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	82 static F77_INT ormlq_lwork (char& side, char& trans,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	83 F77_INT m, F77_INT n, F77_INT k,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	84 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	85 P tau, P c, F77_INT ldc,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	86 P *work, F77_INT lwork, F77_INT& info);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	87
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	88 static F77_INT ormqr_lwork (char& side, char& trans,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	89 F77_INT m, F77_INT n, F77_INT k,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	90 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	91 P tau, P c, F77_INT ldc,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	92 P *work, F77_INT lwork, F77_INT& info);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	93 };
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	94
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	95 #define GEJSV_REAL_QP3_LWORK(f, F) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	96 F77_XFCN (f, F, (m, n, a, lda, jpvt, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	97 tau, work, lwork, info))
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	98
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	99 #define GEJSV_REAL_QR_LWORK(f, F) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	100 F77_XFCN (f, F, (m, n, a, lda, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	101 tau, work, lwork, info))
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	102
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	103 #define GEJSV_REAL_ORM_LWORK(f, F) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	104 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&side, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	105 F77_CONST_CHAR_ARG2 (&trans, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	106 m, n, k, a, lda, tau, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	107 c, ldc, work, lwork, info \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	108 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	109 F77_CHAR_ARG_LEN (1)))
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	110
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	111 // For Matrix
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	112 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	113 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	114 gejsv_lwork<Matrix>::geqp3_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	115 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	116 F77_INT* jpvt, P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	117 F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	118 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	119 GEJSV_REAL_QP3_LWORK (dgeqp3, DGEQP3);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	120 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	121 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	122
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	123 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	124 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	125 gejsv_lwork<Matrix>::geqrf_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	126 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	127 P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	128 F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	129 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	130 GEJSV_REAL_QR_LWORK (dgeqrf, DGEQRF);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	131 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	132 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	133
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	134 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	135 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	136 gejsv_lwork<Matrix>::gelqf_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	137 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	138 P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	139 F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	140 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	141 GEJSV_REAL_QR_LWORK (dgelqf, DGELQF);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	142 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	143 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	144
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	145 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	146 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	147 gejsv_lwork<Matrix>::ormlq_lwork (char& side, char& trans,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	148 F77_INT m, F77_INT n, F77_INT k,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	149 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	150 P tau, P c, F77_INT ldc,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	151 P *work, F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	152 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	153 GEJSV_REAL_ORM_LWORK (dormlq, DORMLQ);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	154 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	155 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	156
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	157 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	158 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	159 gejsv_lwork<Matrix>::ormqr_lwork (char& side, char& trans,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	160 F77_INT m, F77_INT n, F77_INT k,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	161 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	162 P tau, P c, F77_INT ldc,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	163 P *work, F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	164 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	165 GEJSV_REAL_ORM_LWORK (dormqr, DORMQR);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	166 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	167 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	168
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	169 // For FloatMatrix
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	170 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	171 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	172 gejsv_lwork<FloatMatrix>::geqp3_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	173 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	174 F77_INT* jpvt, P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	175 F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	176 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	177 GEJSV_REAL_QP3_LWORK (sgeqp3, SGEQP3);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	178 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	179 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	180
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	181 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	182 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	183 gejsv_lwork<FloatMatrix>::geqrf_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	184 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	185 P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	186 F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	187 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	188 GEJSV_REAL_QR_LWORK (sgeqrf, SGEQRF);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	189 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	190 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	191
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	192 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	193 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	194 gejsv_lwork<FloatMatrix>::gelqf_lwork (F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	195 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	196 P tau, P work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	197 F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	198 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	199 GEJSV_REAL_QR_LWORK (sgelqf, SGELQF);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	200 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	201 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	202
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	203 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	204 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	205 gejsv_lwork<FloatMatrix>::ormlq_lwork (char& side, char& trans,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	206 F77_INT m, F77_INT n, F77_INT k,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	207 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	208 P tau, P c, F77_INT ldc,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	209 P *work, F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	210 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	211 GEJSV_REAL_ORM_LWORK (sormlq, SORMLQ);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	212 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	213 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	214
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	215 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	216 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	217 gejsv_lwork<FloatMatrix>::ormqr_lwork (char& side, char& trans,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	218 F77_INT m, F77_INT n, F77_INT k,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	219 P *a, F77_INT lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	220 P tau, P c, F77_INT ldc,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	221 P *work, F77_INT lwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	222 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	223 GEJSV_REAL_ORM_LWORK (sormqr, SORMQR);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	224 return static_cast<F77_INT> (work[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	225 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	226
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	227 #undef GEJSV_REAL_QP3_LWORK
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	228 #undef GEJSV_REAL_QR_LWORK
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	229 #undef GEJSV_REAL_ORM_LWORK
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	230
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	231 template<typename T>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	232 F77_INT
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	233 gejsv_lwork<T>::optimal (char& joba, char& jobu, char& jobv,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	234 F77_INT m, F77_INT n)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	235 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	236 F77_INT lwork = -1;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	237 std::vector<P> work (2); // dummy work space
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	238
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	239 // variables that mimic running environment of gejsv
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	240 F77_INT lda = std::max<F77_INT> (m, 1);
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	241 F77_INT ierr = 0;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	242 char side = 'L';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	243 char trans = 'N';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	244 std::vector<P> mat_a (1);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	245 P *a = mat_a.data (); // dummy input matrix
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	246 std::vector<F77_INT> vec_jpvt = {0};
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	247 P *tau = work.data ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	248 P *u = work.data ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	249 P *v = work.data ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	250
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	251 bool need_lsvec = jobu == 'U' \|\| jobu == 'F';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	252 bool need_rsvec = jobv == 'V' \|\| jobv == 'J';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	253
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	254 F77_INT lw_pocon = 3 * n; // for [s,d]pocon
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	255 F77_INT lw_geqp3 = geqp3_lwork (m, n, a, lda, vec_jpvt.data (),
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	256 tau, work.data (), -1, ierr);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	257 F77_INT lw_geqrf = geqrf_lwork (m, n, a, lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	258 tau, work.data (), -1, ierr);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	259
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	260 if (! (need_lsvec \|\| need_rsvec) )
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	261 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	262 // only SIGMA is needed
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	263 if (! (joba == 'E' \|\| joba == 'G') )
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	264 lwork = std::max<F77_INT> ({2*m + n, n + lw_geqp3, n + lw_geqrf, 7});
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	265 else
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	266 lwork = std::max<F77_INT> ({2*m + n, n + lw_geqp3, n + lw_geqrf,
3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	267 n + n*n + lw_pocon, 7});
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	268 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	269 else if (need_rsvec && ! need_lsvec)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	270 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	271 // SIGMA and the right singular vectors are needed
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	272 F77_INT lw_gelqf = gelqf_lwork (n, n, a, lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	273 tau, work.data (), -1, ierr);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	274 trans = 'T';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	275 F77_INT lw_ormlq = ormlq_lwork (side, trans, n, n, n, a, lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	276 tau, v, n, work.data (), -1, ierr);
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	277 lwork = std::max<F77_INT> ({2*m + n, n + lw_geqp3, n + lw_pocon,
3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	278 n + lw_gelqf, 2*n + lw_geqrf, n + lw_ormlq});
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	279 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	280 else if (need_lsvec && ! need_rsvec)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	281 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	282 // SIGMA and the left singular vectors are needed
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	283 F77_INT n1 = (jobu == 'U') ? n : m; // size of U is m x n1
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	284 F77_INT lw_ormqr = ormqr_lwork (side, trans, m, n1, n, a, lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	285 tau, u, m, work.data (), -1, ierr);
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	286 lwork = std::max<F77_INT> ({2*m + n, n + lw_geqp3, n + lw_pocon,
3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	287 2*n + lw_geqrf, n + lw_ormqr});
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	288 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	289 else // full SVD is needed
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	290 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	291 if (jobv == 'V')
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	292 lwork = std::max (2m + n, 6n + 2nn);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	293 else if (jobv == 'J')
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	294 lwork = std::max<F77_INT> ({2m + n, 4n + nn, 2n + n*n + 6});
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	295
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	296 F77_INT n1 = (jobu == 'U') ? n : m; // size of U is m x n1
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	297 F77_INT lw_ormqr = ormqr_lwork (side, trans, m, n1, n, a, lda,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	298 tau, u, m, work.data (), -1, ierr);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	299 lwork = std::max (lwork, n + lw_ormqr);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	300 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	301
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	302 return lwork;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	303 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	304
22317 6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22297 diff changeset	305 namespace octave
6ca3acf5fad8 move some new numeric classes to namespace octave::math John W. Eaton <jwe@octave.org> parents: 22297 diff changeset	306 {
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	307 namespace math
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	308 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	309 template <typename T>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	310 T
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	311 svd<T>::left_singular_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	312 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	313 if (m_type == svd::Type::sigma_only)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	314 (*current_liboctave_error_handler)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	315 ("svd: U not computed because type == svd::sigma_only");
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 19697 diff changeset	316
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	317 return left_sm;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	318 }
1543 d6e96e0bc681 [project @ 1995-10-06 05:49:21 by jwe] jwe parents: 1368 diff changeset	319
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	320 template <typename T>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	321 T
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	322 svd<T>::right_singular_matrix (void) const
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	323 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	324 if (m_type == svd::Type::sigma_only)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	325 (*current_liboctave_error_handler)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	326 ("svd: V not computed because type == svd::sigma_only");
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 19697 diff changeset	327
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	328 return right_sm;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	329 }
1543 d6e96e0bc681 [project @ 1995-10-06 05:49:21 by jwe] jwe parents: 1368 diff changeset	330
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	331 // GESVD specializations
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	332
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	333 #define GESVD_REAL_STEP(f, F) \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	334 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&jobu, 1), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	335 F77_CONST_CHAR_ARG2 (&jobv, 1), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	336 m, n, tmp_data, m1, s_vec, u, m1, vt, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	337 nrow_vt1, work.data (), lwork, info \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	338 F77_CHAR_ARG_LEN (1) \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	339 F77_CHAR_ARG_LEN (1)))
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	340
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	341 #define GESVD_COMPLEX_STEP(f, F, CMPLX_ARG) \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	342 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&jobu, 1), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	343 F77_CONST_CHAR_ARG2 (&jobv, 1), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	344 m, n, CMPLX_ARG (tmp_data), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	345 m1, s_vec, CMPLX_ARG (u), m1, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	346 CMPLX_ARG (vt), nrow_vt1, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	347 CMPLX_ARG (work.data ()), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	348 lwork, rwork.data (), info \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	349 F77_CHAR_ARG_LEN (1) \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	350 F77_CHAR_ARG_LEN (1)))
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	351
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	352 // DGESVD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	353 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	354 OCTAVE_API void
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	355 svd<Matrix>::gesvd (char& jobu, char& jobv, F77_INT m, F77_INT n,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	356 double tmp_data, F77_INT m1, double s_vec,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	357 double u, double vt, F77_INT nrow_vt1,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	358 std::vector<double>& work, F77_INT& lwork,
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	359 F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	360 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	361 GESVD_REAL_STEP (dgesvd, DGESVD);
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	362
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	363 lwork = static_cast<F77_INT> (work[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	364 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	365
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	366 GESVD_REAL_STEP (dgesvd, DGESVD);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	367 }
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	368
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	369 // SGESVD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	370 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	371 OCTAVE_API void
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	372 svd<FloatMatrix>::gesvd (char& jobu, char& jobv, F77_INT m, F77_INT n,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	373 float tmp_data, F77_INT m1, float s_vec,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	374 float u, float vt, F77_INT nrow_vt1,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	375 std::vector<float>& work, F77_INT& lwork,
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	376 F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	377 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	378 GESVD_REAL_STEP (sgesvd, SGESVD);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	379
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	380 lwork = static_cast<F77_INT> (work[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	381 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	382
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	383 GESVD_REAL_STEP (sgesvd, SGESVD);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	384 }
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	385
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	386 // ZGESVD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	387 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	388 OCTAVE_API void
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	389 svd<ComplexMatrix>::gesvd (char& jobu, char& jobv, F77_INT m, F77_INT n,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	390 Complex tmp_data, F77_INT m1, double s_vec,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	391 Complex u, Complex vt, F77_INT nrow_vt1,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	392 std::vector<Complex>& work, F77_INT& lwork,
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	393 F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	394 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	395 std::vector<double> rwork (5 * std::max (m, n));
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	396
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	397 GESVD_COMPLEX_STEP (zgesvd, ZGESVD, F77_DBLE_CMPLX_ARG);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	398
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	399 lwork = static_cast<F77_INT> (work[0].real ());
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	400 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	401
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	402 GESVD_COMPLEX_STEP (zgesvd, ZGESVD, F77_DBLE_CMPLX_ARG);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	403 }
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	404
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	405 // CGESVD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	406 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	407 OCTAVE_API void
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	408 svd<FloatComplexMatrix>::gesvd (char& jobu, char& jobv, F77_INT m,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	409 F77_INT n, FloatComplex *tmp_data,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	410 F77_INT m1, float s_vec, FloatComplex u,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	411 FloatComplex *vt, F77_INT nrow_vt1,
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	412 std::vector<FloatComplex>& work,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	413 F77_INT& lwork, F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	414 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	415 std::vector<float> rwork (5 * std::max (m, n));
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	416
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	417 GESVD_COMPLEX_STEP (cgesvd, CGESVD, F77_CMPLX_ARG);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	418
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	419 lwork = static_cast<F77_INT> (work[0].real ());
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	420 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	421
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	422 GESVD_COMPLEX_STEP (cgesvd, CGESVD, F77_CMPLX_ARG);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	423 }
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	424
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	425 #undef GESVD_REAL_STEP
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	426 #undef GESVD_COMPLEX_STEP
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	427
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	428 // GESDD specializations
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	429
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	430 #define GESDD_REAL_STEP(f, F) \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	431 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&jobz, 1), \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	432 m, n, tmp_data, m1, s_vec, u, m1, vt, nrow_vt1, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	433 work.data (), lwork, iwork, info \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	434 F77_CHAR_ARG_LEN (1)))
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	435
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	436 #define GESDD_COMPLEX_STEP(f, F, CMPLX_ARG) \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	437 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&jobz, 1), m, n, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	438 CMPLX_ARG (tmp_data), m1, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	439 s_vec, CMPLX_ARG (u), m1, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	440 CMPLX_ARG (vt), nrow_vt1, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	441 CMPLX_ARG (work.data ()), lwork, \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	442 rwork.data (), iwork, info \
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	443 F77_CHAR_ARG_LEN (1)))
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	444
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	445 // DGESDD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	446 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	447 OCTAVE_API void
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	448 svd<Matrix>::gesdd (char& jobz, F77_INT m, F77_INT n, double *tmp_data,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	449 F77_INT m1, double s_vec, double u, double *vt,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	450 F77_INT nrow_vt1, std::vector<double>& work,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	451 F77_INT& lwork, F77_INT *iwork, F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	452 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	453 GESDD_REAL_STEP (dgesdd, DGESDD);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	454
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	455 lwork = static_cast<F77_INT> (work[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	456 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	457
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	458 GESDD_REAL_STEP (dgesdd, DGESDD);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	459 }
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	460
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	461 // SGESDD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	462 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	463 OCTAVE_API void
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	464 svd<FloatMatrix>::gesdd (char& jobz, F77_INT m, F77_INT n, float *tmp_data,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	465 F77_INT m1, float s_vec, float u, float *vt,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	466 F77_INT nrow_vt1, std::vector<float>& work,
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	467 F77_INT& lwork, F77_INT* iwork, F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	468 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	469 GESDD_REAL_STEP (sgesdd, SGESDD);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	470
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	471 lwork = static_cast<F77_INT> (work[0]);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	472 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	473
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	474 GESDD_REAL_STEP (sgesdd, SGESDD);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	475 }
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	476
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	477 // ZGESDD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	478 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	479 OCTAVE_API void
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	480 svd<ComplexMatrix>::gesdd (char& jobz, F77_INT m, F77_INT n,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	481 Complex tmp_data, F77_INT m1, double s_vec,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	482 Complex u, Complex vt, F77_INT nrow_vt1,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	483 std::vector<Complex>& work, F77_INT& lwork,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	484 F77_INT *iwork, F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	485 {
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	486
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	487 F77_INT min_mn = std::min (m, n);
25349 a8a4b6e6e754 Fix lrwork dimension in svd for complex cases (bug #53823). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	488 F77_INT max_mn = std::max (m, n);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	489
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	490 F77_INT lrwork;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	491 if (jobz == 'N')
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	492 lrwork = 7*min_mn;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	493 else
25349 a8a4b6e6e754 Fix lrwork dimension in svd for complex cases (bug #53823). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	494 lrwork = min_mn * std::max (5min_mn+5, 2max_mn+2*min_mn+1);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	495
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	496 std::vector<double> rwork (lrwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	497
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	498 GESDD_COMPLEX_STEP (zgesdd, ZGESDD, F77_DBLE_CMPLX_ARG);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	499
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	500 lwork = static_cast<F77_INT> (work[0].real ());
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	501 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	502
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	503 GESDD_COMPLEX_STEP (zgesdd, ZGESDD, F77_DBLE_CMPLX_ARG);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	504 }
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	505
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	506 // CGESDD
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	507 template<>
29228 5c14f81e0937 Set API tags in files in liboctave/numeric (patch #8919). Markus Mützel <markus.muetzel@gmx.de> parents: 27923 diff changeset	508 OCTAVE_API void
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	509 svd<FloatComplexMatrix>::gesdd (char& jobz, F77_INT m, F77_INT n,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	510 FloatComplex *tmp_data, F77_INT m1,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	511 float s_vec, FloatComplex u,
c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	512 FloatComplex *vt, F77_INT nrow_vt1,
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	513 std::vector<FloatComplex>& work,
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	514 F77_INT& lwork, F77_INT *iwork,
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	515 F77_INT& info)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	516 {
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	517 F77_INT min_mn = std::min (m, n);
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	518 F77_INT max_mn = std::max (m, n);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	519
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	520 F77_INT lrwork;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	521 if (jobz == 'N')
25349 a8a4b6e6e754 Fix lrwork dimension in svd for complex cases (bug #53823). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	522 lrwork = 7*min_mn;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	523 else
25349 a8a4b6e6e754 Fix lrwork dimension in svd for complex cases (bug #53823). Marco Caliari <marco.caliari@univr.it> parents: 25054 diff changeset	524 lrwork = min_mn * std::max (5min_mn+5, 2max_mn+2*min_mn+1);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	525 std::vector<float> rwork (lrwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	526
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	527 GESDD_COMPLEX_STEP (cgesdd, CGESDD, F77_CMPLX_ARG);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	528
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	529 lwork = static_cast<F77_INT> (work[0].real ());
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	530 work.reserve (lwork);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	531
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	532 GESDD_COMPLEX_STEP (cgesdd, CGESDD, F77_CMPLX_ARG);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	533 }
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	534
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	535 #undef GESDD_REAL_STEP
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	536 #undef GESDD_COMPLEX_STEP
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	537
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	538 // GEJSV specializations
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	539
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	540 #define GEJSV_REAL_STEP(f, F) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	541 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&joba, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	542 F77_CONST_CHAR_ARG2 (&jobu, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	543 F77_CONST_CHAR_ARG2 (&jobv, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	544 F77_CONST_CHAR_ARG2 (&jobr, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	545 F77_CONST_CHAR_ARG2 (&jobt, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	546 F77_CONST_CHAR_ARG2 (&jobp, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	547 m, n, tmp_data, m1, s_vec, u, m1, v, nrow_v1, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	548 work.data (), lwork, iwork.data (), info \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	549 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	550 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	551 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	552 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	553 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	554 F77_CHAR_ARG_LEN (1)))
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	555
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	556 #define GEJSV_COMPLEX_STEP(f, F, CMPLX_ARG) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	557 F77_XFCN (f, F, (F77_CONST_CHAR_ARG2 (&joba, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	558 F77_CONST_CHAR_ARG2 (&jobu, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	559 F77_CONST_CHAR_ARG2 (&jobv, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	560 F77_CONST_CHAR_ARG2 (&jobr, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	561 F77_CONST_CHAR_ARG2 (&jobt, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	562 F77_CONST_CHAR_ARG2 (&jobp, 1), \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	563 m, n, CMPLX_ARG (tmp_data), m1, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	564 s_vec, CMPLX_ARG (u), m1, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	565 CMPLX_ARG (v), nrow_v1, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	566 CMPLX_ARG (work.data ()), lwork, \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	567 rwork.data (), lrwork, iwork.data (), info \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	568 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	569 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	570 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	571 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	572 F77_CHAR_ARG_LEN (1) \
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	573 F77_CHAR_ARG_LEN (1)))
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	574
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	575 // DGEJSV
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	576 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	577 void
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	578 svd<Matrix>::gejsv (char& joba, char& jobu, char& jobv,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	579 char& jobr, char& jobt, char& jobp,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	580 F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	581 P tmp_data, F77_INT m1, DM_P s_vec, P *u,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	582 P *v, F77_INT nrow_v1, std::vector<P>& work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	583 F77_INT& lwork, std::vector<F77_INT>& iwork,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	584 F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	585 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	586 lwork = gejsv_lwork<Matrix>::optimal(joba, jobu, jobv, m, n);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	587 work.reserve (lwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	588
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	589 GEJSV_REAL_STEP (dgejsv, DGEJSV);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	590 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	591
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	592 // SGEJSV
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	593 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	594 void
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	595 svd<FloatMatrix>::gejsv (char& joba, char& jobu, char& jobv,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	596 char& jobr, char& jobt, char& jobp,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	597 F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	598 P tmp_data, F77_INT m1, DM_P s_vec, P *u,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	599 P *v, F77_INT nrow_v1, std::vector<P>& work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	600 F77_INT& lwork, std::vector<F77_INT>& iwork,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	601 F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	602 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	603 lwork = gejsv_lwork<FloatMatrix>::optimal(joba, jobu, jobv, m, n);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	604 work.reserve (lwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	605
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	606 GEJSV_REAL_STEP (sgejsv, SGEJSV);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	607 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	608
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	609 // ZGEJSV
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	610 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	611 void
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	612 svd<ComplexMatrix>::gejsv (char& joba, char& jobu, char& jobv,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	613 char& jobr, char& jobt, char& jobp,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	614 F77_INT m, F77_INT n,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	615 P tmp_data, F77_INT m1, DM_P s_vec, P *u,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	616 P *v, F77_INT nrow_v1, std::vector<P>& work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	617 F77_INT& lwork, std::vector<F77_INT>& iwork,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	618 F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	619 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	620 F77_INT lrwork = -1; // work space size query
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	621 std::vector<double> rwork (1);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	622 work.reserve (2);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	623
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	624 GEJSV_COMPLEX_STEP (zgejsv, ZGEJSV, F77_DBLE_CMPLX_ARG);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	625
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	626 lwork = static_cast<F77_INT> (work[0].real ());
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	627 work.reserve (lwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	628
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	629 lrwork = static_cast<F77_INT> (rwork[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	630 rwork.reserve (lrwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	631
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	632 F77_INT liwork = static_cast<F77_INT> (iwork[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	633 iwork.reserve (liwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	634
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	635 GEJSV_COMPLEX_STEP (zgejsv, ZGEJSV, F77_DBLE_CMPLX_ARG);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	636 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	637
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	638 // CGEJSV
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	639 template<>
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	640 void
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	641 svd<FloatComplexMatrix>::gejsv (char& joba, char& jobu, char& jobv,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	642 char& jobr, char& jobt, char& jobp,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	643 F77_INT m, F77_INT n, P *tmp_data,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	644 F77_INT m1, DM_P s_vec, P u, P *v,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	645 F77_INT nrow_v1, std::vector<P>& work,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	646 F77_INT& lwork,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	647 std::vector<F77_INT>& iwork, F77_INT& info)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	648 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	649 F77_INT lrwork = -1; // work space size query
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	650 std::vector<float> rwork (1);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	651 work.reserve (2);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	652
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	653 GEJSV_COMPLEX_STEP (cgejsv, CGEJSV, F77_CMPLX_ARG);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	654
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	655 lwork = static_cast<F77_INT> (work[0].real ());
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	656 work.reserve (lwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	657
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	658 lrwork = static_cast<F77_INT> (rwork[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	659 rwork.reserve (lrwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	660
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	661 F77_INT liwork = static_cast<F77_INT> (iwork[0]);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	662 iwork.reserve (liwork);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	663
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	664 GEJSV_COMPLEX_STEP (cgejsv, CGEJSV, F77_CMPLX_ARG);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	665 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	666
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	667 #undef GEJSV_REAL_STEP
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	668 #undef GEJSV_COMPLEX_STEP
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	669
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	670 template<typename T>
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	671 svd<T>::svd (const T& a, svd::Type type,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	672 svd::Driver driver)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	673 : m_type (type), m_driver (driver), left_sm (), sigma (), right_sm ()
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	674 {
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	675 F77_INT info;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	676
23696 08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	677 F77_INT m = to_f77_int (a.rows ());
08036a7f3660 remove octave:: namespace tag from symbols used inside octave namespace John W. Eaton <jwe@octave.org> parents: 23475 diff changeset	678 F77_INT n = to_f77_int (a.cols ());
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	679
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	680 if (m == 0 \|\| n == 0)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	681 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	682 switch (m_type)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	683 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	684 case svd::Type::std:
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	685 left_sm = T (m, m, 0);
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	686 for (F77_INT i = 0; i < m; i++)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	687 left_sm.xelem (i, i) = 1;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	688 sigma = DM_T (m, n);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	689 right_sm = T (n, n, 0);
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	690 for (F77_INT i = 0; i < n; i++)
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	691 right_sm.xelem (i, i) = 1;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	692 break;
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	693
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	694 case svd::Type::economy:
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	695 left_sm = T (m, 0, 0);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	696 sigma = DM_T (0, 0);
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	697 right_sm = T (n, 0, 0);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	698 break;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	699
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	700 case svd::Type::sigma_only:
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	701 default:
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	702 sigma = DM_T (0, 1);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	703 break;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	704 }
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	705 return;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	706 }
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	707
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	708 T atmp = a;
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	709 P *tmp_data = atmp.fortran_vec ();
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	710
23450 855122b993da maint: Wrap tertiary operator in parentheses "(COND ? x : y)". Rik <rik@octave.org> parents: 23449 diff changeset	711 F77_INT min_mn = (m < n ? m : n);
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	712
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	713 char jobu = 'A';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	714 char jobv = 'A';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	715
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	716 F77_INT ncol_u = m;
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	717 F77_INT nrow_vt = n;
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	718 F77_INT nrow_s = m;
5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	719 F77_INT ncol_s = n;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	720
22228 4afe3705ea75 svd: do not shadow member variables with local variables. Carnë Draug <carandraug@octave.org> parents: 22204 diff changeset	721 switch (m_type)
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	722 {
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	723 case svd::Type::economy:
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	724 jobu = jobv = 'S';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	725 ncol_u = nrow_vt = nrow_s = ncol_s = min_mn;
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	726 break;
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	727
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	728 case svd::Type::sigma_only:
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	729
23475 d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	730 // Note: for this case, both jobu and jobv should be 'N', but there
d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	731 // seems to be a bug in dgesvd from Lapack V2.0. To demonstrate the
d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	732 // bug, set both jobu and jobv to 'N' and find the singular values of
d691ed308237 maint: Clean up #includes in liboctave/numeric directory. Rik <rik@octave.org> parents: 23450 diff changeset	733 // [eye(3), eye(3)]. The result is [-sqrt(2), -sqrt(2), -sqrt(2)].
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	734 //
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	735 // For Lapack 3.0, this problem seems to be fixed.
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	736
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	737 jobu = jobv = 'N';
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	738 ncol_u = nrow_vt = 1;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	739 break;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	740
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	741 default:
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	742 break;
469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	743 }
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	744
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	745 if (! (jobu == 'N' \|\| jobu == 'O'))
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	746 left_sm.resize (m, ncol_u);
22204 469c817eb256 svd: reduce code duplication with more use of template and macro. Carnë Draug <carandraug@octave.org> parents: 22135 diff changeset	747
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	748 P *u = left_sm.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	749
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	750 sigma.resize (nrow_s, ncol_s);
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	751 DM_P *s_vec = sigma.fortran_vec ();
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	752
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	753 if (! (jobv == 'N' \|\| jobv == 'O'))
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	754 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	755 if (m_driver == svd::Driver::GEJSV)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	756 right_sm.resize (n, nrow_vt);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	757 else
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	758 right_sm.resize (nrow_vt, n);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	759 }
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	760
23449 c763214a8260 maint: Use convention 'int x' for naming pointers. Rik <rik@octave.org>* parents: 23448 diff changeset	761 P *vt = right_sm.fortran_vec ();
537 4ecbfd3c3710 [project @ 1994-07-21 22:30:34 by jwe] jwe parents: 457 diff changeset	762
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	763 // Query _GESVD for the correct dimension of WORK.
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	764
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	765 F77_INT lwork = -1;
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	766
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	767 std::vector<P> work (1);
537 4ecbfd3c3710 [project @ 1994-07-21 22:30:34 by jwe] jwe parents: 457 diff changeset	768
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	769 const F77_INT f77_int_one = static_cast<F77_INT> (1);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	770 F77_INT m1 = std::max (m, f77_int_one);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	771 F77_INT nrow_vt1 = std::max (nrow_vt, f77_int_one);
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	772
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	773 if (m_driver == svd::Driver::GESVD)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	774 gesvd (jobu, jobv, m, n, tmp_data, m1, s_vec, u, vt, nrow_vt1,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	775 work, lwork, info);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	776 else if (m_driver == svd::Driver::GESDD)
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	777 {
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	778 assert (jobu == jobv);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	779 char jobz = jobu;
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	780
22944 5126040e8f49 use F77_INT instead of octave_idx_type for liboctave svd class John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	781 std::vector<F77_INT> iwork (8 * std::min (m, n));
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	782
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	783 gesdd (jobz, m, n, tmp_data, m1, s_vec, u, vt, nrow_vt1,
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	784 work, lwork, iwork.data (), info);
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	785 }
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	786 else if (m_driver == svd::Driver::GEJSV)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	787 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	788 bool transposed = false;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	789 if (n > m)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	790 {
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	791 // GEJSV only accepts m >= n, thus we need to transpose here
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	792 transposed = true;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	793
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	794 std::swap (m, n);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	795 m1 = std::max (m, f77_int_one);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	796 nrow_vt1 = std::max (n, f77_int_one); // we have m > n
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	797 if (m_type == svd::Type::sigma_only)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	798 nrow_vt1 = 1;
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	799 std::swap (jobu, jobv);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	800
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	801 atmp = atmp.hermitian ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	802 tmp_data = atmp.fortran_vec ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	803
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	804 // Swap pointers of U and V.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	805 u = right_sm.fortran_vec ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	806 vt = left_sm.fortran_vec ();
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	807 }
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	808
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	809 // translate jobu and jobv from gesvd to gejsv.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	810 assert ('A' <= jobu && jobu <= 'S' && 'A' <= jobv && jobv <= 'S');
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	811 char job_svd2jsv[1 + 'S' - 'A'][2] = {0};
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	812 job_svd2jsv['A' - 'A'][0] = 'F';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	813 job_svd2jsv['A' - 'A'][1] = 'J';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	814 job_svd2jsv['S' - 'A'][0] = 'U';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	815 job_svd2jsv['S' - 'A'][1] = 'V';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	816 job_svd2jsv['O' - 'A'][0] = 'W';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	817 job_svd2jsv['O' - 'A'][1] = 'W';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	818 job_svd2jsv['N' - 'A'][0] = 'N';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	819 job_svd2jsv['N' - 'A'][1] = 'N';
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	820 jobu = job_svd2jsv[jobu - 'A'][0];
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	821 jobv = job_svd2jsv[jobv - 'A'][1];
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	822
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	823 char joba = 'F'; // 'F': most conservative
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	824 char jobr = 'R'; // 'R' is recommended.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	825 char jobt = 'N'; // or 'T', but that requires U and V appear together
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	826 char jobp = 'N'; // use 'P' if denormal is poorly implemented.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	827
29416 3f45f5b6fa7f svd.cc: Use std::max<F77_INT> for Fortran integers (bug #55727). Markus Mützel <markus.muetzel@gmx.de> parents: 29412 diff changeset	828 std::vector<F77_INT> iwork (std::max<F77_INT> (m + 3*n, 1));
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	829
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	830 gejsv (joba, jobu, jobv, jobr, jobt, jobp, m, n, tmp_data, m1,
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	831 s_vec, u, vt, nrow_vt1, work, lwork, iwork, info);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	832
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	833 if (iwork[2] == 1)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	834 (*current_liboctave_warning_with_id_handler)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	835 ("Octave:convergence", "svd: (driver: GEJSV) "
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	836 "Denormal occured, possible loss of accuracy.");
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	837
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	838 if (info < 0)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	839 (*current_liboctave_error_handler)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	840 ("svd: (driver: GEJSV) Illegal argument at #%d",
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	841 static_cast<int> (-info));
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	842 else if (info > 0)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	843 (*current_liboctave_warning_with_id_handler)
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	844 ("Octave:convergence", "svd: (driver: GEJSV) "
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	845 "Fail to converge within max sweeps, "
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	846 "possible inaccurate result.");
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	847
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	848 if (transposed) // put things that need to transpose back here
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	849 std::swap (m, n);
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	850 }
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	851 else
23448 e1c02e7126a9 avoid calling abort if possible; include cstdlib for declaration if needed John W. Eaton <jwe@octave.org> parents: 23220 diff changeset	852 (*current_liboctave_error_handler) ("svd: unknown driver");
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	853
26731 75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	854 // LAPACK can return -0 which is a small problem (bug #55710).
75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	855 for (octave_idx_type i = 0; i < sigma.diag_length (); i++)
75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	856 {
75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	857 if (! sigma.dgxelem (i))
75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	858 sigma.dgxelem (i) = DM_P (0);
75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	859 }
75d79c39ac92 Replace '-0' in sigma matrix from svd() with '+0' (bug #55710). Rik <rik@octave.org> parents: 26376 diff changeset	860
29412 9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	861 // GESVD and GESDD return VT instead of V, GEJSV return V.
9b6bf68ea663 Add additional svd driver GEJSV for accurate SVD (bug #55727). Eddy Xiao <bewantbe@gmail.com> parents: 29359 diff changeset	862 if (! (jobv == 'N' \|\| jobv == 'O') && (m_driver != svd::Driver::GEJSV))
22592 ee0df00e12d6 fix complex svd (bug #49232) John W. Eaton <jwe@octave.org> parents: 22402 diff changeset	863 right_sm = right_sm.hermitian ();
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	864 }
cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	865
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	866 // Instantiations we need.
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	867
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	868 template class svd<Matrix>;
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	869
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	870 template class svd<FloatMatrix>;
21273 cbced1c09916 better use of templates for svd classes John W. Eaton <jwe@octave.org> parents: 21202 diff changeset	871
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	872 template class svd<ComplexMatrix>;
3336 08ad797989f8 [project @ 1999-11-03 21:41:34 by jwe] jwe parents: 3335 diff changeset	873
22329 7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	874 template class svd<FloatComplexMatrix>;
7f3c7a8bd131 maint: Indent namespaces in liboctave/numeric files. John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	875 }
457 3d4b4f0fa5ba [project @ 1994-06-06 00:33:33 by jwe] jwe parents: diff changeset	876 }

Mercurial > octave

annotate liboctave/numeric/svd.cc @ 29416:3f45f5b6fa7f