Mercurial > hg > octave-nkf
annotate liboctave/dSparse.h @ 9790:a5035bc7fbfb
rewrite dispatch part & slightly improve min,max,cummin,cummax
| author | Jaroslav Hajek <highegg@gmail.com> |
|---|---|
| date | Mon, 09 Nov 2009 13:09:49 +0100 |
| parents | 91d53dc37f79 |
| children | f80c566bc751 |
| rev | line source |
|---|---|
| 5164 | 1 /* |
| 2 | |
| 8920 | 3 Copyright (C) 2004, 2005, 2006, 2007, 2008, 2009 David Bateman |
| 7016 | 4 Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004 Andy Adler |
| 5 | |
| 6 This file is part of Octave. | |
| 5164 | 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. | |
| 5164 | 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/>. | |
| 5164 | 21 |
| 22 */ | |
| 23 | |
| 24 #if !defined (octave_dSparse_h) | |
| 25 #define octave_dSparse_h 1 | |
| 26 | |
| 27 #include "dMatrix.h" | |
| 28 #include "dNDArray.h" | |
| 29 #include "CMatrix.h" | |
| 30 #include "dColVector.h" | |
| 31 #include "CColVector.h" | |
| 32 | |
| 8335 | 33 #include "DET.h" |
| 5164 | 34 #include "MSparse.h" |
| 35 #include "MSparse-defs.h" | |
| 36 #include "Sparse-op-defs.h" | |
| 5785 | 37 #include "MatrixType.h" |
| 5164 | 38 |
|
8968
91d53dc37f79
Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations.
Jason Riedy <jason@acm.org>
parents:
8966
diff
changeset
|
39 class PermMatrix; |
|
8964
f4f4d65faaa0
Implement sparse * diagonal and diagonal * sparse operations, double-prec only.
Jason Riedy <jason@acm.org>
parents:
8920
diff
changeset
|
40 class DiagMatrix; |
| 5164 | 41 class SparseComplexMatrix; |
| 42 class SparseBoolMatrix; | |
| 43 | |
| 44 class | |
| 6108 | 45 OCTAVE_API |
| 5164 | 46 SparseMatrix : public MSparse<double> |
| 47 { | |
| 48 public: | |
| 49 | |
| 50 typedef void (*solve_singularity_handler) (double rcond); | |
| 51 | |
| 52 SparseMatrix (void) : MSparse<double> () { } | |
| 53 | |
| 5275 | 54 SparseMatrix (octave_idx_type r, octave_idx_type c) : MSparse<double> (r, c) { } |
| 5164 | 55 |
| 6823 | 56 SparseMatrix (const dim_vector& dv, octave_idx_type nz = 0) : |
| 57 MSparse<double> (dv, nz) { } | |
| 58 | |
| 5275 | 59 explicit SparseMatrix (octave_idx_type r, octave_idx_type c, double val) |
| 5164 | 60 : MSparse<double> (r, c, val) { } |
| 61 | |
| 62 SparseMatrix (const SparseMatrix& a) : MSparse<double> (a) { } | |
| 63 | |
| 64 SparseMatrix (const SparseMatrix& a, const dim_vector& dv) | |
| 65 : MSparse<double> (a, dv) { } | |
| 66 | |
| 67 SparseMatrix (const MSparse<double>& a) : MSparse<double> (a) { } | |
| 68 | |
| 69 explicit SparseMatrix (const SparseBoolMatrix& a); | |
| 70 | |
| 71 explicit SparseMatrix (const Matrix& a) : MSparse<double> (a) { } | |
| 72 | |
| 73 explicit SparseMatrix (const NDArray& a) : MSparse<double> (a) { } | |
| 74 | |
| 5275 | 75 explicit SparseMatrix (const Array<double> a, const Array<octave_idx_type>& r, |
| 76 const Array<octave_idx_type>& c, octave_idx_type nr = -1, | |
| 77 octave_idx_type nc = -1, bool sum_terms = true) | |
| 5164 | 78 : MSparse<double> (a, r, c, nr, nc, sum_terms) { } |
| 79 | |
| 80 explicit SparseMatrix (const Array<double> a, const Array<double>& r, | |
| 5275 | 81 const Array<double>& c, octave_idx_type nr = -1, |
| 82 octave_idx_type nc = -1, bool sum_terms = true) | |
| 5164 | 83 : MSparse<double> (a, r, c, nr, nc, sum_terms) { } |
| 84 | |
|
8366
8b1a2555c4e2
implement diagonal matrix objects
Jaroslav Hajek <highegg@gmail.com>
parents:
8335
diff
changeset
|
85 explicit SparseMatrix (const DiagMatrix& a); |
|
8b1a2555c4e2
implement diagonal matrix objects
Jaroslav Hajek <highegg@gmail.com>
parents:
8335
diff
changeset
|
86 |
|
8910
6e9f26506804
optimize diag -> sparse and perm -> sparse conversions
Jaroslav Hajek <highegg@gmail.com>
parents:
8366
diff
changeset
|
87 explicit SparseMatrix (const PermMatrix& a); |
|
6e9f26506804
optimize diag -> sparse and perm -> sparse conversions
Jaroslav Hajek <highegg@gmail.com>
parents:
8366
diff
changeset
|
88 |
| 5275 | 89 SparseMatrix (octave_idx_type r, octave_idx_type c, octave_idx_type num_nz) : MSparse<double> (r, c, num_nz) { } |
| 5164 | 90 |
| 91 SparseMatrix& operator = (const SparseMatrix& a) | |
| 92 { | |
| 93 MSparse<double>::operator = (a); | |
| 94 return *this; | |
| 95 } | |
| 96 | |
| 97 bool operator == (const SparseMatrix& a) const; | |
| 98 bool operator != (const SparseMatrix& a) const; | |
| 99 | |
| 100 bool is_symmetric (void) const; | |
| 101 | |
|
9790
a5035bc7fbfb
rewrite dispatch part & slightly improve min,max,cummin,cummax
Jaroslav Hajek <highegg@gmail.com>
parents:
8968
diff
changeset
|
102 SparseMatrix max (int dim = -1) const; |
|
a5035bc7fbfb
rewrite dispatch part & slightly improve min,max,cummin,cummax
Jaroslav Hajek <highegg@gmail.com>
parents:
8968
diff
changeset
|
103 SparseMatrix max (Array<octave_idx_type>& index, int dim = -1) const; |
|
a5035bc7fbfb
rewrite dispatch part & slightly improve min,max,cummin,cummax
Jaroslav Hajek <highegg@gmail.com>
parents:
8968
diff
changeset
|
104 SparseMatrix min (int dim = -1) const; |
|
a5035bc7fbfb
rewrite dispatch part & slightly improve min,max,cummin,cummax
Jaroslav Hajek <highegg@gmail.com>
parents:
8968
diff
changeset
|
105 SparseMatrix min (Array<octave_idx_type>& index, int dim = -1) const; |
| 5164 | 106 |
| 107 // destructive insert/delete/reorder operations | |
| 108 | |
| 5275 | 109 SparseMatrix& insert (const SparseMatrix& a, octave_idx_type r, octave_idx_type c); |
| 5164 | 110 |
| 6823 | 111 SparseMatrix& insert (const SparseMatrix& a, const Array<octave_idx_type>& indx); |
| 112 | |
| 5275 | 113 SparseMatrix concat (const SparseMatrix& rb, const Array<octave_idx_type>& ra_idx); |
| 5164 | 114 SparseComplexMatrix concat (const SparseComplexMatrix& rb, |
| 5275 | 115 const Array<octave_idx_type>& ra_idx); |
| 5164 | 116 |
| 6108 | 117 friend OCTAVE_API SparseMatrix real (const SparseComplexMatrix& a); |
| 118 friend OCTAVE_API SparseMatrix imag (const SparseComplexMatrix& a); | |
| 5164 | 119 |
| 6108 | 120 friend OCTAVE_API SparseMatrix atan2 (const double& x, const SparseMatrix& y); |
| 121 friend OCTAVE_API SparseMatrix atan2 (const SparseMatrix& x, const double& y); | |
| 122 friend OCTAVE_API SparseMatrix atan2 (const SparseMatrix& x, const SparseMatrix& y); | |
| 5164 | 123 |
| 124 SparseMatrix transpose (void) const | |
| 125 { | |
| 126 return MSparse<double>::transpose (); | |
| 127 } | |
| 5506 | 128 SparseMatrix hermitian (void) const { return transpose (); } |
| 5164 | 129 |
|
8303
b11c31849b44
improve norm computation capabilities
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
130 // extract row or column i. |
|
b11c31849b44
improve norm computation capabilities
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
131 |
|
b11c31849b44
improve norm computation capabilities
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
132 RowVector row (octave_idx_type i) const; |
|
b11c31849b44
improve norm computation capabilities
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
133 |
|
b11c31849b44
improve norm computation capabilities
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
134 ColumnVector column (octave_idx_type i) const; |
|
b11c31849b44
improve norm computation capabilities
Jaroslav Hajek <highegg@gmail.com>
parents:
7922
diff
changeset
|
135 |
| 5506 | 136 private: |
| 5785 | 137 SparseMatrix dinverse (MatrixType &mattyp, octave_idx_type& info, |
| 5506 | 138 double& rcond, const bool force = false, |
| 139 const bool calccond = true) const; | |
| 140 | |
| 5785 | 141 SparseMatrix tinverse (MatrixType &mattyp, octave_idx_type& info, |
| 5506 | 142 double& rcond, const bool force = false, |
| 143 const bool calccond = true) const; | |
| 144 | |
| 145 public: | |
| 5164 | 146 SparseMatrix inverse (void) const; |
| 5785 | 147 SparseMatrix inverse (MatrixType& mattype) const; |
| 148 SparseMatrix inverse (MatrixType& mattype, octave_idx_type& info) const; | |
| 149 SparseMatrix inverse (MatrixType& mattype, octave_idx_type& info, | |
| 5506 | 150 double& rcond, int force = 0, int calc_cond = 1) const; |
| 5164 | 151 |
| 152 DET determinant (void) const; | |
| 5275 | 153 DET determinant (octave_idx_type& info) const; |
| 154 DET determinant (octave_idx_type& info, double& rcond, int calc_cond = 1) const; | |
| 5164 | 155 |
| 156 private: | |
| 157 // Diagonal matrix solvers | |
| 5785 | 158 Matrix dsolve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5681 | 159 double& rcond, solve_singularity_handler sing_handler, |
| 160 bool calc_cond = false) const; | |
| 5164 | 161 |
| 5785 | 162 ComplexMatrix dsolve (MatrixType &typ, const ComplexMatrix& b, |
| 5681 | 163 octave_idx_type& info, double& rcond, |
| 164 solve_singularity_handler sing_handler, | |
| 165 bool calc_cond = false) const; | |
| 5164 | 166 |
| 5785 | 167 SparseMatrix dsolve (MatrixType &typ, const SparseMatrix& b, |
| 5681 | 168 octave_idx_type& info, double& rcond, |
| 169 solve_singularity_handler sing_handler, | |
| 170 bool calc_cond = false) const; | |
| 5164 | 171 |
| 5785 | 172 SparseComplexMatrix dsolve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 173 octave_idx_type& info, double& rcond, |
| 5681 | 174 solve_singularity_handler sing_handler, |
| 175 bool calc_cond = false) const; | |
| 5164 | 176 |
| 177 // Upper triangular matrix solvers | |
| 5785 | 178 Matrix utsolve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5681 | 179 double& rcond, solve_singularity_handler sing_handler, |
| 180 bool calc_cond = false) const; | |
| 5164 | 181 |
| 5785 | 182 ComplexMatrix utsolve (MatrixType &typ, const ComplexMatrix& b, |
| 5681 | 183 octave_idx_type& info, double& rcond, |
| 184 solve_singularity_handler sing_handler, | |
| 185 bool calc_cond = false) const; | |
| 5164 | 186 |
| 5785 | 187 SparseMatrix utsolve (MatrixType &typ, const SparseMatrix& b, |
| 5681 | 188 octave_idx_type& info, double& rcond, |
| 189 solve_singularity_handler sing_handler, | |
| 190 bool calc_cond = false) const; | |
| 5164 | 191 |
| 5785 | 192 SparseComplexMatrix utsolve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 193 octave_idx_type& info, double& rcond, |
| 5681 | 194 solve_singularity_handler sing_handler, |
| 195 bool calc_cond = false) const; | |
| 5164 | 196 |
| 197 // Lower triangular matrix solvers | |
| 5785 | 198 Matrix ltsolve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5681 | 199 double& rcond, solve_singularity_handler sing_handler, |
| 200 bool calc_cond = false) const; | |
| 5164 | 201 |
| 5785 | 202 ComplexMatrix ltsolve (MatrixType &typ, const ComplexMatrix& b, |
| 5681 | 203 octave_idx_type& info, double& rcond, |
| 204 solve_singularity_handler sing_handler, | |
| 205 bool calc_cond = false) const; | |
| 5164 | 206 |
| 5785 | 207 SparseMatrix ltsolve (MatrixType &typ, const SparseMatrix& b, |
| 5681 | 208 octave_idx_type& info, double& rcond, |
| 209 solve_singularity_handler sing_handler, | |
| 210 bool calc_cond = false) const; | |
| 5164 | 211 |
| 5785 | 212 SparseComplexMatrix ltsolve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 213 octave_idx_type& info, double& rcond, |
| 5681 | 214 solve_singularity_handler sing_handler, |
| 215 bool calc_cond = false) const; | |
| 5164 | 216 |
| 217 // Tridiagonal matrix solvers | |
| 5785 | 218 Matrix trisolve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5681 | 219 double& rcond, solve_singularity_handler sing_handler, |
| 220 bool calc_cond = false) const; | |
| 5164 | 221 |
| 5785 | 222 ComplexMatrix trisolve (MatrixType &typ, const ComplexMatrix& b, |
| 5681 | 223 octave_idx_type& info, double& rcond, |
| 224 solve_singularity_handler sing_handler, | |
| 225 bool calc_cond = false) const; | |
| 5164 | 226 |
| 5785 | 227 SparseMatrix trisolve (MatrixType &typ, const SparseMatrix& b, |
| 5681 | 228 octave_idx_type& info, double& rcond, |
| 229 solve_singularity_handler sing_handler, | |
| 230 bool calc_cond = false) const; | |
| 5164 | 231 |
| 5785 | 232 SparseComplexMatrix trisolve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 233 octave_idx_type& info, double& rcond, |
| 5681 | 234 solve_singularity_handler sing_handler, |
| 235 bool calc_cond = false) const; | |
| 5164 | 236 |
| 237 // Banded matrix solvers (umfpack/cholesky) | |
| 5785 | 238 Matrix bsolve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5681 | 239 double& rcond, solve_singularity_handler sing_handler, |
| 240 bool calc_cond = false) const; | |
| 5164 | 241 |
| 5785 | 242 ComplexMatrix bsolve (MatrixType &typ, const ComplexMatrix& b, |
| 5681 | 243 octave_idx_type& info, double& rcond, |
| 244 solve_singularity_handler sing_handler, | |
| 245 bool calc_cond = false) const; | |
| 5164 | 246 |
| 5785 | 247 SparseMatrix bsolve (MatrixType &typ, const SparseMatrix& b, |
| 5681 | 248 octave_idx_type& info, double& rcond, |
| 249 solve_singularity_handler sing_handler, | |
| 250 bool calc_cond = false) const; | |
| 5164 | 251 |
| 5785 | 252 SparseComplexMatrix bsolve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 253 octave_idx_type& info, double& rcond, |
| 5681 | 254 solve_singularity_handler sing_handler, |
| 255 bool calc_cond = false) const; | |
| 5164 | 256 |
| 257 // Full matrix solvers (umfpack/cholesky) | |
| 5681 | 258 void * factorize (octave_idx_type& err, double &rcond, Matrix &Control, |
| 259 Matrix &Info, solve_singularity_handler sing_handler, | |
| 260 bool calc_cond = false) const; | |
| 5164 | 261 |
| 5785 | 262 Matrix fsolve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5681 | 263 double& rcond, solve_singularity_handler sing_handler, |
| 264 bool calc_cond = false) const; | |
| 5164 | 265 |
| 5785 | 266 ComplexMatrix fsolve (MatrixType &typ, const ComplexMatrix& b, |
| 5681 | 267 octave_idx_type& info, double& rcond, |
| 268 solve_singularity_handler sing_handler, | |
| 269 bool calc_cond = false) const; | |
| 5164 | 270 |
| 5785 | 271 SparseMatrix fsolve (MatrixType &typ, const SparseMatrix& b, |
| 5681 | 272 octave_idx_type& info, double& rcond, |
| 273 solve_singularity_handler sing_handler, | |
| 274 bool calc_cond = false) const; | |
| 5164 | 275 |
| 5785 | 276 SparseComplexMatrix fsolve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5681 | 277 octave_idx_type& info, double& rcond, |
| 278 solve_singularity_handler sing_handler, | |
| 279 bool calc_cond = false) const; | |
| 5164 | 280 |
| 281 public: | |
| 282 // Generic interface to solver with no probing of type | |
| 5785 | 283 Matrix solve (MatrixType &typ, const Matrix& b) const; |
| 284 Matrix solve (MatrixType &typ, const Matrix& b, octave_idx_type& info) const; | |
| 285 Matrix solve (MatrixType &typ, const Matrix& b, octave_idx_type& info, | |
| 5164 | 286 double& rcond) const; |
| 5785 | 287 Matrix solve (MatrixType &typ, const Matrix& b, octave_idx_type& info, |
| 5697 | 288 double& rcond, solve_singularity_handler sing_handler, |
| 289 bool singular_fallback = true) const; | |
| 5164 | 290 |
| 5785 | 291 ComplexMatrix solve (MatrixType &typ, const ComplexMatrix& b) const; |
| 292 ComplexMatrix solve (MatrixType &typ, const ComplexMatrix& b, | |
| 5275 | 293 octave_idx_type& info) const; |
| 5785 | 294 ComplexMatrix solve (MatrixType &typ, const ComplexMatrix& b, |
| 5697 | 295 octave_idx_type& info, double& rcond) const; |
| 5785 | 296 ComplexMatrix solve (MatrixType &typ, const ComplexMatrix& b, |
| 5697 | 297 octave_idx_type& info, double& rcond, |
| 298 solve_singularity_handler sing_handler, | |
| 299 bool singular_fallback = true) const; | |
| 5164 | 300 |
| 5785 | 301 SparseMatrix solve (MatrixType &typ, const SparseMatrix& b) const; |
| 302 SparseMatrix solve (MatrixType &typ, const SparseMatrix& b, | |
| 5275 | 303 octave_idx_type& info) const; |
| 5785 | 304 SparseMatrix solve (MatrixType &typ, const SparseMatrix& b, |
| 5697 | 305 octave_idx_type& info, double& rcond) const; |
| 5785 | 306 SparseMatrix solve (MatrixType &typ, const SparseMatrix& b, |
| 5697 | 307 octave_idx_type& info, double& rcond, |
| 308 solve_singularity_handler sing_handler, | |
| 309 bool singular_fallback = true) const; | |
| 5164 | 310 |
| 5785 | 311 SparseComplexMatrix solve (MatrixType &typ, |
| 5164 | 312 const SparseComplexMatrix& b) const; |
| 5785 | 313 SparseComplexMatrix solve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 314 octave_idx_type& info) const; |
| 5785 | 315 SparseComplexMatrix solve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5275 | 316 octave_idx_type& info, double& rcond) const; |
| 5785 | 317 SparseComplexMatrix solve (MatrixType &typ, const SparseComplexMatrix& b, |
| 5697 | 318 octave_idx_type& info, double& rcond, |
| 319 solve_singularity_handler sing_handler, | |
| 320 bool singular_fallabck = true) const; | |
| 5164 | 321 |
| 5785 | 322 ColumnVector solve (MatrixType &typ, const ColumnVector& b) const; |
| 323 ColumnVector solve (MatrixType &typ, const ColumnVector& b, | |
| 5275 | 324 octave_idx_type& info) const; |
| 5785 | 325 ColumnVector solve (MatrixType &typ, const ColumnVector& b, |
| 5275 | 326 octave_idx_type& info, double& rcond) const; |
| 5785 | 327 ColumnVector solve (MatrixType &typ, const ColumnVector& b, |
| 5697 | 328 octave_idx_type& info, double& rcond, |
| 329 solve_singularity_handler sing_handler) const; | |
| 5164 | 330 |
| 5785 | 331 ComplexColumnVector solve (MatrixType &typ, |
| 5164 | 332 const ComplexColumnVector& b) const; |
| 5785 | 333 ComplexColumnVector solve (MatrixType &typ, const ComplexColumnVector& b, |
| 5697 | 334 octave_idx_type& info) const; |
| 5785 | 335 ComplexColumnVector solve (MatrixType &typ, const ComplexColumnVector& b, |
| 5275 | 336 octave_idx_type& info, double& rcond) const; |
| 5785 | 337 ComplexColumnVector solve (MatrixType &typ, const ComplexColumnVector& b, |
| 5275 | 338 octave_idx_type& info, double& rcond, |
| 5164 | 339 solve_singularity_handler sing_handler) const; |
| 340 | |
| 341 // Generic interface to solver with probing of type | |
| 342 Matrix solve (const Matrix& b) const; | |
| 5275 | 343 Matrix solve (const Matrix& b, octave_idx_type& info) const; |
| 344 Matrix solve (const Matrix& b, octave_idx_type& info, double& rcond) const; | |
| 345 Matrix solve (const Matrix& b, octave_idx_type& info, double& rcond, | |
| 5164 | 346 solve_singularity_handler sing_handler) const; |
| 347 | |
| 348 ComplexMatrix solve (const ComplexMatrix& b) const; | |
| 5275 | 349 ComplexMatrix solve (const ComplexMatrix& b, octave_idx_type& info) const; |
| 350 ComplexMatrix solve (const ComplexMatrix& b, octave_idx_type& info, | |
| 5164 | 351 double& rcond) const; |
| 5275 | 352 ComplexMatrix solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond, |
| 5164 | 353 solve_singularity_handler sing_handler) const; |
| 354 | |
| 355 SparseMatrix solve (const SparseMatrix& b) const; | |
| 5275 | 356 SparseMatrix solve (const SparseMatrix& b, octave_idx_type& info) const; |
| 357 SparseMatrix solve (const SparseMatrix& b, octave_idx_type& info, | |
| 5164 | 358 double& rcond) const; |
| 5275 | 359 SparseMatrix solve (const SparseMatrix& b, octave_idx_type& info, double& rcond, |
| 5164 | 360 solve_singularity_handler sing_handler) const; |
| 361 | |
| 362 SparseComplexMatrix solve (const SparseComplexMatrix& b) const; | |
| 5275 | 363 SparseComplexMatrix solve (const SparseComplexMatrix& b, octave_idx_type& info) const; |
| 364 SparseComplexMatrix solve (const SparseComplexMatrix& b, octave_idx_type& info, | |
| 5164 | 365 double& rcond) const; |
| 5275 | 366 SparseComplexMatrix solve (const SparseComplexMatrix& b, octave_idx_type& info, |
| 5164 | 367 double& rcond, |
| 368 solve_singularity_handler sing_handler) const; | |
| 369 | |
| 370 ColumnVector solve (const ColumnVector& b) const; | |
| 5275 | 371 ColumnVector solve (const ColumnVector& b, octave_idx_type& info) const; |
| 372 ColumnVector solve (const ColumnVector& b, octave_idx_type& info, double& rcond) const; | |
| 373 ColumnVector solve (const ColumnVector& b, octave_idx_type& info, double& rcond, | |
| 5164 | 374 solve_singularity_handler sing_handler) const; |
| 375 | |
| 376 ComplexColumnVector solve (const ComplexColumnVector& b) const; | |
| 5275 | 377 ComplexColumnVector solve (const ComplexColumnVector& b, octave_idx_type& info) const; |
| 378 ComplexColumnVector solve (const ComplexColumnVector& b, octave_idx_type& info, | |
| 5164 | 379 double& rcond) const; |
| 5275 | 380 ComplexColumnVector solve (const ComplexColumnVector& b, octave_idx_type& info, |
| 5164 | 381 double& rcond, |
| 382 solve_singularity_handler sing_handler) const; | |
| 383 | |
| 384 // other operations | |
| 385 | |
| 386 bool any_element_is_negative (bool = false) const; | |
|
7922
935be827eaf8
error for NaN values in & and | expressions
John W. Eaton <jwe@octave.org>
parents:
7802
diff
changeset
|
387 bool any_element_is_nan (void) const; |
| 5164 | 388 bool any_element_is_inf_or_nan (void) const; |
| 6989 | 389 bool all_elements_are_zero (void) const; |
| 5164 | 390 bool all_elements_are_int_or_inf_or_nan (void) const; |
| 391 bool all_integers (double& max_val, double& min_val) const; | |
| 392 bool too_large_for_float (void) const; | |
| 393 | |
| 394 SparseBoolMatrix operator ! (void) const; | |
| 395 | |
| 396 SparseBoolMatrix all (int dim = -1) const; | |
| 397 SparseBoolMatrix any (int dim = -1) const; | |
| 398 | |
| 399 SparseMatrix cumprod (int dim = -1) const; | |
| 400 SparseMatrix cumsum (int dim = -1) const; | |
| 401 SparseMatrix prod (int dim = -1) const; | |
| 402 SparseMatrix sum (int dim = -1) const; | |
| 403 SparseMatrix sumsq (int dim = -1) const; | |
| 404 SparseMatrix abs (void) const; | |
| 405 | |
| 5275 | 406 SparseMatrix diag (octave_idx_type k = 0) const; |
| 5164 | 407 |
| 408 Matrix matrix_value (void) const; | |
| 409 | |
| 410 SparseMatrix squeeze (void) const; | |
| 411 | |
| 412 SparseMatrix index (idx_vector& i, int resize_ok) const; | |
| 413 | |
| 414 SparseMatrix index (idx_vector& i, idx_vector& j, int resize_ok) const; | |
| 415 | |
| 416 SparseMatrix index (Array<idx_vector>& ra_idx, int resize_ok) const; | |
| 417 | |
| 418 SparseMatrix reshape (const dim_vector& new_dims) const; | |
| 419 | |
| 5275 | 420 SparseMatrix permute (const Array<octave_idx_type>& vec, bool inv = false) const; |
| 5164 | 421 |
| 5275 | 422 SparseMatrix ipermute (const Array<octave_idx_type>& vec) const; |
| 5164 | 423 |
| 424 // i/o | |
| 425 | |
| 6108 | 426 friend OCTAVE_API std::ostream& operator << (std::ostream& os, const SparseMatrix& a); |
| 427 friend OCTAVE_API std::istream& operator >> (std::istream& is, SparseMatrix& a); | |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
428 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
429 typedef double (*dmapper) (double); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
430 typedef Complex (*cmapper) (const Complex&); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
431 typedef bool (*bmapper) (double); |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
432 SparseMatrix map (dmapper fcn) const; |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
433 SparseComplexMatrix map (cmapper fcn) const; |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
434 SparseBoolMatrix map (bmapper fcn) const; |
| 5164 | 435 }; |
| 436 | |
| 5509 | 437 // Publish externally used friend functions. |
| 438 | |
| 6108 | 439 extern OCTAVE_API SparseMatrix real (const SparseComplexMatrix& a); |
| 440 extern OCTAVE_API SparseMatrix imag (const SparseComplexMatrix& a); | |
| 5509 | 441 |
| 442 // Other operators. | |
| 443 | |
| 6108 | 444 extern OCTAVE_API SparseMatrix operator * (const SparseMatrix& a, |
| 5164 | 445 const SparseMatrix& b); |
| 6108 | 446 extern OCTAVE_API Matrix operator * (const Matrix& a, |
| 5429 | 447 const SparseMatrix& b); |
|
7802
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7503
diff
changeset
|
448 extern OCTAVE_API Matrix mul_trans (const Matrix& a, |
|
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7503
diff
changeset
|
449 const SparseMatrix& b); |
| 6108 | 450 extern OCTAVE_API Matrix operator * (const SparseMatrix& a, |
| 5429 | 451 const Matrix& b); |
|
7802
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7503
diff
changeset
|
452 extern OCTAVE_API Matrix trans_mul (const SparseMatrix& a, |
|
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7503
diff
changeset
|
453 const Matrix& b); |
| 5164 | 454 |
|
8964
f4f4d65faaa0
Implement sparse * diagonal and diagonal * sparse operations, double-prec only.
Jason Riedy <jason@acm.org>
parents:
8920
diff
changeset
|
455 extern OCTAVE_API SparseMatrix operator * (const DiagMatrix&, const SparseMatrix&); |
|
f4f4d65faaa0
Implement sparse * diagonal and diagonal * sparse operations, double-prec only.
Jason Riedy <jason@acm.org>
parents:
8920
diff
changeset
|
456 extern OCTAVE_API SparseMatrix operator * (const SparseMatrix&, const DiagMatrix&); |
|
f4f4d65faaa0
Implement sparse * diagonal and diagonal * sparse operations, double-prec only.
Jason Riedy <jason@acm.org>
parents:
8920
diff
changeset
|
457 |
|
8966
1bba53c0a38d
Implement diag + sparse, diag - sparse, sparse + diag, sparse - diag.
Jason Riedy <jason@acm.org>
parents:
8964
diff
changeset
|
458 extern OCTAVE_API SparseMatrix operator + (const DiagMatrix&, const SparseMatrix&); |
|
1bba53c0a38d
Implement diag + sparse, diag - sparse, sparse + diag, sparse - diag.
Jason Riedy <jason@acm.org>
parents:
8964
diff
changeset
|
459 extern OCTAVE_API SparseMatrix operator + (const SparseMatrix&, const DiagMatrix&); |
|
1bba53c0a38d
Implement diag + sparse, diag - sparse, sparse + diag, sparse - diag.
Jason Riedy <jason@acm.org>
parents:
8964
diff
changeset
|
460 extern OCTAVE_API SparseMatrix operator - (const DiagMatrix&, const SparseMatrix&); |
|
1bba53c0a38d
Implement diag + sparse, diag - sparse, sparse + diag, sparse - diag.
Jason Riedy <jason@acm.org>
parents:
8964
diff
changeset
|
461 extern OCTAVE_API SparseMatrix operator - (const SparseMatrix&, const DiagMatrix&); |
|
1bba53c0a38d
Implement diag + sparse, diag - sparse, sparse + diag, sparse - diag.
Jason Riedy <jason@acm.org>
parents:
8964
diff
changeset
|
462 |
|
8968
91d53dc37f79
Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations.
Jason Riedy <jason@acm.org>
parents:
8966
diff
changeset
|
463 extern OCTAVE_API SparseMatrix operator * (const PermMatrix&, const SparseMatrix&); |
|
91d53dc37f79
Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations.
Jason Riedy <jason@acm.org>
parents:
8966
diff
changeset
|
464 extern OCTAVE_API SparseMatrix operator * (const SparseMatrix&, const PermMatrix&); |
|
91d53dc37f79
Add perm * sparse, perm \ sparse, sparse * perm, and sparse / perm operations.
Jason Riedy <jason@acm.org>
parents:
8966
diff
changeset
|
465 |
| 6108 | 466 extern OCTAVE_API SparseMatrix min (double d, const SparseMatrix& m); |
| 467 extern OCTAVE_API SparseMatrix min (const SparseMatrix& m, double d); | |
| 468 extern OCTAVE_API SparseMatrix min (const SparseMatrix& a, const SparseMatrix& b); | |
| 5164 | 469 |
| 6108 | 470 extern OCTAVE_API SparseMatrix max (double d, const SparseMatrix& m); |
| 471 extern OCTAVE_API SparseMatrix max (const SparseMatrix& m, double d); | |
| 472 extern OCTAVE_API SparseMatrix max (const SparseMatrix& a, const SparseMatrix& b); | |
| 5164 | 473 |
| 6708 | 474 SPARSE_SMS_CMP_OP_DECLS (SparseMatrix, double, OCTAVE_API) |
| 475 SPARSE_SMS_BOOL_OP_DECLS (SparseMatrix, double, OCTAVE_API) | |
| 5164 | 476 |
| 6708 | 477 SPARSE_SSM_CMP_OP_DECLS (double, SparseMatrix, OCTAVE_API) |
| 478 SPARSE_SSM_BOOL_OP_DECLS (double, SparseMatrix, OCTAVE_API) | |
| 5164 | 479 |
| 6708 | 480 SPARSE_SMSM_CMP_OP_DECLS (SparseMatrix, SparseMatrix, OCTAVE_API) |
| 481 SPARSE_SMSM_BOOL_OP_DECLS (SparseMatrix, SparseMatrix, OCTAVE_API) | |
| 5164 | 482 |
| 483 SPARSE_FORWARD_DEFS (MSparse, SparseMatrix, Matrix, double) | |
| 484 | |
| 5351 | 485 #ifdef IDX_TYPE_LONG |
| 5322 | 486 #define UMFPACK_DNAME(name) umfpack_dl_ ## name |
| 487 #else | |
| 488 #define UMFPACK_DNAME(name) umfpack_di_ ## name | |
| 489 #endif | |
| 490 | |
| 5164 | 491 #endif |
| 492 | |
| 493 /* | |
| 494 ;;; Local Variables: *** | |
| 495 ;;; mode: C++ *** | |
| 496 ;;; End: *** | |
| 497 */ |