GNU Octave  9.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-s-m.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-scalar.h"
33 #include "ov-float.h"
34 #include "ov-re-mat.h"
35 #include "ov-flt-re-mat.h"
36 #include "ov-typeinfo.h"
37 #include "ops.h"
38 #include "xdiv.h"
39 #include "xpow.h"
40 
42 
43 // scalar by matrix ops.
44 
45 DEFNDBINOP_OP (add, scalar, matrix, scalar, array, +)
46 DEFNDBINOP_OP (sub, scalar, matrix, scalar, array, -)
47 DEFNDBINOP_OP (mul, scalar, matrix, scalar, array, *)
48 
49 DEFBINOP (div, scalar, matrix)
50 {
51  OCTAVE_CAST_BASE_VALUE (const octave_scalar&, v1, a1);
53 
54  Matrix m1 = v1.matrix_value ();
55  Matrix m2 = v2.matrix_value ();
56  MatrixType typ = v2.matrix_type ();
57 
58  Matrix ret = xdiv (m1, m2, typ);
59 
60  v2.matrix_type (typ);
61  return ret;
62 }
63 
64 DEFBINOP_FN (pow, scalar, matrix, xpow)
65 
66 DEFBINOP (ldiv, scalar, matrix)
67 {
68  OCTAVE_CAST_BASE_VALUE (const octave_scalar&, v1, a1);
70 
71  return octave_value (v2.array_value () / v1.double_value ());
72 }
73 
74 DEFNDBINOP_FN (lt, scalar, matrix, scalar, array, mx_el_lt)
75 DEFNDBINOP_FN (le, scalar, matrix, scalar, array, mx_el_le)
76 DEFNDBINOP_FN (eq, scalar, matrix, scalar, array, mx_el_eq)
77 DEFNDBINOP_FN (ge, scalar, matrix, scalar, array, mx_el_ge)
78 DEFNDBINOP_FN (gt, scalar, matrix, scalar, array, mx_el_gt)
79 DEFNDBINOP_FN (ne, scalar, matrix, scalar, array, mx_el_ne)
80 
81 DEFNDBINOP_OP (el_mul, scalar, matrix, scalar, array, *)
82 DEFNDBINOP_FN (el_div, scalar, matrix, scalar, array, elem_xdiv)
83 DEFNDBINOP_FN (el_pow, scalar, matrix, scalar, array, elem_xpow)
84 
85 DEFBINOP (el_ldiv, scalar, matrix)
86 {
87  OCTAVE_CAST_BASE_VALUE (const octave_scalar&, v1, a1);
89 
90  return octave_value (v2.array_value () / v1.double_value ());
91 }
92 
93 DEFNDBINOP_FN (el_and, scalar, matrix, scalar, array, mx_el_and)
94 DEFNDBINOP_FN (el_or, scalar, matrix, scalar, array, mx_el_or)
95 
96 DEFNDCATOP_FN (s_m, scalar, matrix, array, array, concat)
97 
98 DEFCONV (matrix_conv, scalar, matrix)
99 {
100  OCTAVE_CAST_BASE_VALUE (const octave_scalar&, v, a);
101 
102  return new octave_matrix (v.matrix_value ());
103 }
104 
105 void
106 install_s_m_ops (octave::type_info& ti)
107 {
126 
128 
132 
133  INSTALL_WIDENOP_TI (ti, octave_scalar, octave_matrix, matrix_conv);
134 }
135 
136 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
Definition: dMatrix.h:42
MatrixType matrix_type() const
Definition: ov-base-mat.h:139
Matrix matrix_value(bool=false) const
Definition: ov-ch-mat.h:114
NDArray array_value(bool=false) const
Definition: ov-ch-mat.h:120
OCTAVE_BEGIN_NAMESPACE(octave) static octave_value daspk_fcn
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
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
void install_s_m_ops(octave::type_info &ti)
Definition: op-s-m.cc:106
#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 DEFCONV(name, a_dummy, b_dummy)
Definition: ops.h:187
#define DEFBINOP(name, t1, t2)
Definition: ops.h:245
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:305
#define INSTALL_ASSIGNCONV_TI(ti, t1, t2, tr)
Definition: ops.h:88
#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_WIDENOP_TI(ti, t1, t2, f)
Definition: ops.h:92
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:75
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 xdiv(const Matrix &a, const SparseMatrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:137
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