[Octave-bug-tracker] [bug #60273] Wrong result in GSVD

[Leonardo]

Follow-up Comment #3, bug #60273 (project octave):

From

https://www.netlib.org/lapack/lug/node36.html

it seems that matrix X is just the matrix Q because

>> inv(C)*U'*A*X
ans =

   2.6026e+00  -3.5966e+00  -2.3962e+00
  -5.3213e-16  -3.8168e+00   5.0481e-01
  -7.7716e-16  -8.8818e-16   1.6511e+00

is upper triangular.

Also, matrix X is orthogonal (in general, it should not be)

>> X'*X
ans =

   1.0000e+00  -1.8041e-16  -5.5511e-17
  -1.8041e-16   1.0000e+00   1.2490e-16
  -5.5511e-17   1.2490e-16   1.0000e+00




    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?60273>

_______________________________________________
  Mensagem enviada pelo Savannah
  https://savannah.gnu.org/