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