octave: liboctave/numeric/base-lu.cc annotate

annotate liboctave/numeric/base-lu.cc @ 17769:49a5a4be04a1

maint: Use GNU style coding conventions for code in liboctave/ * liboctave/array/Array-C.cc, liboctave/array/Array-b.cc, liboctave/array/Array-ch.cc, liboctave/array/Array-d.cc, liboctave/array/Array-f.cc, liboctave/array/Array-fC.cc, liboctave/array/Array-util.cc, liboctave/array/Array-util.h, liboctave/array/Array.cc, liboctave/array/Array.h, liboctave/array/Array3.h, liboctave/array/CColVector.cc, liboctave/array/CColVector.h, liboctave/array/CDiagMatrix.cc, liboctave/array/CDiagMatrix.h, liboctave/array/CMatrix.cc, liboctave/array/CMatrix.h, liboctave/array/CNDArray.cc, liboctave/array/CNDArray.h, liboctave/array/CRowVector.cc, liboctave/array/CRowVector.h, liboctave/array/CSparse.cc, liboctave/array/CSparse.h, liboctave/array/DiagArray2.h, liboctave/array/MArray.cc, liboctave/array/MArray.h, liboctave/array/MDiagArray2.cc, liboctave/array/MDiagArray2.h, liboctave/array/MSparse.cc, liboctave/array/MSparse.h, liboctave/array/MatrixType.cc, liboctave/array/MatrixType.h, liboctave/array/PermMatrix.h, liboctave/array/Range.cc, liboctave/array/Range.h, liboctave/array/Sparse.cc, liboctave/array/Sparse.h, liboctave/array/boolMatrix.cc, liboctave/array/boolMatrix.h, liboctave/array/boolNDArray.cc, liboctave/array/boolNDArray.h, liboctave/array/boolSparse.cc, liboctave/array/boolSparse.h, liboctave/array/chMatrix.cc, liboctave/array/chMatrix.h, liboctave/array/chNDArray.cc, liboctave/array/chNDArray.h, liboctave/array/dColVector.h, liboctave/array/dDiagMatrix.cc, liboctave/array/dDiagMatrix.h, liboctave/array/dMatrix.cc, liboctave/array/dMatrix.h, liboctave/array/dNDArray.cc, liboctave/array/dNDArray.h, liboctave/array/dRowVector.h, liboctave/array/dSparse.cc, liboctave/array/dSparse.h, liboctave/array/dim-vector.cc, liboctave/array/dim-vector.h, liboctave/array/fCColVector.cc, liboctave/array/fCColVector.h, liboctave/array/fCDiagMatrix.cc, liboctave/array/fCDiagMatrix.h, liboctave/array/fCMatrix.cc, liboctave/array/fCMatrix.h, liboctave/array/fCNDArray.cc, liboctave/array/fCNDArray.h, liboctave/array/fCRowVector.cc, liboctave/array/fCRowVector.h, liboctave/array/fColVector.h, liboctave/array/fDiagMatrix.cc, liboctave/array/fDiagMatrix.h, liboctave/array/fMatrix.cc, liboctave/array/fMatrix.h, liboctave/array/fNDArray.cc, liboctave/array/fNDArray.h, liboctave/array/fRowVector.h, liboctave/array/idx-vector.cc, liboctave/array/idx-vector.h, liboctave/array/intNDArray.cc, liboctave/array/intNDArray.h, liboctave/cruft/misc/blaswrap.c, liboctave/cruft/misc/quit.cc, liboctave/numeric/CmplxCHOL.cc, liboctave/numeric/CmplxCHOL.h, liboctave/numeric/CmplxGEPBAL.cc, liboctave/numeric/CmplxGEPBAL.h, liboctave/numeric/CmplxHESS.h, liboctave/numeric/CmplxLU.cc, liboctave/numeric/CmplxLU.h, liboctave/numeric/CmplxQR.cc, liboctave/numeric/CmplxQRP.cc, liboctave/numeric/CmplxQRP.h, liboctave/numeric/CmplxSCHUR.h, liboctave/numeric/CmplxSVD.cc, liboctave/numeric/CmplxSVD.h, liboctave/numeric/CollocWt.h, liboctave/numeric/DAE.h, liboctave/numeric/DAEFunc.h, liboctave/numeric/DAERT.h, liboctave/numeric/DAERTFunc.h, liboctave/numeric/DASPK.cc, liboctave/numeric/DASRT.cc, liboctave/numeric/DASRT.h, liboctave/numeric/DASSL.cc, liboctave/numeric/DET.h, liboctave/numeric/EIG.cc, liboctave/numeric/EIG.h, liboctave/numeric/LSODE.cc, liboctave/numeric/ODE.h, liboctave/numeric/ODEFunc.h, liboctave/numeric/ODES.h, liboctave/numeric/ODESFunc.h, liboctave/numeric/Quad.cc, liboctave/numeric/Quad.h, liboctave/numeric/SparseCmplxCHOL.h, liboctave/numeric/SparseCmplxLU.cc, liboctave/numeric/SparseCmplxLU.h, liboctave/numeric/SparseCmplxQR.cc, liboctave/numeric/SparseCmplxQR.h, liboctave/numeric/SparseQR.cc, liboctave/numeric/SparseQR.h, liboctave/numeric/SparsedbleCHOL.h, liboctave/numeric/SparsedbleLU.cc, liboctave/numeric/SparsedbleLU.h, liboctave/numeric/base-aepbal.h, liboctave/numeric/base-dae.h, liboctave/numeric/base-de.h, liboctave/numeric/base-lu.cc, liboctave/numeric/base-lu.h, liboctave/numeric/base-min.h, liboctave/numeric/base-qr.h, liboctave/numeric/bsxfun.h, liboctave/numeric/dbleCHOL.cc, liboctave/numeric/dbleCHOL.h, liboctave/numeric/dbleGEPBAL.h, liboctave/numeric/dbleHESS.h, liboctave/numeric/dbleLU.cc, liboctave/numeric/dbleLU.h, liboctave/numeric/dbleQR.cc, liboctave/numeric/dbleQRP.cc, liboctave/numeric/dbleQRP.h, liboctave/numeric/dbleSCHUR.cc, liboctave/numeric/dbleSCHUR.h, liboctave/numeric/dbleSVD.cc, liboctave/numeric/dbleSVD.h, liboctave/numeric/eigs-base.cc, liboctave/numeric/fCmplxAEPBAL.cc, liboctave/numeric/fCmplxAEPBAL.h, liboctave/numeric/fCmplxCHOL.cc, liboctave/numeric/fCmplxCHOL.h, liboctave/numeric/fCmplxGEPBAL.cc, liboctave/numeric/fCmplxGEPBAL.h, liboctave/numeric/fCmplxHESS.h, liboctave/numeric/fCmplxLU.cc, liboctave/numeric/fCmplxLU.h, liboctave/numeric/fCmplxQR.cc, liboctave/numeric/fCmplxQR.h, liboctave/numeric/fCmplxQRP.cc, liboctave/numeric/fCmplxQRP.h, liboctave/numeric/fCmplxSCHUR.cc, liboctave/numeric/fCmplxSCHUR.h, liboctave/numeric/fCmplxSVD.h, liboctave/numeric/fEIG.cc, liboctave/numeric/fEIG.h, liboctave/numeric/floatCHOL.cc, liboctave/numeric/floatCHOL.h, liboctave/numeric/floatGEPBAL.cc, liboctave/numeric/floatGEPBAL.h, liboctave/numeric/floatHESS.h, liboctave/numeric/floatLU.cc, liboctave/numeric/floatLU.h, liboctave/numeric/floatQR.cc, liboctave/numeric/floatQRP.cc, liboctave/numeric/floatQRP.h, liboctave/numeric/floatSCHUR.cc, liboctave/numeric/floatSCHUR.h, liboctave/numeric/floatSVD.cc, liboctave/numeric/floatSVD.h, liboctave/numeric/lo-mappers.cc, liboctave/numeric/lo-mappers.h, liboctave/numeric/lo-specfun.cc, liboctave/numeric/lo-specfun.h, liboctave/numeric/oct-convn.cc, liboctave/numeric/oct-fftw.cc, liboctave/numeric/oct-fftw.h, liboctave/numeric/oct-norm.cc, liboctave/numeric/oct-rand.cc, liboctave/numeric/oct-rand.h, liboctave/numeric/randgamma.c, liboctave/numeric/randgamma.h, liboctave/numeric/randmtzig.c, liboctave/numeric/randpoisson.c, liboctave/numeric/randpoisson.h, liboctave/numeric/sparse-base-chol.h, liboctave/numeric/sparse-base-lu.h, liboctave/numeric/sparse-dmsolve.cc, liboctave/operators/Sparse-diag-op-defs.h, liboctave/operators/Sparse-op-defs.h, liboctave/operators/mx-inlines.cc, liboctave/system/dir-ops.h, liboctave/system/file-ops.cc, liboctave/system/file-stat.cc, liboctave/system/file-stat.h, liboctave/system/lo-sysdep.cc, liboctave/system/lo-sysdep.h, liboctave/system/mach-info.cc, liboctave/system/mach-info.h, liboctave/system/oct-env.cc, liboctave/system/oct-group.cc, liboctave/system/oct-syscalls.cc, liboctave/system/oct-syscalls.h, liboctave/system/oct-time.h, liboctave/system/tempname.c, liboctave/util/action-container.h, liboctave/util/base-list.h, liboctave/util/cmd-edit.cc, liboctave/util/cmd-edit.h, liboctave/util/cmd-hist.cc, liboctave/util/cmd-hist.h, liboctave/util/data-conv.cc, liboctave/util/data-conv.h, liboctave/util/kpse.cc, liboctave/util/lo-array-gripes.cc, liboctave/util/lo-cieee.c, liboctave/util/lo-regexp.cc, liboctave/util/lo-utils.cc, liboctave/util/oct-alloc.cc, liboctave/util/oct-base64.cc, liboctave/util/oct-binmap.h, liboctave/util/oct-cmplx.h, liboctave/util/oct-glob.cc, liboctave/util/oct-inttypes.cc, liboctave/util/oct-inttypes.h, liboctave/util/oct-locbuf.cc, liboctave/util/oct-locbuf.h, liboctave/util/oct-mem.h, liboctave/util/oct-mutex.cc, liboctave/util/oct-refcount.h, liboctave/util/oct-shlib.cc, liboctave/util/oct-shlib.h, liboctave/util/oct-sort.cc, liboctave/util/oct-sort.h, liboctave/util/pathsearch.cc, liboctave/util/pathsearch.h, liboctave/util/sparse-util.cc, liboctave/util/str-vec.cc, liboctave/util/str-vec.h, liboctave/util/unwind-prot.h, liboctave/util/url-transfer.cc, liboctave/util/url-transfer.h: Use GNU style coding conventions.

