GNU Octave  4.0.0
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-cs.cc
Go to the documentation of this file.
1 /*
2 
3 Copyright (C) 1996-2015 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 #ifdef HAVE_CONFIG_H
24 #include <config.h>
25 #endif
26 
27 #include "Array-util.h"
28 
29 #include "gripes.h"
30 #include "oct-obj.h"
31 #include "ov.h"
32 #include "ov-complex.h"
33 #include "ov-cx-mat.h"
34 #include "ov-flt-cx-mat.h"
35 #include "ov-typeinfo.h"
36 #include "ov-null-mat.h"
37 #include "ops.h"
38 #include "xdiv.h"
39 #include "xpow.h"
40 
41 // unary complex scalar ops.
42 
43 DEFUNOP (not, complex)
44 {
46  Complex x = v.complex_value ();
47  if (xisnan (x))
49  return octave_value (x == 0.0);
50 }
51 
52 DEFUNOP_OP (uplus, complex, /* no-op */)
53 DEFUNOP_OP (uminus, complex, -)
54 DEFUNOP_OP (transpose, complex, /* no-op */)
55 
56 DEFUNOP (hermitian, complex)
57 {
59 
60  return octave_value (conj (v.complex_value ()));
61 }
62 
63 DEFNCUNOP_METHOD (incr, complex, increment)
64 DEFNCUNOP_METHOD (decr, complex, decrement)
65 
66 // complex scalar by complex scalar ops.
67 
68 DEFBINOP_OP (add, complex, complex, +)
69 DEFBINOP_OP (sub, complex, complex, -)
70 DEFBINOP_OP (mul, complex, complex, *)
71 
72 DEFBINOP (div, complex, complex)
73 {
75 
77 
78  if (d == 0.0)
80 
81  return octave_value (v1.complex_value () / d);
82 }
83 
84 DEFBINOP_FN (pow, complex, complex, xpow)
85 
86 DEFBINOP (ldiv, complex, complex)
87 {
89 
91 
92  if (d == 0.0)
94 
95  return octave_value (v2.complex_value () / d);
96 }
97 
98 DEFCMPLXCMPOP_OP (lt, complex, complex, <)
99 DEFCMPLXCMPOP_OP (le, complex, complex, <=)
100 DEFCMPLXCMPOP_OP (eq, complex, complex, ==)
101 DEFCMPLXCMPOP_OP (ge, complex, complex, >=)
102 DEFCMPLXCMPOP_OP (gt, complex, complex, >)
103 DEFCMPLXCMPOP_OP (ne, complex, complex, !=)
104 
105 DEFBINOP_OP (el_mul, complex, complex, *)
106 
107 DEFBINOP (el_div, complex, complex)
108 {
110 
111  Complex d = v2.complex_value ();
112 
113  if (d == 0.0)
115 
116  return octave_value (v1.complex_value () / d);
117 }
118 
119 DEFBINOP_FN (el_pow, complex, complex, xpow)
120 
121 DEFBINOP (el_ldiv, complex, complex)
122 {
124 
125  Complex d = v1.complex_value ();
126 
127  if (d == 0.0)
129 
130  return octave_value (v2.complex_value () / d);
131 }
132 
133 DEFBINOP (el_and, complex, complex)
134 {
136 
137  return v1.complex_value () != 0.0 && v2.complex_value () != 0.0;
138 }
139 
140 DEFBINOP (el_or, complex, complex)
141 {
143 
144  return v1.complex_value () != 0.0 || v2.complex_value () != 0.0;
145 }
146 
147 DEFNDCATOP_FN (cs_cs, complex, complex, complex_array, complex_array, concat)
148 
149 CONVDECL (complex_to_float_complex)
150 {
151  CAST_CONV_ARG (const octave_complex&);
152 
154  static_cast<FloatComplex>
155  (v.complex_value ())));
156 }
157 
158 void
160 {
166 
167  INSTALL_NCUNOP (op_incr, octave_complex, incr);
168  INSTALL_NCUNOP (op_decr, octave_complex, decr);
169 
188 
190 
192 
198 
200  complex_to_float_complex);
201 }
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:59
octave_value op_uplus(const octave_value &a)
Definition: ov.h:1294
#define DEFBINOP(name, t1, t2)
Definition: ops.h:279
#define DEFUNOP(name, t)
Definition: ops.h:231
octave_value op_el_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1342
#define INSTALL_NCUNOP(op, t, f)
Definition: ops.h:42
bool xisnan(double x)
Definition: lo-mappers.cc:144
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1335
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1343
static void transpose(octave_idx_type N, const octave_idx_type *ridx, const octave_idx_type *cidx, octave_idx_type *ridx2, octave_idx_type *cidx2)
Definition: symrcm.cc:382
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:664
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1328
const octave_base_value const Array< octave_idx_type > &ra_idx octave_int16_scalar & v1
#define DEFUNOP_OP(name, t, op)
Definition: ops.h:234
F77_RET_T const double const double double * d
#define CAST_BINOP_ARGS(t1, t2)
Definition: ops.h:79
#define INSTALL_CONVOP(t1, t2, f)
Definition: ops.h:68
#define DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:282
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1326
virtual Complex complex_value(bool=false) const
Definition: ov-base.cc:574
octave_value op_el_or(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1345
ComplexColumnVector conj(const ComplexColumnVector &a)
Definition: CColVector.cc:244
#define INSTALL_BINOP(op, t1, t2, f)
Definition: ops.h:46
#define CAST_CONV_ARG(t)
Definition: ops.h:83
octave_value op_not(const octave_value &a)
Definition: ov.h:1293
octave_value op_transpose(const octave_value &a)
Definition: ov.h:1297
void gripe_nan_to_logical_conversion(void)
octave_value op_el_and(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1344
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
octave_value op_le(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1334
octave_value op_lt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1333
#define INSTALL_CATOP(t1, t2, f)
Definition: ops.h:51
octave_value op_el_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1341
void install_cs_cs_ops(void)
Definition: op-cs-cs.cc:159
const octave_char_matrix & v2
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1338
#define DEFCMPLXCMPOP_OP(name, t1, t2, op)
Definition: ops.h:290
octave_value op_add(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1323
#define INSTALL_UNOP(op, t, f)
Definition: ops.h:38
octave_value op_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1329
octave_value op_sub(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1324
Complex complex_value(bool=false) const
Definition: ov-ch-mat.cc:91
#define CAST_UNOP_ARG(t)
Definition: ops.h:76
octave_value op_el_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1340
#define DEFBINOP_FN(name, t1, t2, f)
Definition: ops.h:323
std::complex< double > Complex
Definition: oct-cmplx.h:29
octave_value op_hermitian(const octave_value &a)
Definition: ov.h:1298
#define INSTALL_ASSIGNCONV(t1, t2, tr)
Definition: ops.h:64
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1336
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:364
octave_value op_uminus(const octave_value &a)
Definition: ov.h:1295
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1325
void gripe_divide_by_zero(void)
Definition: gripes.cc:195
octave_value op_gt(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1337
#define CONVDECL(name)
Definition: ops.h:166
return octave_value(v1.char_array_value().concat(v2.char_array_value(), ra_idx),((a1.is_sq_string()||a2.is_sq_string())? '\'': '"'))
F77_RET_T const double * x
#define DEFNCUNOP_METHOD(name, t, method)
Definition: ops.h:264