[bug #68434] Function gsvd crashes when matrices have 10 or more columns
Ian Proudler <[email protected]> Tue, 16 Jun 2026 10:25:51 -0400 (EDT)
| Newsgroups | gmane.comp.gnu.octave.bugs |
|---|---|
| Message-ID | <[email protected]> |
Please use the bug tracker to post updates to a bug report. The mailing list is intended as a read-only notification stream. Info posted to this mailing list address won't appear in the tracker database where it is most useful.
Follow-up Comment #18, bug #68434 (group octave):
I patched my version of 11.3.0 and ran some tests on random matrices. It all
looks good with one minor issue.
From a Linear Algebra point of view C and S are supposed to have R columns
where R = rank([A;B]). In fact they have the same number of columns as A & B.
I can believe that this is not how OpenBLAS works since it would require a
determination of rank. I cannot find my way around the OpenBLAS documentation
so can't check.
One consequence of this is that C'*C + S'*S != eye(columns(A)) as stated in
the Octave documentation. There are some zeros appended if R < columns(A).
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?68434>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
signature.asc
(application/pgp-signature, 228 B)
-----BEGIN PGP SIGNATURE----- iHUEABYIAB0WIQQk97aszIMMAvLLwm6qLAuaBUf3TgUCajFc7wAKCRCqLAuaBUf3 Tv1dAQDyqulRmO3HE/U/WPVfmpSTTT0Cv5MEELSXHh+oNYeAbgEAoH1berCzxbQN FfLC716CQzEUu1CtWdIOtFIRVC9Q9wc= =R6XT -----END PGP SIGNATURE-----