[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-----