GNU Octave  4.4.1
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-cs-m.cc
Go to the documentation of this file.
1 /*
2 
3 Copyright (C) 1996-2018 John W. Eaton
4 
5 This file is part of Octave.
6 
7 Octave is free software: you can redistribute it and/or modify it
8 under the terms of the GNU General Public License as published by
9 the Free Software Foundation, either version 3 of the License, or
10 (at your option) any later version.
11 
12 Octave is distributed in the hope that it will be useful, but
13 WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU General Public License for more details.
16 
17 You should have received a copy of the GNU General Public License
18 along with Octave; see the file COPYING. If not, see
19 <https://www.gnu.org/licenses/>.
20 
21 */
22 
23 #if defined (HAVE_CONFIG_H)
24 # include "config.h"
25 #endif
26 
27 #include "mx-cs-nda.h"
28 #include "mx-nda-cs.h"
29 #include "mx-cs-nda.h"
30 #include "mx-nda-cs.h"
31 
32 #include "errwarn.h"
33 #include "ovl.h"
34 #include "ov.h"
35 #include "ov-complex.h"
36 #include "ov-cx-mat.h"
37 #include "ov-re-mat.h"
38 #include "ov-typeinfo.h"
39 #include "ops.h"
40 #include "xdiv.h"
41 #include "xpow.h"
42 
43 // complex scalar by matrix ops.
44 
45 DEFNDBINOP_OP (add, complex, matrix, complex, array, +)
46 DEFNDBINOP_OP (sub, complex, matrix, complex, array, -)
47 DEFNDBINOP_OP (mul, complex, matrix, complex, array, *)
48 
49 DEFBINOP (div, complex, matrix)
50 {
51  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
52  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
53 
55  Matrix m2 = v2.matrix_value ();
56  MatrixType typ = v2.matrix_type ();
57 
58  ComplexMatrix ret = xdiv (m1, m2, typ);
59 
60  v2.matrix_type (typ);
61  return ret;
62 }
63 
64 DEFBINOP_FN (pow, complex, matrix, xpow)
65 
66 DEFBINOP (ldiv, complex, matrix)
67 {
68  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
69  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
70 
71  Complex d = v1.complex_value ();
72 
73  if (d == 0.0)
75 
76  return octave_value (v2.array_value () / d);
77 }
78 
79 DEFNDCMPLXCMPOP_FN (lt, complex, matrix, complex, array, mx_el_lt)
80 DEFNDCMPLXCMPOP_FN (le, complex, matrix, complex, array, mx_el_le)
81 DEFNDCMPLXCMPOP_FN (eq, complex, matrix, complex, array, mx_el_eq)
82 DEFNDCMPLXCMPOP_FN (ge, complex, matrix, complex, array, mx_el_ge)
83 DEFNDCMPLXCMPOP_FN (gt, complex, matrix, complex, array, mx_el_gt)
84 DEFNDCMPLXCMPOP_FN (ne, complex, matrix, complex, array, mx_el_ne)
85 
86 DEFNDBINOP_OP (el_mul, complex, matrix, complex, array, *)
87 DEFNDBINOP_FN (el_div, complex, matrix, complex, array, x_el_div)
88 DEFNDBINOP_FN (el_pow, complex, matrix, complex, array, elem_xpow)
89 
90 DEFBINOP (el_ldiv, complex, matrix)
91 {
92  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
93  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
94 
95  Complex d = v1.complex_value ();
96 
97  if (d == 0.0)
99 
100  return octave_value (v2.array_value () / d);
101 }
102 
103 DEFNDBINOP_FN (el_and, complex, matrix, complex, array, mx_el_and)
104 DEFNDBINOP_FN (el_or, complex, matrix, complex, array, mx_el_or)
105 
106 DEFNDCATOP_FN (cs_m, complex, matrix, complex_array, array, concat)
107 
108 void
109 install_cs_m_ops (octave::type_info& ti)
110 {
129 
131 
134 }
ComplexMatrix complex_matrix_value(bool=false) const
Definition: ov-complex.cc:237
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:58
#define DEFBINOP(name, t1, t2)
Definition: ops.h:264
octave_value op_el_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1620
#define INSTALL_BINOP_TI(ti, op, t1, t2, f)
Definition: ops.h:53
characters Given a string matrix
Definition: hex2num.cc:155
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1613
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1621
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:653
boolMatrix mx_el_le(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:89
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1608
#define DEFNDBINOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:335
F77_RET_T const F77_REAL const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
const octave_base_value & a2
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1606
then the function must return scalars which will be concatenated into the return array(s). If code
Definition: cellfun.cc:400
octave_value op_el_or(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1623
boolMatrix mx_el_ge(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:89
Matrix x_el_div(double a, const SparseMatrix &b)
Definition: sparse-xdiv.cc:371
boolMatrix mx_el_gt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:89
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:58
octave_value elem_xpow(double a, const SparseMatrix &b)
Definition: sparse-xpow.cc:247
octave_value op_el_and(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1622
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
void install_cs_m_ops(octave::type_info &ti)
Definition: op-cs-m.cc:109
boolMatrix mx_el_ne(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:89
void warn_divide_by_zero(void)
Definition: errwarn.cc:326
Definition: dMatrix.h:36
octave_value op_le(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1612
octave_value op_lt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1611
octave_value op_el_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1619
#define INSTALL_ASSIGNCONV_TI(ti, t1, t2, tr)
Definition: ops.h:71
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
boolMatrix mx_el_or(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:86
NDArray array_value(bool=false) const
Definition: ov-ch-mat.h:118
Complex complex_value(bool=false) const
Definition: ov-complex.cc:225
Matrix xdiv(const Matrix &a, const SparseMatrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:133
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1616
boolMatrix mx_el_and(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:86
MatrixType matrix_type(void) const
Definition: ov-base-mat.h:124
octave_value op_add(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1603
octave_value op_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1609
Matrix matrix_value(bool=false) const
Definition: ov-ch-mat.h:112
#define DEFNDCMPLXCMPOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:346
octave_value op_sub(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1604
boolMatrix mx_el_lt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:89
octave_value op_el_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1618
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:324
std::complex< double > Complex
Definition: oct-cmplx.h:31
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1614
#define DEFNDBINOP_OP(name, t1, t2, e1, e2, op)
Definition: ops.h:310
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:382
boolMatrix mx_el_eq(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:89
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1605
octave_value op_gt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1615