GNU Octave  9.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
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)