GNU Octave  8.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-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_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 (void) : 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  OCTAVE_API bool operator == (const SparseComplexMatrix& a) const;
109  OCTAVE_API bool operator != (const SparseComplexMatrix& a) const;
110 
111  OCTAVE_API bool ishermitian (void) const;
112 
113  OCTAVE_API SparseComplexMatrix max (int dim = -1) const;
114  OCTAVE_API SparseComplexMatrix
115  max (Array<octave_idx_type>& index, int dim = -1) const;
116  OCTAVE_API SparseComplexMatrix min (int dim = -1) const;
117  OCTAVE_API SparseComplexMatrix
118  min (Array<octave_idx_type>& index, int dim = -1) const;
119 
120  OCTAVE_API SparseComplexMatrix&
121  insert (const SparseComplexMatrix& a, octave_idx_type r, octave_idx_type c);
122  OCTAVE_API SparseComplexMatrix&
123  insert (const SparseMatrix& a, octave_idx_type r, octave_idx_type c);
124  OCTAVE_API SparseComplexMatrix&
125  insert (const SparseComplexMatrix& a, const Array<octave_idx_type>& indx);
126  OCTAVE_API SparseComplexMatrix&
127  insert (const SparseMatrix& a, const Array<octave_idx_type>& indx);
128 
129  OCTAVE_API SparseComplexMatrix
130  concat (const SparseComplexMatrix& rb, const Array<octave_idx_type>& ra_idx);
131  OCTAVE_API SparseComplexMatrix
132  concat (const SparseMatrix& rb, const Array<octave_idx_type>& ra_idx);
133 
134  OCTAVE_API ComplexMatrix matrix_value (void) const;
135 
136  OCTAVE_API SparseComplexMatrix hermitian (void) const; // complex conjugate transpose
137  SparseComplexMatrix transpose (void) const
138  { return MSparse<Complex>::transpose (); }
139 
140  friend OCTAVE_API SparseComplexMatrix conj (const SparseComplexMatrix& a);
141 
142  // extract row or column i.
143 
144  OCTAVE_API ComplexRowVector row (octave_idx_type i) const;
145 
146  OCTAVE_API ComplexColumnVector column (octave_idx_type i) const;
147 
148 private:
149  SparseComplexMatrix dinverse (MatrixType& mattype, octave_idx_type& info,
150  double& rcond, const bool force = false,
151  const bool calccond = true) const;
152 
153  SparseComplexMatrix tinverse (MatrixType& mattype, octave_idx_type& info,
154  double& rcond, const bool force = false,
155  const bool calccond = true) const;
156 
157 public:
158  OCTAVE_API SparseComplexMatrix inverse (void) const;
159  OCTAVE_API SparseComplexMatrix inverse (MatrixType& mattype) const;
160  OCTAVE_API SparseComplexMatrix
161  inverse (MatrixType& mattype, octave_idx_type& info) const;
162  OCTAVE_API SparseComplexMatrix
163  inverse (MatrixType& mattype, octave_idx_type& info, double& rcond,
164  bool force = false, bool calc_cond = true) const;
165 
166  OCTAVE_API ComplexDET determinant (void) const;
167  OCTAVE_API ComplexDET determinant (octave_idx_type& info) const;
168  OCTAVE_API ComplexDET
169  determinant (octave_idx_type& info, double& rcond,
170  bool calc_cond = true) const;
171 
172 private:
173  // Diagonal matrix solvers
174  ComplexMatrix dsolve (MatrixType& mattype, const Matrix& b,
175  octave_idx_type& info, double& rcond,
176  solve_singularity_handler sing_handler,
177  bool calc_cond = false) const;
178 
179  ComplexMatrix dsolve (MatrixType& mattype, const ComplexMatrix& b,
180  octave_idx_type& info, double& rcond,
181  solve_singularity_handler sing_handler,
182  bool calc_cond = false) const;
183 
184  SparseComplexMatrix dsolve (MatrixType& mattype, const SparseMatrix& b,
185  octave_idx_type& info, double& rcond,
186  solve_singularity_handler sing_handler,
187  bool calc_cond = false) const;
188 
189  SparseComplexMatrix dsolve (MatrixType& mattype, const SparseComplexMatrix& b,
190  octave_idx_type& info, double& rcond,
191  solve_singularity_handler sing_handler,
192  bool calc_cond = false) const;
193 
194  // Upper triangular matrix solvers
195  ComplexMatrix utsolve (MatrixType& mattype, const Matrix& b,
196  octave_idx_type& info, double& rcond,
197  solve_singularity_handler sing_handler,
198  bool calc_cond = false) const;
199 
200  ComplexMatrix utsolve (MatrixType& mattype, const ComplexMatrix& b,
201  octave_idx_type& info, double& rcond,
202  solve_singularity_handler sing_handler,
203  bool calc_cond = false) const;
204 
205  SparseComplexMatrix utsolve (MatrixType& mattype, const SparseMatrix& b,
206  octave_idx_type& info, double& rcond,
207  solve_singularity_handler sing_handler,
208  bool calc_cond = false) const;
209 
210  SparseComplexMatrix utsolve (MatrixType& mattype, const SparseComplexMatrix& b,
211  octave_idx_type& info, double& rcond,
212  solve_singularity_handler sing_handler,
213  bool calc_cond = false) const;
214 
215  // Lower triangular matrix solvers
216  ComplexMatrix ltsolve (MatrixType& mattype, const Matrix& b,
217  octave_idx_type& info, double& rcond,
218  solve_singularity_handler sing_handler,
219  bool calc_cond = false) const;
220 
221  ComplexMatrix ltsolve (MatrixType& mattype, const ComplexMatrix& b,
222  octave_idx_type& info, double& rcond,
223  solve_singularity_handler sing_handler,
224  bool calc_cond = false) const;
225 
226  SparseComplexMatrix ltsolve (MatrixType& mattype, const SparseMatrix& b,
227  octave_idx_type& info, double& rcond,
228  solve_singularity_handler sing_handler,
229  bool calc_cond = false) const;
230 
231  SparseComplexMatrix ltsolve (MatrixType& mattype, const SparseComplexMatrix& b,
232  octave_idx_type& info, double& rcond,
233  solve_singularity_handler sing_handler,
234  bool calc_cond = false) const;
235 
236  // Tridiagonal matrix solvers
237  ComplexMatrix trisolve (MatrixType& mattype, const Matrix& b,
238  octave_idx_type& info, double& rcond,
239  solve_singularity_handler sing_handler,
240  bool calc_cond = false) const;
241 
242  ComplexMatrix trisolve (MatrixType& mattype, const ComplexMatrix& b,
243  octave_idx_type& info, double& rcond,
244  solve_singularity_handler sing_handler,
245  bool calc_cond = false) const;
246 
247  SparseComplexMatrix trisolve (MatrixType& mattype, const SparseMatrix& b,
248  octave_idx_type& info, double& rcond,
249  solve_singularity_handler sing_handler,
250  bool calc_cond = false) const;
251 
252  SparseComplexMatrix trisolve (MatrixType& mattype, const SparseComplexMatrix& b,
253  octave_idx_type& info, double& rcond,
254  solve_singularity_handler sing_handler,
255  bool calc_cond = false) const;
256 
257  // Banded matrix solvers (umfpack/cholesky)
258  ComplexMatrix bsolve (MatrixType& mattype, const Matrix& b,
259  octave_idx_type& info, double& rcond,
260  solve_singularity_handler sing_handler,
261  bool calc_cond = false) const;
262 
263  ComplexMatrix bsolve (MatrixType& mattype, const ComplexMatrix& b,
264  octave_idx_type& info, double& rcond,
265  solve_singularity_handler sing_handler,
266  bool calc_cond = false) const;
267 
268  SparseComplexMatrix bsolve (MatrixType& mattype, const SparseMatrix& b,
269  octave_idx_type& info, double& rcond,
270  solve_singularity_handler sing_handler,
271  bool calc_cond = false) const;
272 
273  SparseComplexMatrix bsolve (MatrixType& mattype, const SparseComplexMatrix& b,
274  octave_idx_type& info, double& rcond,
275  solve_singularity_handler sing_handler,
276  bool calc_cond = false) const;
277 
278  // Full matrix solvers (umfpack/cholesky)
279  void * factorize (octave_idx_type& err, double& rcond, Matrix& Control,
280  Matrix& Info, solve_singularity_handler sing_handler,
281  bool calc_cond) const;
282 
283  ComplexMatrix fsolve (MatrixType& mattype, const Matrix& b,
284  octave_idx_type& info, double& rcond,
285  solve_singularity_handler sing_handler,
286  bool calc_cond = false) const;
287 
288  ComplexMatrix fsolve (MatrixType& mattype, const ComplexMatrix& b,
289  octave_idx_type& info, double& rcond,
290  solve_singularity_handler sing_handler,
291  bool calc_cond = false) const;
292 
293  SparseComplexMatrix fsolve (MatrixType& mattype, const SparseMatrix& b,
294  octave_idx_type& info, double& rcond,
295  solve_singularity_handler sing_handler,
296  bool calc_cond = false) const;
297 
298  SparseComplexMatrix fsolve (MatrixType& mattype, const SparseComplexMatrix& b,
299  octave_idx_type& info, double& rcond,
300  solve_singularity_handler sing_handler,
301  bool calc_cond = false) const;
302 
303 public:
304  // Generic interface to solver with no probing of type
305  OCTAVE_API ComplexMatrix solve (MatrixType& mattype, const Matrix& b) const;
306  OCTAVE_API ComplexMatrix
307  solve (MatrixType& mattype, const Matrix& b, octave_idx_type& info) const;
308  OCTAVE_API ComplexMatrix
309  solve (MatrixType& mattype, const Matrix& b, octave_idx_type& info,
310  double& rcond) const;
311  OCTAVE_API ComplexMatrix
312  solve (MatrixType& mattype, const Matrix& b, octave_idx_type& info,
313  double& rcond, solve_singularity_handler sing_handler,
314  bool singular_fallback = true) const;
315 
316  OCTAVE_API ComplexMatrix
317  solve (MatrixType& mattype, const ComplexMatrix& b) const;
318  OCTAVE_API ComplexMatrix
319  solve (MatrixType& mattype, const ComplexMatrix& b,
320  octave_idx_type& info) const;
321  OCTAVE_API ComplexMatrix
322  solve (MatrixType& mattype, const ComplexMatrix& b,
323  octave_idx_type& info, double& rcond) const;
324  OCTAVE_API ComplexMatrix
325  solve (MatrixType& mattype, const ComplexMatrix& b,
326  octave_idx_type& info, double& rcond,
327  solve_singularity_handler sing_handler,
328  bool singular_fallback = true) const;
329 
330  OCTAVE_API SparseComplexMatrix
331  solve (MatrixType& mattype, const SparseMatrix& b) const;
332  OCTAVE_API SparseComplexMatrix
333  solve (MatrixType& mattype, const SparseMatrix& b,
334  octave_idx_type& info) const;
335  OCTAVE_API SparseComplexMatrix
336  solve (MatrixType& mattype, const SparseMatrix& b, octave_idx_type& info,
337  double& rcond) const;
338  OCTAVE_API SparseComplexMatrix
339  solve (MatrixType& mattype, const SparseMatrix& b, octave_idx_type& info,
340  double& rcond, solve_singularity_handler sing_handler,
341  bool singular_fallback = true) const;
342 
343  OCTAVE_API SparseComplexMatrix
344  solve (MatrixType& mattype, const SparseComplexMatrix& b) const;
345  OCTAVE_API SparseComplexMatrix
346  solve (MatrixType& mattype, const SparseComplexMatrix& b,
347  octave_idx_type& info) const;
348  OCTAVE_API SparseComplexMatrix
349  solve (MatrixType& mattype, const SparseComplexMatrix& b,
350  octave_idx_type& info, double& rcond) const;
351  OCTAVE_API SparseComplexMatrix
352  solve (MatrixType& mattype, const SparseComplexMatrix& b,
353  octave_idx_type& info, double& rcond,
354  solve_singularity_handler sing_handler,
355  bool singular_fallback = true) const;
356 
357  OCTAVE_API ComplexColumnVector
358  solve (MatrixType& mattype, const ColumnVector& b) const;
359  OCTAVE_API ComplexColumnVector
360  solve (MatrixType& mattype, const ColumnVector& b,
361  octave_idx_type& info) const;
362  OCTAVE_API ComplexColumnVector
363  solve (MatrixType& mattype, const ColumnVector& b,
364  octave_idx_type& info, double& rcond) const;
365  OCTAVE_API ComplexColumnVector
366  solve (MatrixType& mattype, const ColumnVector& b,
367  octave_idx_type& info, double& rcond,
368  solve_singularity_handler sing_handler) const;
369 
370  OCTAVE_API ComplexColumnVector
371  solve (MatrixType& mattype, const ComplexColumnVector& b) const;
372  OCTAVE_API ComplexColumnVector
373  solve (MatrixType& mattype, const ComplexColumnVector& b,
374  octave_idx_type& info) const;
375  OCTAVE_API ComplexColumnVector
376  solve (MatrixType& mattype, const ComplexColumnVector& b,
377  octave_idx_type& info, double& rcond) const;
378  OCTAVE_API ComplexColumnVector
379  solve (MatrixType& mattype, const ComplexColumnVector& b,
380  octave_idx_type& info, double& rcond,
381  solve_singularity_handler sing_handler) const;
382 
383  // Generic interface to solver with probing of type
384  OCTAVE_API ComplexMatrix solve (const Matrix& b) const;
385  OCTAVE_API ComplexMatrix
386  solve (const Matrix& b, octave_idx_type& info) const;
387  OCTAVE_API ComplexMatrix
388  solve (const Matrix& b, octave_idx_type& info, double& rcond) const;
389  OCTAVE_API ComplexMatrix
390  solve (const Matrix& b, octave_idx_type& info, double& rcond,
391  solve_singularity_handler sing_handler) const;
392 
393  OCTAVE_API ComplexMatrix solve (const ComplexMatrix& b) const;
394  OCTAVE_API ComplexMatrix
395  solve (const ComplexMatrix& b, octave_idx_type& info) const;
396  OCTAVE_API ComplexMatrix
397  solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond) const;
398  OCTAVE_API ComplexMatrix
399  solve (const ComplexMatrix& b, octave_idx_type& info, double& rcond,
400  solve_singularity_handler sing_handler) const;
401 
402  OCTAVE_API SparseComplexMatrix solve (const SparseMatrix& b) const;
403  OCTAVE_API SparseComplexMatrix
404  solve (const SparseMatrix& b, octave_idx_type& info) const;
405  OCTAVE_API SparseComplexMatrix
406  solve (const SparseMatrix& b, octave_idx_type& info, double& rcond) const;
407  OCTAVE_API SparseComplexMatrix
408  solve (const SparseMatrix& b, octave_idx_type& info, double& rcond,
409  solve_singularity_handler sing_handler) const;
410 
411  OCTAVE_API SparseComplexMatrix solve (const SparseComplexMatrix& b) const;
412  OCTAVE_API SparseComplexMatrix
413  solve (const SparseComplexMatrix& b, octave_idx_type& info) const;
414  OCTAVE_API SparseComplexMatrix
415  solve (const SparseComplexMatrix& b, octave_idx_type& info,
416  double& rcond) const;
417  OCTAVE_API SparseComplexMatrix
418  solve (const SparseComplexMatrix& b, octave_idx_type& info, double& rcond,
419  solve_singularity_handler sing_handler) const;
420 
421  OCTAVE_API ComplexColumnVector solve (const ColumnVector& b) const;
422  OCTAVE_API ComplexColumnVector
423  solve (const ColumnVector& b, octave_idx_type& info) const;
424  OCTAVE_API ComplexColumnVector
425  solve (const ColumnVector& b, octave_idx_type& info, double& rcond) const;
426  OCTAVE_API ComplexColumnVector
427  solve (const ColumnVector& b, octave_idx_type& info, double& rcond,
428  solve_singularity_handler sing_handler) const;
429 
430  OCTAVE_API ComplexColumnVector solve (const ComplexColumnVector& b) const;
431  OCTAVE_API ComplexColumnVector
432  solve (const ComplexColumnVector& b, octave_idx_type& info) const;
433  OCTAVE_API ComplexColumnVector
434  solve (const ComplexColumnVector& b, octave_idx_type& info,
435  double& rcond) const;
436  OCTAVE_API ComplexColumnVector
437  solve (const ComplexColumnVector& b, octave_idx_type& info, double& rcond,
438  solve_singularity_handler sing_handler) const;
439 
440  OCTAVE_API SparseComplexMatrix squeeze (void) const;
441 
442  OCTAVE_API SparseComplexMatrix reshape (const dim_vector& new_dims) const;
443 
444  OCTAVE_API SparseComplexMatrix
445  permute (const Array<octave_idx_type>& vec, bool inv = false) const;
446 
447  OCTAVE_API SparseComplexMatrix
448  ipermute (const Array<octave_idx_type>& vec) const;
449 
450  OCTAVE_API bool any_element_is_nan (void) const;
451  OCTAVE_API bool any_element_is_inf_or_nan (void) const;
452  OCTAVE_API bool all_elements_are_real (void) const;
453  OCTAVE_API bool all_integers (double& max_val, double& min_val) const;
454  OCTAVE_API bool too_large_for_float (void) const;
455 
456  OCTAVE_API SparseBoolMatrix operator ! (void) const;
457 
458  OCTAVE_API SparseBoolMatrix all (int dim = -1) const;
459  OCTAVE_API SparseBoolMatrix any (int dim = -1) const;
460 
461  OCTAVE_API SparseComplexMatrix cumprod (int dim = -1) const;
462  OCTAVE_API SparseComplexMatrix cumsum (int dim = -1) const;
463  OCTAVE_API SparseComplexMatrix prod (int dim = -1) const;
464  OCTAVE_API SparseComplexMatrix sum (int dim = -1) const;
465  OCTAVE_API SparseComplexMatrix sumsq (int dim = -1) const;
466  OCTAVE_API SparseMatrix abs (void) const;
467 
468  OCTAVE_API SparseComplexMatrix diag (octave_idx_type k = 0) const;
469 
470  // i/o
471  friend OCTAVE_API std::ostream&
472  operator << (std::ostream& os, const SparseComplexMatrix& a);
473  friend OCTAVE_API std::istream&
474  operator >> (std::istream& is, SparseComplexMatrix& a);
475 };
476 
477 extern OCTAVE_API SparseComplexMatrix
478 operator * (const SparseMatrix&, const SparseComplexMatrix&);
479 extern OCTAVE_API SparseComplexMatrix
480 operator * (const SparseComplexMatrix&, const SparseMatrix&);
481 extern OCTAVE_API SparseComplexMatrix
482 operator * (const SparseComplexMatrix&, const SparseComplexMatrix&);
483 
484 extern OCTAVE_API ComplexMatrix
485 operator * (const Matrix&, const SparseComplexMatrix&);
486 extern OCTAVE_API ComplexMatrix
487 operator * (const ComplexMatrix&, const SparseMatrix&);
488 extern OCTAVE_API ComplexMatrix
489 operator * (const ComplexMatrix&, const SparseComplexMatrix&);
490 extern OCTAVE_API ComplexMatrix
491 mul_trans (const ComplexMatrix&, const SparseComplexMatrix&);
492 extern OCTAVE_API ComplexMatrix
493 mul_herm (const ComplexMatrix&, const SparseComplexMatrix&);
494 
495 extern OCTAVE_API ComplexMatrix
496 operator * (const SparseMatrix&, const ComplexMatrix&);
497 extern OCTAVE_API ComplexMatrix
498 operator * (const SparseComplexMatrix&, const Matrix&);
499 extern OCTAVE_API ComplexMatrix
500 operator * (const SparseComplexMatrix&, const ComplexMatrix&);
501 extern OCTAVE_API ComplexMatrix
502 trans_mul (const SparseComplexMatrix&, const ComplexMatrix&);
503 extern OCTAVE_API ComplexMatrix
504 herm_mul (const SparseComplexMatrix&, const ComplexMatrix&);
505 
506 extern OCTAVE_API SparseComplexMatrix
507 operator * (const DiagMatrix&, const SparseComplexMatrix&);
508 extern OCTAVE_API SparseComplexMatrix
509 operator * (const SparseComplexMatrix&, const DiagMatrix&);
510 
511 extern OCTAVE_API SparseComplexMatrix
512 operator * (const ComplexDiagMatrix&, const SparseMatrix&);
513 extern OCTAVE_API SparseComplexMatrix
514 operator * (const SparseMatrix&, const ComplexDiagMatrix&);
515 
516 extern OCTAVE_API SparseComplexMatrix
517 operator * (const ComplexDiagMatrix&, const SparseComplexMatrix&);
518 extern OCTAVE_API SparseComplexMatrix
519 operator * (const SparseComplexMatrix&, const ComplexDiagMatrix&);
520 
521 extern OCTAVE_API SparseComplexMatrix
522 operator + (const ComplexDiagMatrix&, const SparseMatrix&);
523 extern OCTAVE_API SparseComplexMatrix
524 operator + (const DiagMatrix&, const SparseComplexMatrix&);
525 extern OCTAVE_API SparseComplexMatrix
526 operator + (const ComplexDiagMatrix&, const SparseComplexMatrix&);
527 extern OCTAVE_API SparseComplexMatrix
528 operator + (const SparseMatrix&, const ComplexDiagMatrix&);
529 extern OCTAVE_API SparseComplexMatrix
530 operator + (const SparseComplexMatrix&, const DiagMatrix&);
531 extern OCTAVE_API SparseComplexMatrix
532 operator + (const SparseComplexMatrix&, const ComplexDiagMatrix&);
533 
534 extern OCTAVE_API SparseComplexMatrix
535 operator - (const ComplexDiagMatrix&, const SparseMatrix&);
536 extern OCTAVE_API SparseComplexMatrix
537 operator - (const DiagMatrix&, const SparseComplexMatrix&);
538 extern OCTAVE_API SparseComplexMatrix
539 operator - (const ComplexDiagMatrix&, const SparseComplexMatrix&);
540 extern OCTAVE_API SparseComplexMatrix
541 operator - (const SparseMatrix&, const ComplexDiagMatrix&);
542 extern OCTAVE_API SparseComplexMatrix
543 operator - (const SparseComplexMatrix&, const DiagMatrix&);
544 extern OCTAVE_API SparseComplexMatrix
545 operator - (const SparseComplexMatrix&, const ComplexDiagMatrix&);
546 
547 extern OCTAVE_API SparseComplexMatrix
548 operator * (const PermMatrix&, const SparseComplexMatrix&);
549 extern OCTAVE_API SparseComplexMatrix
550 operator * (const SparseComplexMatrix&, const PermMatrix&);
551 
552 extern OCTAVE_API SparseComplexMatrix
553 min (const Complex& c, const SparseComplexMatrix& m);
554 extern OCTAVE_API SparseComplexMatrix
555 min (const SparseComplexMatrix& m, const Complex& c);
556 extern OCTAVE_API SparseComplexMatrix
557 min (const SparseComplexMatrix& a, const SparseComplexMatrix& b);
558 
559 extern OCTAVE_API SparseComplexMatrix
560 max (const Complex& c, const SparseComplexMatrix& m);
561 extern OCTAVE_API SparseComplexMatrix
562 max (const SparseComplexMatrix& m, const Complex& c);
563 extern OCTAVE_API SparseComplexMatrix
564 max (const SparseComplexMatrix& a, const SparseComplexMatrix& b);
565 
566 SPARSE_SMS_CMP_OP_DECLS (SparseComplexMatrix, Complex, OCTAVE_API)
567 SPARSE_SMS_BOOL_OP_DECLS (SparseComplexMatrix, Complex, OCTAVE_API)
568 
569 SPARSE_SSM_CMP_OP_DECLS (Complex, SparseComplexMatrix, OCTAVE_API)
570 SPARSE_SSM_BOOL_OP_DECLS (Complex, SparseComplexMatrix, OCTAVE_API)
571 
572 SPARSE_SMSM_CMP_OP_DECLS (SparseComplexMatrix, SparseComplexMatrix, OCTAVE_API)
573 SPARSE_SMSM_BOOL_OP_DECLS (SparseComplexMatrix, SparseComplexMatrix, OCTAVE_API)
574 
575 SPARSE_FORWARD_DEFS (MSparse, SparseComplexMatrix, ComplexMatrix, Complex)
576 
577 #endif
operator<<
template OCTAVE_API 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
mul_trans
OCTAVE_API ComplexMatrix mul_trans(const ComplexMatrix &, const SparseComplexMatrix &)
Definition: CSparse.cc:7511
operator+
OCTAVE_API SparseComplexMatrix operator+(const ComplexDiagMatrix &, const SparseMatrix &)
Definition: CSparse.cc:7587
operator*
OCTAVE_API SparseComplexMatrix operator*(const SparseMatrix &, const SparseComplexMatrix &)
Definition: CSparse.cc:7481
mul_herm
OCTAVE_API ComplexMatrix mul_herm(const ComplexMatrix &, const SparseComplexMatrix &)
Definition: CSparse.cc:7517
operator-
OCTAVE_API SparseComplexMatrix operator-(const ComplexDiagMatrix &, const SparseMatrix &)
Definition: CSparse.cc:7618
trans_mul
OCTAVE_API ComplexMatrix trans_mul(const SparseComplexMatrix &, const ComplexMatrix &)
Definition: CSparse.cc:7541
min
OCTAVE_API SparseComplexMatrix min(const Complex &c, const SparseComplexMatrix &m)
Definition: CSparse.cc:7669
herm_mul
OCTAVE_API ComplexMatrix herm_mul(const SparseComplexMatrix &, const ComplexMatrix &)
Definition: CSparse.cc:7547
max
OCTAVE_API SparseComplexMatrix max(const Complex &c, const SparseComplexMatrix &m)
Definition: CSparse.cc:7786
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(void) 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::transpose
SparseComplexMatrix transpose(void) const
Definition: CSparse.h:137
SparseComplexMatrix::SparseComplexMatrix
SparseComplexMatrix(void)
Definition: CSparse.h:55
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(octave_idx_type r, octave_idx_type c, octave_idx_type num_nz)
Definition: CSparse.h:98
SparseComplexMatrix::solve
OCTAVE_API ComplexMatrix solve(const ComplexMatrix &b) const
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
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
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
SparseComplexMatrix
class OCTAVE_API SparseComplexMatrix
Definition: mx-fwd.h:56
m
T octave_idx_type m
Definition: mx-inlines.cc:773
r
T * r
Definition: mx-inlines.cc:773
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:1838
abs
static T abs(T x)
Definition: pr-output.cc:1678