Mercurial > octave-nkf

comparison liboctave/CMatrix.cc @ 3233:98d0ee053ba4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1999-01-27 20:23:40 by jwe]
author jwe
date Wed, 27 Jan 1999 20:23:46 +0000
parents 7aae2c3636a7
children 68259f410026
comparison
equal deleted inserted replaced
3232:3f6ae91e86b0 3233:98d0ee053ba4
1693 operator * (const ComplexColumnVector& v, const ComplexRowVector& a) 1693 operator * (const ComplexColumnVector& v, const ComplexRowVector& a)
1694 { 1694 {
1695 ComplexMatrix retval; 1695 ComplexMatrix retval;
1696 1696
1697 int len = v.length (); 1697 int len = v.length ();
1698 int a_len = a.length (); 1698
1699 1699 if (len != 0)
1700 if (len != a_len) 1700 {
1701 gripe_nonconformant ("operator *", len, 1, 1, a_len); 1701 int a_len = a.length ();
1702 else 1702
1703 { 1703 retval.resize (len, a_len);
1704 if (len != 0) 1704 Complex *c = retval.fortran_vec ();
1705 { 1705
1706 retval.resize (len, a_len); 1706 F77_XFCN (zgemm, ZGEMM, ("N", "N", len, a_len, 1, 1.0,
1707 Complex *c = retval.fortran_vec (); 1707 v.data (), len, a.data (), 1, 0.0,
1708 1708 c, len, 1L, 1L));
1709 F77_XFCN (zgemm, ZGEMM, ("N", "N", len, a_len, 1, 1.0, 1709
1710 v.data (), len, a.data (), 1, 0.0, 1710 if (f77_exception_encountered)
1711 c, len, 1L, 1L)); 1711 (*current_liboctave_error_handler)
1712 1712 ("unrecoverable error in zgemm");
1713 if (f77_exception_encountered)
1714 (*current_liboctave_error_handler)
1715 ("unrecoverable error in zgemm");
1716 }
1717 } 1713 }
1718 1714
1719 return retval; 1715 return retval;
1720 } 1716 }
1721 1717