Mercurial > hg > octave-avbm
changeset 13829:8e25d6d1db10
ensure that complex matrix multiplication results are fully initialized
* CMatrix.cc (ComplexMatrix::xgemm): Ensure that all elements of
retval are initialized.
* fCMatrix.cc (FloatComplexMatrix::xgemm): Likewise.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Sat, 05 Nov 2011 02:02:54 -0400 |
parents | f1b023fd098d |
children | 462b5f556346 |
files | liboctave/CMatrix.cc liboctave/fCMatrix.cc |
diffstat | 2 files changed, 2 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/CMatrix.cc +++ b/liboctave/CMatrix.cc @@ -3797,7 +3797,7 @@ octave_idx_type lda = a.rows (), tda = a.cols (); octave_idx_type ldb = b.rows (), tdb = b.cols (); - retval = ComplexMatrix (a_nr, b_nc); + retval = ComplexMatrix (a_nr, b_nc, 0.0); Complex *c = retval.fortran_vec (); if (b_nc == 1 && a_nr == 1)
--- a/liboctave/fCMatrix.cc +++ b/liboctave/fCMatrix.cc @@ -3793,7 +3793,7 @@ octave_idx_type lda = a.rows (), tda = a.cols (); octave_idx_type ldb = b.rows (), tdb = b.cols (); - retval = FloatComplexMatrix (a_nr, b_nc); + retval = FloatComplexMatrix (a_nr, b_nc, 0.0); FloatComplex *c = retval.fortran_vec (); if (b_nc == 1 && a_nr == 1)