author	Rik <rik@octave.org>
date	Sat, 26 Oct 2013 18:57:05 -0700
parents	d63878346099
children	aa9ca67f09fb

rev	line source
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	1 /*
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	2
17744 d63878346099 maint: Update copyright notices for release. John W. Eaton <jwe@octave.org> parents: 15271 diff changeset	3 Copyright (C) 1996-2013 John W. Eaton
10521 4d1fc073fbb7 add some missing copyright stmts Jaroslav Hajek <highegg@gmail.com> parents: 10350 diff changeset	4 Copyright (C) 2009 VZLU Prague
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	5
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	7
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	11 option) any later version.
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	12
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	16 for more details.
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	17
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	18 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5307 diff changeset	20 <http://www.gnu.org/licenses/>.
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	21
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	22 */
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	23
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	24 #ifdef HAVE_CONFIG_H
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	25 #include <config.h>
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	26 #endif
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	27
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	28 #include "base-lu.h"
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	29
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	30 template <class lu_type>
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	31 base_lu<lu_type>::base_lu (const lu_type& l, const lu_type& u,
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	32 const PermMatrix& p)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	33 : a_fact (u), l_fact (l), ipvt (p.pvec ())
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	34 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	35 if (l.columns () != u.rows ())
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	36 (*current_liboctave_error_handler) ("lu: dimension mismatch");
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	37 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	38
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	39 template <class lu_type>
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	40 bool
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	41 base_lu <lu_type> :: packed (void) const
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	42 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	43 return l_fact.dims () == dim_vector ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	44 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	45
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	46 template <class lu_type>
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	47 void
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	48 base_lu <lu_type> :: unpack (void)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	49 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	50 if (packed ())
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	51 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	52 l_fact = L ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	53 a_fact = U (); // FIXME: sub-optimal
9697 51c17bd18563 omission in base-lu::unpack Jaroslav Hajek <highegg@gmail.com> parents: 9694 diff changeset	54 ipvt = getp ();
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	55 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	56 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	57
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	58 template <class lu_type>
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	59 lu_type
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	60 base_lu <lu_type> :: L (void) const
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	61 {
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	62 if (packed ())
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	63 {
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	64 octave_idx_type a_nr = a_fact.rows ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	65 octave_idx_type a_nc = a_fact.cols ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	66 octave_idx_type mn = (a_nr < a_nc ? a_nr : a_nc);
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	67
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	68 lu_type l (a_nr, mn, lu_elt_type (0.0));
4811 198f3712c692 [project @ 2004-03-03 19:24:33 by jwe] jwe parents: 4329 diff changeset	69
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	70 for (octave_idx_type i = 0; i < a_nr; i++)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	71 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	72 if (i < a_nc)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	73 l.xelem (i, i) = 1.0;
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	74
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	75 for (octave_idx_type j = 0; j < (i < a_nc ? i : a_nc); j++)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	76 l.xelem (i, j) = a_fact.xelem (i, j);
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	77 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	78
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	79 return l;
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	80 }
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	81 else
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	82 return l_fact;
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	83 }
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	84
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	85 template <class lu_type>
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	86 lu_type
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	87 base_lu <lu_type> :: U (void) const
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	88 {
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	89 if (packed ())
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	90 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	91 octave_idx_type a_nr = a_fact.rows ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	92 octave_idx_type a_nc = a_fact.cols ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	93 octave_idx_type mn = (a_nr < a_nc ? a_nr : a_nc);
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	94
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	95 lu_type u (mn, a_nc, lu_elt_type (0.0));
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	96
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	97 for (octave_idx_type i = 0; i < mn; i++)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	98 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	99 for (octave_idx_type j = i; j < a_nc; j++)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	100 u.xelem (i, j) = a_fact.xelem (i, j);
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	101 }
4329 d53c33d93440 [project @ 2003-02-18 20:00:48 by jwe] jwe parents: 4192 diff changeset	102
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	103 return u;
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	104 }
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	105 else
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	106 return a_fact;
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	107 }
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	108
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	109 template <class lu_type>
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	110 lu_type
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	111 base_lu <lu_type> :: Y (void) const
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	112 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	113 if (! packed ())
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	114 (*current_liboctave_error_handler)
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	115 ("lu: Y () not implemented for unpacked form");
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	116 return a_fact;
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	117 }
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	118
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	119 template <class lu_type>
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	120 Array<octave_idx_type>
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	121 base_lu <lu_type> :: getp (void) const
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	122 {
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	123 if (packed ())
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	124 {
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	125 octave_idx_type a_nr = a_fact.rows ();
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	126
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	127 Array<octave_idx_type> pvt (dim_vector (a_nr, 1));
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	128
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	129 for (octave_idx_type i = 0; i < a_nr; i++)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	130 pvt.xelem (i) = i;
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	131
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	132 for (octave_idx_type i = 0; i < ipvt.length (); i++)
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	133 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	134 octave_idx_type k = ipvt.xelem (i);
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	135
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	136 if (k != i)
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	137 {
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	138 octave_idx_type tmp = pvt.xelem (k);
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	139 pvt.xelem (k) = pvt.xelem (i);
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	140 pvt.xelem (i) = tmp;
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	141 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	142 }
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	143
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	144 return pvt;
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	145 }
9694 50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	146 else
50db3c5175b5 allow unpacked form of LU Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	147 return ipvt;
1991 413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	148 }
413f6a81868f [project @ 1996-03-03 00:40:53 by jwe] jwe parents: diff changeset	149
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	150 template <class lu_type>
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	151 PermMatrix
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	152 base_lu <lu_type> :: P (void) const
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	153 {
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	154 return PermMatrix (getp (), false);
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	155 }
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	156
445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	157 template <class lu_type>
7515 f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	158 ColumnVector
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	159 base_lu <lu_type> :: P_vec (void) const
7515 f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	160 {
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	161 octave_idx_type a_nr = a_fact.rows ();
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	162
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	163 ColumnVector p (a_nr);
7515 f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	164
8367 445d27d79f4e support permutation matrix objects Jaroslav Hajek <highegg@gmail.com> parents: 7515 diff changeset	165 Array<octave_idx_type> pvt = getp ();
7515 f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	166
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	167 for (octave_idx_type i = 0; i < a_nr; i++)
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	168 p.xelem (i) = static_cast<double> (pvt.xelem (i) + 1);
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	169
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	170 return p;
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	171 }
f3c00dc0912b Eliminate the rest of the dispatched sparse functions David Bateman <dbateman@free.fr> parents: 7017 diff changeset	172
9715 9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	173 template <class lu_type>
9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	174 bool
9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	175 base_lu<lu_type>::regular (void) const
9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	176 {
13217 08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	177 bool retval = true;
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	178
9715 9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	179 octave_idx_type k = std::min (a_fact.rows (), a_fact.columns ());
13217 08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	180
9715 9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	181 for (octave_idx_type i = 0; i < k; i++)
13217 08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	182 {
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	183 if (a_fact(i, i) == lu_elt_type ())
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	184 {
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	185 retval = false;
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	186 break;
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	187 }
08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	188 }
9715 9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	189
13217 08ba694723f5 fix warnings for unused but set variables John W. Eaton <jwe@octave.org> parents: 11590 diff changeset	190 return retval;
9715 9f27172fbd1e auto-set MatrixType from certain functions Jaroslav Hajek <highegg@gmail.com> parents: 9697 diff changeset	191 }

Mercurial > octave

annotate liboctave/numeric/base-lu.cc @ 17769:49a5a4be04a1