Mercurial > octave-nkf

changeset 13829:8e25d6d1db10

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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	Sat Nov 05 01:59:40 2011 -0400
+++ b/liboctave/CMatrix.cc	Sat Nov 05 02:02:54 2011 -0400
@@ -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	Sat Nov 05 01:59:40 2011 -0400
+++ b/liboctave/fCMatrix.cc	Sat Nov 05 02:02:54 2011 -0400
@@ -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)