Mercurial > octave

comparison liboctave/dMatrix.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
1525 operator * (const ColumnVector& v, const RowVector& a) 1525 operator * (const ColumnVector& v, const RowVector& a)
1526 { 1526 {
1527 Matrix retval; 1527 Matrix retval;
1528 1528
1529 int len = v.length (); 1529 int len = v.length ();
1530 int a_len = a.length (); 1530
1531 1531 if (len != 0)
1532 if (len != a_len) 1532 {
1533 gripe_nonconformant ("operator *", len, 1, 1, a_len); 1533 int a_len = a.length ();
1534 else 1534
1535 { 1535 retval.resize (len, a_len);
1536 if (len != 0) 1536 double *c = retval.fortran_vec ();
1537 { 1537
1538 retval.resize (len, a_len); 1538 F77_XFCN (dgemm, DGEMM, ("N", "N", len, a_len, 1, 1.0,
1539 double *c = retval.fortran_vec (); 1539 v.data (), len, a.data (), 1, 0.0,
1540 1540 c, len, 1L, 1L));
1541 F77_XFCN (dgemm, DGEMM, ("N", "N", len, a_len, 1, 1.0, 1541
1542 v.data (), len, a.data (), 1, 0.0, 1542 if (f77_exception_encountered)
1543 c, len, 1L, 1L)); 1543 (*current_liboctave_error_handler)
1544 1544 ("unrecoverable error in dgemm");
1545 if (f77_exception_encountered)
1546 (*current_liboctave_error_handler)
1547 ("unrecoverable error in dgemm");
1548 }
1549 } 1545 }
1550 1546
1551 return retval; 1547 return retval;
1552 } 1548 }
1553 1549