Mercurial > hg > octave-nkf

comparison liboctave/fCMatrix.cc @ 11596:a9cf422ed849

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
avoid apparent bug in ATLAS versions of CHERK/ZHERK
author John W. Eaton <jwe@octave.org>
date Fri, 21 Jan 2011 11:05:21 -0500
parents 12df7854fa7c
children 353c71c76f22
comparison
equal deleted inserted replaced
11595:5ec6aa05638d 11596:a9cf422ed849
3748 retval = FloatComplexMatrix (a_nr, b_nc, 0.0); 3748 retval = FloatComplexMatrix (a_nr, b_nc, 0.0);
3749 else if (a.data () == b.data () && a_nr == b_nc && tra != trb) 3749 else if (a.data () == b.data () && a_nr == b_nc && tra != trb)
3750 { 3750 {
3751 octave_idx_type lda = a.rows (); 3751 octave_idx_type lda = a.rows ();
3752 3752
3753 retval = FloatComplexMatrix (a_nr, b_nc); 3753 // FIXME -- looking at the reference BLAS, it appears that it
3754 // should not be necessary to initialize the output matrix if
3755 // BETA is 0 in the call to CHERK, but ATLAS appears to
3756 // use the result matrix before zeroing the elements.
3757
3758 retval = FloatComplexMatrix (a_nr, b_nc, 0.0);
3754 FloatComplex *c = retval.fortran_vec (); 3759 FloatComplex *c = retval.fortran_vec ();
3755 3760
3756 const char ctra = get_blas_trans_arg (tra, cja); 3761 const char ctra = get_blas_trans_arg (tra, cja);
3757 if (cja || cjb) 3762 if (cja || cjb)
3758 { 3763 {