Mercurial > octave
annotate liboctave/numeric/qrp.cc @ 30564:796f54d4ddbf stable
update Octave Project Developers copyright for the new year
In files that have the "Octave Project Developers" copyright notice,
update for 2021.
In all .txi and .texi files except gpl.txi and gpl.texi in the
doc/liboctave and doc/interpreter directories, change the copyright
to "Octave Project Developers", the same as used for other source
files. Update copyright notices for 2022 (not done since 2019). For
gpl.txi and gpl.texi, change the copyright notice to be "Free Software
Foundation, Inc." and leave the date at 2007 only because this file
only contains the text of the GPL, not anything created by the Octave
Project Developers.
Add Paul Thomas to contributors.in.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 28 Dec 2021 18:22:40 -0500 |
parents | abf2e77cca42 |
children | e88a07dec498 |
rev | line source |
---|---|
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
1 //////////////////////////////////////////////////////////////////////// |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
2 // |
30564
796f54d4ddbf
update Octave Project Developers copyright for the new year
John W. Eaton <jwe@octave.org>
parents:
30069
diff
changeset
|
3 // Copyright (C) 1994-2022 The Octave Project Developers |
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
4 // |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
5 // See the file COPYRIGHT.md in the top-level directory of this |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
6 // distribution or <https://octave.org/copyright/>. |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
7 // |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
8 // This file is part of Octave. |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
9 // |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
10 // Octave is free software: you can redistribute it and/or modify it |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
11 // under the terms of the GNU General Public License as published by |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
12 // the Free Software Foundation, either version 3 of the License, or |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
13 // (at your option) any later version. |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
14 // |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
15 // Octave is distributed in the hope that it will be useful, but |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
16 // WITHOUT ANY WARRANTY; without even the implied warranty of |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
18 // GNU General Public License for more details. |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
19 // |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
20 // You should have received a copy of the GNU General Public License |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
21 // along with Octave; see the file COPYING. If not, see |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
22 // <https://www.gnu.org/licenses/>. |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
23 // |
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
24 //////////////////////////////////////////////////////////////////////// |
538 | 25 |
21724
aba2e6293dd8
use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
26 #if defined (HAVE_CONFIG_H) |
21301
40de9f8f23a6
Use '#include "config.h"' rather than <config.h>.
Rik <rik@octave.org>
parents:
21280
diff
changeset
|
27 # include "config.h" |
538 | 28 #endif |
29 | |
1367 | 30 #include <cassert> |
538 | 31 |
23475
d691ed308237
maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents:
23450
diff
changeset
|
32 #include <algorithm> |
d691ed308237
maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents:
23450
diff
changeset
|
33 |
d691ed308237
maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents:
23450
diff
changeset
|
34 #include "Array.h" |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
35 #include "CMatrix.h" |
23475
d691ed308237
maint: Clean up #includes in liboctave/numeric directory.
Rik <rik@octave.org>
parents:
23450
diff
changeset
|
36 #include "MArray.h" |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
37 #include "dMatrix.h" |
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
38 #include "dRowVector.h" |
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
39 #include "fCMatrix.h" |
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
40 #include "fMatrix.h" |
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
41 #include "fRowVector.h" |
22322
93b3cdd36854
move most f77 function decls to separate header files
John W. Eaton <jwe@octave.org>
parents:
22317
diff
changeset
|
42 #include "lo-lapack-proto.h" |
8597
c86718093c1b
improve & fix QR classes
Jaroslav Hajek <highegg@gmail.com>
parents:
8375
diff
changeset
|
43 #include "oct-locbuf.h" |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
44 #include "qrp.h" |
538 | 45 |
22317
6ca3acf5fad8
move some new numeric classes to namespace octave::math
John W. Eaton <jwe@octave.org>
parents:
22135
diff
changeset
|
46 namespace octave |
6ca3acf5fad8
move some new numeric classes to namespace octave::math
John W. Eaton <jwe@octave.org>
parents:
22135
diff
changeset
|
47 { |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
48 namespace math |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
49 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
50 // Specialization. |
538 | 51 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
52 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
53 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
54 void |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
55 qrp<Matrix>::init (const Matrix& a, type qr_type) |
8597
c86718093c1b
improve & fix QR classes
Jaroslav Hajek <highegg@gmail.com>
parents:
8375
diff
changeset
|
56 { |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
57 assert (qr_type != qr<Matrix>::raw); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
58 |
23696
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
59 F77_INT m = to_f77_int (a.rows ()); |
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
60 F77_INT n = to_f77_int (a.cols ()); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
61 |
23450
855122b993da
maint: Wrap tertiary operator in parentheses "(COND ? x : y)".
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
62 F77_INT min_mn = (m < n ? m : n); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
63 OCTAVE_LOCAL_BUFFER (double, tau, min_mn); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
64 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
65 F77_INT info = 0; |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
66 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
67 Matrix afact = a; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
68 if (m > n && qr_type == qr<Matrix>::std) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
69 afact.resize (m, m); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
70 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
71 MArray<F77_INT> jpvt (dim_vector (n, 1), 0); |
8368
c72c1c9bccdc
call blocked permuted qr factorization routines from LAPACK
Jaroslav Hajek <highegg@gmail.com>
parents:
8367
diff
changeset
|
72 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
73 if (m > 0) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
74 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
75 // workspace query. |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
76 double rlwork; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
77 F77_XFCN (dgeqp3, DGEQP3, (m, n, afact.fortran_vec (), |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
78 m, jpvt.fortran_vec (), tau, |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
79 &rlwork, -1, info)); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
80 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
81 // allocate buffer and do the job. |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
82 F77_INT lwork = static_cast<F77_INT> (rlwork); |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
83 lwork = std::max (lwork, static_cast<F77_INT> (1)); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
84 OCTAVE_LOCAL_BUFFER (double, work, lwork); |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
85 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
86 F77_XFCN (dgeqp3, DGEQP3, (m, n, afact.fortran_vec (), |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
87 m, jpvt.fortran_vec (), tau, |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
88 work, lwork, info)); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
89 } |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
90 else |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
91 { |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
92 for (F77_INT i = 0; i < n; i++) |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
93 jpvt(i) = i+1; |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
94 } |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
95 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
96 // Form Permutation matrix (if economy is requested, return the |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
97 // indices only!) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
98 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
99 jpvt -= static_cast<F77_INT> (1); |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
100 m_p = PermMatrix (jpvt, true); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
101 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
102 form (n, afact, tau, qr_type); |
8597
c86718093c1b
improve & fix QR classes
Jaroslav Hajek <highegg@gmail.com>
parents:
8375
diff
changeset
|
103 } |
1922 | 104 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
105 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
106 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
107 qrp<Matrix>::qrp (const Matrix& a, type qr_type) |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
108 : qr<Matrix> (), m_p () |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
109 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
110 init (a, qr_type); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
111 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
112 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
113 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
114 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
115 RowVector |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
116 qrp<Matrix>::Pvec (void) const |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
117 { |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
118 Array<double> pa (m_p.col_perm_vec ()); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
119 RowVector pv (MArray<double> (pa) + 1.0); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
120 return pv; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
121 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
122 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
123 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
124 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
125 void |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
126 qrp<FloatMatrix>::init (const FloatMatrix& a, type qr_type) |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
127 { |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
128 assert (qr_type != qr<FloatMatrix>::raw); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
129 |
23696
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
130 F77_INT m = to_f77_int (a.rows ()); |
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
131 F77_INT n = to_f77_int (a.cols ()); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
132 |
23450
855122b993da
maint: Wrap tertiary operator in parentheses "(COND ? x : y)".
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
133 F77_INT min_mn = (m < n ? m : n); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
134 OCTAVE_LOCAL_BUFFER (float, tau, min_mn); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
135 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
136 F77_INT info = 0; |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
137 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
138 FloatMatrix afact = a; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
139 if (m > n && qr_type == qr<FloatMatrix>::std) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
140 afact.resize (m, m); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
141 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
142 MArray<F77_INT> jpvt (dim_vector (n, 1), 0); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
143 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
144 if (m > 0) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
145 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
146 // workspace query. |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
147 float rlwork; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
148 F77_XFCN (sgeqp3, SGEQP3, (m, n, afact.fortran_vec (), |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
149 m, jpvt.fortran_vec (), tau, |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
150 &rlwork, -1, info)); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
151 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
152 // allocate buffer and do the job. |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
153 F77_INT lwork = static_cast<F77_INT> (rlwork); |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
154 lwork = std::max (lwork, static_cast<F77_INT> (1)); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
155 OCTAVE_LOCAL_BUFFER (float, work, lwork); |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
156 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
157 F77_XFCN (sgeqp3, SGEQP3, (m, n, afact.fortran_vec (), |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
158 m, jpvt.fortran_vec (), tau, |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
159 work, lwork, info)); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
160 } |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
161 else |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
162 { |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
163 for (F77_INT i = 0; i < n; i++) |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
164 jpvt(i) = i+1; |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
165 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
166 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
167 // Form Permutation matrix (if economy is requested, return the |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
168 // indices only!) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
169 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
170 jpvt -= static_cast<F77_INT> (1); |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
171 m_p = PermMatrix (jpvt, true); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
172 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
173 form (n, afact, tau, qr_type); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
174 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
175 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
176 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
177 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
178 qrp<FloatMatrix>::qrp (const FloatMatrix& a, type qr_type) |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
179 : qr<FloatMatrix> (), m_p () |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
180 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
181 init (a, qr_type); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
182 } |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
183 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
184 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
185 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
186 FloatRowVector |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
187 qrp<FloatMatrix>::Pvec (void) const |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
188 { |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
189 Array<float> pa (m_p.col_perm_vec ()); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
190 FloatRowVector pv (MArray<float> (pa) + 1.0f); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
191 return pv; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
192 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
193 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
194 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
195 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
196 void |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
197 qrp<ComplexMatrix>::init (const ComplexMatrix& a, type qr_type) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
198 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
199 assert (qr_type != qr<ComplexMatrix>::raw); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
200 |
23696
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
201 F77_INT m = to_f77_int (a.rows ()); |
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
202 F77_INT n = to_f77_int (a.cols ()); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
203 |
23450
855122b993da
maint: Wrap tertiary operator in parentheses "(COND ? x : y)".
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
204 F77_INT min_mn = (m < n ? m : n); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
205 OCTAVE_LOCAL_BUFFER (Complex, tau, min_mn); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
206 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
207 F77_INT info = 0; |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
208 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
209 ComplexMatrix afact = a; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
210 if (m > n && qr_type == qr<ComplexMatrix>::std) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
211 afact.resize (m, m); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
212 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
213 MArray<F77_INT> jpvt (dim_vector (n, 1), 0); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
214 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
215 if (m > 0) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
216 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
217 OCTAVE_LOCAL_BUFFER (double, rwork, 2*n); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
218 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
219 // workspace query. |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
220 Complex clwork; |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
221 F77_XFCN (zgeqp3, ZGEQP3, (m, n, |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
222 F77_DBLE_CMPLX_ARG (afact.fortran_vec ()), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
223 m, jpvt.fortran_vec (), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
224 F77_DBLE_CMPLX_ARG (tau), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
225 F77_DBLE_CMPLX_ARG (&clwork), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
226 -1, rwork, info)); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
227 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
228 // allocate buffer and do the job. |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
229 F77_INT lwork = static_cast<F77_INT> (clwork.real ()); |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
230 lwork = std::max (lwork, static_cast<F77_INT> (1)); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
231 OCTAVE_LOCAL_BUFFER (Complex, work, lwork); |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
232 |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
233 F77_XFCN (zgeqp3, ZGEQP3, (m, n, |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
234 F77_DBLE_CMPLX_ARG (afact.fortran_vec ()), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
235 m, jpvt.fortran_vec (), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
236 F77_DBLE_CMPLX_ARG (tau), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
237 F77_DBLE_CMPLX_ARG (work), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
238 lwork, rwork, info)); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
239 } |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
240 else |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
241 { |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
242 for (F77_INT i = 0; i < n; i++) |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
243 jpvt(i) = i+1; |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
244 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
245 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
246 // Form Permutation matrix (if economy is requested, return the |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
247 // indices only!) |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
248 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
249 jpvt -= static_cast<F77_INT> (1); |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
250 m_p = PermMatrix (jpvt, true); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
251 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
252 form (n, afact, tau, qr_type); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
253 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
254 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
255 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
256 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
257 qrp<ComplexMatrix>::qrp (const ComplexMatrix& a, type qr_type) |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
258 : qr<ComplexMatrix> (), m_p () |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
259 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
260 init (a, qr_type); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
261 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
262 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
263 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
264 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
265 RowVector |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
266 qrp<ComplexMatrix>::Pvec (void) const |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
267 { |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
268 Array<double> pa (m_p.col_perm_vec ()); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
269 RowVector pv (MArray<double> (pa) + 1.0); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
270 return pv; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
271 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
272 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
273 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
274 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
275 void |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
276 qrp<FloatComplexMatrix>::init (const FloatComplexMatrix& a, type qr_type) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
277 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
278 assert (qr_type != qr<FloatComplexMatrix>::raw); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
279 |
23696
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
280 F77_INT m = to_f77_int (a.rows ()); |
08036a7f3660
remove octave:: namespace tag from symbols used inside octave namespace
John W. Eaton <jwe@octave.org>
parents:
23475
diff
changeset
|
281 F77_INT n = to_f77_int (a.cols ()); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
282 |
23450
855122b993da
maint: Wrap tertiary operator in parentheses "(COND ? x : y)".
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
283 F77_INT min_mn = (m < n ? m : n); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
284 OCTAVE_LOCAL_BUFFER (FloatComplex, tau, min_mn); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
285 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
286 F77_INT info = 0; |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
287 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
288 FloatComplexMatrix afact = a; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
289 if (m > n && qr_type == qr<FloatComplexMatrix>::std) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
290 afact.resize (m, m); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
291 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
292 MArray<F77_INT> jpvt (dim_vector (n, 1), 0); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
293 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
294 if (m > 0) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
295 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
296 OCTAVE_LOCAL_BUFFER (float, rwork, 2*n); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
297 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
298 // workspace query. |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
299 FloatComplex clwork; |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
300 F77_XFCN (cgeqp3, CGEQP3, (m, n, |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
301 F77_CMPLX_ARG (afact.fortran_vec ()), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
302 m, jpvt.fortran_vec (), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
303 F77_CMPLX_ARG (tau), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
304 F77_CMPLX_ARG (&clwork), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
305 -1, rwork, info)); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
306 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
307 // allocate buffer and do the job. |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
308 F77_INT lwork = static_cast<F77_INT> (clwork.real ()); |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
309 lwork = std::max (lwork, static_cast<F77_INT> (1)); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
310 OCTAVE_LOCAL_BUFFER (FloatComplex, work, lwork); |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
311 |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
312 F77_XFCN (cgeqp3, CGEQP3, (m, n, |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
313 F77_CMPLX_ARG (afact.fortran_vec ()), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
314 m, jpvt.fortran_vec (), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
315 F77_CMPLX_ARG (tau), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
316 F77_CMPLX_ARG (work), |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
317 lwork, rwork, info)); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
318 } |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
319 else |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
320 { |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
321 for (F77_INT i = 0; i < n; i++) |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
322 jpvt(i) = i+1; |
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
323 } |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
324 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
325 // Form Permutation matrix (if economy is requested, return the |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
326 // indices only!) |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
327 |
22951
2e6a327e7903
use F77_INT instead of octave_idx_type for liboctave QRP functions
John W. Eaton <jwe@octave.org>
parents:
22755
diff
changeset
|
328 jpvt -= static_cast<F77_INT> (1); |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
329 m_p = PermMatrix (jpvt, true); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
330 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
331 form (n, afact, tau, qr_type); |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
332 } |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
333 |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
334 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
335 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
336 qrp<FloatComplexMatrix>::qrp (const FloatComplexMatrix& a, type qr_type) |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
337 : qr<FloatComplexMatrix> (), m_p () |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
338 { |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
339 init (a, qr_type); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
340 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
341 |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
342 template <> |
29228
5c14f81e0937
Set API tags in files in liboctave/numeric (patch #8919).
Markus Mützel <markus.muetzel@gmx.de>
parents:
27923
diff
changeset
|
343 OCTAVE_API |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
344 FloatRowVector |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
345 qrp<FloatComplexMatrix>::Pvec (void) const |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
346 { |
30069
abf2e77cca42
maint: use "m_" prefix for member variables in class qrp.
Rik <rik@octave.org>
parents:
29359
diff
changeset
|
347 Array<float> pa (m_p.col_perm_vec ()); |
22329
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
348 FloatRowVector pv (MArray<float> (pa) + 1.0f); |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
349 return pv; |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
350 } |
7f3c7a8bd131
maint: Indent namespaces in liboctave/numeric files.
John W. Eaton <jwe@octave.org>
parents:
22323
diff
changeset
|
351 } |
21280
ebdf74c15722
better use of templates for qrp classes
John W. Eaton <jwe@octave.org>
parents:
21279
diff
changeset
|
352 } |