GNU Octave 10.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
 
Loading...
Searching...
No Matches
dSparse.h
Go to the documentation of this file.
1////////////////////////////////////////////////////////////////////////
2//
3// Copyright (C) 1998-2025 The Octave Project Developers
4//
5// See the file COPYRIGHT.md in the top-level directory of this
6// distribution or <https://octave.org/copyright/>.
7//
8// This file is part of Octave.
9//
10// Octave is free software: you can redistribute it and/or modify it
11// under the terms of the GNU General Public License as published by
12// the Free Software Foundation, either version 3 of the License, or
13// (at your option) any later version.
14//
15// Octave is distributed in the hope that it will be useful, but
16// WITHOUT ANY WARRANTY; without even the implied warranty of
17// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18// GNU General Public License for more details.
19//
20// You should have received a copy of the GNU General Public License
21// along with Octave; see the file COPYING. If not, see
22// <https://www.gnu.org/licenses/>.
23//
24////////////////////////////////////////////////////////////////////////
25
26#if ! defined (octave_dSparse_h)
27#define octave_dSparse_h 1
28
29#include "octave-config.h"
30
31#include "mx-fwd.h"
32
33#include "CColVector.h"
34#include "CMatrix.h"
35#include "DET.h"
36#include "MSparse.h"
37#include "MatrixType.h"
38#include "Sparse-op-decls.h"
39#include "dColVector.h"
40#include "dMatrix.h"
41#include "dNDArray.h"
42
43class SparseMatrix : public MSparse<double>
44{
45public:
46
47 // Corresponding dense matrix type for this sparse matrix type.
48 typedef Matrix dense_matrix_type;
49
50 typedef void (*solve_singularity_handler) (double rcond);
51
52 SparseMatrix () : MSparse<double> () { }
53
56
57 SparseMatrix (const dim_vector& dv, octave_idx_type nz = 0) :
58 MSparse<double> (dv, nz) { }
59
60 explicit SparseMatrix (octave_idx_type r, octave_idx_type c, double val)
61 : MSparse<double> (r, c, val) { }
62
63 SparseMatrix (const SparseMatrix& a) : MSparse<double> (a) { }
64
65 SparseMatrix (const SparseMatrix& a, const dim_vector& dv)
66 : MSparse<double> (a, dv) { }
67
68 SparseMatrix (const MSparse<double>& a) : MSparse<double> (a) { }
69
70 SparseMatrix (const Sparse<double>& a) : MSparse<double> (a) { }
71
72 explicit OCTAVE_API SparseMatrix (const SparseBoolMatrix& a);
73
74 explicit SparseMatrix (const Matrix& a) : MSparse<double> (a) { }
75
76 explicit SparseMatrix (const NDArray& a) : MSparse<double> (a) { }
77
78 SparseMatrix (const Array<double>& a, const octave::idx_vector& r,
79 const octave::idx_vector& c, octave_idx_type nr = -1,
80 octave_idx_type nc = -1, bool sum_terms = true,
81 octave_idx_type nzm = -1)
82 : MSparse<double> (a, r, c, nr, nc, sum_terms, nzm) { }
83
84 explicit OCTAVE_API SparseMatrix (const DiagMatrix& a);
85
86 explicit SparseMatrix (const PermMatrix& a) : MSparse<double> (a) { }
87
90
91 SparseMatrix& operator = (const SparseMatrix& a)
92 {
93 MSparse<double>::operator = (a);
94 return *this;
95 }
96
97 ~SparseMatrix () = default;
98
99 OCTAVE_API bool operator == (const SparseMatrix& a) const;
100 OCTAVE_API bool operator != (const SparseMatrix& a) const;
101
102 OCTAVE_API bool issymmetric () const;
103
104 OCTAVE_API SparseMatrix max (int dim = -1) const;
105 OCTAVE_API SparseMatrix max (Array<octave_idx_type>& index, int dim = -1) const;
106 OCTAVE_API SparseMatrix min (int dim = -1) const;
107 OCTAVE_API SparseMatrix min (Array<octave_idx_type>& index, int dim = -1) const;
108
109 // destructive insert/delete/reorder operations
110
111 OCTAVE_API SparseMatrix&
112 insert (const SparseMatrix& a, octave_idx_type r, octave_idx_type c);
113
114 OCTAVE_API SparseMatrix&
115 insert (const SparseMatrix& a, const Array<octave_idx_type>& indx);
116
117 OCTAVE_API SparseMatrix
118 concat (const SparseMatrix& rb, const Array<octave_idx_type>& ra_idx);
119 OCTAVE_API SparseComplexMatrix
120 concat (const SparseComplexMatrix& rb, const Array<octave_idx_type>& ra_idx);
121
122 friend OCTAVE_API SparseMatrix real (const SparseComplexMatrix& a);
123 friend OCTAVE_API SparseMatrix imag (const SparseComplexMatrix& a);
124
125 SparseMatrix transpose () const
126 {
127 return MSparse<double>::transpose ();
128 }
129 SparseMatrix hermitian () const { return transpose (); }
130
131 // extract row or column i.
132
133 OCTAVE_API RowVector row (octave_idx_type i) const;
134
135 OCTAVE_API ColumnVector column (octave_idx_type i) const;
136
137private:
138 OCTAVE_API SparseMatrix
139 dinverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
140 const bool force = false, const bool calccond = true) const;
141
142 OCTAVE_API SparseMatrix
143 tinverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
144 const bool force = false, const bool calccond = true) const;
145
146public:
147 OCTAVE_API SparseMatrix inverse () const;
148 OCTAVE_API SparseMatrix inverse (MatrixType& mattype) const;
149 OCTAVE_API SparseMatrix
150 inverse (MatrixType& mattype, octave_idx_type& info) const;
151 OCTAVE_API SparseMatrix
152 inverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
153 bool force = false, bool calc_cond = true) const;
154
155 OCTAVE_API DET determinant () const;
156 OCTAVE_API DET determinant (octave_idx_type& info) const;
157 OCTAVE_API DET determinant (octave_idx_type& info, double& rcond,
158 bool calc_cond = true) const;
159
160private:
161 // Diagonal matrix solvers
162 OCTAVE_API Matrix
163 dsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
164 double& rcond, solve_singularity_handler sing_handler,
165 bool calc_cond = false) const;
166
167 OCTAVE_API ComplexMatrix
168 dsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
169 double& rcond, solve_singularity_handler sing_handler,
170 bool calc_cond = false) const;
171
172 OCTAVE_API SparseMatrix
173 dsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
174 double& rcond, solve_singularity_handler sing_handler,
175 bool calc_cond = false) const;
176
177 OCTAVE_API SparseComplexMatrix
178 dsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
179 double& rcond, solve_singularity_handler sing_handler,
180 bool calc_cond = false) const;
181
182 // Upper triangular matrix solvers
183 OCTAVE_API Matrix
184 utsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
185 double& rcond, solve_singularity_handler sing_handler,
186 bool calc_cond = false) const;
187
188 OCTAVE_API ComplexMatrix
189 utsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
190 double& rcond, solve_singularity_handler sing_handler,
191 bool calc_cond = false) const;
192
193 OCTAVE_API SparseMatrix
194 utsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
195 double& rcond, solve_singularity_handler sing_handler,
196 bool calc_cond = false) const;
197
198 OCTAVE_API SparseComplexMatrix
199 utsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
200 double& rcond, solve_singularity_handler sing_handler,
201 bool calc_cond = false) const;
202
203 // Lower triangular matrix solvers
204 OCTAVE_API Matrix
205 ltsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
206 double& rcond, solve_singularity_handler sing_handler,
207 bool calc_cond = false) const;
208
209 OCTAVE_API ComplexMatrix
210 ltsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
211 double& rcond, solve_singularity_handler sing_handler,
212 bool calc_cond = false) const;
213
214 OCTAVE_API SparseMatrix
215 ltsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
216 double& rcond, solve_singularity_handler sing_handler,
217 bool calc_cond = false) const;
218
219 OCTAVE_API SparseComplexMatrix
220 ltsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
221 double& rcond, solve_singularity_handler sing_handler,
222 bool calc_cond = false) const;
223
224 // Tridiagonal matrix solvers
225 OCTAVE_API Matrix
226 trisolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
227 double& rcond, solve_singularity_handler sing_handler,
228 bool calc_cond = false) const;
229
230 OCTAVE_API ComplexMatrix
231 trisolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
232 double& rcond, solve_singularity_handler sing_handler,
233 bool calc_cond = false) const;
234
235 OCTAVE_API SparseMatrix
236 trisolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
237 double& rcond, solve_singularity_handler sing_handler,
238 bool calc_cond = false) const;
239
240 OCTAVE_API SparseComplexMatrix
241 trisolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
242 double& rcond, solve_singularity_handler sing_handler,
243 bool calc_cond = false) const;
244
245 // Banded matrix solvers (umfpack/cholesky)
246 OCTAVE_API Matrix
247 bsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
248 double& rcond, solve_singularity_handler sing_handler,
249 bool calc_cond = false) const;
250
251 OCTAVE_API ComplexMatrix
252 bsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
253 double& rcond, solve_singularity_handler sing_handler,
254 bool calc_cond = false) const;
255
256 OCTAVE_API SparseMatrix
257 bsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
258 double& rcond, solve_singularity_handler sing_handler,
259 bool calc_cond = false) const;
260
261 OCTAVE_API SparseComplexMatrix
262 bsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
263 double& rcond, solve_singularity_handler sing_handler,
264 bool calc_cond = false) const;
265
266 // Full matrix solvers (umfpack/cholesky)
267 OCTAVE_API void *
268 factorize (octave_idx_type& err, double& rcond, Matrix& Control,
269 Matrix& Info, solve_singularity_handler sing_handler,
270 bool calc_cond = false) const;
271
272 OCTAVE_API Matrix
273 fsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
274 double& rcond, solve_singularity_handler sing_handler,
275 bool calc_cond = false) const;
276
277 OCTAVE_API ComplexMatrix
278 fsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
279 double& rcond, solve_singularity_handler sing_handler,
280 bool calc_cond = false) const;
281
282 OCTAVE_API SparseMatrix
283 fsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
284 double& rcond, solve_singularity_handler sing_handler,
285 bool calc_cond = false) const;
286
287 OCTAVE_API SparseComplexMatrix
288 fsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
289 double& rcond, solve_singularity_handler sing_handler,
290 bool calc_cond = false) const;
291
292public:
293 // Generic interface to solver with no probing of type
294 OCTAVE_API Matrix solve (MatrixType& typ, const Matrix& b) const;
295 OCTAVE_API Matrix
296 solve (MatrixType& typ, const Matrix& b, octave_idx_type& info) const;
297 OCTAVE_API Matrix
298 solve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
299 double& rcond) const;
300 OCTAVE_API Matrix
301 solve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
302 double& rcond, solve_singularity_handler sing_handler,
303 bool singular_fallback = true) const;
304
305 OCTAVE_API ComplexMatrix
306 solve (MatrixType& typ, const ComplexMatrix& b) const;
307 OCTAVE_API ComplexMatrix
308 solve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info) const;
309 OCTAVE_API ComplexMatrix
310 solve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
311 double& rcond) const;
312 OCTAVE_API ComplexMatrix
313 solve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
314 double& rcond, solve_singularity_handler sing_handler,
315 bool singular_fallback = true) const;
316
317 OCTAVE_API SparseMatrix solve (MatrixType& typ, const SparseMatrix& b) const;
318 OCTAVE_API SparseMatrix
319 solve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info) const;
320 OCTAVE_API SparseMatrix
321 solve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
322 double& rcond) const;
323 OCTAVE_API SparseMatrix
324 solve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
325 double& rcond, solve_singularity_handler sing_handler,
326 bool singular_fallback = true) const;
327
328 OCTAVE_API SparseComplexMatrix
329 solve (MatrixType& typ, const SparseComplexMatrix& b) const;
330 OCTAVE_API SparseComplexMatrix
331 solve (MatrixType& typ, const SparseComplexMatrix& b,
332 octave_idx_type& info) const;
333 OCTAVE_API SparseComplexMatrix
334 solve (MatrixType& typ, const SparseComplexMatrix& b,
335 octave_idx_type& info, double& rcond) const;
336 OCTAVE_API SparseComplexMatrix
337 solve (MatrixType& typ, const SparseComplexMatrix& b,
338 octave_idx_type& info, double& rcond,
339 solve_singularity_handler sing_handler,
340 bool singular_fallabck = true) const;
341
342 OCTAVE_API ColumnVector solve (MatrixType& typ, const ColumnVector& b) const;
343 OCTAVE_API ColumnVector
344 solve (MatrixType& typ, const ColumnVector& b, octave_idx_type& info) const;
345 OCTAVE_API ColumnVector
346 solve (MatrixType& typ, const ColumnVector& b, octave_idx_type& info,
347 double& rcond) const;
348 OCTAVE_API ColumnVector
349 solve (MatrixType& typ, const ColumnVector& b, octave_idx_type& info,
350 double& rcond, solve_singularity_handler sing_handler) const;
351
352 OCTAVE_API ComplexColumnVector
353 solve (MatrixType& typ, const ComplexColumnVector& b) const;
354 OCTAVE_API ComplexColumnVector
355 solve (MatrixType& typ, const ComplexColumnVector& b,
356 octave_idx_type& info) const;
357 OCTAVE_API ComplexColumnVector
358 solve (MatrixType& typ, const ComplexColumnVector& b,
359 octave_idx_type& info, double& rcond) const;
360 OCTAVE_API ComplexColumnVector
361 solve (MatrixType& typ, const ComplexColumnVector& b,
362 octave_idx_type& info, double& rcond,
363 solve_singularity_handler sing_handler) const;
364
365 // Generic interface to solver with probing of type
366 OCTAVE_API Matrix solve (const Matrix& b) const;
367 OCTAVE_API Matrix solve (const Matrix& b, octave_idx_type& info) const;
368 OCTAVE_API Matrix
369 solve (const Matrix& b, octave_idx_type& info, double& rcond) const;
370 OCTAVE_API Matrix
371 solve (const Matrix& b, octave_idx_type& info, double& rcond,
372 solve_singularity_handler sing_handler) const;
373
374 OCTAVE_API ComplexMatrix solve (const ComplexMatrix& b) const;
375 OCTAVE_API ComplexMatrix
376 solve (const ComplexMatrix& b, octave_idx_type& info) const;
377 OCTAVE_API ComplexMatrix
378 solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond) const;
379 OCTAVE_API ComplexMatrix
380 solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond,
381 solve_singularity_handler sing_handler) const;
382
383 OCTAVE_API SparseMatrix solve (const SparseMatrix& b) const;
384 OCTAVE_API SparseMatrix
385 solve (const SparseMatrix& b, octave_idx_type& info) const;
386 OCTAVE_API SparseMatrix
387 solve (const SparseMatrix& b, octave_idx_type& info, double& rcond) const;
388 OCTAVE_API SparseMatrix
389 solve (const SparseMatrix& b, octave_idx_type& info, double& rcond,
390 solve_singularity_handler sing_handler) const;
391
392 OCTAVE_API SparseComplexMatrix solve (const SparseComplexMatrix& b) const;
393 OCTAVE_API SparseComplexMatrix
394 solve (const SparseComplexMatrix& b, octave_idx_type& info) const;
395 OCTAVE_API SparseComplexMatrix
396 solve (const SparseComplexMatrix& b, octave_idx_type& info,
397 double& rcond) const;
398 OCTAVE_API SparseComplexMatrix
399 solve (const SparseComplexMatrix& b, octave_idx_type& info,
400 double& rcond, solve_singularity_handler sing_handler) const;
401
402 OCTAVE_API ColumnVector solve (const ColumnVector& b) const;
403 OCTAVE_API ColumnVector
404 solve (const ColumnVector& b, octave_idx_type& info) const;
405 OCTAVE_API ColumnVector
406 solve (const ColumnVector& b, octave_idx_type& info, double& rcond) const;
407 OCTAVE_API ColumnVector
408 solve (const ColumnVector& b, octave_idx_type& info, double& rcond,
409 solve_singularity_handler sing_handler) const;
410
411 OCTAVE_API ComplexColumnVector solve (const ComplexColumnVector& b) const;
412 OCTAVE_API ComplexColumnVector
413 solve (const ComplexColumnVector& b, octave_idx_type& info) const;
414 OCTAVE_API ComplexColumnVector
415 solve (const ComplexColumnVector& b, octave_idx_type& info,
416 double& rcond) const;
417 OCTAVE_API ComplexColumnVector
418 solve (const ComplexColumnVector& b, octave_idx_type& info,
419 double& rcond, solve_singularity_handler sing_handler) const;
420
421 // other operations
422
423 OCTAVE_API bool any_element_is_negative (bool = false) const;
424 OCTAVE_API bool any_element_is_nan () const;
425 OCTAVE_API bool any_element_is_inf_or_nan () const;
426 OCTAVE_API bool any_element_not_one_or_zero () const;
427 OCTAVE_API bool all_elements_are_zero () const;
428 OCTAVE_API bool all_elements_are_int_or_inf_or_nan () const;
429 OCTAVE_API bool all_integers (double& max_val, double& min_val) const;
430 OCTAVE_API bool too_large_for_float () const;
431
432 OCTAVE_API SparseBoolMatrix operator ! () const;
433
434 OCTAVE_API SparseBoolMatrix all (int dim = -1) const;
435 OCTAVE_API SparseBoolMatrix any (int dim = -1) const;
436
437 OCTAVE_API SparseMatrix cumprod (int dim = -1) const;
438 OCTAVE_API SparseMatrix cumsum (int dim = -1) const;
439 OCTAVE_API SparseMatrix prod (int dim = -1) const;
440 OCTAVE_API SparseMatrix sum (int dim = -1) const;
441 OCTAVE_API SparseMatrix sumsq (int dim = -1) const;
442 OCTAVE_API SparseMatrix abs () const;
443
444 OCTAVE_API SparseMatrix diag (octave_idx_type k = 0) const;
445
446 OCTAVE_API Matrix matrix_value () const;
447
448 OCTAVE_API SparseMatrix squeeze () const;
449
450 OCTAVE_API SparseMatrix reshape (const dim_vector& new_dims) const;
451
452 OCTAVE_API SparseMatrix
453 permute (const Array<octave_idx_type>& vec, bool inv = false) const;
454
455 OCTAVE_API SparseMatrix ipermute (const Array<octave_idx_type>& vec) const;
456
457 // i/o
458
459 friend OCTAVE_API std::ostream& operator << (std::ostream& os,
460 const SparseMatrix& a);
461 friend OCTAVE_API std::istream& operator >> (std::istream& is,
462 SparseMatrix& a);
463
464};
465
466// Publish externally used friend functions.
467
468extern OCTAVE_API SparseMatrix real (const SparseComplexMatrix& a);
469extern OCTAVE_API SparseMatrix imag (const SparseComplexMatrix& a);
470
471// Other operators.
472
473extern OCTAVE_API SparseMatrix operator * (const SparseMatrix& a,
474 const SparseMatrix& b);
475extern OCTAVE_API Matrix operator * (const Matrix& a,
476 const SparseMatrix& b);
477extern OCTAVE_API Matrix mul_trans (const Matrix& a,
478 const SparseMatrix& b);
479extern OCTAVE_API Matrix operator * (const SparseMatrix& a,
480 const Matrix& b);
481extern OCTAVE_API Matrix trans_mul (const SparseMatrix& a,
482 const Matrix& b);
483
484extern OCTAVE_API SparseMatrix operator * (const DiagMatrix&,
485 const SparseMatrix&);
486extern OCTAVE_API SparseMatrix operator * (const SparseMatrix&,
487 const DiagMatrix&);
488
489extern OCTAVE_API SparseMatrix operator + (const DiagMatrix&,
490 const SparseMatrix&);
491extern OCTAVE_API SparseMatrix operator + (const SparseMatrix&,
492 const DiagMatrix&);
493extern OCTAVE_API SparseMatrix operator - (const DiagMatrix&,
494 const SparseMatrix&);
495extern OCTAVE_API SparseMatrix operator - (const SparseMatrix&,
496 const DiagMatrix&);
497
498extern OCTAVE_API SparseMatrix operator * (const PermMatrix&,
499 const SparseMatrix&);
500extern OCTAVE_API SparseMatrix operator * (const SparseMatrix&,
501 const PermMatrix&);
502
503extern OCTAVE_API SparseMatrix min (double d, const SparseMatrix& m);
504extern OCTAVE_API SparseMatrix min (const SparseMatrix& m, double d);
505extern OCTAVE_API SparseMatrix min (const SparseMatrix& a,
506 const SparseMatrix& b);
507
508extern OCTAVE_API SparseMatrix max (double d, const SparseMatrix& m);
509extern OCTAVE_API SparseMatrix max (const SparseMatrix& m, double d);
510extern OCTAVE_API SparseMatrix max (const SparseMatrix& a,
511 const SparseMatrix& b);
512
513SPARSE_SMS_CMP_OP_DECLS (SparseMatrix, double, OCTAVE_API)
514SPARSE_SMS_BOOL_OP_DECLS (SparseMatrix, double, OCTAVE_API)
515
516SPARSE_SSM_CMP_OP_DECLS (double, SparseMatrix, OCTAVE_API)
517SPARSE_SSM_BOOL_OP_DECLS (double, SparseMatrix, OCTAVE_API)
518
519SPARSE_SMSM_CMP_OP_DECLS (SparseMatrix, SparseMatrix, OCTAVE_API)
520SPARSE_SMSM_BOOL_OP_DECLS (SparseMatrix, SparseMatrix, OCTAVE_API)
521
522SPARSE_FORWARD_DEFS (MSparse, SparseMatrix, Matrix, double)
523
524#endif
SPARSE_FORWARD_DEFS
#define SPARSE_FORWARD_DEFS(B, R, F, T)
Definition MSparse.h:192
SPARSE_SMSM_CMP_OP_DECLS
#define SPARSE_SMSM_CMP_OP_DECLS(M1, M2, API)
Definition Sparse-op-decls.h:108
SPARSE_SMSM_BOOL_OP_DECLS
#define SPARSE_SMSM_BOOL_OP_DECLS(M1, M2, API)
Definition Sparse-op-decls.h:120
SPARSE_SSM_CMP_OP_DECLS
#define SPARSE_SSM_CMP_OP_DECLS(S, M, API)
Definition Sparse-op-decls.h:79
SPARSE_SMS_BOOL_OP_DECLS
#define SPARSE_SMS_BOOL_OP_DECLS(M, S, API)
Definition Sparse-op-decls.h:62
SPARSE_SSM_BOOL_OP_DECLS
#define SPARSE_SSM_BOOL_OP_DECLS(S, M, API)
Definition Sparse-op-decls.h:91
SPARSE_SMS_CMP_OP_DECLS
#define SPARSE_SMS_CMP_OP_DECLS(M, S, API)
Definition Sparse-op-decls.h:50
Array
N Dimensional Array with copy-on-write semantics.
Definition Array.h:130
MSparse::transpose
MSparse< T > transpose() const
Definition MSparse.h:96
MSparse::operator=
MSparse< T > & operator=(const MSparse< T > &a)
Definition MSparse.h:78
SparseMatrix::insert
SparseMatrix & insert(const SparseMatrix &a, octave_idx_type r, octave_idx_type c)
Definition dSparse.cc:168
SparseMatrix::squeeze
SparseMatrix squeeze() const
Definition dSparse.cc:7433
SparseMatrix::issymmetric
bool issymmetric() const
Definition dSparse.cc:128
SparseMatrix::any_element_not_one_or_zero
bool any_element_not_one_or_zero() const
Definition dSparse.cc:7202
SparseMatrix::reshape
SparseMatrix reshape(const dim_vector &new_dims) const
Definition dSparse.cc:7439
SparseMatrix::sum
SparseMatrix sum(int dim=-1) const
Definition dSparse.cc:7355
SparseMatrix::any
SparseBoolMatrix any(int dim=-1) const
Definition dSparse.cc:7325
SparseMatrix::any_element_is_inf_or_nan
bool any_element_is_inf_or_nan() const
Definition dSparse.cc:7187
SparseMatrix::all_integers
bool all_integers(double &max_val, double &min_val) const
Definition dSparse.cc:7249
SparseMatrix::any_element_is_nan
bool any_element_is_nan() const
Definition dSparse.cc:7172
SparseMatrix::too_large_for_float
bool too_large_for_float() const
Definition dSparse.cc:7277
SparseMatrix::operator==
bool operator==(const SparseMatrix &a) const
Definition dSparse.cc:98
SparseMatrix::transpose
SparseMatrix transpose() const
Definition dSparse.h:125
SparseMatrix::SparseMatrix
SparseMatrix(const Sparse< double > &a)
Definition dSparse.h:70
SparseMatrix::determinant
DET determinant() const
Definition dSparse.cc:1017
SparseMatrix::abs
SparseMatrix abs() const
Definition dSparse.cc:7379
SparseMatrix::inverse
SparseMatrix inverse() const
Definition dSparse.cc:608
SparseMatrix::operator!=
bool operator!=(const SparseMatrix &a) const
Definition dSparse.cc:122
SparseMatrix::row
RowVector row(octave_idx_type i) const
Definition dSparse.cc:500
SparseMatrix::solve
ComplexMatrix solve(const ComplexMatrix &b) const
SparseMatrix::cumprod
SparseMatrix cumprod(int dim=-1) const
Definition dSparse.cc:7331
SparseMatrix::real
friend SparseMatrix real(const SparseComplexMatrix &a)
Definition dSparse.cc:551
SparseMatrix::diag
SparseMatrix diag(octave_idx_type k=0) const
Definition dSparse.cc:7392
SparseMatrix::prod
SparseMatrix prod(int dim=-1) const
Definition dSparse.cc:7343
SparseMatrix::SparseMatrix
SparseMatrix(const SparseMatrix &a, const dim_vector &dv)
Definition dSparse.h:65
SparseMatrix::all_elements_are_zero
bool all_elements_are_zero() const
Definition dSparse.cc:7217
SparseMatrix::~SparseMatrix
~SparseMatrix()=default
SparseMatrix::all
SparseBoolMatrix all(int dim=-1) const
Definition dSparse.cc:7319
SparseMatrix::SparseMatrix
SparseMatrix(const SparseMatrix &a)
Definition dSparse.h:63
SparseMatrix::min
SparseMatrix min(int dim=-1) const
Definition dSparse.cc:334
SparseMatrix::permute
SparseMatrix permute(const Array< octave_idx_type > &vec, bool inv=false) const
Definition dSparse.cc:7445
SparseMatrix::cumsum
SparseMatrix cumsum(int dim=-1) const
Definition dSparse.cc:7337
SparseMatrix::SparseMatrix
SparseMatrix(octave_idx_type r, octave_idx_type c)
Definition dSparse.h:54
SparseMatrix::column
ColumnVector column(octave_idx_type i) const
Definition dSparse.cc:519
SparseMatrix::concat
SparseMatrix concat(const SparseMatrix &rb, const Array< octave_idx_type > &ra_idx)
Definition dSparse.cc:531
SparseMatrix::SparseMatrix
SparseMatrix(const MSparse< double > &a)
Definition dSparse.h:68
SparseMatrix::solve_singularity_handler
void(* solve_singularity_handler)(double rcond)
Definition dSparse.h:50
SparseMatrix::ipermute
SparseMatrix ipermute(const Array< octave_idx_type > &vec) const
Definition dSparse.cc:7451
SparseMatrix::dense_matrix_type
Matrix dense_matrix_type
Definition dSparse.h:48
SparseMatrix::operator!
SparseBoolMatrix operator!() const
Definition dSparse.cc:7283
SparseMatrix::sumsq
SparseMatrix sumsq(int dim=-1) const
Definition dSparse.cc:7361
SparseMatrix::SparseMatrix
SparseMatrix(const Array< double > &a, const octave::idx_vector &r, const octave::idx_vector &c, octave_idx_type nr=-1, octave_idx_type nc=-1, bool sum_terms=true, octave_idx_type nzm=-1)
Definition dSparse.h:78
SparseMatrix::hermitian
SparseMatrix hermitian() const
Definition dSparse.h:129
SparseMatrix::SparseMatrix
SparseMatrix(octave_idx_type r, octave_idx_type c, octave_idx_type num_nz)
Definition dSparse.h:88
SparseMatrix::SparseMatrix
SparseMatrix(const NDArray &a)
Definition dSparse.h:76
SparseMatrix::all_elements_are_int_or_inf_or_nan
bool all_elements_are_int_or_inf_or_nan() const
Definition dSparse.cc:7229
SparseMatrix::solve
Matrix solve(MatrixType &typ, const Matrix &b) const
Definition dSparse.cc:6661
SparseMatrix::SparseMatrix
SparseMatrix(const PermMatrix &a)
Definition dSparse.h:86
SparseMatrix::SparseMatrix
SparseMatrix(const Matrix &a)
Definition dSparse.h:74
SparseMatrix::imag
friend SparseMatrix imag(const SparseComplexMatrix &a)
Definition dSparse.cc:572
SparseMatrix::operator<<
friend std::ostream & operator<<(std::ostream &os, const SparseMatrix &a)
Definition dSparse.cc:7404
SparseMatrix::SparseMatrix
SparseMatrix(const dim_vector &dv, octave_idx_type nz=0)
Definition dSparse.h:57
SparseMatrix::operator=
SparseMatrix & operator=(const SparseMatrix &a)
Definition dSparse.h:91
SparseMatrix::max
SparseMatrix max(int dim=-1) const
Definition dSparse.cc:183
SparseMatrix::any_element_is_negative
bool any_element_is_negative(bool=false) const
Definition dSparse.cc:7151
SparseMatrix::matrix_value
Matrix matrix_value() const
Definition dSparse.cc:7398
SparseMatrix::operator>>
friend std::istream & operator>>(std::istream &is, SparseMatrix &a)
Definition dSparse.cc:7425
SparseMatrix::SparseMatrix
SparseMatrix(octave_idx_type r, octave_idx_type c, double val)
Definition dSparse.h:60
Sparse::index
Sparse< T, Alloc > index(const octave::idx_vector &i, bool resize_ok=false) const
Definition Sparse.cc:1433
base_det
Definition DET.h:38
dim_vector
Vector representing the dimensions (size) of an Array.
Definition dim-vector.h:90
octave_idx_type
operator+
SparseMatrix operator+(const DiagMatrix &, const SparseMatrix &)
Definition dSparse.cc:7503
trans_mul
Matrix trans_mul(const SparseMatrix &a, const Matrix &b)
Definition dSparse.cc:7483
operator*
SparseMatrix operator*(const SparseMatrix &a, const SparseMatrix &b)
Definition dSparse.cc:7459
real
SparseMatrix real(const SparseComplexMatrix &a)
Definition dSparse.cc:551
min
SparseMatrix min(double d, const SparseMatrix &m)
Definition dSparse.cc:7547
operator-
SparseMatrix operator-(const DiagMatrix &, const SparseMatrix &)
Definition dSparse.cc:7509
max
SparseMatrix max(double d, const SparseMatrix &m)
Definition dSparse.cc:7697
mul_trans
Matrix mul_trans(const Matrix &a, const SparseMatrix &b)
Definition dSparse.cc:7471
imag
SparseMatrix imag(const SparseComplexMatrix &a)
Definition dSparse.cc:572
d
F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
Definition lo-slatec-proto.h:39
OCTAVE_API
#define OCTAVE_API
Definition main.in.cc:55
ra_idx
const octave_base_value const Array< octave_idx_type > & ra_idx
Definition op-int-concat.cc:170