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-fcs-fcs.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 "gripes.h"
28 #include "oct-obj.h"
29 #include "ov.h"
30 #include "ov-complex.h"
31 #include "ov-flt-complex.h"
32 #include "ov-flt-cx-mat.h"
33 #include "ov-typeinfo.h"
34 #include "ov-null-mat.h"
35 #include "ops.h"
36 #include "xdiv.h"
37 #include "xpow.h"
38 
39 // unary complex scalar ops.
40 
41 DEFUNOP (not, float_complex)
42 {
44  FloatComplex x = v.float_complex_value ();
45  if (xisnan (x))
47  return octave_value (x == 0.0f);
48 }
49 
50 DEFUNOP_OP (uplus, float_complex, /* no-op */)
51 DEFUNOP_OP (uminus, float_complex, -)
52 DEFUNOP_OP (transpose, float_complex, /* no-op */)
53 
54 DEFUNOP (hermitian, float_complex)
55 {
57 
58  return octave_value (conj (v.float_complex_value ()));
59 }
60 
61 DEFNCUNOP_METHOD (incr, float_complex, increment)
62 DEFNCUNOP_METHOD (decr, float_complex, decrement)
63 
64 // complex scalar by complex scalar ops.
65 
66 DEFBINOP_OP (add, float_complex, float_complex, +)
67 DEFBINOP_OP (sub, float_complex, float_complex, -)
68 DEFBINOP_OP (mul, float_complex, float_complex, *)
69 
70 DEFBINOP (div, float_complex, float_complex)
71 {
73 
75 
76  if (d == 0.0f)
78 
79  return octave_value (v1.float_complex_value () / d);
80 }
81 
82 DEFBINOP_FN (pow, float_complex, float_complex, xpow)
83 
84 DEFBINOP (ldiv, float_complex, float_complex)
85 {
87 
89 
90  if (d == 0.0f)
92 
93  return octave_value (v2.float_complex_value () / d);
94 }
95 
96 DEFCMPLXCMPOP_OP (lt, float_complex, float_complex, <)
97 DEFCMPLXCMPOP_OP (le, float_complex, float_complex, <=)
98 DEFCMPLXCMPOP_OP (eq, float_complex, float_complex, ==)
99 DEFCMPLXCMPOP_OP (ge, float_complex, float_complex, >=)
100 DEFCMPLXCMPOP_OP (gt, float_complex, float_complex, >)
101 DEFCMPLXCMPOP_OP (ne, float_complex, float_complex, !=)
102 
103 DEFBINOP_OP (el_mul, float_complex, float_complex, *)
104 
105 DEFBINOP (el_div, float_complex, float_complex)
106 {
108 
110 
111  if (d == 0.0f)
113 
114  return octave_value (v1.float_complex_value () / d);
115 }
116 
117 DEFBINOP_FN (el_pow, float_complex, float_complex, xpow)
118 
119 DEFBINOP (el_ldiv, float_complex, float_complex)
120 {
122 
124 
125  if (d == 0.0f)
127 
128  return octave_value (v2.float_complex_value () / d);
129 }
130 
131 DEFBINOP (el_and, float_complex, float_complex)
132 {
134 
135  return (v1.float_complex_value () != 0.0f
136  && v2.float_complex_value () != 0.0f);
137 }
138 
139 DEFBINOP (el_or, float_complex, float_complex)
140 {
142 
143  return (v1.float_complex_value () != 0.0f
144  || v2.float_complex_value () != 0.0f);
145 }
146 
147 DEFNDCATOP_FN (fcs_fcs, float_complex, float_complex, float_complex_array,
148  float_complex_array, concat)
149 
150 DEFNDCATOP_FN (cs_fcs, complex, float_complex, float_complex_array,
151  float_complex_array, concat)
152 
153 DEFNDCATOP_FN (fcs_cs, float_complex, complex, float_complex_array,
154  float_complex_array, concat)
155 
156 CONVDECL (float_complex_to_complex)
157 {
159 
160  return new octave_complex_matrix
161  (ComplexMatrix (1, 1,
162  static_cast<Complex>(v.float_complex_value ())));
163 }
164 
165 void
167 {
173 
174  INSTALL_NCUNOP (op_incr, octave_float_complex, incr);
175  INSTALL_NCUNOP (op_decr, octave_float_complex, decr);
176 
193  el_ldiv);
196 
200 
203 
206 
213 
215  float_complex_to_complex);
216 }
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
FloatComplex float_complex_value(bool=false) const
Definition: ov-ch-mat.cc:111
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
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
F77_RET_T const double const double * f
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)
void install_fcs_fcs_ops(void)
Definition: op-fcs-fcs.cc:166
virtual FloatComplex float_complex_value(bool=false) const
Definition: ov-base.cc:583
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
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
std::complex< float > FloatComplex
Definition: oct-cmplx.h:30
#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
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