GNU Octave  6.2.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-sm-cs.cc
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 //
3 // Copyright (C) 1998-2021 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 (HAVE_CONFIG_H)
27 # include "config.h"
28 #endif
29 
30 #include "ovl.h"
31 #include "ov.h"
32 #include "ov-typeinfo.h"
33 #include "ov-complex.h"
34 #include "ops.h"
35 #include "xpow.h"
36 
37 #include "sparse-xpow.h"
38 #include "sparse-xdiv.h"
39 #include "ov-re-sparse.h"
40 #include "ov-cx-sparse.h"
41 #include "smx-sm-cs.h"
42 #include "smx-cs-sm.h"
43 
44 // sparse matrix by scalar ops.
45 
46 DEFBINOP_OP (add, sparse_matrix, complex, +)
47 DEFBINOP_OP (sub, sparse_matrix, complex, -)
48 DEFBINOP_OP (mul, sparse_matrix, complex, *)
49 
50 DEFBINOP (div, sparse_matrix, complex)
51 {
52  const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1);
53  const octave_complex& v2 = dynamic_cast<const octave_complex&> (a2);
54 
56 }
57 
58 DEFBINOP (pow, sparse_matrix, complex)
59 {
60  const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1);
61  const octave_complex& v2 = dynamic_cast<const octave_complex&> (a2);
62  return xpow (v1.matrix_value (), v2.complex_value ());
63 }
64 
65 DEFBINOP (ldiv, sparse_matrix, complex)
66 {
67  const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1);
68  const octave_complex& v2 = dynamic_cast<const octave_complex&> (a2);
69 
70  if (v1.rows () == 1 && v1.columns () == 1)
72  / v1.scalar_value ()));
73  else
74  {
75  MatrixType typ = v1.matrix_type ();
78  ComplexMatrix ret = xleftdiv (m1, m2, typ);
79  v1.matrix_type (typ);
80  return ret;
81  }
82 }
83 
84 DEFBINOP_FN (lt, sparse_matrix, complex, mx_el_lt)
85 DEFBINOP_FN (le, sparse_matrix, complex, mx_el_le)
86 DEFBINOP_FN (eq, sparse_matrix, complex, mx_el_eq)
87 DEFBINOP_FN (ge, sparse_matrix, complex, mx_el_ge)
88 DEFBINOP_FN (gt, sparse_matrix, complex, mx_el_gt)
89 DEFBINOP_FN (ne, sparse_matrix, complex, mx_el_ne)
90 
91 DEFBINOP_OP (el_mul, sparse_matrix, complex, *)
92 
93 DEFBINOP (el_div, sparse_matrix, complex)
94 {
95  const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1);
96  const octave_complex& v2 = dynamic_cast<const octave_complex&> (a2);
97 
99 }
100 
101 DEFBINOP_FN (el_pow, sparse_matrix, complex, elem_xpow)
102 
103 DEFBINOP (el_ldiv, sparse_matrix, complex)
104 {
105  const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1);
106  const octave_complex& v2 = dynamic_cast<const octave_complex&> (a2);
107 
109  v1.sparse_matrix_value ()));
110 }
111 
112 DEFBINOP_FN (el_and, sparse_matrix, complex, mx_el_and)
113 DEFBINOP_FN (el_or, sparse_matrix, complex, mx_el_or)
114 
115 DEFCATOP (sm_cs, sparse_matrix, complex)
116 {
117  octave_sparse_matrix& v1 = dynamic_cast<octave_sparse_matrix&> (a1);
118  const octave_complex& v2 = dynamic_cast<const octave_complex&> (a2);
119  SparseComplexMatrix tmp (1, 1, v2.complex_value ());
120  return octave_value (v1.sparse_matrix_value (). concat (tmp, ra_idx));
121 }
122 
123 void
125 {
132 
143  el_ldiv);
146 
148 
151 }
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:418
boolMatrix mx_el_le(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:91
boolMatrix mx_el_ne(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:91
boolMatrix mx_el_or(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:88
boolMatrix mx_el_and(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:88
boolMatrix mx_el_gt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:91
boolMatrix mx_el_eq(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:91
boolMatrix mx_el_lt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:91
boolMatrix mx_el_ge(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:91
MatrixType matrix_type(void) const
octave_idx_type columns(void) const
Definition: ov-base.h:325
octave_idx_type rows(void) const
Definition: ov-base.h:318
Complex complex_value(bool=false) const
Definition: ov-ch-mat.cc:127
SparseMatrix sparse_matrix_value(bool=false) const
Definition: ov-re-sparse.h:127
double scalar_value(bool frc_str_conv=false) const
Definition: ov-re-sparse.h:110
Matrix matrix_value(bool=false) const
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
const octave_base_value & a2
const octave_base_value const Array< octave_idx_type > & ra_idx
return octave_value(v1.char_array_value() . concat(v2.char_array_value(), ra_idx),((a1.is_sq_string()||a2.is_sq_string()) ? '\'' :'"'))
const octave_char_matrix & v2
void install_sm_cs_ops(octave::type_info &ti)
Definition: op-sm-cs.cc:124
#define DEFBINOP(name, t1, t2)
Definition: ops.h:230
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:290
#define DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:235
#define INSTALL_ASSIGNCONV_TI(ti, t1, t2, tr)
Definition: ops.h:73
#define INSTALL_BINOP_TI(ti, op, t1, t2, f)
Definition: ops.h:55
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:60
#define DEFCATOP(name, t1, t2)
Definition: ops.h:328
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1579
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1584
octave_value op_add(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1566
octave_value op_gt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1578
octave_value op_el_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1582
octave_value op_el_or(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1586
octave_value op_sub(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1567
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1568
octave_value op_le(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1575
octave_value op_el_and(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1585
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1576
octave_value op_el_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1583
octave_value op_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1572
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1571
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1577
octave_value op_el_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1581
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1569
octave_value op_lt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1574
Matrix x_el_div(double a, const SparseMatrix &b)
Definition: sparse-xdiv.cc:373
Matrix xleftdiv(const SparseMatrix &a, const Matrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:466
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:60
octave_value elem_xpow(double a, const SparseMatrix &b)
Definition: sparse-xpow.cc:249