Mercurial > octave
annotate liboctave/CmplxQR.h @ 7554:40574114c514
implement Cholesky factorization updating
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Tue, 04 Mar 2008 22:25:50 -0500 |
| parents | 56be6f31dd4e |
| children | efccca5f2ad7 |
| rev | line source |
|---|---|
| 457 | 1 /* |
| 2 | |
| 7017 | 3 Copyright (C) 1994, 1995, 1996, 1997, 2000, 2002, 2004, 2005, 2006, |
| 4 2007 John W. Eaton | |
| 457 | 5 |
| 6 This file is part of Octave. | |
| 7 | |
| 8 Octave is free software; you can redistribute it and/or modify it | |
| 9 under the terms of the GNU General Public License as published by the | |
| 7016 | 10 Free Software Foundation; either version 3 of the License, or (at your |
| 11 option) any later version. | |
| 457 | 12 |
| 13 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 16 for more details. | |
| 17 | |
| 18 You should have received a copy of the GNU General Public License | |
| 7016 | 19 along with Octave; see the file COPYING. If not, see |
| 20 <http://www.gnu.org/licenses/>. | |
| 457 | 21 |
| 22 */ | |
| 23 | |
|
7554
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7553
diff
changeset
|
24 // updating/downdating by Jaroslav Hajek 2008 |
|
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7553
diff
changeset
|
25 |
| 457 | 26 #if !defined (octave_ComplexQR_h) |
| 27 #define octave_ComplexQR_h 1 | |
| 28 | |
| 3503 | 29 #include <iostream> |
| 457 | 30 |
| 31 #include "CMatrix.h" | |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
32 #include "CColVector.h" |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
33 #include "CRowVector.h" |
| 539 | 34 #include "dbleQR.h" |
| 457 | 35 |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
36 |
| 1881 | 37 class |
| 6108 | 38 OCTAVE_API |
| 1881 | 39 ComplexQR |
| 457 | 40 { |
| 41 public: | |
| 42 | |
| 1927 | 43 ComplexQR (void) : q (), r () { } |
| 457 | 44 |
| 2763 | 45 ComplexQR (const ComplexMatrix&, QR::type = QR::std); |
| 457 | 46 |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
47 ComplexQR (const ComplexMatrix& q, const ComplexMatrix& r); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
48 |
| 1881 | 49 ComplexQR (const ComplexQR& a) : q (a.q), r (a.r) { } |
| 457 | 50 |
| 1528 | 51 ComplexQR& operator = (const ComplexQR& a) |
| 52 { | |
| 1881 | 53 if (this != &a) |
| 54 { | |
| 55 q = a.q; | |
| 56 r = a.r; | |
| 57 } | |
| 1528 | 58 return *this; |
| 59 } | |
| 60 | |
| 1927 | 61 ~ComplexQR (void) { } |
| 1921 | 62 |
| 2763 | 63 void init (const ComplexMatrix&, QR::type = QR::std); |
| 64 | |
| 1528 | 65 ComplexMatrix Q (void) const { return q; } |
| 2763 | 66 |
| 1528 | 67 ComplexMatrix R (void) const { return r; } |
| 457 | 68 |
|
7553
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
69 void update (const ComplexMatrix& u, const ComplexMatrix& v); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
70 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
71 void insert_col (const ComplexMatrix& u, octave_idx_type j); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
72 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
73 void delete_col (octave_idx_type j); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
74 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
75 void insert_row (const ComplexMatrix& u, octave_idx_type j); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
76 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
77 void delete_row (octave_idx_type j); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
78 |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
79 void economize(); |
|
56be6f31dd4e
implementation of QR factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
80 |
| 3504 | 81 friend std::ostream& operator << (std::ostream&, const ComplexQR&); |
| 457 | 82 |
| 539 | 83 protected: |
| 457 | 84 |
| 85 ComplexMatrix q; | |
| 86 ComplexMatrix r; | |
| 87 }; | |
| 88 | |
| 89 #endif | |
| 90 | |
| 91 /* | |
| 92 ;;; Local Variables: *** | |
| 93 ;;; mode: C++ *** | |
| 94 ;;; End: *** | |
| 95 */ |