Mercurial > hg > octave-max
changeset 12063:4c7b51b2541d release-3-2-x
fix typos in complex xgemm
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Sun, 23 Aug 2009 11:10:35 +0200 |
| parents | 26fb50dc7401 |
| children | 8e8da2c36e7e |
| files | liboctave/CMatrix.cc liboctave/ChangeLog liboctave/fCMatrix.cc |
| diffstat | 3 files changed, 9 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/CMatrix.cc +++ b/liboctave/CMatrix.cc @@ -3746,6 +3746,7 @@ %!assert([1+i 2+i 3+i] * [ 4+i ; 5+i ; 6+i], 29+21i, 1e-14) %!assert([1+i 2+i ; 3+i 4+i ] * [5+i ; 6+i], [15 + 14i ; 37 + 18i], 1e-14) %!assert([1+i 2+i ; 3+i 4+i ] * [5+i 6+i ; 7+i 8+i], [17 + 15i 20 + 17i; 41 + 19i 48 + 21i], 1e-14) +%!assert([1 i]*[i 0]', -i); */ /* Test some simple identities @@ -3839,7 +3840,7 @@ F77_FUNC (xzdotu, XZDOTU) (a_nc, a.data (), 1, b.data (), 1, *c); if (conja) *c = std::conj (*c); } - else if (conjb) + else if (conja) F77_FUNC (xzdotc, XZDOTC) (a_nc, a.data (), 1, b.data (), 1, *c); else F77_FUNC (xzdotc, XZDOTC) (a_nc, b.data (), 1, a.data (), 1, *c);
--- a/liboctave/ChangeLog +++ b/liboctave/ChangeLog @@ -1,3 +1,8 @@ +2009-08-16 Jaroslav Hajek <highegg@gmail.com> + + * CMatrix.cc (xgemm): Fix typo. Add test. + * fCMatrix.cc (xgemm): Ditto. + 2009-07-20 John W. Eaton <jwe@octave.org> * lo-ieee.cc (octave_ieee_init) [__NetBSD__]: Call nan to
--- a/liboctave/fCMatrix.cc +++ b/liboctave/fCMatrix.cc @@ -3739,6 +3739,7 @@ %!assert([1+i 2+i 3+i] * [ 4+i ; 5+i ; 6+i], 29+21i, 1e-14) %!assert([1+i 2+i ; 3+i 4+i ] * [5+i ; 6+i], [15 + 14i ; 37 + 18i], 1e-14) %!assert([1+i 2+i ; 3+i 4+i ] * [5+i 6+i ; 7+i 8+i], [17 + 15i 20 + 17i; 41 + 19i 48 + 21i], 1e-14) +%!assert([1 i]*[i 0]', -i); */ /* Test some simple identities @@ -3832,7 +3833,7 @@ F77_FUNC (xcdotu, XCDOTU) (a_nc, a.data (), 1, b.data (), 1, *c); if (conja) *c = std::conj (*c); } - else if (conjb) + else if (conja) F77_FUNC (xcdotc, XCDOTC) (a_nc, a.data (), 1, b.data (), 1, *c); else F77_FUNC (xcdotc, XCDOTC) (a_nc, b.data (), 1, a.data (), 1, *c);