GNU Octave  9.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-fcm-fcs.cc
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 //
3 // Copyright (C) 1996-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 (HAVE_CONFIG_H)
27 # include "config.h"
28 #endif
29 
30 #include "ovl.h"
31 #include "ov.h"
32 #include "ov-cx-mat.h"
33 #include "ov-flt-cx-mat.h"
34 #include "ov-flt-complex.h"
35 #include "ov-complex.h"
36 #include "ov-typeinfo.h"
37 #include "ops.h"
38 #include "xdiv.h"
39 #include "xpow.h"
40 
42 
43 // complex matrix by complex scalar ops.
44 
45 DEFNDBINOP_OP (add, float_complex_matrix, float_complex,
46  float_complex_array, float_complex, +)
47 DEFNDBINOP_OP (sub, float_complex_matrix, float_complex,
48  float_complex_array, float_complex, -)
49 DEFNDBINOP_OP (mul, float_complex_matrix, float_complex,
50  float_complex_array, float_complex, *)
51 
52 DEFBINOP (div, float_complex_matrix, float_complex)
53 {
56 
57  return octave_value (v1.float_complex_array_value () / v2.float_complex_value ());
58 }
59 
60 DEFBINOP_FN (pow, float_complex_matrix, float_complex, xpow)
61 
62 DEFBINOP (ldiv, float_complex_matrix, float_complex)
63 {
66 
67  FloatComplexMatrix m1 = v1.float_complex_matrix_value ();
69  MatrixType typ = v1.matrix_type ();
70 
71  FloatComplexMatrix ret = xleftdiv (m1, m2, typ);
72  v1.matrix_type (typ);
73  return ret;
74 }
75 
76 DEFNDCMPLXCMPOP_FN (lt, float_complex_matrix, float_complex,
77  float_complex_array, float_complex, mx_el_lt)
78 DEFNDCMPLXCMPOP_FN (le, float_complex_matrix, float_complex,
79  float_complex_array, float_complex, mx_el_le)
80 DEFNDCMPLXCMPOP_FN (eq, float_complex_matrix, float_complex,
81  float_complex_array, float_complex, mx_el_eq)
82 DEFNDCMPLXCMPOP_FN (ge, float_complex_matrix, float_complex,
83  float_complex_array, float_complex, mx_el_ge)
84 DEFNDCMPLXCMPOP_FN (gt, float_complex_matrix, float_complex,
85  float_complex_array, float_complex, mx_el_gt)
86 DEFNDCMPLXCMPOP_FN (ne, float_complex_matrix, float_complex,
87  float_complex_array, float_complex, mx_el_ne)
88 
89 DEFNDBINOP_OP (el_mul, float_complex_matrix, float_complex,
90  float_complex_array, float_complex, *)
91 
92 DEFBINOP (el_div, float_complex_matrix, float_complex)
93 {
96 
97  return octave_value (v1.float_complex_array_value () / v2.float_complex_value ());
98 }
99 
100 DEFNDBINOP_FN (el_pow, float_complex_matrix, float_complex,
101  float_complex_array, float_complex, elem_xpow)
102 
103 DEFBINOP (el_ldiv, float_complex_matrix, float_complex)
104 {
107 
108  return elem_xdiv (v2.float_complex_value (), v1.float_complex_array_value ());
109 }
110 
111 DEFNDBINOP_FN (el_and, float_complex_matrix, float_complex,
112  float_complex_array, float_complex, mx_el_and)
113 DEFNDBINOP_FN (el_or, float_complex_matrix, float_complex,
114  float_complex_array, float_complex, mx_el_or)
115 
116 DEFNDCATOP_FN (fcm_fcs, float_complex_matrix, float_complex,
117  float_complex_array, float_complex_array, concat)
118 
119 DEFNDCATOP_FN (cm_fcs, complex_matrix, float_complex,
120  float_complex_array, float_complex_array, concat)
121 
122 DEFNDCATOP_FN (fcm_cs, float_complex_matrix, complex,
123  float_complex_array, float_complex_array, concat)
124 
125 DEFNDASSIGNOP_FN (assign, float_complex_matrix, float_complex,
126  float_complex, assign)
127 DEFNDASSIGNOP_FN (dbl_assign, complex_matrix, float_complex,
128  complex, assign)
129 
130 DEFNDASSIGNOP_OP (assign_add, float_complex_matrix, float_complex_scalar,
131  float_complex, +=)
132 DEFNDASSIGNOP_OP (assign_sub, float_complex_matrix, float_complex_scalar,
133  float_complex, -=)
134 DEFNDASSIGNOP_OP (assign_mul, float_complex_matrix, float_complex_scalar,
135  float_complex, *=)
136 DEFNDASSIGNOP_OP (assign_div, float_complex_matrix, float_complex_scalar,
137  float_complex, /=)
138 
139 void
141 {
143  octave_float_complex, add);
145  octave_float_complex, sub);
147  octave_float_complex, mul);
149  octave_float_complex, div);
153  octave_float_complex, ldiv);
155  lt);
157  le);
159  eq);
161  ge);
163  gt);
165  ne);
167  octave_float_complex, el_mul);
169  octave_float_complex, el_div);
171  octave_float_complex, el_pow);
173  octave_float_complex, el_ldiv);
175  octave_float_complex, el_and);
177  octave_float_complex, el_or);
178 
180  fcm_fcs);
183 
185  octave_float_complex, assign);
187  octave_float_complex, dbl_assign);
188 
190  octave_float_complex_scalar, assign_add);
192  octave_float_complex_scalar, assign_sub);
194  octave_float_complex_scalar, assign_mul);
196  octave_float_complex_scalar, assign_div);
197 }
198 
199 OCTAVE_END_NAMESPACE(octave)
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
FloatComplex float_complex_value(bool=false) const
Definition: ov-ch-mat.cc:139
FloatComplexMatrix float_complex_matrix_value(bool=false) const
Definition: ov-ch-mat.h:133
OCTAVE_BEGIN_NAMESPACE(octave) static octave_value daspk_fcn
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
void install_fcm_fcs_ops(octave::type_info &ti)
Definition: op-fcm-fcs.cc:140
const octave_base_value & a2
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
#define DEFNDASSIGNOP_FN(name, t1, t2, e, f)
Definition: ops.h:127
#define DEFNDBINOP_OP(name, t1, t2, e1, e2, op)
Definition: ops.h:291
#define DEFNDBINOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:316
#define DEFNDCMPLXCMPOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:327
#define DEFBINOP(name, t1, t2)
Definition: ops.h:245
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:305
#define OCTAVE_CAST_BASE_VALUE(T, T_VAL, BASE_VAL)
Definition: ops.h:52
#define INSTALL_BINOP_TI(ti, op, t1, t2, f)
Definition: ops.h:70
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:364
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:75
#define DEFNDASSIGNOP_OP(name, t1, t2, f, op)
Definition: ops.h:141
#define INSTALL_ASSIGNOP_TI(ti, op, t1, t2, f)
Definition: ops.h:79
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1676
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1681
octave_value op_add(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1663
octave_value op_gt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1675
octave_value op_el_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1679
octave_value op_el_or(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1683
octave_value op_sub(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1664
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1665
octave_value op_le(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1672
octave_value op_el_and(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1682
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1673
octave_value op_el_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1680
octave_value op_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1669
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1668
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1674
octave_value op_el_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1678
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1666
octave_value op_lt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1671
Matrix elem_xdiv(double a, const SparseMatrix &b)
Definition: sparse-xdiv.cc:376
Matrix xleftdiv(const SparseMatrix &a, const Matrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:469
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:62
octave_value elem_xpow(double a, const SparseMatrix &b)
Definition: sparse-xpow.cc:347