GNU Octave  9.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
dSparse.h
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 //
3 // Copyright (C) 1998-2024 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 
43 class
44 SparseMatrix : public MSparse<double>
45 {
46 public:
47 
48  // Corresponding dense matrix type for this sparse matrix type.
49  typedef Matrix dense_matrix_type;
50 
51  typedef void (*solve_singularity_handler) (double rcond);
52 
53  SparseMatrix () : MSparse<double> () { }
54 
55  SparseMatrix (octave_idx_type r, octave_idx_type c)
56  : MSparse<double> (r, c) { }
57 
58  SparseMatrix (const dim_vector& dv, octave_idx_type nz = 0) :
59  MSparse<double> (dv, nz) { }
60 
61  explicit SparseMatrix (octave_idx_type r, octave_idx_type c, double val)
62  : MSparse<double> (r, c, val) { }
63 
64  SparseMatrix (const SparseMatrix& a) : MSparse<double> (a) { }
65 
66  SparseMatrix (const SparseMatrix& a, const dim_vector& dv)
67  : MSparse<double> (a, dv) { }
68 
69  SparseMatrix (const MSparse<double>& a) : MSparse<double> (a) { }
70 
71  SparseMatrix (const Sparse<double>& a) : MSparse<double> (a) { }
72 
73  explicit OCTAVE_API SparseMatrix (const SparseBoolMatrix& a);
74 
75  explicit SparseMatrix (const Matrix& a) : MSparse<double> (a) { }
76 
77  explicit SparseMatrix (const NDArray& a) : MSparse<double> (a) { }
78 
79  SparseMatrix (const Array<double>& a, const octave::idx_vector& r,
80  const octave::idx_vector& c, octave_idx_type nr = -1,
81  octave_idx_type nc = -1, bool sum_terms = true,
82  octave_idx_type nzm = -1)
83  : MSparse<double> (a, r, c, nr, nc, sum_terms, nzm) { }
84 
85  explicit OCTAVE_API SparseMatrix (const DiagMatrix& a);
86 
87  explicit SparseMatrix (const PermMatrix& a) : MSparse<double> (a) { }
88 
89  SparseMatrix (octave_idx_type r, octave_idx_type c,
90  octave_idx_type num_nz) : MSparse<double> (r, c, num_nz) { }
91 
92  SparseMatrix& operator = (const SparseMatrix& a)
93  {
94  MSparse<double>::operator = (a);
95  return *this;
96  }
97 
98  ~SparseMatrix () = default;
99 
100  OCTAVE_API bool operator == (const SparseMatrix& a) const;
101  OCTAVE_API bool operator != (const SparseMatrix& a) const;
102 
103  OCTAVE_API bool issymmetric () const;
104 
105  OCTAVE_API SparseMatrix max (int dim = -1) const;
106  OCTAVE_API SparseMatrix max (Array<octave_idx_type>& index, int dim = -1) const;
107  OCTAVE_API SparseMatrix min (int dim = -1) const;
108  OCTAVE_API SparseMatrix min (Array<octave_idx_type>& index, int dim = -1) const;
109 
110  // destructive insert/delete/reorder operations
111 
112  OCTAVE_API SparseMatrix&
113  insert (const SparseMatrix& a, octave_idx_type r, octave_idx_type c);
114 
115  OCTAVE_API SparseMatrix&
116  insert (const SparseMatrix& a, const Array<octave_idx_type>& indx);
117 
118  OCTAVE_API SparseMatrix
119  concat (const SparseMatrix& rb, const Array<octave_idx_type>& ra_idx);
120  OCTAVE_API SparseComplexMatrix
121  concat (const SparseComplexMatrix& rb, const Array<octave_idx_type>& ra_idx);
122 
123  friend OCTAVE_API SparseMatrix real (const SparseComplexMatrix& a);
124  friend OCTAVE_API SparseMatrix imag (const SparseComplexMatrix& a);
125 
126  SparseMatrix transpose () const
127  {
128  return MSparse<double>::transpose ();
129  }
130  SparseMatrix hermitian () const { return transpose (); }
131 
132  // extract row or column i.
133 
134  OCTAVE_API RowVector row (octave_idx_type i) const;
135 
136  OCTAVE_API ColumnVector column (octave_idx_type i) const;
137 
138 private:
139  OCTAVE_API SparseMatrix
140  dinverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
141  const bool force = false, const bool calccond = true) const;
142 
143  OCTAVE_API SparseMatrix
144  tinverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
145  const bool force = false, const bool calccond = true) const;
146 
147 public:
148  OCTAVE_API SparseMatrix inverse () const;
149  OCTAVE_API SparseMatrix inverse (MatrixType& mattype) const;
150  OCTAVE_API SparseMatrix
151  inverse (MatrixType& mattype, octave_idx_type& info) const;
152  OCTAVE_API SparseMatrix
153  inverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
154  bool force = false, bool calc_cond = true) const;
155 
156  OCTAVE_API DET determinant () const;
157  OCTAVE_API DET determinant (octave_idx_type& info) const;
158  OCTAVE_API DET determinant (octave_idx_type& info, double& rcond,
159  bool calc_cond = true) const;
160 
161 private:
162  // Diagonal matrix solvers
163  OCTAVE_API Matrix
164  dsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
165  double& rcond, solve_singularity_handler sing_handler,
166  bool calc_cond = false) const;
167 
168  OCTAVE_API ComplexMatrix
169  dsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
170  double& rcond, solve_singularity_handler sing_handler,
171  bool calc_cond = false) const;
172 
173  OCTAVE_API SparseMatrix
174  dsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
175  double& rcond, solve_singularity_handler sing_handler,
176  bool calc_cond = false) const;
177 
178  OCTAVE_API SparseComplexMatrix
179  dsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
180  double& rcond, solve_singularity_handler sing_handler,
181  bool calc_cond = false) const;
182 
183  // Upper triangular matrix solvers
184  OCTAVE_API Matrix
185  utsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
186  double& rcond, solve_singularity_handler sing_handler,
187  bool calc_cond = false) const;
188 
189  OCTAVE_API ComplexMatrix
190  utsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
191  double& rcond, solve_singularity_handler sing_handler,
192  bool calc_cond = false) const;
193 
194  OCTAVE_API SparseMatrix
195  utsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
196  double& rcond, solve_singularity_handler sing_handler,
197  bool calc_cond = false) const;
198 
199  OCTAVE_API SparseComplexMatrix
200  utsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
201  double& rcond, solve_singularity_handler sing_handler,
202  bool calc_cond = false) const;
203 
204  // Lower triangular matrix solvers
205  OCTAVE_API Matrix
206  ltsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
207  double& rcond, solve_singularity_handler sing_handler,
208  bool calc_cond = false) const;
209 
210  OCTAVE_API ComplexMatrix
211  ltsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
212  double& rcond, solve_singularity_handler sing_handler,
213  bool calc_cond = false) const;
214 
215  OCTAVE_API SparseMatrix
216  ltsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
217  double& rcond, solve_singularity_handler sing_handler,
218  bool calc_cond = false) const;
219 
220  OCTAVE_API SparseComplexMatrix
221  ltsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
222  double& rcond, solve_singularity_handler sing_handler,
223  bool calc_cond = false) const;
224 
225  // Tridiagonal matrix solvers
226  OCTAVE_API Matrix
227  trisolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
228  double& rcond, solve_singularity_handler sing_handler,
229  bool calc_cond = false) const;
230 
231  OCTAVE_API ComplexMatrix
232  trisolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
233  double& rcond, solve_singularity_handler sing_handler,
234  bool calc_cond = false) const;
235 
236  OCTAVE_API SparseMatrix
237  trisolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
238  double& rcond, solve_singularity_handler sing_handler,
239  bool calc_cond = false) const;
240 
241  OCTAVE_API SparseComplexMatrix
242  trisolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
243  double& rcond, solve_singularity_handler sing_handler,
244  bool calc_cond = false) const;
245 
246  // Banded matrix solvers (umfpack/cholesky)
247  OCTAVE_API Matrix
248  bsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
249  double& rcond, solve_singularity_handler sing_handler,
250  bool calc_cond = false) const;
251 
252  OCTAVE_API ComplexMatrix
253  bsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
254  double& rcond, solve_singularity_handler sing_handler,
255  bool calc_cond = false) const;
256 
257  OCTAVE_API SparseMatrix
258  bsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
259  double& rcond, solve_singularity_handler sing_handler,
260  bool calc_cond = false) const;
261 
262  OCTAVE_API SparseComplexMatrix
263  bsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
264  double& rcond, solve_singularity_handler sing_handler,
265  bool calc_cond = false) const;
266 
267  // Full matrix solvers (umfpack/cholesky)
268  OCTAVE_API void *
269  factorize (octave_idx_type& err, double& rcond, Matrix& Control,
270  Matrix& Info, solve_singularity_handler sing_handler,
271  bool calc_cond = false) const;
272 
273  OCTAVE_API Matrix
274  fsolve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
275  double& rcond, solve_singularity_handler sing_handler,
276  bool calc_cond = false) const;
277 
278  OCTAVE_API ComplexMatrix
279  fsolve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
280  double& rcond, solve_singularity_handler sing_handler,
281  bool calc_cond = false) const;
282 
283  OCTAVE_API SparseMatrix
284  fsolve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
285  double& rcond, solve_singularity_handler sing_handler,
286  bool calc_cond = false) const;
287 
288  OCTAVE_API SparseComplexMatrix
289  fsolve (MatrixType& typ, const SparseComplexMatrix& b, octave_idx_type& info,
290  double& rcond, solve_singularity_handler sing_handler,
291  bool calc_cond = false) const;
292 
293 public:
294  // Generic interface to solver with no probing of type
295  OCTAVE_API Matrix solve (MatrixType& typ, const Matrix& b) const;
296  OCTAVE_API Matrix
297  solve (MatrixType& typ, const Matrix& b, octave_idx_type& info) const;
298  OCTAVE_API Matrix
299  solve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
300  double& rcond) const;
301  OCTAVE_API Matrix
302  solve (MatrixType& typ, const Matrix& b, octave_idx_type& info,
303  double& rcond, solve_singularity_handler sing_handler,
304  bool singular_fallback = true) const;
305 
306  OCTAVE_API ComplexMatrix
307  solve (MatrixType& typ, const ComplexMatrix& b) const;
308  OCTAVE_API ComplexMatrix
309  solve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info) const;
310  OCTAVE_API ComplexMatrix
311  solve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
312  double& rcond) const;
313  OCTAVE_API ComplexMatrix
314  solve (MatrixType& typ, const ComplexMatrix& b, octave_idx_type& info,
315  double& rcond, solve_singularity_handler sing_handler,
316  bool singular_fallback = true) const;
317 
318  OCTAVE_API SparseMatrix solve (MatrixType& typ, const SparseMatrix& b) const;
319  OCTAVE_API SparseMatrix
320  solve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info) const;
321  OCTAVE_API SparseMatrix
322  solve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
323  double& rcond) const;
324  OCTAVE_API SparseMatrix
325  solve (MatrixType& typ, const SparseMatrix& b, octave_idx_type& info,
326  double& rcond, solve_singularity_handler sing_handler,
327  bool singular_fallback = true) const;
328 
329  OCTAVE_API SparseComplexMatrix
330  solve (MatrixType& typ, const SparseComplexMatrix& b) const;
331  OCTAVE_API SparseComplexMatrix
332  solve (MatrixType& typ, const SparseComplexMatrix& b,
333  octave_idx_type& info) const;
334  OCTAVE_API SparseComplexMatrix
335  solve (MatrixType& typ, const SparseComplexMatrix& b,
336  octave_idx_type& info, double& rcond) const;
337  OCTAVE_API SparseComplexMatrix
338  solve (MatrixType& typ, const SparseComplexMatrix& b,
339  octave_idx_type& info, double& rcond,
340  solve_singularity_handler sing_handler,
341  bool singular_fallabck = true) const;
342 
343  OCTAVE_API ColumnVector solve (MatrixType& typ, const ColumnVector& b) const;
344  OCTAVE_API ColumnVector
345  solve (MatrixType& typ, const ColumnVector& b, octave_idx_type& info) const;
346  OCTAVE_API ColumnVector
347  solve (MatrixType& typ, const ColumnVector& b, octave_idx_type& info,
348  double& rcond) const;
349  OCTAVE_API ColumnVector
350  solve (MatrixType& typ, const ColumnVector& b, octave_idx_type& info,
351  double& rcond, solve_singularity_handler sing_handler) const;
352 
353  OCTAVE_API ComplexColumnVector
354  solve (MatrixType& typ, const ComplexColumnVector& b) const;
355  OCTAVE_API ComplexColumnVector
356  solve (MatrixType& typ, const ComplexColumnVector& b,
357  octave_idx_type& info) const;
358  OCTAVE_API ComplexColumnVector
359  solve (MatrixType& typ, const ComplexColumnVector& b,
360  octave_idx_type& info, double& rcond) const;
361  OCTAVE_API ComplexColumnVector
362  solve (MatrixType& typ, const ComplexColumnVector& b,
363  octave_idx_type& info, double& rcond,
364  solve_singularity_handler sing_handler) const;
365 
366  // Generic interface to solver with probing of type
367  OCTAVE_API Matrix solve (const Matrix& b) const;
368  OCTAVE_API Matrix solve (const Matrix& b, octave_idx_type& info) const;
369  OCTAVE_API Matrix
370  solve (const Matrix& b, octave_idx_type& info, double& rcond) const;
371  OCTAVE_API Matrix
372  solve (const Matrix& b, octave_idx_type& info, double& rcond,
373  solve_singularity_handler sing_handler) const;
374 
375  OCTAVE_API ComplexMatrix solve (const ComplexMatrix& b) const;
376  OCTAVE_API ComplexMatrix
377  solve (const ComplexMatrix& b, octave_idx_type& info) const;
378  OCTAVE_API ComplexMatrix
379  solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond) const;
380  OCTAVE_API ComplexMatrix
381  solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond,
382  solve_singularity_handler sing_handler) const;
383 
384  OCTAVE_API SparseMatrix solve (const SparseMatrix& b) const;
385  OCTAVE_API SparseMatrix
386  solve (const SparseMatrix& b, octave_idx_type& info) const;
387  OCTAVE_API SparseMatrix
388  solve (const SparseMatrix& b, octave_idx_type& info, double& rcond) const;
389  OCTAVE_API SparseMatrix
390  solve (const SparseMatrix& b, octave_idx_type& info, double& rcond,
391  solve_singularity_handler sing_handler) const;
392 
393  OCTAVE_API SparseComplexMatrix solve (const SparseComplexMatrix& b) const;
394  OCTAVE_API SparseComplexMatrix
395  solve (const SparseComplexMatrix& b, octave_idx_type& info) const;
396  OCTAVE_API SparseComplexMatrix
397  solve (const SparseComplexMatrix& b, octave_idx_type& info,
398  double& rcond) const;
399  OCTAVE_API SparseComplexMatrix
400  solve (const SparseComplexMatrix& b, octave_idx_type& info,
401  double& rcond, solve_singularity_handler sing_handler) const;
402 
403  OCTAVE_API ColumnVector solve (const ColumnVector& b) const;
404  OCTAVE_API ColumnVector
405  solve (const ColumnVector& b, octave_idx_type& info) const;
406  OCTAVE_API ColumnVector
407  solve (const ColumnVector& b, octave_idx_type& info, double& rcond) const;
408  OCTAVE_API ColumnVector
409  solve (const ColumnVector& b, octave_idx_type& info, double& rcond,
410  solve_singularity_handler sing_handler) const;
411 
412  OCTAVE_API ComplexColumnVector solve (const ComplexColumnVector& b) const;
413  OCTAVE_API ComplexColumnVector
414  solve (const ComplexColumnVector& b, octave_idx_type& info) const;
415  OCTAVE_API ComplexColumnVector
416  solve (const ComplexColumnVector& b, octave_idx_type& info,
417  double& rcond) const;
418  OCTAVE_API ComplexColumnVector
419  solve (const ComplexColumnVector& b, octave_idx_type& info,
420  double& rcond, solve_singularity_handler sing_handler) const;
421 
422  // other operations
423 
424  OCTAVE_API bool any_element_is_negative (bool = false) const;
425  OCTAVE_API bool any_element_is_nan () const;
426  OCTAVE_API bool any_element_is_inf_or_nan () const;
427  OCTAVE_API bool any_element_not_one_or_zero () const;
428  OCTAVE_API bool all_elements_are_zero () const;
429  OCTAVE_API bool all_elements_are_int_or_inf_or_nan () const;
430  OCTAVE_API bool all_integers (double& max_val, double& min_val) const;
431  OCTAVE_API bool too_large_for_float () const;
432 
433  OCTAVE_API SparseBoolMatrix operator ! () const;
434 
435  OCTAVE_API SparseBoolMatrix all (int dim = -1) const;
436  OCTAVE_API SparseBoolMatrix any (int dim = -1) const;
437 
438  OCTAVE_API SparseMatrix cumprod (int dim = -1) const;
439  OCTAVE_API SparseMatrix cumsum (int dim = -1) const;
440  OCTAVE_API SparseMatrix prod (int dim = -1) const;
441  OCTAVE_API SparseMatrix sum (int dim = -1) const;
442  OCTAVE_API SparseMatrix sumsq (int dim = -1) const;
443  OCTAVE_API SparseMatrix abs () const;
444 
445  OCTAVE_API SparseMatrix diag (octave_idx_type k = 0) const;
446 
447  OCTAVE_API Matrix matrix_value () const;
448 
449  OCTAVE_API SparseMatrix squeeze () const;
450 
451  OCTAVE_API SparseMatrix reshape (const dim_vector& new_dims) const;
452 
453  OCTAVE_API SparseMatrix
454  permute (const Array<octave_idx_type>& vec, bool inv = false) const;
455 
456  OCTAVE_API SparseMatrix ipermute (const Array<octave_idx_type>& vec) const;
457 
458  // i/o
459 
460  friend OCTAVE_API std::ostream& operator << (std::ostream& os,
461  const SparseMatrix& a);
462  friend OCTAVE_API std::istream& operator >> (std::istream& is,
463  SparseMatrix& a);
464 
465 };
466 
467 // Publish externally used friend functions.
468 
469 extern OCTAVE_API SparseMatrix real (const SparseComplexMatrix& a);
470 extern OCTAVE_API SparseMatrix imag (const SparseComplexMatrix& a);
471 
472 // Other operators.
473 
474 extern OCTAVE_API SparseMatrix operator * (const SparseMatrix& a,
475  const SparseMatrix& b);
476 extern OCTAVE_API Matrix operator * (const Matrix& a,
477  const SparseMatrix& b);
478 extern OCTAVE_API Matrix mul_trans (const Matrix& a,
479  const SparseMatrix& b);
480 extern OCTAVE_API Matrix operator * (const SparseMatrix& a,
481  const Matrix& b);
482 extern OCTAVE_API Matrix trans_mul (const SparseMatrix& a,
483  const Matrix& b);
484 
485 extern OCTAVE_API SparseMatrix operator * (const DiagMatrix&,
486  const SparseMatrix&);
487 extern OCTAVE_API SparseMatrix operator * (const SparseMatrix&,
488  const DiagMatrix&);
489 
490 extern OCTAVE_API SparseMatrix operator + (const DiagMatrix&,
491  const SparseMatrix&);
492 extern OCTAVE_API SparseMatrix operator + (const SparseMatrix&,
493  const DiagMatrix&);
494 extern OCTAVE_API SparseMatrix operator - (const DiagMatrix&,
495  const SparseMatrix&);
496 extern OCTAVE_API SparseMatrix operator - (const SparseMatrix&,
497  const DiagMatrix&);
498 
499 extern OCTAVE_API SparseMatrix operator * (const PermMatrix&,
500  const SparseMatrix&);
501 extern OCTAVE_API SparseMatrix operator * (const SparseMatrix&,
502  const PermMatrix&);
503 
504 extern OCTAVE_API SparseMatrix min (double d, const SparseMatrix& m);
505 extern OCTAVE_API SparseMatrix min (const SparseMatrix& m, double d);
506 extern OCTAVE_API SparseMatrix min (const SparseMatrix& a,
507  const SparseMatrix& b);
508 
509 extern OCTAVE_API SparseMatrix max (double d, const SparseMatrix& m);
510 extern OCTAVE_API SparseMatrix max (const SparseMatrix& m, double d);
511 extern OCTAVE_API SparseMatrix max (const SparseMatrix& a,
512  const SparseMatrix& b);
513 
514 SPARSE_SMS_CMP_OP_DECLS (SparseMatrix, double, OCTAVE_API)
515 SPARSE_SMS_BOOL_OP_DECLS (SparseMatrix, double, OCTAVE_API)
516 
517 SPARSE_SSM_CMP_OP_DECLS (double, SparseMatrix, OCTAVE_API)
518 SPARSE_SSM_BOOL_OP_DECLS (double, SparseMatrix, OCTAVE_API)
519 
520 SPARSE_SMSM_CMP_OP_DECLS (SparseMatrix, SparseMatrix, OCTAVE_API)
521 SPARSE_SMSM_BOOL_OP_DECLS (SparseMatrix, SparseMatrix, OCTAVE_API)
522 
523 SPARSE_FORWARD_DEFS (MSparse, SparseMatrix, Matrix, double)
524 
525 #endif
operator<<
template std::ostream & operator<<(std::ostream &, const Array< bool > &)
concat
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:418
SPARSE_FORWARD_DEFS
#define SPARSE_FORWARD_DEFS(B, R, F, T)
Definition: MSparse.h:194
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
operator>>
std::istream & operator>>(std::istream &is, SparseBoolMatrix &a)
Definition: boolSparse.cc:279
MSparse::transpose
MSparse< T > transpose() const
Definition: MSparse.h:98
MSparse::operator=
MSparse< T > & operator=(const MSparse< T > &a)
Definition: MSparse.h:80
Matrix
Definition: dMatrix.h:42
SparseMatrix::transpose
SparseMatrix transpose() const
Definition: dSparse.h:126
SparseMatrix::SparseMatrix
SparseMatrix(const Sparse< double > &a)
Definition: dSparse.h:71
SparseMatrix::solve
ComplexMatrix solve(const ComplexMatrix &b) const
SparseMatrix::SparseMatrix
SparseMatrix(const SparseMatrix &a, const dim_vector &dv)
Definition: dSparse.h:66
SparseMatrix::~SparseMatrix
~SparseMatrix()=default
SparseMatrix::SparseMatrix
SparseMatrix(const SparseMatrix &a)
Definition: dSparse.h:64
SparseMatrix::SparseMatrix
SparseMatrix(octave_idx_type r, octave_idx_type c)
Definition: dSparse.h:55
SparseMatrix::SparseMatrix
SparseMatrix(const MSparse< double > &a)
Definition: dSparse.h:69
SparseMatrix::dense_matrix_type
Matrix dense_matrix_type
Definition: dSparse.h:49
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:79
SparseMatrix::hermitian
SparseMatrix hermitian() const
Definition: dSparse.h:130
SparseMatrix::SparseMatrix
SparseMatrix(octave_idx_type r, octave_idx_type c, octave_idx_type num_nz)
Definition: dSparse.h:89
SparseMatrix::SparseMatrix
SparseMatrix()
Definition: dSparse.h:53
SparseMatrix::SparseMatrix
SparseMatrix(const NDArray &a)
Definition: dSparse.h:77
SparseMatrix::SparseMatrix
SparseMatrix(const PermMatrix &a)
Definition: dSparse.h:87
SparseMatrix::SparseMatrix
SparseMatrix(const Matrix &a)
Definition: dSparse.h:75
SparseMatrix::SparseMatrix
SparseMatrix(const dim_vector &dv, octave_idx_type nz=0)
Definition: dSparse.h:58
SparseMatrix::SparseMatrix
SparseMatrix(octave_idx_type r, octave_idx_type c, double val)
Definition: dSparse.h:61
Sparse
Definition: Sparse.h:49
base_det
Definition: DET.h:39
dim_vector
Vector representing the dimensions (size) of an Array.
Definition: dim-vector.h:94
octave_idx_type
operator+
SparseMatrix operator+(const DiagMatrix &, const SparseMatrix &)
Definition: dSparse.cc:7555
trans_mul
Matrix trans_mul(const SparseMatrix &a, const Matrix &b)
Definition: dSparse.cc:7535
operator*
SparseMatrix operator*(const SparseMatrix &a, const SparseMatrix &b)
Definition: dSparse.cc:7511
real
SparseMatrix real(const SparseComplexMatrix &a)
Definition: dSparse.cc:554
min
SparseMatrix min(double d, const SparseMatrix &m)
Definition: dSparse.cc:7599
operator-
SparseMatrix operator-(const DiagMatrix &, const SparseMatrix &)
Definition: dSparse.cc:7561
max
SparseMatrix max(double d, const SparseMatrix &m)
Definition: dSparse.cc:7749
mul_trans
Matrix mul_trans(const Matrix &a, const SparseMatrix &b)
Definition: dSparse.cc:7523
imag
SparseMatrix imag(const SparseComplexMatrix &a)
Definition: dSparse.cc:575
operator!=
bool operator!=(const dim_vector &a, const dim_vector &b)
Definition: dim-vector.h:536
operator==
bool operator==(const dim_vector &a, const dim_vector &b)
Definition: dim-vector.h:520
idx_vector
octave::idx_vector idx_vector
Definition: idx-vector.h:1022
d
F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
Definition: lo-slatec-proto.h:39
too_large_for_float
bool too_large_for_float(double x)
Definition: lo-utils.cc:56
OCTAVE_API
#define OCTAVE_API
Definition: main.cc:55
m
T octave_idx_type m
Definition: mx-inlines.cc:781
r
T * r
Definition: mx-inlines.cc:781
ra_idx
const octave_base_value const Array< octave_idx_type > & ra_idx
Definition: op-int-concat.cc:170
operator!
octave_value operator!(const octave_value &a)
Definition: ov.h:1633
abs
template int8_t abs(int8_t)