[Octave-bug-tracker] [bug #50463] gsvd implementation uses functions removed from lapack >= 3.6.0

[Markus Mützel]

Follow-up Comment #33, bug #50463 (project octave):

Carlo, thank you for testing the patch on macOS. Could you run this test (one
of the BISTs) and see whether it gives any errors?

 A0 = randn (5, 3);
 B0 = diag ([1 2 4]);
 A = A0;
 B = B0;

## A (5x3) and B (3x3) are full rank
%!test <48807>
 [U, V, X, C, S, R] = gsvd (A, B);
 D1 = zeros (5, 3);  D1(1:3, 1:3) = C;
 D2 = S;
 assert (norm (diag (C).^2 + diag (S).^2 - ones (3, 1)) <= 1e-6);
 assert (norm ((U'*A*X) - D1*R) <= 1e-6);
 assert (norm ((V'*B*X) - D2*R) <= 1e-6);



    _______________________________________________________

Reply to this item at:

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

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