GNU Octave  9.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-m-scm.cc
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////
2 //
3 // Copyright (C) 1998-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-typeinfo.h"
33 #include "ov-re-mat.h"
34 #include "ov-cx-mat.h"
35 #include "ops.h"
36 #include "xdiv.h"
37 
38 #include "sparse-xpow.h"
39 #include "sparse-xdiv.h"
40 #include "smx-scm-m.h"
41 #include "smx-m-scm.h"
42 #include "ov-cx-sparse.h"
43 
45 
46 // matrix by sparse complex matrix ops.
47 
48 DEFBINOP_OP (add, matrix, sparse_complex_matrix, +)
49 DEFBINOP_OP (sub, matrix, sparse_complex_matrix, -)
50 
51 DEFBINOP_OP (mul, matrix, sparse_complex_matrix, *)
52 
53 DEFBINOP (div, matrix, sparse_complex_matrix)
54 {
55  OCTAVE_CAST_BASE_VALUE (const octave_matrix&, v1, a1);
57 
58  if (v2.rows () == 1 && v2.columns () == 1)
59  return octave_value (v1.array_value () / v2.complex_value ());
60  else
61  {
62  MatrixType typ = v2.matrix_type ();
63 
64  ComplexMatrix ret = xdiv (v1.matrix_value (),
66 
67  v2.matrix_type (typ);
68  return ret;
69  }
70 }
71 
72 DEFBINOPX (pow, matrix, sparse_complex_matrix)
73 {
74  error ("can't do A ^ B for A and B both matrices");
75 }
76 
77 DEFBINOP (ldiv, matrix, sparse_complex_matrix)
78 {
79  OCTAVE_CAST_BASE_VALUE (const octave_matrix&, v1, a1);
81  MatrixType typ = v1.matrix_type ();
82 
83  ComplexMatrix ret = xleftdiv (v1.matrix_value (),
84  v2.complex_matrix_value (), typ);
85 
86  v1.matrix_type (typ);
87  return ret;
88 }
89 
90 DEFBINOP_FN (lt, matrix, sparse_complex_matrix, mx_el_lt)
91 DEFBINOP_FN (le, matrix, sparse_complex_matrix, mx_el_le)
92 DEFBINOP_FN (eq, matrix, sparse_complex_matrix, mx_el_eq)
93 DEFBINOP_FN (ge, matrix, sparse_complex_matrix, mx_el_ge)
94 DEFBINOP_FN (gt, matrix, sparse_complex_matrix, mx_el_gt)
95 DEFBINOP_FN (ne, matrix, sparse_complex_matrix, mx_el_ne)
96 
97 DEFBINOP_FN (el_mul, matrix, sparse_complex_matrix, product)
98 DEFBINOP_FN (el_div, matrix, sparse_complex_matrix, quotient)
99 
100 DEFBINOP (el_pow, matrix, sparse_complex_matrix)
101 {
102  OCTAVE_CAST_BASE_VALUE (const octave_matrix&, v1, a1);
104 
105  return octave_value
106  (elem_xpow (SparseMatrix (v1.matrix_value ()),
108 }
109 
110 DEFBINOP (el_ldiv, matrix, sparse_complex_matrix)
111 {
112  OCTAVE_CAST_BASE_VALUE (const octave_matrix&, v1, a1);
114  return octave_value
116 }
117 
118 DEFBINOP_FN (el_and, matrix, sparse_complex_matrix, mx_el_and)
119 DEFBINOP_FN (el_or, matrix, sparse_complex_matrix, mx_el_or)
120 
121 DEFCATOP (m_scm, matrix, sparse_complex_matrix)
122 {
123  OCTAVE_CAST_BASE_VALUE (const octave_matrix&, v1, a1);
125  SparseMatrix tmp (v1.matrix_value ());
127  ra_idx));
128 }
129 
130 DEFCONV (sparse_complex_matrix_conv, matrix, sparse_complex_matrix)
131 {
132  OCTAVE_CAST_BASE_VALUE (const octave_matrix&, v, a);
134  (SparseComplexMatrix (v.complex_matrix_value ()));
135 }
136 
137 void
138 install_m_scm_ops (octave::type_info& ti)
139 {
146  ldiv);
154  el_mul);
156  el_div);
158  el_pow);
160  el_ldiv);
162  el_and);
164  el_or);
165 
167 
170 
172  sparse_complex_matrix_conv);
173 }
174 
175 OCTAVE_END_NAMESPACE(octave)
ComplexColumnVector quotient(const ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:158
ComplexColumnVector product(const ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:158
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
ComplexMatrix matrix_value() const
Definition: CSparse.cc:600
MatrixType matrix_type() const
Definition: ov-base-mat.h:139
octave_idx_type rows() const
Definition: ov-base.h:374
octave_idx_type columns() const
Definition: ov-base.h:381
virtual SparseComplexMatrix sparse_complex_matrix_value(bool=false) const
Definition: ov-base.cc:701
ComplexMatrix complex_matrix_value(bool=false) const
Definition: ov-ch-mat.h:130
Complex complex_value(bool=false) const
Definition: ov-ch-mat.cc:127
OCTAVE_BEGIN_NAMESPACE(octave) static octave_value daspk_fcn
void() error(const char *fmt,...)
Definition: error.cc:988
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_m_scm_ops(octave::type_info &ti)
Definition: op-m-scm.cc:138
#define DEFBINOPX(name, t1, t2)
Definition: ops.h:240
#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 DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:250
#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 INSTALL_WIDENOP_TI(ti, t1, t2, f)
Definition: ops.h:92
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:75
#define DEFCATOP(name, t1, t2)
Definition: ops.h:344
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 xdiv(const Matrix &a, const SparseMatrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:137
Matrix xleftdiv(const SparseMatrix &a, const Matrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:469
octave_value elem_xpow(double a, const SparseMatrix &b)
Definition: sparse-xpow.cc:347