Mercurial > octave-nkf
annotate liboctave/dbleSVD.h @ 10843:229675bb7647 ss-3-3-52
version is now 3.3.52
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Sun, 01 Aug 2010 11:49:45 -0400 |
| parents | 3ce0c530a9c9 |
| children | 367bfee35ba0 |
| rev | line source |
|---|---|
| 457 | 1 /* |
| 2 | |
| 7017 | 3 Copyright (C) 1994, 1995, 1996, 1997, 2000, 2002, 2003, 2004, 2005, |
| 9245 | 4 2006, 2007, 2009 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 | |
| 24 #if !defined (octave_SVD_h) | |
| 25 #define octave_SVD_h 1 | |
| 26 | |
|
8950
d865363208d6
include <iosfwd> instead of <iostream> in header files
John W. Eaton <jwe@octave.org>
parents:
7017
diff
changeset
|
27 #include <iosfwd> |
| 457 | 28 |
| 29 #include "dDiagMatrix.h" | |
| 30 #include "dMatrix.h" | |
| 31 | |
| 1880 | 32 class |
| 6108 | 33 OCTAVE_API |
| 1880 | 34 SVD |
| 457 | 35 { |
| 36 public: | |
| 37 | |
| 537 | 38 enum type |
| 39 { | |
| 40 std, | |
| 41 economy, | |
| 2802 | 42 sigma_only |
| 537 | 43 }; |
| 44 | |
| 10601 | 45 enum driver |
| 46 { | |
| 47 GESVD, | |
| 48 GESDD | |
| 49 }; | |
| 50 | |
| 1881 | 51 SVD (void) : sigma (), left_sm (), right_sm () { } |
| 1528 | 52 |
| 10601 | 53 SVD (const Matrix& a, |
| 54 type svd_type = SVD::std, driver svd_driver = SVD::GESVD) | |
| 55 { init (a, svd_type, svd_driver); } | |
| 457 | 56 |
| 10601 | 57 SVD (const Matrix& a, octave_idx_type& info, |
| 58 type svd_type = SVD::std, driver svd_driver = SVD::GESVD) | |
| 1528 | 59 { |
| 10601 | 60 info = init (a, svd_type, svd_driver); |
| 1528 | 61 } |
| 457 | 62 |
| 1528 | 63 SVD (const SVD& a) |
| 4374 | 64 : type_computed (a.type_computed), |
| 65 sigma (a.sigma), left_sm (a.left_sm), right_sm (a.right_sm) { } | |
| 457 | 66 |
| 1528 | 67 SVD& operator = (const SVD& a) |
| 68 { | |
| 1880 | 69 if (this != &a) |
|
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
70 { |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
71 type_computed = a.type_computed; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
72 sigma = a.sigma; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
73 left_sm = a.left_sm; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
74 right_sm = a.right_sm; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
75 } |
| 1528 | 76 |
| 77 return *this; | |
| 78 } | |
| 79 | |
| 1930 | 80 ~SVD (void) { } |
| 81 | |
| 1528 | 82 DiagMatrix singular_values (void) const { return sigma; } |
| 83 | |
| 1544 | 84 Matrix left_singular_matrix (void) const; |
| 1528 | 85 |
| 1544 | 86 Matrix right_singular_matrix (void) const; |
| 457 | 87 |
| 3504 | 88 friend std::ostream& operator << (std::ostream& os, const SVD& a); |
| 457 | 89 |
| 90 private: | |
| 91 | |
| 1544 | 92 SVD::type type_computed; |
| 93 | |
| 457 | 94 DiagMatrix sigma; |
| 95 Matrix left_sm; | |
| 96 Matrix right_sm; | |
| 1880 | 97 |
| 10601 | 98 octave_idx_type init (const Matrix& a, |
| 99 type svd_type = std, driver svd_driver = GESVD); | |
| 457 | 100 }; |
| 101 | |
| 102 #endif |