Mercurial > octave-nkf
comparison liboctave/numeric/fEIG.cc @ 20525:ff904ae0285b
eig: Return correct solution for a pair of hermitian matrices (bug #45511)
* liboctave/numeric/EIG.cc (EIG::init): Use correct form of hermitian_init.
* liboctave/numeric/fEIG.cc (FloatEIG::init): Likewise.
* libinterp/corefcn/eig.cc: Add %!test cases.
Thanks to Marco Caliari for identifying the fix.
| author | Mike Miller <mtmiller@octave.org> |
|---|---|
| date | Tue, 15 Sep 2015 08:48:35 -0400 |
| parents | 4197fc428c7d |
| children |
comparison
equal
deleted
inserted
replaced
| 20524:ba4088aee342 | 20525:ff904ae0285b |
|---|---|
| 710 info | 710 info |
| 711 F77_CHAR_ARG_LEN (1) | 711 F77_CHAR_ARG_LEN (1) |
| 712 F77_CHAR_ARG_LEN (1))); | 712 F77_CHAR_ARG_LEN (1))); |
| 713 | 713 |
| 714 if (a.is_hermitian () && b.is_hermitian () && info == 0) | 714 if (a.is_hermitian () && b.is_hermitian () && info == 0) |
| 715 return hermitian_init (a, calc_ev); | 715 return hermitian_init (a, b, calc_ev); |
| 716 | 716 |
| 717 FloatComplexMatrix atmp = a; | 717 FloatComplexMatrix atmp = a; |
| 718 FloatComplex *atmp_data = atmp.fortran_vec (); | 718 FloatComplex *atmp_data = atmp.fortran_vec (); |
| 719 | 719 |
| 720 FloatComplexMatrix btmp = b; | 720 FloatComplexMatrix btmp = b; |