# HG changeset patch
# User Jaroslav Hajek <highegg@gmail.com>
# Date 1250398515 -7200
# Node ID f3ebc728ffd7f23ae8c46a0dc2ea4ed25b124971
# Parent  691545147aca93c32d91dfcb5a85d2113daf8951
fix typos in complex xgemm

diff --git a/liboctave/CMatrix.cc b/liboctave/CMatrix.cc
--- 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);
diff --git a/liboctave/ChangeLog b/liboctave/ChangeLog
--- 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-08-14  John W. Eaton  <jwe@octave.org>
 
 	* Makefile.in (LINK_DEPS): Include BLAS_LIBS in the list.
diff --git a/liboctave/fCMatrix.cc b/liboctave/fCMatrix.cc
--- 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);