[bug #48807] gsvd function needs matlab compatible interface

"Dmitri A. Sergatskov" <[email protected]> Wed, 24 Jun 2026 06:20:28 -0400 (EDT)
Newsgroups gmane.comp.gnu.octave.bugs
Message-ID <[email protected]>
--===============1561491888686833974==
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline

UGxlYXNlIHVzZSB0aGUgYnVnIHRyYWNrZXIgdG8gcG9zdCB1cGRhdGVzIHRvIGEgYnVnIHJlcG9y
dC4gIFRoZSBtYWlsaW5nIGxpc3QgaXMgaW50ZW5kZWQgYXMgYSByZWFkLW9ubHkgbm90aWZpY2F0
aW9uIHN0cmVhbS4gIEluZm8gcG9zdGVkIHRvIHRoaXMgbWFpbGluZyBsaXN0IGFkZHJlc3Mgd29u
J3QgYXBwZWFyIGluIHRoZSB0cmFja2VyIGRhdGFiYXNlIHdoZXJlIGl0IGlzIG1vc3QgdXNlZnVs
LgoK

--===============1561491888686833974==
Content-Type: MULTIPART/signed; micalg=pgp-sha256; protocol="application/pgp-signature"; BOUNDARY="8323329-1726956429-1782296428=:1050108"
Content-Disposition: inline

--8323329-1726956429-1782296428=:1050108
Content-Type: TEXT/plain; CHARSET=utf-8
Content-Transfer-Encoding: QUOTED-PRINTABLE
Content-Disposition: inline

Follow-up Comment #30, bug #48807 (group octave):=0A=0AOne potential Matlab=
 compatibility issue. With posted patch I see:=0A=0AM =3D 4; P =3D 2; N =3D=
 8;=0AA =3D [1:8; 1:8; 1:8; 1:8];=0AB =3D [8:-1:1; 8:-1:1];=0A[U, V, X, C, =
S] =3D gsvd (A, B)=0A=0A=0AU =3D=0A=0A  -8.6603e-01  -5.0000e-01  -5.6156e-=
17  -2.5375e-17=0A   2.8868e-01  -5.0000e-01   8.1650e-01  -9.0181e-17=0A  =
 2.8868e-01  -5.0000e-01  -4.0825e-01  -7.0711e-01=0A   2.8868e-01  -5.0000=
e-01  -4.0825e-01   7.0711e-01=0A=0AV =3D=0A=0A  -0.7071  -0.7071=0A  -0.70=
71   0.7071=0A=0AX =3D=0A=0A         0         0         0         0       =
  0         0  -11.3137=0A-2.0000=0A         0         0         0         =
0         0         0   -9.8995=0A-4.0000=0A         0         0         0 =
        0         0         0   -8.4853=0A-6.0000=0A         0         0   =
      0         0         0         0   -7.0711=0A-8.0000=0A         0     =
    0         0         0         0         0   -5.6569=0A-10.0000=0A      =
   0         0         0         0         0         0   -4.2426=0A-12.0000=
=0A         0         0         0         0         0         0   -2.8284=
=0A-14.0000=0A         0         0         0         0         0         0 =
  -1.4142=0A-16.0000=0A=0AC =3D=0A=0A        0        0        0        0  =
      0        0   0.0000        0=0A        0        0        0        0  =
      0        0        0   1.0000=0A        0        0        0        0  =
      0        0        0        0=0A        0        0        0        0  =
      0        0        0        0=0A=0AS =3D=0A=0A   0   0   0   0   0   0=
   1   0=0A   0   0   0   0   0   0   0   0=0A=0A=0A=0ASo identity C^2+S^2 =
=3D I does not hold:=0A=0A=0Aoctave:8> C'*C + S'*S=0Aans =3D=0A=0A   0   0 =
  0   0   0   0   0   0=0A   0   0   0   0   0   0   0   0=0A   0   0   0  =
 0   0   0   0   0=0A   0   0   0   0   0   0   0   0=0A   0   0   0   0   =
0   0   0   0=0A   0   0   0   0   0   0   0   0=0A   0   0   0   0   0   0=
   1   0=0A   0   0   0   0   0   0   0   1=0A=0A=0A=0Abut Matlab may "chea=
t" and fill it up. Can anyone confirm one way or another?=0A=0ADmitri.=0A--=
 =0A=0A=0A=0A=0A    _______________________________________________________=
=0A=0AReply to this item at:=0A=0A  <https://savannah.gnu.org/bugs/?48807>=
=0A=0A_______________________________________________=0AMessage sent via Sa=
vannah=0Ahttps://savannah.gnu.org/=0A
--8323329-1726956429-1782296428=:1050108
Content-Type: APPLICATION/pgp-signature; name=signature.asc

-----BEGIN PGP SIGNATURE-----

iHUEABYIAB0WIQQk97aszIMMAvLLwm6qLAuaBUf3TgUCajuvbAAKCRCqLAuaBUf3
ToPnAP9Jx9RG0Cld5asFGpoPIaFxjEOq3CwmVcAb2X7R8UMd9wD/RO21BTYtg8mD
3kT8tNSQK52qG62CR8Lh9shvg9M3iw4=
=7snJ
-----END PGP SIGNATURE-----

--8323329-1726956429-1782296428=:1050108--



--===============1561491888686833974==--