Mercurial > hg > octave-max
comparison liboctave/CMatrix.cc @ 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 | a9d20c890190 |
children | 899122410ea5 |
comparison
equal
deleted
inserted
replaced
12062:26fb50dc7401 | 12063:4c7b51b2541d |
---|---|
3744 | 3744 |
3745 /* Simple Dot Product, Matrix-Vector and Matrix-Matrix Unit tests | 3745 /* Simple Dot Product, Matrix-Vector and Matrix-Matrix Unit tests |
3746 %!assert([1+i 2+i 3+i] * [ 4+i ; 5+i ; 6+i], 29+21i, 1e-14) | 3746 %!assert([1+i 2+i 3+i] * [ 4+i ; 5+i ; 6+i], 29+21i, 1e-14) |
3747 %!assert([1+i 2+i ; 3+i 4+i ] * [5+i ; 6+i], [15 + 14i ; 37 + 18i], 1e-14) | 3747 %!assert([1+i 2+i ; 3+i 4+i ] * [5+i ; 6+i], [15 + 14i ; 37 + 18i], 1e-14) |
3748 %!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) | 3748 %!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) |
3749 %!assert([1 i]*[i 0]', -i); | |
3749 */ | 3750 */ |
3750 | 3751 |
3751 /* Test some simple identities | 3752 /* Test some simple identities |
3752 %!shared M, cv, rv | 3753 %!shared M, cv, rv |
3753 %! M = randn(10,10)+i*rand(10,10); | 3754 %! M = randn(10,10)+i*rand(10,10); |
3837 if (conja == conjb) | 3838 if (conja == conjb) |
3838 { | 3839 { |
3839 F77_FUNC (xzdotu, XZDOTU) (a_nc, a.data (), 1, b.data (), 1, *c); | 3840 F77_FUNC (xzdotu, XZDOTU) (a_nc, a.data (), 1, b.data (), 1, *c); |
3840 if (conja) *c = std::conj (*c); | 3841 if (conja) *c = std::conj (*c); |
3841 } | 3842 } |
3842 else if (conjb) | 3843 else if (conja) |
3843 F77_FUNC (xzdotc, XZDOTC) (a_nc, a.data (), 1, b.data (), 1, *c); | 3844 F77_FUNC (xzdotc, XZDOTC) (a_nc, a.data (), 1, b.data (), 1, *c); |
3844 else | 3845 else |
3845 F77_FUNC (xzdotc, XZDOTC) (a_nc, b.data (), 1, a.data (), 1, *c); | 3846 F77_FUNC (xzdotc, XZDOTC) (a_nc, b.data (), 1, a.data (), 1, *c); |
3846 } | 3847 } |
3847 else if (b_nc == 1 && ! conjb) | 3848 else if (b_nc == 1 && ! conjb) |