Mercurial > hg > octave-lyh
annotate liboctave/dbleQR.h @ 9551:19d298e6f7e5
make ! operator check for NaNs, simplify implementations in liboctave
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Fri, 21 Aug 2009 08:18:16 +0200 |
| parents | d865363208d6 |
| children | 7918eb15040c |
| rev | line source |
|---|---|
| 457 | 1 /* |
| 2 | |
| 7017 | 3 Copyright (C) 1994, 1995, 1996, 1997, 2000, 2002, 2004, 2005, 2006, |
| 4 2007 John W. Eaton | |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
5 Copyright (C) 2008, 2009 Jaroslav Hajek |
| 457 | 6 |
| 7 This file is part of Octave. | |
| 8 | |
| 9 Octave is free software; you can redistribute it and/or modify it | |
| 10 under the terms of the GNU General Public License as published by the | |
| 7016 | 11 Free Software Foundation; either version 3 of the License, or (at your |
| 12 option) any later version. | |
| 457 | 13 |
| 14 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 16 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 17 for more details. | |
| 18 | |
| 19 You should have received a copy of the GNU General Public License | |
| 7016 | 20 along with Octave; see the file COPYING. If not, see |
| 21 <http://www.gnu.org/licenses/>. | |
| 457 | 22 |
| 23 */ | |
| 24 | |
| 25 #if !defined (octave_QR_h) | |
| 26 #define octave_QR_h 1 | |
| 27 | |
|
8950
d865363208d6
include <iosfwd> instead of <iostream> in header files
John W. Eaton <jwe@octave.org>
parents:
8597
diff
changeset
|
28 #include <iosfwd> |
| 457 | 29 |
| 30 #include "dMatrix.h" | |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
31 #include "dColVector.h" |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
32 #include "dRowVector.h" |
| 457 | 33 |
| 1881 | 34 class |
| 6108 | 35 OCTAVE_API |
| 1881 | 36 QR |
| 457 | 37 { |
| 38 public: | |
| 39 | |
| 539 | 40 enum type |
| 41 { | |
| 42 std, | |
| 43 raw, | |
| 2802 | 44 economy |
| 539 | 45 }; |
| 46 | |
| 1927 | 47 QR (void) : q (), r () { } |
| 457 | 48 |
| 2763 | 49 QR (const Matrix&, QR::type = QR::std); |
| 457 | 50 |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
51 QR (const Matrix& q, const Matrix& r); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
52 |
| 1881 | 53 QR (const QR& a) : q (a.q), r (a.r) { } |
| 457 | 54 |
| 1528 | 55 QR& operator = (const QR& a) |
| 56 { | |
| 1881 | 57 if (this != &a) |
| 58 { | |
| 59 q = a.q; | |
| 60 r = a.r; | |
| 61 } | |
| 1528 | 62 return *this; |
| 63 } | |
| 64 | |
| 1927 | 65 ~QR (void) { } |
| 1921 | 66 |
| 2763 | 67 void init (const Matrix&, QR::type); |
| 68 | |
| 1528 | 69 Matrix Q (void) const { return q; } |
| 70 | |
| 71 Matrix R (void) const { return r; } | |
| 457 | 72 |
|
8562
a6edd5c23cb5
use replacement methods if qrupdate is not available
Jaroslav Hajek <highegg@gmail.com>
parents:
8547
diff
changeset
|
73 QR::type get_type (void) const; |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
74 |
|
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
75 void update (const ColumnVector& u, const ColumnVector& v); |
|
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
76 |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
77 void update (const Matrix& u, const Matrix& v); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
78 |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
79 void insert_col (const ColumnVector& u, octave_idx_type j); |
|
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
80 |
|
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
81 void insert_col (const Matrix& u, const Array<octave_idx_type>& j); |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
82 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
83 void delete_col (octave_idx_type j); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
84 |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
85 void delete_col (const Array<octave_idx_type>& j); |
|
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
86 |
|
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
87 void insert_row (const RowVector& u, octave_idx_type j); |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
88 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
89 void delete_row (octave_idx_type j); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
90 |
|
7700
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7553
diff
changeset
|
91 void shift_cols (octave_idx_type i, octave_idx_type j); |
|
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7553
diff
changeset
|
92 |
| 3504 | 93 friend std::ostream& operator << (std::ostream&, const QR&); |
| 457 | 94 |
| 539 | 95 protected: |
| 457 | 96 |
|
8597
c86718093c1b
improve & fix QR classes
Jaroslav Hajek <highegg@gmail.com>
parents:
8562
diff
changeset
|
97 void form (octave_idx_type n, Matrix& afact, |
|
c86718093c1b
improve & fix QR classes
Jaroslav Hajek <highegg@gmail.com>
parents:
8562
diff
changeset
|
98 double *tau, QR::type qr_type); |
|
c86718093c1b
improve & fix QR classes
Jaroslav Hajek <highegg@gmail.com>
parents:
8562
diff
changeset
|
99 |
| 457 | 100 Matrix q; |
| 101 Matrix r; | |
| 102 }; | |
| 103 | |
|
8562
a6edd5c23cb5
use replacement methods if qrupdate is not available
Jaroslav Hajek <highegg@gmail.com>
parents:
8547
diff
changeset
|
104 #ifndef HAVE_QRUPDATE |
|
a6edd5c23cb5
use replacement methods if qrupdate is not available
Jaroslav Hajek <highegg@gmail.com>
parents:
8547
diff
changeset
|
105 void warn_qrupdate_once (void); |
|
a6edd5c23cb5
use replacement methods if qrupdate is not available
Jaroslav Hajek <highegg@gmail.com>
parents:
8547
diff
changeset
|
106 #endif |
|
a6edd5c23cb5
use replacement methods if qrupdate is not available
Jaroslav Hajek <highegg@gmail.com>
parents:
8547
diff
changeset
|
107 |
| 457 | 108 #endif |
| 109 | |
| 110 /* | |
| 111 ;;; Local Variables: *** | |
| 112 ;;; mode: C++ *** | |
| 113 ;;; End: *** | |
| 114 */ |