octave-bug-tracker
[Top][All Lists]

## [Octave-bug-tracker] [bug #48807] new gsvd function is incompatible w/Ma

 From: Nir Krakauer Subject: [Octave-bug-tracker] [bug #48807] new gsvd function is incompatible w/Matlab Date: Sun, 11 Sep 2016 20:39:13 +0000 (UTC) User-agent: Mozilla/5.0 (Macintosh; Intel Mac OS X 10_11_6) AppleWebKit/601.7.7 (KHTML, like Gecko) Version/9.1.2 Safari/601.7.7

```Follow-up Comment #6, bug #48807 (project octave):

Can someone please run the following test cases in Matlab?

A = [4 2 6; 3 5 7; 8 7 8; 7 4 1; 8 2 9];
A1 = [4 2 6; 3 5 1; 8 7 9; 7 4 10; 8 2 14];
[U, V, X, C, S] = gsvd(A, diag([1 2 4]))
[U, V, X, C, S] = gsvd(A, diag([1 0 4]))
[U, V, X, C, S] = gsvd(A1, diag([1 2 4]))
[U, V, X, C, S] = gsvd(A1, diag([1 0 4]))
[U, V, X, C, S] = gsvd(A', diag([1 2 4 8 16]))
[U, V, X, C, S] = gsvd(A', diag([1 0 4 8 16]))
[U, V, X, C, S] = gsvd(A1', diag([1 2 4 8 16]))
[U, V, X, C, S] = gsvd(A1', diag([1 0 4 8 16]))
Z = A + i*[8 8 1; 9 9 4; 1 7 5; 8 7 9; 3 1 5];
Z1 = A1 + i*[8 8 8; 9 9 9; 1 7 -5; 8 7 9; 3 1 5];
[U, V, X, C, S] = gsvd(Z, diag([1 2 4])+i*diag([4 -2 -1]))
[U, V, X, C, S] = gsvd(Z, diag([1 0 4])+i*diag([4 0 -1]))
[U, V, X, C, S] = gsvd(Z1, diag([1 2 4])+i*diag([4 -2 -1]))
[U, V, X, C, S] = gsvd(Z1, diag([1 0 4])+i*diag([4 0 -1]))
[U, V, X, C, S] = gsvd(Z', diag([1 2 4 8 16])+i*diag([-5 4 -3 2 -1]))
[U, V, X, C, S] = gsvd(Z', diag([1 0 4 8 16])+i*diag([-5 0 -3 2 -1]))
[U, V, X, C, S] = gsvd(Z1', diag([1 2 4 8 16])+i*diag([-5 4 -3 2 -1]))
[U, V, X, C, S] = gsvd(Z1', diag([1 0 4 8 16])+i*diag([-5 0 -3 2 -1]))

_______________________________________________________

<http://savannah.gnu.org/bugs/?48807>

_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/

```