gsvd test error

[Thomas D. Dean]

>> ver()
----------------------------------------------------------------------
GNU Octave Version: 10.0.0 (hg id: 3a0d7c21384c+)
GNU Octave License: GNU General Public License
Operating System: Linux 6.5.0-17-generic #17~22.04.1-Ubuntu SMP 
PREEMPT_DYNAMIC Tue Jan 16 14:32:32 UTC 2 x86_64
----------------------------------------------------------------------
>> version("-hgid")
>> A = reshape (1:15,5,3);
B = magic (3);
[U,V,X,C,S] = gsvd (A,B);
assert (size (U), [5, 5]);
assert (size (V), [3, 3]);
assert (size (X), [3, 3]);
assert (size (C), [5, 3]);
assert (C(4:5, :), zeros (2,3));
assert (size (S), [3, 3]);
assert (U*C*X', A, 50*eps);
assert (V*S*X', B, 50*eps);
S0 = gsvd (A, B);
assert (size (S0), [3, 1]);
S1 = sort (svd (A / B));
assert (S0, S1, 10*eps)
error: ASSERT errors for:  assert (S0,S1,10 * eps)

  Location  |  Observed  |  Expected  |  Reason
    (3)         5.0123       5.0123      Abs err 2.6645e-15 exceeds tol 
2.2204e-15 by 4e-16

Tom Dean