GNU Octave  6.2.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-fm-fs.cc
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 //
3 // Copyright (C) 1996-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-flt-re-mat.h"
33 #include "ov-float.h"
34 #include "ov-scalar.h"
35 #include "ov-typeinfo.h"
36 #include "ops.h"
37 #include "xdiv.h"
38 #include "xpow.h"
39 
40 // matrix by scalar ops.
41 
42 DEFNDBINOP_OP (add, float_matrix, float_scalar, float_array, float_scalar, +)
43 DEFNDBINOP_OP (sub, float_matrix, float_scalar, float_array, float_scalar, -)
44 DEFNDBINOP_OP (mul, float_matrix, float_scalar, float_array, float_scalar, *)
45 
46 DEFBINOP (div, float_matrix, float)
47 {
48  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
49  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
50 
51  return octave_value (v1.float_array_value () / v2.float_value ());
52 }
53 
54 DEFBINOP_FN (pow, float_matrix, float_scalar, xpow)
55 
56 DEFBINOP (ldiv, float_matrix, float)
57 {
58  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
59  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
60 
63  MatrixType typ = v1.matrix_type ();
64 
65  FloatMatrix ret = xleftdiv (m1, m2, typ);
66 
67  v1.matrix_type (typ);
68  return ret;
69 }
70 
71 DEFNDBINOP_FN (lt, float_matrix, float_scalar, float_array,
72  float_scalar, mx_el_lt)
73 DEFNDBINOP_FN (le, float_matrix, float_scalar, float_array,
74  float_scalar, mx_el_le)
75 DEFNDBINOP_FN (eq, float_matrix, float_scalar, float_array,
76  float_scalar, mx_el_eq)
77 DEFNDBINOP_FN (ge, float_matrix, float_scalar, float_array,
78  float_scalar, mx_el_ge)
79 DEFNDBINOP_FN (gt, float_matrix, float_scalar, float_array,
80  float_scalar, mx_el_gt)
81 DEFNDBINOP_FN (ne, float_matrix, float_scalar, float_array,
82  float_scalar, mx_el_ne)
83 
84 DEFNDBINOP_OP (el_mul, float_matrix, float_scalar, float_array, float_scalar, *)
85 
86 DEFBINOP (el_div, float_matrix, float)
87 {
88  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
89  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
90 
91  return octave_value (v1.float_array_value () / v2.float_value ());
92 }
93 
94 DEFNDBINOP_FN (el_pow, float_matrix, float_scalar, float_array,
95  float_scalar, elem_xpow)
96 
97 DEFBINOP (el_ldiv, float_matrix, float)
98 {
99  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
100  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
101 
102  return x_el_div (v2.float_value (), v1.float_array_value ());
103 }
104 
105 DEFNDBINOP_FN (el_and, float_matrix, float_scalar, float_array,
106  float_scalar, mx_el_and)
107 DEFNDBINOP_FN (el_or, float_matrix, float_scalar, float_array,
108  float_scalar, mx_el_or)
109 
110 DEFNDCATOP_FN (fm_fs, float_matrix, float_scalar, float_array,
111  float_array, concat)
112 
113 DEFNDCATOP_FN (m_fs, matrix, float_scalar, float_array, float_array, concat)
114 
115 DEFNDCATOP_FN (fm_s, float_matrix, scalar, float_array, float_array, concat)
116 
117 DEFNDASSIGNOP_FN (assign, float_matrix, float_scalar, float_scalar, assign)
118 DEFNDASSIGNOP_FN (dbl_assign, matrix, float_scalar, scalar, assign)
119 
120 DEFNDASSIGNOP_OP (assign_add, float_matrix, float_scalar, float_scalar, +=)
121 DEFNDASSIGNOP_OP (assign_sub, float_matrix, float_scalar, float_scalar, -=)
122 DEFNDASSIGNOP_OP (assign_mul, float_matrix, float_scalar, float_scalar, *=)
123 DEFNDASSIGNOP_OP (assign_div, float_matrix, float_scalar, float_scalar, /=)
124 
125 void
126 install_fm_fs_ops (octave::type_info& ti)
127 {
141  el_mul);
143  el_div);
145  el_pow);
147  el_ldiv);
149  el_and);
151  el_or);
152 
156 
158  assign);
160  dbl_assign);
161 
163  assign_add);
165  assign_sub);
167  assign_mul);
169  assign_div);
170 }
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
Definition: ov-base-mat.h:131
float float_value(bool=false) const
Definition: ov-ch-mat.cc:83
FloatMatrix float_matrix_value(bool=false) const
Definition: ov-ch-mat.h:117
FloatNDArray float_array_value(bool=false) const
FloatMatrix float_matrix_value(bool=false) const
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
void install_fm_fs_ops(octave::type_info &ti)
Definition: op-fm-fs.cc:126
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:112
#define DEFNDBINOP_OP(name, t1, t2, e1, e2, op)
Definition: ops.h:276
#define DEFNDBINOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:301
#define DEFBINOP(name, t1, t2)
Definition: ops.h:230
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:290
#define INSTALL_BINOP_TI(ti, op, t1, t2, f)
Definition: ops.h:55
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:348
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:60
#define DEFNDASSIGNOP_OP(name, t1, t2, f, op)
Definition: ops.h:126
#define INSTALL_ASSIGNOP_TI(ti, op, t1, t2, f)
Definition: ops.h:64
static bool scalar(const dim_vector &dims)
Definition: ov-struct.cc:669
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