Mercurial > hg > octave-lyh
diff liboctave/CMatrix.cc @ 9528:ec066ba012c8
more fixes & tests for matrix multiply
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Sun, 16 Aug 2009 08:25:24 +0200 |
parents | f3ebc728ffd7 |
children | 3d6a9aea2aea |
line wrap: on
line diff
--- a/liboctave/CMatrix.cc +++ b/liboctave/CMatrix.cc @@ -3755,7 +3755,11 @@ %! cv = randn(10,1)+i*rand(10,1); %! rv = randn(1,10)+i*rand(1,10); %!assert([M*cv,M*cv],M*[cv,cv],1e-14) +%!assert([M.'*cv,M.'*cv],M.'*[cv,cv],1e-14) +%!assert([M'*cv,M'*cv],M'*[cv,cv],1e-14) %!assert([rv*M;rv*M],[rv;rv]*M,1e-14) +%!assert([rv*M.';rv*M.'],[rv;rv]*M.',1e-14) +%!assert([rv*M';rv*M'],[rv;rv]*M',1e-14) %!assert(2*rv*cv,[rv,rv]*[cv;cv],1e-14) */ @@ -3853,7 +3857,7 @@ b.data (), 1, 0.0, c, 1 F77_CHAR_ARG_LEN (1))); } - else if (a_nr == 1 && ! conja) + else if (a_nr == 1 && ! conja && ! conjb) { const char *crevtransb = get_blas_trans_arg (! transb, conjb); F77_XFCN (zgemv, ZGEMV, (F77_CONST_CHAR_ARG2 (crevtransb, 1),