Mercurial > hg > octave-nkf
annotate liboctave/numeric/dbleCHOL.h @ 16583:e74ef19d2268
use octave_idx_type instead of int
* ls-mat-ascii.cc (looks_like_mat_ascii_file): USe octave_idx_type
instead of int for temporary variables.
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Sun, 28 Apr 2013 23:21:43 -0400 |
| parents | 648dabbb4c6b |
| children | d63878346099 |
| rev | line source |
|---|---|
| 457 | 1 /* |
| 2 | |
|
14138
72c96de7a403
maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents:
11586
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_CHOL_h) | |
| 25 #define octave_CHOL_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 "dMatrix.h" | |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
30 #include "dColVector.h" |
| 457 | 31 |
| 1881 | 32 class |
| 6108 | 33 OCTAVE_API |
| 1881 | 34 CHOL |
| 457 | 35 { |
| 36 public: | |
| 37 | |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
38 CHOL (void) : chol_mat (), xrcond (0) { } |
| 1528 | 39 |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
40 CHOL (const Matrix& 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 } |
| 457 | 45 |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
46 CHOL (const Matrix& a, octave_idx_type& info, bool calc_cond = false) |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
47 : chol_mat (), xrcond (0) |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
48 { |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
49 info = init (a, calc_cond); |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10312
diff
changeset
|
50 } |
| 1528 | 51 |
| 6486 | 52 CHOL (const CHOL& a) : chol_mat (a.chol_mat), xrcond (a.xrcond) { } |
| 457 | 53 |
| 1528 | 54 CHOL& operator = (const CHOL& a) |
| 55 { | |
| 1881 | 56 if (this != &a) |
|
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
57 { |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
58 chol_mat = a.chol_mat; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
59 xrcond = a.xrcond; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
60 } |
| 1528 | 61 return *this; |
| 62 } | |
| 63 | |
| 64 Matrix chol_matrix (void) const { return chol_mat; } | |
| 65 | |
| 6486 | 66 double rcond (void) const { return xrcond; } |
| 67 | |
| 5340 | 68 // Compute the inverse of a matrix using the Cholesky factorization. |
| 69 Matrix inverse (void) const; | |
| 70 | |
|
7554
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
71 void set (const Matrix& R); |
|
40574114c514
implement Cholesky factorization updating
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
72 |
|
8547
d66c9b6e506a
imported patch qrupdate.diff
Jaroslav Hajek <highegg@gmail.com>
parents:
7700
diff
changeset
|
73 void update (const ColumnVector& u); |
|
7554
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 octave_idx_type downdate (const ColumnVector& 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 insert_sym (const ColumnVector& u, octave_idx_type j); |
|
7700
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
78 |
|
efccca5f2ad7
more QR & Cholesky updating functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7554
diff
changeset
|
79 void delete_sym (octave_idx_type j); |
|
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 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
|
82 |
| 6108 | 83 friend OCTAVE_API std::ostream& operator << (std::ostream& os, const CHOL& a); |
| 457 | 84 |
| 85 private: | |
| 86 | |
| 1881 | 87 Matrix chol_mat; |
| 457 | 88 |
| 6486 | 89 double xrcond; |
| 90 | |
| 91 octave_idx_type init (const Matrix& a, bool calc_cond); | |
| 457 | 92 }; |
| 93 | |
| 6108 | 94 Matrix OCTAVE_API chol2inv (const Matrix& r); |
| 5340 | 95 |
| 457 | 96 #endif |