GNU Octave  4.4.1
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-cs-s.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 "errwarn.h"
28 #include "ovl.h"
29 #include "ov.h"
30 #include "ov-complex.h"
31 #include "ov-cx-mat.h"
32 #include "ov-scalar.h"
33 #include "ov-typeinfo.h"
34 #include "ops.h"
35 #include "xdiv.h"
36 #include "xpow.h"
37 
38 // complex scalar by scalar ops.
39 
40 DEFBINOP_OP (add, complex, scalar, +)
41 DEFBINOP_OP (sub, complex, scalar, -)
42 DEFBINOP_OP (mul, complex, scalar, *)
43 
44 DEFBINOP (div, complex, scalar)
45 {
46  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
47  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
48 
49  double d = v2.double_value ();
50 
51  if (d == 0.0)
53 
54  return octave_value (v1.complex_value () / d);
55 }
56 
57 DEFBINOP_FN (pow, complex, scalar, xpow)
58 
59 DEFBINOP (ldiv, complex, scalar)
60 {
61  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
62  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
63 
64  Complex d = v1.complex_value ();
65 
66  if (d == 0.0)
68 
69  return octave_value (v2.double_value () / d);
70 }
71 
72 DEFCMPLXCMPOP_OP (lt, complex, scalar, <)
73 DEFCMPLXCMPOP_OP (le, complex, scalar, <=)
74 DEFCMPLXCMPOP_OP (eq, complex, scalar, ==)
75 DEFCMPLXCMPOP_OP (ge, complex, scalar, >=)
76 DEFCMPLXCMPOP_OP (gt, complex, scalar, >)
77 DEFCMPLXCMPOP_OP (ne, complex, scalar, !=)
78 
79 DEFBINOP_OP (el_mul, complex, scalar, *)
80 
81 DEFBINOP (el_div, complex, scalar)
82 {
83  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
84  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
85 
86  double d = v2.double_value ();
87 
88  if (d == 0.0)
90 
91  return octave_value (v1.complex_value () / d);
92 }
93 
94 DEFBINOP_FN (el_pow, complex, scalar, xpow)
95 
96 DEFBINOP (el_ldiv, complex, scalar)
97 {
98  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
99  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
100 
101  Complex d = v1.complex_value ();
102 
103  if (d == 0.0)
105 
106  return octave_value (v2.double_value () / d);
107 }
108 
109 DEFBINOP (el_and, complex, scalar)
110 {
111  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
112  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
113 
114  return v1.complex_value () != 0.0 && v2.double_value ();
115 }
116 
117 DEFBINOP (el_or, complex, scalar)
118 {
119  const octave_complex& v1 = dynamic_cast<const octave_complex&> (a1);
120  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
121 
122  return v1.complex_value () != 0.0 || v2.double_value ();
123 }
124 
125 DEFNDCATOP_FN (cs_s, complex, scalar, complex_array, array, concat)
126 
127 void
129 {
148 
150 
153 }
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
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1613
nd group nd example oindent but is performed more efficiently If only and it is a scalar
Definition: data.cc:5264
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
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1608
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
#define DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:269
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
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:58
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_s_ops(octave::type_info &ti)
Definition: op-cs-s.cc:128
double double_value(bool=false) const
Definition: ov-ch-mat.cc:66
void warn_divide_by_zero(void)
Definition: errwarn.cc:326
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
Complex complex_value(bool=false) const
Definition: ov-complex.cc:225
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1616
#define DEFCMPLXCMPOP_OP(name, t1, t2, op)
Definition: ops.h:281
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
octave_value op_sub(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1604
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 DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:382
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