Mercurial > octave
annotate libinterp/corefcn/gsvd.cc @ 29958:32c3a5805893
move DEFUN and DEFMETHOD functions inside octave namespace
* mk-opts.pl: Surround emitted function definitions with
OCTAVE_NAMESPACE_BEGIN and OCTAVE_NAMESPACE_END tags.
* mk-builtins.pl: Surround emitted function declarations with
OCTAVE_NAMESPACE_BEGIN and OCTAVE_NAMESPACE_END tags. Emit deprecated
global inline functions.
* __betainc__.cc, __contourc__.cc, __dsearchn__.cc, __eigs__.cc,
__expint__.cc, __ftp__.cc, __gammainc__.cc, __ichol__.cc, __ilu__.cc,
__lin_interpn__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc,
amd.cc, balance.cc, besselj.cc, bitfcns.cc, bsxfun.cc, call-stack.cc,
ccolamd.cc, cellfun.cc, chol.cc, colamd.cc, colloc.cc, conv2.cc,
daspk.cc, dasrt.cc, dassl.cc, data.cc, debug.cc, defaults.cc,
defun.cc, det.cc, dirfns.cc, display.cc, dlmread.cc, dmperm.cc,
dot.cc, eig.cc, ellipj.cc, environment.cc, error.cc, event-manager.cc,
fcn-info.cc, fft.cc, fft2.cc, fftn.cc, file-io.cc, filter.cc, find.cc,
gcd.cc, getgrent.cc, getpwent.cc, getrusage.cc, givens.cc,
graphics.cc, gsvd.cc, hash.cc, help.cc, hess.cc, hex2num.cc, input.cc,
interpreter.cc, inv.cc, jsondecode.cc, jsonencode.cc, kron.cc,
load-path.cc, load-save.cc, lookup.cc, ls-oct-text.cc, lsode.cc,
lu.cc, mappers.cc, matrix_type.cc, max.cc, mgorth.cc, nproc.cc,
oct-hist.cc, ordqz.cc, ordschur.cc, pager.cc, pinv.cc, pr-flt-fmt.cc,
pr-output.cc, psi.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc,
rcond.cc, regexp.cc, schur.cc, settings.cc, sighandlers.cc, sparse.cc,
spparms.cc, sqrtm.cc, stream-euler.cc, strfind.cc, strfns.cc,
sub2ind.cc, svd.cc, sylvester.cc, symbfact.cc, symrcm.cc, symtab.cc,
syscalls.cc, sysdep.cc, time.cc, toplev.cc, tril.cc, tsearch.cc,
typecast.cc, urlwrite.cc, utils.cc, variables.cc, __delaunayn__.cc,
__fltk_uigetfile__.cc, __glpk__.cc, __init_gnuplot__.cc, __ode15__.cc,
__voronoi__.cc, audiodevinfo.cc, audioread.cc, convhulln.cc, fftw.cc,
gzip.cc, ov-base.cc, ov-bool-mat.cc, ov-cell.cc, ov-class.cc,
ov-classdef.cc, ov-fcn-handle.cc, ov-java.cc, ov-null-mat.cc,
ov-oncleanup.cc, ov-struct.cc, ov-typeinfo.cc, ov-usr-fcn.cc, ov.cc,
octave.cc, lex.ll, oct-parse.yy, profiler.cc, pt-eval.cc: Surround
DEFUN and DEFMETHOD function defnitions with OCTAVE_NAMESPACE_BEGIN
and OCTAVE_NAMESPACE_END tags.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 13 Aug 2021 21:53:51 -0400 |
parents | 7854d5752dd2 |
children | 7d6709900da7 |
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) 1997-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 //////////////////////////////////////////////////////////////////////// |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
25 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
26 #ifdef HAVE_CONFIG_H |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
27 # include <config.h> |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
28 #endif |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
29 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
30 #include "dMatrix.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
31 #include "CMatrix.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
32 #include "dDiagMatrix.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
33 #include "gsvd.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
34 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
35 #include "defun.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
36 #include "defun-int.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
37 #include "error.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
38 #include "errwarn.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
39 #include "utils.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
40 #include "ovl.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
41 #include "ov.h" |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
42 |
29958
32c3a5805893
move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
29359
diff
changeset
|
43 OCTAVE_NAMESPACE_BEGIN |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
44 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
45 template <typename T> |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
46 static typename octave::math::gsvd<T>::Type |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
47 gsvd_type (int nargout) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
48 { |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
49 return ((nargout == 0 || nargout == 1) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
50 ? octave::math::gsvd<T>::Type::sigma_only |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
51 : (nargout > 5) ? octave::math::gsvd<T>::Type::std |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
52 : octave::math::gsvd<T>::Type::economy); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
53 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
54 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
55 // Named like this to avoid conflicts with the gsvd class. |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
56 template <typename T> |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
57 static octave_value_list |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
58 do_gsvd (const T& A, const T& B, const octave_idx_type nargout, |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
59 bool is_single = false) |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
60 { |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
61 octave::math::gsvd<T> result (A, B, gsvd_type<T> (nargout)); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
62 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
63 octave_value_list retval (nargout); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
64 if (nargout < 2) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
65 { |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
66 if (is_single) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
67 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
68 FloatDiagMatrix sigA = result.singular_values_A (); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
69 FloatDiagMatrix sigB = result.singular_values_B (); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
70 for (int i = sigA.rows () - 1; i >= 0; i--) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
71 sigA.dgxelem(i) /= sigB.dgxelem(i); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
72 retval(0) = sigA.diag (); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
73 } |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
74 else |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
75 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
76 DiagMatrix sigA = result.singular_values_A (); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
77 DiagMatrix sigB = result.singular_values_B (); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
78 for (int i = sigA.rows () - 1; i >= 0; i--) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
79 sigA.dgxelem(i) /= sigB.dgxelem(i); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
80 retval(0) = sigA.diag (); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
81 } |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
82 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
83 else |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
84 { |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
85 retval(0) = result.left_singular_matrix_A (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
86 retval(1) = result.left_singular_matrix_B (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
87 if (nargout > 2) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
88 retval(2) = result.right_singular_matrix (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
89 if (nargout > 3) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
90 retval(3) = result.singular_values_A (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
91 if (nargout > 4) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
92 retval(4) = result.singular_values_B (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
93 if (nargout > 5) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
94 retval(5) = result.R_matrix (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
95 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
96 return retval; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
97 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
98 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
99 DEFUN (gsvd, args, nargout, |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
100 doc: /* -*- texinfo -*- |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
101 @deftypefn {} {@var{S} =} gsvd (@var{A}, @var{B}) |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
102 @deftypefnx {} {[@var{U}, @var{V}, @var{X}, @var{C}, @var{S}] =} gsvd (@var{A}, @var{B}) |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
103 @deftypefnx {} {[@var{U}, @var{V}, @var{X}, @var{C}, @var{S}] =} gsvd (@var{A}, @var{B}, 0) |
25106
d7ad543255c5
doc: Shorten very long first sentences of docstrings (bug #53388).
Rik <rik@octave.org>
parents:
25054
diff
changeset
|
104 Compute the generalized singular value decomposition of (@var{A}, @var{B}). |
d7ad543255c5
doc: Shorten very long first sentences of docstrings (bug #53388).
Rik <rik@octave.org>
parents:
25054
diff
changeset
|
105 |
d7ad543255c5
doc: Shorten very long first sentences of docstrings (bug #53388).
Rik <rik@octave.org>
parents:
25054
diff
changeset
|
106 The generalized singular value decomposition is defined by the following |
d7ad543255c5
doc: Shorten very long first sentences of docstrings (bug #53388).
Rik <rik@octave.org>
parents:
25054
diff
changeset
|
107 relations: |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
108 |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
109 @tex |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
110 $$ A = U C X^\dagger $$ |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
111 $$ B = V S X^\dagger $$ |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
112 $$ C^\dagger C + S^\dagger S = eye (columns (A)) $$ |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
113 @end tex |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
114 @ifnottex |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
115 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
116 @example |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
117 @group |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
118 A = U*C*X' |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
119 B = V*S*X' |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
120 C'*C + S'*S = eye (columns (A)) |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
121 @end group |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
122 @end example |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
123 |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
124 @end ifnottex |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
125 |
23565
3a730821e4a2
doc: Peridoc grammarcheck of documentation.
Rik <rik@octave.org>
parents:
23221
diff
changeset
|
126 The function @code{gsvd} normally returns just the vector of generalized |
3a730821e4a2
doc: Peridoc grammarcheck of documentation.
Rik <rik@octave.org>
parents:
23221
diff
changeset
|
127 singular values |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
128 @tex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
129 $$ \sqrt{{{diag (C^\dagger C)} \over {diag (S^\dagger S)}}} $$ |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
130 @end tex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
131 @ifnottex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
132 @code{sqrt (diag (C'*C) ./ diag (S'*S))}. |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
133 @end ifnottex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
134 If asked for five return values, it also computes |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
135 @tex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
136 $U$, $V$, $X$, and $C$. |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
137 @end tex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
138 @ifnottex |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
139 U, V, X, and C. |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
140 @end ifnottex |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
141 |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
142 If the optional third input is present, @code{gsvd} constructs the |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
143 "economy-sized" decomposition where the number of columns of @var{U}, @var{V} |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
144 and the number of rows of @var{C}, @var{S} is less than or equal to the number |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
145 of columns of @var{A}. This option is not yet implemented. |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
146 |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
147 Programming Note: the code is a wrapper to the corresponding @sc{lapack} dggsvd |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
148 and zggsvd routines. |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
149 |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
150 @seealso{svd} |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
151 @end deftypefn */) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
152 { |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
153 int nargin = args.length (); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
154 |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
155 if (nargin < 2 || nargin > 3) |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
156 print_usage (); |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
157 else if (nargin == 3) |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
158 warning ("gsvd: economy-sized decomposition is not yet implemented, returning full decomposition"); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
159 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
160 octave_value_list retval; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
161 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
162 octave_value argA = args(0); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
163 octave_value argB = args(1); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
164 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
165 octave_idx_type nr = argA.rows (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
166 octave_idx_type nc = argA.columns (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
167 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
168 octave_idx_type np = argB.columns (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
169 |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
170 // FIXME: This "special" case should be handled in the gsvd class, not here |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
171 if (nr == 0 || nc == 0) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
172 { |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
173 retval = octave_value_list (nargout); |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
174 if (nargout < 2) // S = gsvd (A, B) |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
175 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
176 if (argA.is_single_type () || argB.is_single_type ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
177 retval(0) = FloatMatrix (0, 1); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
178 else |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
179 retval(0) = Matrix (0, 1); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
180 } |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
181 else // [U, V, X, C, S, R] = gsvd (A, B) |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
182 { |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
183 if (argA.is_single_type () || argB.is_single_type ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
184 { |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
185 retval(0) = octave::float_identity_matrix (nc, nc); |
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
186 retval(1) = octave::float_identity_matrix (nc, nc); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
187 if (nargout > 2) |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
188 retval(2) = octave::float_identity_matrix (nr, nr); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
189 if (nargout > 3) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
190 retval(3) = FloatMatrix (nr, nc); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
191 if (nargout > 4) |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
192 retval(4) = octave::float_identity_matrix (nr, nr); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
193 if (nargout > 5) |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
194 retval(5) = octave::float_identity_matrix (nr, nr); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
195 } |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
196 else |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
197 { |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
198 retval(0) = octave::identity_matrix (nc, nc); |
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
199 retval(1) = octave::identity_matrix (nc, nc); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
200 if (nargout > 2) |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
201 retval(2) = octave::identity_matrix (nr, nr); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
202 if (nargout > 3) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
203 retval(3) = Matrix (nr, nc); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
204 if (nargout > 4) |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
205 retval(4) = octave::identity_matrix (nr, nr); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
206 if (nargout > 5) |
25646
4d565baa475e
move libinterp/utils functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
25106
diff
changeset
|
207 retval(5) = octave::identity_matrix (nr, nr); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
208 } |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
209 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
210 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
211 else |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
212 { |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
213 if (nc != np) |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
214 print_usage (); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
215 |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
216 if (argA.is_single_type () || argB.is_single_type ()) |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
217 { |
23582
0cc2011d800e
maint: Deprecate is_real_type and replace with isreal.
Rik <rik@octave.org>
parents:
23581
diff
changeset
|
218 if (argA.isreal () && argB.isreal ()) |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
219 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
220 FloatMatrix tmpA = argA.xfloat_matrix_value ("gsvd: A must be a real or complex matrix"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
221 FloatMatrix tmpB = argB.xfloat_matrix_value ("gsvd: B must be a real or complex matrix"); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
222 |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
223 if (tmpA.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
224 error ("gsvd: A cannot have Inf or NaN values"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
225 if (tmpB.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
226 error ("gsvd: B cannot have Inf or NaN values"); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
227 |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
228 retval = do_gsvd (tmpA, tmpB, nargout, true); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
229 } |
23581
c3075ae020e1
maint: Deprecate is_complex_type and replace with iscomplex.
Rik <rik@octave.org>
parents:
23565
diff
changeset
|
230 else if (argA.iscomplex () || argB.iscomplex ()) |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
231 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
232 FloatComplexMatrix ctmpA = argA.xfloat_complex_matrix_value ("gsvd: A must be a real or complex matrix"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
233 FloatComplexMatrix ctmpB = argB.xfloat_complex_matrix_value ("gsvd: B must be a real or complex matrix"); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
234 |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
235 if (ctmpA.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
236 error ("gsvd: A cannot have Inf or NaN values"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
237 if (ctmpB.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
238 error ("gsvd: B cannot have Inf or NaN values"); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
239 |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
240 retval = do_gsvd (ctmpA, ctmpB, nargout, true); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
241 } |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
242 else |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
243 error ("gsvd: A and B must be real or complex matrices"); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
244 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
245 else |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
246 { |
23582
0cc2011d800e
maint: Deprecate is_real_type and replace with isreal.
Rik <rik@octave.org>
parents:
23581
diff
changeset
|
247 if (argA.isreal () && argB.isreal ()) |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
248 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
249 Matrix tmpA = argA.xmatrix_value ("gsvd: A must be a real or complex matrix"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
250 Matrix tmpB = argB.xmatrix_value ("gsvd: B must be a real or complex matrix"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
251 |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
252 if (tmpA.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
253 error ("gsvd: A cannot have Inf or NaN values"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
254 if (tmpB.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
255 error ("gsvd: B cannot have Inf or NaN values"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
256 |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
257 retval = do_gsvd (tmpA, tmpB, nargout); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
258 } |
23581
c3075ae020e1
maint: Deprecate is_complex_type and replace with iscomplex.
Rik <rik@octave.org>
parents:
23565
diff
changeset
|
259 else if (argA.iscomplex () || argB.iscomplex ()) |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
260 { |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
261 ComplexMatrix ctmpA = argA.xcomplex_matrix_value ("gsvd: A must be a real or complex matrix"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
262 ComplexMatrix ctmpB = argB.xcomplex_matrix_value ("gsvd: B must be a real or complex matrix"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
263 |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
264 if (ctmpA.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
265 error ("gsvd: A cannot have Inf or NaN values"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
266 if (ctmpB.any_element_is_inf_or_nan ()) |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
267 error ("gsvd: B cannot have Inf or NaN values"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
268 |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
269 retval = do_gsvd (ctmpA, ctmpB, nargout); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
270 } |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
271 else |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
272 error ("gsvd: A and B must be real or complex matrices"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
273 } |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
274 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
275 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
276 return retval; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
277 } |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
278 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
279 /* |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
280 |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
281 ## Basic test of decomposition |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
282 %!test <48807> |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
283 %! A = reshape (1:15,5,3); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
284 %! B = magic (3); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
285 %! [U,V,X,C,S] = gsvd (A,B); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
286 %! assert (U*C*X', A, 50*eps); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
287 %! assert (V*S*X', B, 50*eps); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
288 %! S0 = gsvd (A, B); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
289 %! S1 = svd (A / B); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
290 %! assert (S0, S1, 10*eps); |
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
291 |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
292 ## a few tests for gsvd.m |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
293 %!shared A, A0, B, B0, U, V, C, S, X, R, D1, D2 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
294 %! A0 = randn (5, 3); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
295 %! B0 = diag ([1 2 4]); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
296 %! A = A0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
297 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
298 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
299 ## A (5x3) and B (3x3) are full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
300 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
301 %! [U, V, X, C, S, R] = gsvd (A, B); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
302 %! D1 = zeros (5, 3); D1(1:3, 1:3) = C; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
303 %! D2 = S; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
304 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
305 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
306 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
307 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
308 ## A: 5x3 full rank, B: 3x3 rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
309 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
310 %! B(2, 2) = 0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
311 %! [U, V, X, C, S, R] = gsvd (A, B); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
312 %! D1 = zeros (5, 3); D1(1, 1) = 1; D1(2:3, 2:3) = C; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
313 %! D2 = [zeros(2, 1) S; zeros(1, 3)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
314 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (2, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
315 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
316 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
317 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
318 ## A: 5x3 rank deficient, B: 3x3 full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
319 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
320 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
321 %! A(:, 3) = 2*A(:, 1) - A(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
322 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
323 %! D1 = zeros (5, 3); D1(1:3, 1:3) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
324 %! D2 = S; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
325 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
326 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
327 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
328 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
329 ## A and B are both rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
330 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
331 %! B(:, 3) = 2*B(:, 1) - B(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
332 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
333 %! D1 = zeros (5, 2); D1(1:2, 1:2) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
334 %! D2 = [S; zeros(1, 2)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
335 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (2, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
336 %! assert (norm ((U'*A*X) - D1*[zeros(2, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
337 %! assert (norm ((V'*B*X) - D2*[zeros(2, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
338 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
339 ## A (now 3x5) and B (now 5x5) are full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
340 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
341 %! A = A0.'; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
342 %! B0 = diag ([1 2 4 8 16]); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
343 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
344 %! [U, V, X, C, S, R] = gsvd (A, B); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
345 %! D1 = [C zeros(3,2)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
346 %! D2 = [S zeros(3,2); zeros(2, 3) eye(2)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
347 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
348 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
349 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
350 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
351 ## A: 3x5 full rank, B: 5x5 rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
352 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
353 %! B(2, 2) = 0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
354 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
355 %! D1 = zeros (3, 5); D1(1, 1) = 1; D1(2:3, 2:3) = C; |
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
356 %! D2 = zeros (5, 5); D2(1:2, 2:3) = S; D2(3:4, 4:5) = eye (2); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
357 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (2, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
358 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
359 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
360 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
361 ## A: 3x5 rank deficient, B: 5x5 full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
362 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
363 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
364 %! A(3, :) = 2*A(1, :) - A(2, :); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
365 %! [U, V, X, C, S, R] = gsvd (A, B); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
366 %! D1 = zeros (3, 5); D1(1:3, 1:3) = C; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
367 %! D2 = zeros (5, 5); D2(1:3, 1:3) = S; D2(4:5, 4:5) = eye (2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
368 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
369 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
370 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
371 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
372 ## A and B are both rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
373 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
374 %! A = A0.'; B = B0.'; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
375 %! A(:, 3) = 2*A(:, 1) - A(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
376 %! B(:, 3) = 2*B(:, 1) - B(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
377 %! [U, V, X, C, S, R]=gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
378 %! D1 = zeros (3, 4); D1(1:3, 1:3) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
379 %! D2 = eye (4); D2(1:3, 1:3) = S; D2(5,:) = 0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
380 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
381 %! assert (norm ((U'*A*X) - D1*[zeros(4, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
382 %! assert (norm ((V'*B*X) - D2*[zeros(4, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
383 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
384 ## A: 5x3 complex full rank, B: 3x3 complex full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
385 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
386 %! A0 = A0 + j*randn (5, 3); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
387 %! B0 = diag ([1 2 4]) + j*diag ([4 -2 -1]); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
388 %! A = A0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
389 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
390 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
391 %! D1 = zeros (5, 3); D1(1:3, 1:3) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
392 %! D2 = S; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
393 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
394 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
395 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
396 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
397 ## A: 5x3 complex full rank, B: 3x3 complex rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
398 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
399 %! B(2, 2) = 0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
400 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
401 %! D1 = zeros (5, 3); D1(1, 1) = 1; D1(2:3, 2:3) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
402 %! D2 = [zeros(2, 1) S; zeros(1, 3)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
403 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (2, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
404 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
405 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
406 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
407 ## A: 5x3 complex rank deficient, B: 3x3 complex full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
408 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
409 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
410 %! A(:, 3) = 2*A(:, 1) - A(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
411 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
412 %! D1 = zeros (5, 3); D1(1:3, 1:3) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
413 %! D2 = S; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
414 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
415 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
416 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
417 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
418 ## A (5x3) and B (3x3) are both complex rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
419 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
420 %! B(:, 3) = 2*B(:, 1) - B(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
421 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
422 %! D1 = zeros (5, 2); D1(1:2, 1:2) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
423 %! D2 = [S; zeros(1, 2)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
424 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (2, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
425 %! assert (norm ((U'*A*X) - D1*[zeros(2, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
426 %! assert (norm ((V'*B*X) - D2*[zeros(2, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
427 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
428 ## A (now 3x5) complex and B (now 5x5) complex are full rank |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
429 ## now, A is 3x5 |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
430 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
431 %! A = A0.'; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
432 %! B0 = diag ([1 2 4 8 16]) + j*diag ([-5 4 -3 2 -1]); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
433 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
434 %! [U, V, X, C, S, R] = gsvd (A, B); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
435 %! D1 = [C zeros(3,2)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
436 %! D2 = [S zeros(3,2); zeros(2, 3) eye(2)]; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
437 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
438 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
439 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
440 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
441 ## A: 3x5 complex full rank, B: 5x5 complex rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
442 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
443 %! B(2, 2) = 0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
444 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
445 %! D1 = zeros (3, 5); D1(1, 1) = 1; D1(2:3, 2:3) = C; |
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
446 %! D2 = zeros (5,5); D2(1:2, 2:3) = S; D2(3:4, 4:5) = eye (2); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
447 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (2, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
448 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
449 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
450 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
451 ## A: 3x5 complex rank deficient, B: 5x5 complex full rank |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
452 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
453 %! B = B0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
454 %! A(3, :) = 2*A(1, :) - A(2, :); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
455 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
456 %! D1 = zeros (3, 5); D1(1:3, 1:3) = C; |
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
457 %! D2 = zeros (5,5); D2(1:3, 1:3) = S; D2(4:5, 4:5) = eye (2); |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
458 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
459 %! assert (norm ((U'*A*X) - D1*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
460 %! assert (norm ((V'*B*X) - D2*R) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
461 |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
462 ## A and B are both complex rank deficient |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
463 %!test <48807> |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
464 %! A = A0.'; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
465 %! B = B0.'; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
466 %! A(:, 3) = 2*A(:, 1) - A(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
467 %! B(:, 3) = 2*B(:, 1) - B(:, 2); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
468 %! [U, V, X, C, S, R] = gsvd (A, B); |
28915
c40a367a84c0
maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents:
27923
diff
changeset
|
469 %! D1 = zeros (3, 4); D1(1:3, 1:3) = C; |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
470 %! D2 = eye (4); D2(1:3, 1:3) = S; D2(5,:) = 0; |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
471 %! assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
472 %! assert (norm ((U'*A*X) - D1*[zeros(4, 1) R]) <= 1e-6); |
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
473 %! assert (norm ((V'*B*X) - D2*[zeros(4, 1) R]) <= 1e-6); |
22580
98eeed41f372
gsvd.cc: Update gsvd implementation in preparation for Matlab compatibility (bug #48807).
Rik <rik@octave.org>
parents:
22572
diff
changeset
|
474 |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
475 ## Test that single inputs produce single outputs |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
476 %!test |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
477 %! s = gsvd (single (ones (0,1)), B); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
478 %! assert (class (s), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
479 %! s = gsvd (single (ones (1,0)), B); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
480 %! assert (class (s), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
481 %! s = gsvd (single (ones (1,0)), B); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
482 %! [U,V,X,C,S,R] = gsvd (single ([]), B); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
483 %! assert (class (U), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
484 %! assert (class (V), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
485 %! assert (class (X), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
486 %! assert (class (C), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
487 %! assert (class (S), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
488 %! assert (class (R), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
489 %! |
22709
5c04055aa767
maint: Strip trailing whitespace from source files.
John W. Eaton <jwe@octave.org>
parents:
22664
diff
changeset
|
490 %! s = gsvd (single (A), B); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
491 %! assert (class (s), "single"); |
22709
5c04055aa767
maint: Strip trailing whitespace from source files.
John W. Eaton <jwe@octave.org>
parents:
22664
diff
changeset
|
492 %! [U,V,X,C,S,R] = gsvd (single (A), B); |
22664
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
493 %! assert (class (U), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
494 %! assert (class (V), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
495 %! assert (class (X), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
496 %! assert (class (C), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
497 %! assert (class (S), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
498 %! assert (class (R), "single"); |
f1bb2f0bcfec
Connect C++ gsvd to float version in liboctave (bug #48807).
Rik <rik@octave.org>
parents:
22580
diff
changeset
|
499 |
22572
d65ed83dfd13
gsvd: backout cset 73a85c6cacd1.
Carnë Draug <carandraug@octave.org>
parents:
diff
changeset
|
500 */ |
29958
32c3a5805893
move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
29359
diff
changeset
|
501 |
32c3a5805893
move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
29359
diff
changeset
|
502 OCTAVE_NAMESPACE_END |