Mercurial > octave-nkf

comparison liboctave/array/CDiagMatrix.cc @ 18084:8e056300994b

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Follow coding convention of defining and initializing only 1 variable per line in liboctave. * liboctave/array/Array-b.cc, liboctave/array/Array-util.cc, liboctave/array/Array.cc, liboctave/array/CDiagMatrix.cc, liboctave/array/CMatrix.cc, liboctave/array/CSparse.cc, liboctave/array/MDiagArray2.cc, liboctave/array/MatrixType.cc, liboctave/array/PermMatrix.cc, liboctave/array/Sparse.cc, liboctave/array/Sparse.h, liboctave/array/boolSparse.cc, liboctave/array/dDiagMatrix.cc, liboctave/array/dMatrix.cc, liboctave/array/dSparse.cc, liboctave/array/dim-vector.cc, liboctave/array/fCDiagMatrix.cc, liboctave/array/fCMatrix.cc, liboctave/array/fDiagMatrix.cc, liboctave/array/fMatrix.cc, liboctave/array/idx-vector.cc, liboctave/array/idx-vector.h, liboctave/numeric/CmplxLU.cc, liboctave/numeric/CmplxQR.cc, liboctave/numeric/base-qr.cc, liboctave/numeric/bsxfun-defs.cc, liboctave/numeric/bsxfun.h, liboctave/numeric/dbleLU.cc, liboctave/numeric/dbleQR.cc, liboctave/numeric/fCmplxLU.cc, liboctave/numeric/fCmplxQR.cc, liboctave/numeric/floatLU.cc, liboctave/numeric/floatQR.cc, liboctave/numeric/lo-specfun.cc, liboctave/numeric/oct-convn.cc, liboctave/numeric/oct-norm.cc, liboctave/numeric/sparse-dmsolve.cc, liboctave/operators/mx-inlines.cc, liboctave/operators/mx-op-defs.h, liboctave/util/caseless-str.h, liboctave/util/kpse.cc, liboctave/util/lo-utils.cc, liboctave/util/oct-binmap.h, liboctave/util/oct-cmplx.h, liboctave/util/oct-inttypes.cc, liboctave/util/oct-inttypes.h, liboctave/util/oct-sort.cc: Follow coding convention of defining and initializing only 1 variable per line in liboctave.
author Rik <rik@octave.org>
date Wed, 04 Dec 2013 22:13:18 -0800
parents 49a5a4be04a1
children c5a101de2d88
comparison
equal deleted inserted replaced
18083:938f01339043 18084:8e056300994b
446 if (a_nc != b_nr) 446 if (a_nc != b_nr)
447 gripe_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc); 447 gripe_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc);
448 448
449 ComplexDiagMatrix c (a_nr, b_nc); 449 ComplexDiagMatrix c (a_nr, b_nc);
450 450
451 octave_idx_type len = c.length (), lenm = len < a_nc ? len : a_nc; 451 octave_idx_type len = c.length ();
452 octave_idx_type lenm = len < a_nc ? len : a_nc;
452 453
453 for (octave_idx_type i = 0; i < lenm; i++) 454 for (octave_idx_type i = 0; i < lenm; i++)
454 c.dgxelem (i) = a.dgelem (i) * b.dgelem (i); 455 c.dgxelem (i) = a.dgelem (i) * b.dgelem (i);
455 for (octave_idx_type i = lenm; i < len; i++) 456 for (octave_idx_type i = lenm; i < len; i++)
456 c.dgxelem (i) = 0.0; 457 c.dgxelem (i) = 0.0;
547 548
548 double 549 double
549 ComplexDiagMatrix::rcond (void) const 550 ComplexDiagMatrix::rcond (void) const
550 { 551 {
551 ColumnVector av = extract_diag (0).map<double> (std::abs); 552 ColumnVector av = extract_diag (0).map<double> (std::abs);
552 double amx = av.max (), amn = av.min (); 553 double amx = av.max ();
554 double amn = av.min ();
553 return amx == 0 ? 0.0 : amn / amx; 555 return amx == 0 ? 0.0 : amn / amx;
554 } 556 }
555 557
556 // i/o 558 // i/o
557 559