GNU Octave  4.2.1
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
op-cs-s.cc
Go to the documentation of this file.
1 /*
2 
3 Copyright (C) 1996-2017 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 the
9 Free Software Foundation; either version 3 of the License, or (at your
10 option) any later version.
11 
12 Octave is distributed in the hope that it will be useful, but WITHOUT
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
15 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 <http://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 
152 }
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:58
#define DEFBINOP(name, t1, t2)
Definition: ops.h:223
octave_value op_el_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1531
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1524
nd group nd example oindent but is performed more efficiently If only and it is a scalar
Definition: data.cc:5342
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1532
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:655
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1519
F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &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 F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T F77_REAL F77_REAL &F77_RET_T F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
#define DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:228
const octave_base_value & a2
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1517
octave_value op_el_or(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1534
#define INSTALL_BINOP(op, t1, t2, f)
Definition: ops.h:48
then the function must return scalars which will be concatenated into the return array(s).If code
Definition: cellfun.cc:398
octave_value op_el_and(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1533
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
void warn_divide_by_zero(void)
Definition: errwarn.cc:320
octave_value op_le(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1523
octave_value op_lt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1522
#define INSTALL_CATOP(t1, t2, f)
Definition: ops.h:53
octave_value op_el_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1530
double double_value(bool=false) const
Definition: ov-scalar.h:145
const octave_char_matrix & v2
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1527
#define DEFCMPLXCMPOP_OP(name, t1, t2, op)
Definition: ops.h:240
octave_value op_add(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1514
octave_value op_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1520
octave_value op_sub(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1515
octave_value op_el_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1529
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:283
std::complex< double > Complex
Definition: oct-cmplx.h:31
#define INSTALL_ASSIGNCONV(t1, t2, tr)
Definition: ops.h:66
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1525
void install_cs_s_ops(void)
Definition: op-cs-s.cc:128
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:341
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1516
Complex complex_value(bool=false) const
Definition: ov-complex.cc:225
octave_value op_gt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1526
return octave_value(v1.char_array_value().concat(v2.char_array_value(), ra_idx),((a1.is_sq_string()||a2.is_sq_string())? '\'': '"'))