GNU Octave  4.4.1
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
op-fs-fs.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 "Array-util.h"
28 
29 #include "errwarn.h"
30 #include "ovl.h"
31 #include "ov.h"
32 #include "ov-scalar.h"
33 #include "ov-float.h"
34 #include "ov-flt-re-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 // scalar unary ops.
42 
43 DEFUNOP (not, float_scalar)
44 {
45  const octave_float_scalar& v = dynamic_cast<const octave_float_scalar&> (a);
46  float x = v.float_value ();
47  if (octave::math::isnan (x))
49 
50  return octave_value (x == 0.0f);
51 }
52 
53 DEFUNOP_OP (uplus, float_scalar, /* no-op */)
54 DEFUNOP_OP (uminus, float_scalar, -)
55 DEFUNOP_OP (transpose, float_scalar, /* no-op */)
56 DEFUNOP_OP (hermitian, float_scalar, /* no-op */)
57 
58 DEFNCUNOP_METHOD (incr, float_scalar, increment)
59 DEFNCUNOP_METHOD (decr, float_scalar, decrement)
60 
61 // float by float ops.
62 
63 DEFBINOP_OP (add, float_scalar, float_scalar, +)
64 DEFBINOP_OP (sub, float_scalar, float_scalar, -)
65 DEFBINOP_OP (mul, float_scalar, float_scalar, *)
66 
67 DEFBINOP (div, float_scalar, float_scalar)
68 {
69  const octave_float_scalar& v1 = dynamic_cast<const octave_float_scalar&> (a1);
70  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
71 
72  float d = v2.float_value ();
73 
74  if (d == 0.0)
76 
77  return octave_value (v1.float_value () / d);
78 }
79 
80 DEFBINOP_FN (pow, float_scalar, float_scalar, xpow)
81 
82 DEFBINOP (ldiv, float_scalar, float_scalar)
83 {
84  const octave_float_scalar& v1 = dynamic_cast<const octave_float_scalar&> (a1);
85  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
86 
87  float d = v1.float_value ();
88 
89  if (d == 0.0)
91 
92  return octave_value (v2.float_value () / d);
93 }
94 
95 DEFBINOP_OP (lt, float_scalar, float_scalar, <)
96 DEFBINOP_OP (le, float_scalar, float_scalar, <=)
97 DEFBINOP_OP (eq, float_scalar, float_scalar, ==)
98 DEFBINOP_OP (ge, float_scalar, float_scalar, >=)
99 DEFBINOP_OP (gt, float_scalar, float_scalar, >)
100 DEFBINOP_OP (ne, float_scalar, float_scalar, !=)
101 
102 DEFBINOP_OP (el_mul, float_scalar, float_scalar, *)
103 
104 DEFBINOP (el_div, float_scalar, float_scalar)
105 {
106  const octave_float_scalar& v1 = dynamic_cast<const octave_float_scalar&> (a1);
107  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
108 
109  float d = v2.float_value ();
110 
111  if (d == 0.0)
113 
114  return octave_value (v1.float_value () / d);
115 }
116 
117 DEFBINOP_FN (el_pow, float_scalar, float_scalar, xpow)
118 
119 DEFBINOP (el_ldiv, float_scalar, float_scalar)
120 {
121  const octave_float_scalar& v1 = dynamic_cast<const octave_float_scalar&> (a1);
122  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
123 
124  float d = v1.float_value ();
125 
126  if (d == 0.0)
128 
129  return octave_value (v2.float_value () / d);
130 }
131 
132 DEFSCALARBOOLOP_OP (el_and, float_scalar, float_scalar, &&)
133 DEFSCALARBOOLOP_OP (el_or, float_scalar, float_scalar, ||)
134 
135 DEFNDCATOP_FN (fs_fs, float_scalar, float_scalar, float_array, float_array,
137 DEFNDCATOP_FN (s_fs, scalar, float_scalar, float_array, float_array, concat)
138 DEFNDCATOP_FN (fs_s, float_scalar, scalar, float_array, float_array, concat)
139 
140 void
141 install_fs_fs_ops (octave::type_info& ti)
142 {
148 
149  INSTALL_NCUNOP_TI (ti, op_incr, octave_float_scalar, incr);
150  INSTALL_NCUNOP_TI (ti, op_decr, octave_float_scalar, decr);
151 
165  el_mul);
167  el_div);
169  el_pow);
171  el_ldiv);
173  el_and);
175  el_or);
176 
180 
184 
191 }
octave_value xpow(const SparseMatrix &a, double b)
Definition: sparse-xpow.cc:58
octave_value op_uplus(const octave_value &a)
Definition: ov.h:1574
#define DEFBINOP(name, t1, t2)
Definition: ops.h:264
#define DEFUNOP(name, t)
Definition: ops.h:213
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
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 * f
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1621
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:386
bool isnan(bool)
Definition: lo-mappers.h:187
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
#define INSTALL_NCUNOP_TI(ti, op, t, f)
Definition: ops.h:49
#define DEFUNOP_OP(name, t, op)
Definition: ops.h:217
void err_nan_to_logical_conversion(void)
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
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
Definition: cellfun.cc:400
#define INSTALL_UNOP_TI(ti, op, t, f)
Definition: ops.h:45
octave_value op_div(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1606
octave_value op_el_or(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1623
octave_value op_not(const octave_value &a)
Definition: ov.h:1573
#define INSTALL_CATOP_TI(ti, t1, t2, f)
Definition: ops.h:58
octave_value op_transpose(const octave_value &a)
Definition: ov.h:1577
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_fs_fs_ops(octave::type_info &ti)
Definition: op-fs-fs.cc:141
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 DEFSCALARBOOLOP_OP(name, t1, t2, op)
Definition: ops.h:295
#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
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1616
float float_value(bool=false) const
Definition: ov-float.h:147
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
octave_value op_hermitian(const octave_value &a)
Definition: ov.h:1578
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1614
float float_value(bool=false) const
Definition: ov-ch-mat.cc:78
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:382
octave_value op_uminus(const octave_value &a)
Definition: ov.h:1575
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
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 * x
#define DEFNCUNOP_METHOD(name, t, method)
Definition: ops.h:251