GNU Octave  9.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
CSparse.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_CSparse_h)
27 #define octave_CSparse_h 1
28 
29 #include "octave-config.h"
30 
31 #include "CColVector.h"
32 #include "CMatrix.h"
33 #include "CNDArray.h"
34 #include "DET.h"
35 #include "MSparse.h"
36 #include "MatrixType.h"
37 #include "Sparse-op-decls.h"
38 #include "dColVector.h"
39 #include "dMatrix.h"
40 #include "dNDArray.h"
41 #include "oct-cmplx.h"
42 #include "mx-fwd.h"
43 
44 class
45 OCTAVE_API
46 SparseComplexMatrix : public MSparse<Complex>
47 {
48 public:
49 
50  // Corresponding dense matrix type for this sparse matrix type.
51  typedef ComplexMatrix dense_matrix_type;
52 
53  typedef void (*solve_singularity_handler) (double rcond);
54 
55  SparseComplexMatrix () : MSparse<Complex> () { }
56 
57  SparseComplexMatrix (octave_idx_type r,
58  octave_idx_type c) : MSparse<Complex> (r, c) { }
59 
60  SparseComplexMatrix (const dim_vector& dv, octave_idx_type nz = 0)
61  : MSparse<Complex> (dv, nz) { }
62 
63  explicit SparseComplexMatrix (octave_idx_type r, octave_idx_type c,
64  Complex val)
65  : MSparse<Complex> (r, c, val) { }
66 
67  SparseComplexMatrix (octave_idx_type r, octave_idx_type c, double val)
68  : MSparse<Complex> (r, c, Complex (val)) { }
69 
70  SparseComplexMatrix (const SparseComplexMatrix& a)
71  : MSparse<Complex> (a) { }
72 
73  SparseComplexMatrix (const SparseComplexMatrix& a, const dim_vector& dv)
74  : MSparse<Complex> (a, dv) { }
75 
76  SparseComplexMatrix (const MSparse<Complex>& a) : MSparse<Complex> (a) { }
77 
78  SparseComplexMatrix (const Sparse<Complex>& a) : MSparse<Complex> (a) { }
79 
80  explicit SparseComplexMatrix (const ComplexMatrix& a)
81  : MSparse<Complex> (a) { }
82 
83  explicit SparseComplexMatrix (const ComplexNDArray& a)
84  : MSparse<Complex> (a) { }
85 
86  SparseComplexMatrix (const Array<Complex>& a, const octave::idx_vector& r,
87  const octave::idx_vector& c, octave_idx_type nr = -1,
88  octave_idx_type nc = -1, bool sum_terms = true,
89  octave_idx_type nzm = -1)
90  : MSparse<Complex> (a, r, c, nr, nc, sum_terms, nzm) { }
91 
92  explicit OCTAVE_API SparseComplexMatrix (const SparseMatrix& a);
93 
94  explicit OCTAVE_API SparseComplexMatrix (const SparseBoolMatrix& a);
95 
96  explicit OCTAVE_API SparseComplexMatrix (const ComplexDiagMatrix& a);
97 
98  SparseComplexMatrix (octave_idx_type r, octave_idx_type c,
99  octave_idx_type num_nz)
100  : MSparse<Complex> (r, c, num_nz) { }
101 
102  SparseComplexMatrix& operator = (const SparseComplexMatrix& a)
103  {
104  MSparse<Complex>::operator = (a);
105  return *this;
106  }
107 
108  ~SparseComplexMatrix () = default;
109 
110  OCTAVE_API bool operator == (const SparseComplexMatrix& a) const;
111  OCTAVE_API bool operator != (const SparseComplexMatrix& a) const;
112 
113  OCTAVE_API bool ishermitian () const;
114 
115  OCTAVE_API SparseComplexMatrix max (int dim = -1) const;
116  OCTAVE_API SparseComplexMatrix
117  max (Array<octave_idx_type>& index, int dim = -1) const;
118  OCTAVE_API SparseComplexMatrix min (int dim = -1) const;
119  OCTAVE_API SparseComplexMatrix
120  min (Array<octave_idx_type>& index, int dim = -1) const;
121 
122  OCTAVE_API SparseComplexMatrix&
123  insert (const SparseComplexMatrix& a, octave_idx_type r, octave_idx_type c);
124  OCTAVE_API SparseComplexMatrix&
125  insert (const SparseMatrix& a, octave_idx_type r, octave_idx_type c);
126  OCTAVE_API SparseComplexMatrix&
127  insert (const SparseComplexMatrix& a, const Array<octave_idx_type>& indx);
128  OCTAVE_API SparseComplexMatrix&
129  insert (const SparseMatrix& a, const Array<octave_idx_type>& indx);
130 
131  OCTAVE_API SparseComplexMatrix
132  concat (const SparseComplexMatrix& rb, const Array<octave_idx_type>& ra_idx);
133  OCTAVE_API SparseComplexMatrix
134  concat (const SparseMatrix& rb, const Array<octave_idx_type>& ra_idx);
135 
136  OCTAVE_API ComplexMatrix matrix_value () const;
137 
138  OCTAVE_API SparseComplexMatrix hermitian () const; // complex conjugate transpose
139  SparseComplexMatrix transpose () const
140  { return MSparse<Complex>::transpose (); }
141 
142  friend OCTAVE_API SparseComplexMatrix conj (const SparseComplexMatrix& a);
143 
144  // extract row or column i.
145 
146  OCTAVE_API ComplexRowVector row (octave_idx_type i) const;
147 
148  OCTAVE_API ComplexColumnVector column (octave_idx_type i) const;
149 
150 private:
151  SparseComplexMatrix dinverse (MatrixType& mattype, octave_idx_type& info,
152  double& rcond, const bool force = false,
153  const bool calccond = true) const;
154 
155  SparseComplexMatrix tinverse (MatrixType& mattype, octave_idx_type& info,
156  double& rcond, const bool force = false,
157  const bool calccond = true) const;
158 
159 public:
160  OCTAVE_API SparseComplexMatrix inverse () const;
161  OCTAVE_API SparseComplexMatrix inverse (MatrixType& mattype) const;
162  OCTAVE_API SparseComplexMatrix
163  inverse (MatrixType& mattype, octave_idx_type& info) const;
164  OCTAVE_API SparseComplexMatrix
165  inverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
166  bool force = false, bool calc_cond = true) const;
167 
168  OCTAVE_API ComplexDET determinant () const;
169  OCTAVE_API ComplexDET determinant (octave_idx_type& info) const;
170  OCTAVE_API ComplexDET
171  determinant (octave_idx_type& info, double& rcond,
172  bool calc_cond = true) const;
173 
174 private:
175  // Diagonal matrix solvers
176  ComplexMatrix dsolve (MatrixType& mattype, const Matrix& b,
177  octave_idx_type& info, double& rcond,
178  solve_singularity_handler sing_handler,
179  bool calc_cond = false) const;
180 
181  ComplexMatrix dsolve (MatrixType& mattype, const ComplexMatrix& b,
182  octave_idx_type& info, double& rcond,
183  solve_singularity_handler sing_handler,
184  bool calc_cond = false) const;
185 
186  SparseComplexMatrix dsolve (MatrixType& mattype, const SparseMatrix& b,
187  octave_idx_type& info, double& rcond,
188  solve_singularity_handler sing_handler,
189  bool calc_cond = false) const;
190 
191  SparseComplexMatrix dsolve (MatrixType& mattype, const SparseComplexMatrix& b,
192  octave_idx_type& info, double& rcond,
193  solve_singularity_handler sing_handler,
194  bool calc_cond = false) const;
195 
196  // Upper triangular matrix solvers
197  ComplexMatrix utsolve (MatrixType& mattype, const Matrix& b,
198  octave_idx_type& info, double& rcond,
199  solve_singularity_handler sing_handler,
200  bool calc_cond = false) const;
201 
202  ComplexMatrix utsolve (MatrixType& mattype, const ComplexMatrix& b,
203  octave_idx_type& info, double& rcond,
204  solve_singularity_handler sing_handler,
205  bool calc_cond = false) const;
206 
207  SparseComplexMatrix utsolve (MatrixType& mattype, const SparseMatrix& b,
208  octave_idx_type& info, double& rcond,
209  solve_singularity_handler sing_handler,
210  bool calc_cond = false) const;
211 
212  SparseComplexMatrix utsolve (MatrixType& mattype, const SparseComplexMatrix& b,
213  octave_idx_type& info, double& rcond,
214  solve_singularity_handler sing_handler,
215  bool calc_cond = false) const;
216 
217  // Lower triangular matrix solvers
218  ComplexMatrix ltsolve (MatrixType& mattype, const Matrix& b,
219  octave_idx_type& info, double& rcond,
220  solve_singularity_handler sing_handler,
221  bool calc_cond = false) const;
222 
223  ComplexMatrix ltsolve (MatrixType& mattype, const ComplexMatrix& b,
224  octave_idx_type& info, double& rcond,
225  solve_singularity_handler sing_handler,
226  bool calc_cond = false) const;
227 
228  SparseComplexMatrix ltsolve (MatrixType& mattype, const SparseMatrix& b,
229  octave_idx_type& info, double& rcond,
230  solve_singularity_handler sing_handler,
231  bool calc_cond = false) const;
232 
233  SparseComplexMatrix ltsolve (MatrixType& mattype, const SparseComplexMatrix& b,
234  octave_idx_type& info, double& rcond,
235  solve_singularity_handler sing_handler,
236  bool calc_cond = false) const;
237 
238  // Tridiagonal matrix solvers
239  ComplexMatrix trisolve (MatrixType& mattype, const Matrix& b,
240  octave_idx_type& info, double& rcond,
241  solve_singularity_handler sing_handler,
242  bool calc_cond = false) const;
243 
244  ComplexMatrix trisolve (MatrixType& mattype, const ComplexMatrix& b,
245  octave_idx_type& info, double& rcond,
246  solve_singularity_handler sing_handler,
247  bool calc_cond = false) const;
248 
249  SparseComplexMatrix trisolve (MatrixType& mattype, const SparseMatrix& b,
250  octave_idx_type& info, double& rcond,
251  solve_singularity_handler sing_handler,
252  bool calc_cond = false) const;
253 
254  SparseComplexMatrix trisolve (MatrixType& mattype, const SparseComplexMatrix& b,
255  octave_idx_type& info, double& rcond,
256  solve_singularity_handler sing_handler,
257  bool calc_cond = false) const;
258 
259  // Banded matrix solvers (umfpack/cholesky)
260  ComplexMatrix bsolve (MatrixType& mattype, const Matrix& b,
261  octave_idx_type& info, double& rcond,
262  solve_singularity_handler sing_handler,
263  bool calc_cond = false) const;
264 
265  ComplexMatrix bsolve (MatrixType& mattype, const ComplexMatrix& b,
266  octave_idx_type& info, double& rcond,
267  solve_singularity_handler sing_handler,
268  bool calc_cond = false) const;
269 
270  SparseComplexMatrix bsolve (MatrixType& mattype, const SparseMatrix& b,
271  octave_idx_type& info, double& rcond,
272  solve_singularity_handler sing_handler,
273  bool calc_cond = false) const;
274 
275  SparseComplexMatrix bsolve (MatrixType& mattype, const SparseComplexMatrix& b,
276  octave_idx_type& info, double& rcond,
277  solve_singularity_handler sing_handler,
278  bool calc_cond = false) const;
279 
280  // Full matrix solvers (umfpack/cholesky)
281  void * factorize (octave_idx_type& err, double& rcond, Matrix& Control,
282  Matrix& Info, solve_singularity_handler sing_handler,
283  bool calc_cond) const;
284 
285  ComplexMatrix fsolve (MatrixType& mattype, const Matrix& b,
286  octave_idx_type& info, double& rcond,
287  solve_singularity_handler sing_handler,
288  bool calc_cond = false) const;
289 
290  ComplexMatrix fsolve (MatrixType& mattype, const ComplexMatrix& b,
291  octave_idx_type& info, double& rcond,
292  solve_singularity_handler sing_handler,
293  bool calc_cond = false) const;
294 
295  SparseComplexMatrix fsolve (MatrixType& mattype, const SparseMatrix& b,
296  octave_idx_type& info, double& rcond,
297  solve_singularity_handler sing_handler,
298  bool calc_cond = false) const;
299 
300  SparseComplexMatrix fsolve (MatrixType& mattype, const SparseComplexMatrix& b,
301  octave_idx_type& info, double& rcond,
302  solve_singularity_handler sing_handler,
303  bool calc_cond = false) const;
304 
305 public:
306  // Generic interface to solver with no probing of type
307  OCTAVE_API ComplexMatrix solve (MatrixType& mattype, const Matrix& b) const;
308  OCTAVE_API ComplexMatrix
309  solve (MatrixType& mattype, const Matrix& b, octave_idx_type& info) const;
310  OCTAVE_API ComplexMatrix
311  solve (MatrixType& mattype, const Matrix& b, octave_idx_type& info,
312  double& rcond) const;
313  OCTAVE_API ComplexMatrix
314  solve (MatrixType& mattype, const Matrix& b, octave_idx_type& info,
315  double& rcond, solve_singularity_handler sing_handler,
316  bool singular_fallback = true) const;
317 
318  OCTAVE_API ComplexMatrix
319  solve (MatrixType& mattype, const ComplexMatrix& b) const;
320  OCTAVE_API ComplexMatrix
321  solve (MatrixType& mattype, const ComplexMatrix& b,
322  octave_idx_type& info) const;
323  OCTAVE_API ComplexMatrix
324  solve (MatrixType& mattype, const ComplexMatrix& b,
325  octave_idx_type& info, double& rcond) const;
326  OCTAVE_API ComplexMatrix
327  solve (MatrixType& mattype, const ComplexMatrix& b,
328  octave_idx_type& info, double& rcond,
329  solve_singularity_handler sing_handler,
330  bool singular_fallback = true) const;
331 
332  OCTAVE_API SparseComplexMatrix
333  solve (MatrixType& mattype, const SparseMatrix& b) const;
334  OCTAVE_API SparseComplexMatrix
335  solve (MatrixType& mattype, const SparseMatrix& b,
336  octave_idx_type& info) const;
337  OCTAVE_API SparseComplexMatrix
338  solve (MatrixType& mattype, const SparseMatrix& b, octave_idx_type& info,
339  double& rcond) const;
340  OCTAVE_API SparseComplexMatrix
341  solve (MatrixType& mattype, const SparseMatrix& b, octave_idx_type& info,
342  double& rcond, solve_singularity_handler sing_handler,
343  bool singular_fallback = true) const;
344 
345  OCTAVE_API SparseComplexMatrix
346  solve (MatrixType& mattype, const SparseComplexMatrix& b) const;
347  OCTAVE_API SparseComplexMatrix
348  solve (MatrixType& mattype, const SparseComplexMatrix& b,
349  octave_idx_type& info) const;
350  OCTAVE_API SparseComplexMatrix
351  solve (MatrixType& mattype, const SparseComplexMatrix& b,
352  octave_idx_type& info, double& rcond) const;
353  OCTAVE_API SparseComplexMatrix
354  solve (MatrixType& mattype, const SparseComplexMatrix& b,
355  octave_idx_type& info, double& rcond,
356  solve_singularity_handler sing_handler,
357  bool singular_fallback = true) const;
358 
359  OCTAVE_API ComplexColumnVector
360  solve (MatrixType& mattype, const ColumnVector& b) const;
361  OCTAVE_API ComplexColumnVector
362  solve (MatrixType& mattype, const ColumnVector& b,
363  octave_idx_type& info) const;
364  OCTAVE_API ComplexColumnVector
365  solve (MatrixType& mattype, const ColumnVector& b,
366  octave_idx_type& info, double& rcond) const;
367  OCTAVE_API ComplexColumnVector
368  solve (MatrixType& mattype, const ColumnVector& b,
369  octave_idx_type& info, double& rcond,
370  solve_singularity_handler sing_handler) const;
371 
372  OCTAVE_API ComplexColumnVector
373  solve (MatrixType& mattype, const ComplexColumnVector& b) const;
374  OCTAVE_API ComplexColumnVector
375  solve (MatrixType& mattype, const ComplexColumnVector& b,
376  octave_idx_type& info) const;
377  OCTAVE_API ComplexColumnVector
378  solve (MatrixType& mattype, const ComplexColumnVector& b,
379  octave_idx_type& info, double& rcond) const;
380  OCTAVE_API ComplexColumnVector
381  solve (MatrixType& mattype, const ComplexColumnVector& b,
382  octave_idx_type& info, double& rcond,
383  solve_singularity_handler sing_handler) const;
384 
385  // Generic interface to solver with probing of type
386  OCTAVE_API ComplexMatrix solve (const Matrix& b) const;
387  OCTAVE_API ComplexMatrix
388  solve (const Matrix& b, octave_idx_type& info) const;
389  OCTAVE_API ComplexMatrix
390  solve (const Matrix& b, octave_idx_type& info, double& rcond) const;
391  OCTAVE_API ComplexMatrix
392  solve (const Matrix& b, octave_idx_type& info, double& rcond,
393  solve_singularity_handler sing_handler) const;
394 
395  OCTAVE_API ComplexMatrix solve (const ComplexMatrix& b) const;
396  OCTAVE_API ComplexMatrix
397  solve (const ComplexMatrix& b, octave_idx_type& info) const;
398  OCTAVE_API ComplexMatrix
399  solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond) const;
400  OCTAVE_API ComplexMatrix
401  solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond,
402  solve_singularity_handler sing_handler) const;
403 
404  OCTAVE_API SparseComplexMatrix solve (const SparseMatrix& b) const;
405  OCTAVE_API SparseComplexMatrix
406  solve (const SparseMatrix& b, octave_idx_type& info) const;
407  OCTAVE_API SparseComplexMatrix
408  solve (const SparseMatrix& b, octave_idx_type& info, double& rcond) const;
409  OCTAVE_API SparseComplexMatrix
410  solve (const SparseMatrix& b, octave_idx_type& info, double& rcond,
411  solve_singularity_handler sing_handler) const;
412 
413  OCTAVE_API SparseComplexMatrix solve (const SparseComplexMatrix& b) const;
414  OCTAVE_API SparseComplexMatrix
415  solve (const SparseComplexMatrix& b, octave_idx_type& info) const;
416  OCTAVE_API SparseComplexMatrix
417  solve (const SparseComplexMatrix& b, octave_idx_type& info,
418  double& rcond) const;
419  OCTAVE_API SparseComplexMatrix
420  solve (const SparseComplexMatrix& b, octave_idx_type& info, double& rcond,
421  solve_singularity_handler sing_handler) const;
422 
423  OCTAVE_API ComplexColumnVector solve (const ColumnVector& b) const;
424  OCTAVE_API ComplexColumnVector
425  solve (const ColumnVector& b, octave_idx_type& info) const;
426  OCTAVE_API ComplexColumnVector
427  solve (const ColumnVector& b, octave_idx_type& info, double& rcond) const;
428  OCTAVE_API ComplexColumnVector
429  solve (const ColumnVector& b, octave_idx_type& info, double& rcond,
430  solve_singularity_handler sing_handler) const;
431 
432  OCTAVE_API ComplexColumnVector solve (const ComplexColumnVector& b) const;
433  OCTAVE_API ComplexColumnVector
434  solve (const ComplexColumnVector& b, octave_idx_type& info) const;
435  OCTAVE_API ComplexColumnVector
436  solve (const ComplexColumnVector& b, octave_idx_type& info,
437  double& rcond) const;
438  OCTAVE_API ComplexColumnVector
439  solve (const ComplexColumnVector& b, octave_idx_type& info, double& rcond,
440  solve_singularity_handler sing_handler) const;
441 
442  OCTAVE_API SparseComplexMatrix squeeze () const;
443 
444  OCTAVE_API SparseComplexMatrix reshape (const dim_vector& new_dims) const;
445 
446  OCTAVE_API SparseComplexMatrix
447  permute (const Array<octave_idx_type>& vec, bool inv = false) const;
448 
449  OCTAVE_API SparseComplexMatrix
450  ipermute (const Array<octave_idx_type>& vec) const;
451 
452  OCTAVE_API bool any_element_is_nan () const;
453  OCTAVE_API bool any_element_is_inf_or_nan () const;
454  OCTAVE_API bool all_elements_are_real () const;
455  OCTAVE_API bool all_integers (double& max_val, double& min_val) const;
456  OCTAVE_API bool too_large_for_float () const;
457 
458  OCTAVE_API SparseBoolMatrix operator ! () const;
459 
460  OCTAVE_API SparseBoolMatrix all (int dim = -1) const;
461  OCTAVE_API SparseBoolMatrix any (int dim = -1) const;
462 
463  OCTAVE_API SparseComplexMatrix cumprod (int dim = -1) const;
464  OCTAVE_API SparseComplexMatrix cumsum (int dim = -1) const;
465  OCTAVE_API SparseComplexMatrix prod (int dim = -1) const;
466  OCTAVE_API SparseComplexMatrix sum (int dim = -1) const;
467  OCTAVE_API SparseComplexMatrix sumsq (int dim = -1) const;
468  OCTAVE_API SparseMatrix abs () const;
469 
470  OCTAVE_API SparseComplexMatrix diag (octave_idx_type k = 0) const;
471 
472  // i/o
473  friend OCTAVE_API std::ostream&
474  operator << (std::ostream& os, const SparseComplexMatrix& a);
475  friend OCTAVE_API std::istream&
476  operator >> (std::istream& is, SparseComplexMatrix& a);
477 };
478 
479 extern OCTAVE_API SparseComplexMatrix
480 operator * (const SparseMatrix&, const SparseComplexMatrix&);
481 extern OCTAVE_API SparseComplexMatrix
482 operator * (const SparseComplexMatrix&, const SparseMatrix&);
483 extern OCTAVE_API SparseComplexMatrix
484 operator * (const SparseComplexMatrix&, const SparseComplexMatrix&);
485 
486 extern OCTAVE_API ComplexMatrix
487 operator * (const Matrix&, const SparseComplexMatrix&);
488 extern OCTAVE_API ComplexMatrix
489 operator * (const ComplexMatrix&, const SparseMatrix&);
490 extern OCTAVE_API ComplexMatrix
491 operator * (const ComplexMatrix&, const SparseComplexMatrix&);
492 extern OCTAVE_API ComplexMatrix
493 mul_trans (const ComplexMatrix&, const SparseComplexMatrix&);
494 extern OCTAVE_API ComplexMatrix
495 mul_herm (const ComplexMatrix&, const SparseComplexMatrix&);
496 
497 extern OCTAVE_API ComplexMatrix
498 operator * (const SparseMatrix&, const ComplexMatrix&);
499 extern OCTAVE_API ComplexMatrix
500 operator * (const SparseComplexMatrix&, const Matrix&);
501 extern OCTAVE_API ComplexMatrix
502 operator * (const SparseComplexMatrix&, const ComplexMatrix&);
503 extern OCTAVE_API ComplexMatrix
504 trans_mul (const SparseComplexMatrix&, const ComplexMatrix&);
505 extern OCTAVE_API ComplexMatrix
506 herm_mul (const SparseComplexMatrix&, const ComplexMatrix&);
507 
508 extern OCTAVE_API SparseComplexMatrix
509 operator * (const DiagMatrix&, const SparseComplexMatrix&);
510 extern OCTAVE_API SparseComplexMatrix
511 operator * (const SparseComplexMatrix&, const DiagMatrix&);
512 
513 extern OCTAVE_API SparseComplexMatrix
514 operator * (const ComplexDiagMatrix&, const SparseMatrix&);
515 extern OCTAVE_API SparseComplexMatrix
516 operator * (const SparseMatrix&, const ComplexDiagMatrix&);
517 
518 extern OCTAVE_API SparseComplexMatrix
519 operator * (const ComplexDiagMatrix&, const SparseComplexMatrix&);
520 extern OCTAVE_API SparseComplexMatrix
521 operator * (const SparseComplexMatrix&, const ComplexDiagMatrix&);
522 
523 extern OCTAVE_API SparseComplexMatrix
524 operator + (const ComplexDiagMatrix&, const SparseMatrix&);
525 extern OCTAVE_API SparseComplexMatrix
526 operator + (const DiagMatrix&, const SparseComplexMatrix&);
527 extern OCTAVE_API SparseComplexMatrix
528 operator + (const ComplexDiagMatrix&, const SparseComplexMatrix&);
529 extern OCTAVE_API SparseComplexMatrix
530 operator + (const SparseMatrix&, const ComplexDiagMatrix&);
531 extern OCTAVE_API SparseComplexMatrix
532 operator + (const SparseComplexMatrix&, const DiagMatrix&);
533 extern OCTAVE_API SparseComplexMatrix
534 operator + (const SparseComplexMatrix&, const ComplexDiagMatrix&);
535 
536 extern OCTAVE_API SparseComplexMatrix
537 operator - (const ComplexDiagMatrix&, const SparseMatrix&);
538 extern OCTAVE_API SparseComplexMatrix
539 operator - (const DiagMatrix&, const SparseComplexMatrix&);
540 extern OCTAVE_API SparseComplexMatrix
541 operator - (const ComplexDiagMatrix&, const SparseComplexMatrix&);
542 extern OCTAVE_API SparseComplexMatrix
543 operator - (const SparseMatrix&, const ComplexDiagMatrix&);
544 extern OCTAVE_API SparseComplexMatrix
545 operator - (const SparseComplexMatrix&, const DiagMatrix&);
546 extern OCTAVE_API SparseComplexMatrix
547 operator - (const SparseComplexMatrix&, const ComplexDiagMatrix&);
548 
549 extern OCTAVE_API SparseComplexMatrix
550 operator * (const PermMatrix&, const SparseComplexMatrix&);
551 extern OCTAVE_API SparseComplexMatrix
552 operator * (const SparseComplexMatrix&, const PermMatrix&);
553 
554 extern OCTAVE_API SparseComplexMatrix
555 min (const Complex& c, const SparseComplexMatrix& m);
556 extern OCTAVE_API SparseComplexMatrix
557 min (const SparseComplexMatrix& m, const Complex& c);
558 extern OCTAVE_API SparseComplexMatrix
559 min (const SparseComplexMatrix& a, const SparseComplexMatrix& b);
560 
561 extern OCTAVE_API SparseComplexMatrix
562 max (const Complex& c, const SparseComplexMatrix& m);
563 extern OCTAVE_API SparseComplexMatrix
564 max (const SparseComplexMatrix& m, const Complex& c);
565 extern OCTAVE_API SparseComplexMatrix
566 max (const SparseComplexMatrix& a, const SparseComplexMatrix& b);
567 
568 SPARSE_SMS_CMP_OP_DECLS (SparseComplexMatrix, Complex, OCTAVE_API)
569 SPARSE_SMS_BOOL_OP_DECLS (SparseComplexMatrix, Complex, OCTAVE_API)
570 
571 SPARSE_SSM_CMP_OP_DECLS (Complex, SparseComplexMatrix, OCTAVE_API)
572 SPARSE_SSM_BOOL_OP_DECLS (Complex, SparseComplexMatrix, OCTAVE_API)
573 
574 SPARSE_SMSM_CMP_OP_DECLS (SparseComplexMatrix, SparseComplexMatrix, OCTAVE_API)
575 SPARSE_SMSM_BOOL_OP_DECLS (SparseComplexMatrix, SparseComplexMatrix, OCTAVE_API)
576 
577 SPARSE_FORWARD_DEFS (MSparse, SparseComplexMatrix, ComplexMatrix, Complex)
578 
579 #endif
operator<<
template std::ostream & operator<<(std::ostream &, const Array< bool > &)
conj
ComplexColumnVector conj(const ComplexColumnVector &a)
Definition: CColVector.cc:217
concat
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:418
max
SparseComplexMatrix max(const Complex &c, const SparseComplexMatrix &m)
Definition: CSparse.cc:7791
herm_mul
ComplexMatrix herm_mul(const SparseComplexMatrix &, const ComplexMatrix &)
Definition: CSparse.cc:7552
operator+
SparseComplexMatrix operator+(const ComplexDiagMatrix &, const SparseMatrix &)
Definition: CSparse.cc:7592
operator*
SparseComplexMatrix operator*(const SparseMatrix &, const SparseComplexMatrix &)
Definition: CSparse.cc:7486
mul_herm
ComplexMatrix mul_herm(const ComplexMatrix &, const SparseComplexMatrix &)
Definition: CSparse.cc:7522
mul_trans
ComplexMatrix mul_trans(const ComplexMatrix &, const SparseComplexMatrix &)
Definition: CSparse.cc:7516
min
SparseComplexMatrix min(const Complex &c, const SparseComplexMatrix &m)
Definition: CSparse.cc:7674
trans_mul
ComplexMatrix trans_mul(const SparseComplexMatrix &, const ComplexMatrix &)
Definition: CSparse.cc:7546
operator-
SparseComplexMatrix operator-(const ComplexDiagMatrix &, const SparseMatrix &)
Definition: CSparse.cc:7623
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
Array
N Dimensional Array with copy-on-write semantics.
Definition: Array.h:130
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
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const ComplexNDArray &a)
Definition: CSparse.h:83
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const ComplexMatrix &a)
Definition: CSparse.h:80
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const Sparse< Complex > &a)
Definition: CSparse.h:78
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(octave_idx_type r, octave_idx_type c, double val)
Definition: CSparse.h:67
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(octave_idx_type r, octave_idx_type c)
Definition: CSparse.h:57
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const SparseComplexMatrix &a)
Definition: CSparse.h:70
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const dim_vector &dv, octave_idx_type nz=0)
Definition: CSparse.h:60
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const MSparse< Complex > &a)
Definition: CSparse.h:76
SparseComplexMatrix::solve
ComplexMatrix solve(const ComplexMatrix &b) const
SparseComplexMatrix::dense_matrix_type
ComplexMatrix dense_matrix_type
Definition: CSparse.h:51
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(octave_idx_type r, octave_idx_type c, Complex val)
Definition: CSparse.h:63
SparseComplexMatrix::~SparseComplexMatrix
~SparseComplexMatrix()=default
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(octave_idx_type r, octave_idx_type c, octave_idx_type num_nz)
Definition: CSparse.h:98
SparseComplexMatrix::transpose
SparseComplexMatrix transpose() const
Definition: CSparse.h:139
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const SparseComplexMatrix &a, const dim_vector &dv)
Definition: CSparse.h:73
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(const Array< Complex > &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: CSparse.h:86
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!=
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
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
Complex
std::complex< double > Complex
Definition: oct-cmplx.h:33
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)