Mercurial > hg > octave-nkf
annotate liboctave/numeric/CmplxSVD.h @ 15436:32fd31378052
Additional changes to repmat.m (e1f59fd57756).
Bug # 37390.
* scripts/general/repmat.m: Changeset e1f59fd57756 did not handle empties
in all situations correctly. Also, ML allows the inputs to be matrices.
| author | Ben Abbott <bpabbott@mac.com> |
|---|---|
| date | Sat, 22 Sep 2012 16:48:19 -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 |
| 457 | 4 |
| 5 This file is part of Octave. | |
| 6 | |
| 7 Octave is free software; you can redistribute it and/or modify it | |
| 8 under the terms of the GNU General Public License as published by the | |
| 7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
| 10 option) any later version. | |
| 457 | 11 |
| 12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 15 for more details. | |
| 16 | |
| 17 You should have received a copy of the GNU General Public License | |
| 7016 | 18 along with Octave; see the file COPYING. If not, see |
| 19 <http://www.gnu.org/licenses/>. | |
| 457 | 20 |
| 21 */ | |
| 22 | |
| 23 #if !defined (octave_ComplexSVD_h) | |
| 24 #define octave_ComplexSVD_h 1 | |
| 25 | |
|
8950
d865363208d6
include <iosfwd> instead of <iostream> in header files
John W. Eaton <jwe@octave.org>
parents:
7017
diff
changeset
|
26 #include <iosfwd> |
| 457 | 27 |
| 28 #include "dDiagMatrix.h" | |
| 29 #include "CMatrix.h" | |
| 537 | 30 #include "dbleSVD.h" |
| 457 | 31 |
| 1881 | 32 class |
| 6108 | 33 OCTAVE_API |
| 1881 | 34 ComplexSVD |
| 457 | 35 { |
| 36 public: | |
| 37 | |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
38 ComplexSVD (void) |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
39 : type_computed (), sigma (), left_sm (), right_sm () |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
40 { } |
| 1528 | 41 |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
42 ComplexSVD (const ComplexMatrix& a, SVD::type svd_type = SVD::std, |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
43 SVD::driver svd_driver = SVD::GESVD) |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
44 : type_computed (), sigma (), left_sm (), right_sm () |
| 1528 | 45 { |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
46 init (a, svd_type, svd_driver); |
| 1528 | 47 } |
| 48 | |
| 5275 | 49 ComplexSVD (const ComplexMatrix& a, octave_idx_type& info, |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
50 SVD::type svd_type = SVD::std, |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
51 SVD::driver svd_driver = SVD::GESVD) |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
52 : type_computed (), sigma (), left_sm (), right_sm () |
| 1528 | 53 { |
| 10601 | 54 info = init (a, svd_type, svd_driver); |
| 1528 | 55 } |
| 457 | 56 |
| 1528 | 57 ComplexSVD (const ComplexSVD& a) |
|
11498
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
58 : type_computed (a.type_computed), sigma (a.sigma), |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
59 left_sm (a.left_sm), right_sm (a.right_sm) |
|
367bfee35ba0
data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
10601
diff
changeset
|
60 { } |
| 457 | 61 |
| 1528 | 62 ComplexSVD& operator = (const ComplexSVD& a) |
| 63 { | |
| 1881 | 64 if (this != &a) |
|
10312
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
65 { |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
66 type_computed = a.type_computed; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
67 sigma = a.sigma; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
68 left_sm = a.left_sm; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
69 right_sm = a.right_sm; |
|
cbc402e64d83
untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents:
10158
diff
changeset
|
70 } |
| 1528 | 71 return *this; |
| 72 } | |
| 73 | |
| 1930 | 74 ~ComplexSVD (void) { } |
| 75 | |
| 1528 | 76 DiagMatrix singular_values (void) const { return sigma; } |
| 77 | |
| 1544 | 78 ComplexMatrix left_singular_matrix (void) const; |
| 1528 | 79 |
| 1544 | 80 ComplexMatrix right_singular_matrix (void) const; |
| 457 | 81 |
| 3504 | 82 friend std::ostream& operator << (std::ostream& os, const ComplexSVD& a); |
| 457 | 83 |
| 84 private: | |
| 85 | |
| 1544 | 86 SVD::type type_computed; |
| 87 | |
| 457 | 88 DiagMatrix sigma; |
| 89 ComplexMatrix left_sm; | |
| 90 ComplexMatrix right_sm; | |
| 1881 | 91 |
| 10601 | 92 octave_idx_type init (const ComplexMatrix& a, |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
93 SVD::type svd_type = SVD::std, |
| 10601 | 94 SVD::driver svd_driver = SVD::GESVD); |
| 457 | 95 }; |
| 96 | |
| 97 #endif |