Mercurial > octave-nkf
annotate liboctave/CmplxCHOL.h @ 14164:a737b3fb9c4d stable ss-3-5-92
snapshot 3.5.92
* configure.ac (AC_INIT): Version is now 3.5.92.
(OCTAVE_RELEASE_DATE): Now 2012-01-06.
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Fri, 06 Jan 2012 14:34:06 -0500 |
| parents | 72c96de7a403 |
| children |
| rev | line source |
|---|---|
| 457 | 1 /* |
| 2 | |
|
14138
72c96de7a403
maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
3 Copyright (C) 1994-2012 John W. Eaton |
| 11523 | 4 Copyright (C) 2008-2009 Jaroslav Hajek |
| 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 | |
| 24 #if !defined (octave_ComplexCHOL_h) | |
| 25 #define octave_ComplexCHOL_h 1 | |
| 26 | |
|
8950
d865363208d6
include <iosfwd> instead of <iostream> in header files
John W. Eaton <jwe@octave.org>
parents:
8562
diff
changeset
|
27 #include <iosfwd> |
| 457 | 28 |
| 29 #include "CMatrix.h" | |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
30 #include "CColVector.h" |
| 457 | 31 |
| 1881 | 32 class |
| 6108 | 33 OCTAVE_API |
| 1881 | 34 ComplexCHOL |
| 457 | 35 { |
| 36 public: | |
| 37 | |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
38 ComplexCHOL (void) : chol_mat (), xrcond (0) { } |
| 1528 | 39 |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
40 ComplexCHOL (const ComplexMatrix& a, bool calc_cond = false) |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
41 : chol_mat (), xrcond (0) |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
42 { |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
43 init (a, calc_cond); |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
44 } |
| 1528 | 45 |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
46 ComplexCHOL (const ComplexMatrix& a, octave_idx_type& info, |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
47 bool calc_cond = false) |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
48 : chol_mat (), xrcond (0) |
| 1528 | 49 { |
| 6486 | 50 info = init (a, calc_cond); |
| 1528 | 51 } |
| 52 | |
| 4055 | 53 ComplexCHOL (const ComplexCHOL& a) |
| 6486 | 54 : chol_mat (a.chol_mat), xrcond (a.xrcond) { } |
| 1528 | 55 |
| 4055 | 56 ComplexCHOL& operator = (const ComplexCHOL& a) |
| 1528 | 57 { |
| 1881 | 58 if (this != &a) |
|
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
59 { |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
60 chol_mat = a.chol_mat; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
61 xrcond = a.xrcond; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
62 } |
| 1528 | 63 |
| 64 return *this; | |
| 65 } | |
| 66 | |
| 5340 | 67 ComplexMatrix chol_matrix (void) const { return chol_mat; } |
| 68 | |
| 6486 | 69 double rcond (void) const { return xrcond; } |
| 70 | |
| 5340 | 71 ComplexMatrix inverse (void) const; |
| 457 | 72 |
|
7554
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
73 void set (const ComplexMatrix& R); |
|
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
74 |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
75 void update (const ComplexColumnVector& u); |
|
7554
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
76 |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
77 octave_idx_type downdate (const ComplexColumnVector& u); |
|
7554
40574114c514
implement Cholesky 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 octave_idx_type insert_sym (const ComplexColumnVector& u, octave_idx_type j); |
|
7700
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
80 |
|
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
81 void delete_sym (octave_idx_type j); |
|
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
82 |
|
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
83 void shift_sym (octave_idx_type i, octave_idx_type j); |
|
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
84 |
| 6108 | 85 friend OCTAVE_API std::ostream& operator << (std::ostream& os, const ComplexCHOL& a); |
| 457 | 86 |
| 87 private: | |
| 88 | |
| 1881 | 89 ComplexMatrix chol_mat; |
| 457 | 90 |
| 6486 | 91 double xrcond; |
| 92 | |
| 93 octave_idx_type init (const ComplexMatrix& a, bool calc_cond); | |
| 457 | 94 }; |
| 95 | |
| 6108 | 96 ComplexMatrix OCTAVE_API chol2inv (const ComplexMatrix& r); |
| 5340 | 97 |
| 457 | 98 #endif |