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-fm-fs.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-flt-re-mat.h"
31 #include "ov-float.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 // matrix by scalar ops.
39 
40 DEFNDBINOP_OP (add, float_matrix, float_scalar, float_array, float_scalar, +)
41 DEFNDBINOP_OP (sub, float_matrix, float_scalar, float_array, float_scalar, -)
42 DEFNDBINOP_OP (mul, float_matrix, float_scalar, float_array, float_scalar, *)
43 
44 DEFBINOP (div, float_matrix, float)
45 {
46  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
47  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
48 
49  float d = v2.float_value ();
50 
51  if (d == 0.0)
53 
54  return octave_value (v1.float_array_value () / d);
55 }
56 
57 DEFBINOP_FN (pow, float_matrix, float_scalar, xpow)
58 
59 DEFBINOP (ldiv, float_matrix, float)
60 {
61  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
62  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
63 
66  MatrixType typ = v1.matrix_type ();
67 
68  FloatMatrix ret = xleftdiv (m1, m2, typ);
69 
70  v1.matrix_type (typ);
71  return ret;
72 }
73 
74 DEFNDBINOP_FN (lt, float_matrix, float_scalar, float_array,
75  float_scalar, mx_el_lt)
76 DEFNDBINOP_FN (le, float_matrix, float_scalar, float_array,
77  float_scalar, mx_el_le)
78 DEFNDBINOP_FN (eq, float_matrix, float_scalar, float_array,
79  float_scalar, mx_el_eq)
80 DEFNDBINOP_FN (ge, float_matrix, float_scalar, float_array,
81  float_scalar, mx_el_ge)
82 DEFNDBINOP_FN (gt, float_matrix, float_scalar, float_array,
83  float_scalar, mx_el_gt)
84 DEFNDBINOP_FN (ne, float_matrix, float_scalar, float_array,
85  float_scalar, mx_el_ne)
86 
87 DEFNDBINOP_OP (el_mul, float_matrix, float_scalar, float_array, float_scalar, *)
88 
89 DEFBINOP (el_div, float_matrix, float)
90 {
91  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
92  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
93 
94  float d = v2.float_value ();
95 
96  if (d == 0.0)
98 
99  return octave_value (v1.float_array_value () / d);
100 }
101 
102 DEFNDBINOP_FN (el_pow, float_matrix, float_scalar, float_array,
103  float_scalar, elem_xpow)
104 
105 DEFBINOP (el_ldiv, float_matrix, float)
106 {
107  const octave_float_matrix& v1 = dynamic_cast<const octave_float_matrix&> (a1);
108  const octave_float_scalar& v2 = dynamic_cast<const octave_float_scalar&> (a2);
109 
110  return x_el_div (v2.float_value (), v1.float_array_value ());
111 }
112 
113 DEFNDBINOP_FN (el_and, float_matrix, float_scalar, float_array,
114  float_scalar, mx_el_and)
115 DEFNDBINOP_FN (el_or, float_matrix, float_scalar, float_array,
116  float_scalar, mx_el_or)
117 
118 DEFNDCATOP_FN (fm_fs, float_matrix, float_scalar, float_array,
119  float_array, concat)
120 
121 DEFNDCATOP_FN (m_fs, matrix, float_scalar, float_array, float_array, concat)
122 
123 DEFNDCATOP_FN (fm_s, float_matrix, scalar, float_array, float_array, concat)
124 
125 DEFNDASSIGNOP_FN (assign, float_matrix, float_scalar, float_scalar, assign)
126 DEFNDASSIGNOP_FN (dbl_assign, matrix, float_scalar, scalar, assign)
127 
128 DEFNDASSIGNOP_OP (assign_add, float_matrix, float_scalar, float_scalar, +=)
129 DEFNDASSIGNOP_OP (assign_sub, float_matrix, float_scalar, float_scalar, -=)
130 DEFNDASSIGNOP_OP (assign_mul, float_matrix, float_scalar, float_scalar, *=)
131 DEFNDASSIGNOP_OP (assign_div, float_matrix, float_scalar, float_scalar, /=)
132 
133 void
135 {
154 
158 
160  assign);
161  INSTALL_ASSIGNOP (op_asn_eq, octave_matrix, octave_float_scalar, dbl_assign);
162 
164  assign_add);
166  assign_sub);
168  assign_mul);
170  assign_div);
171 }
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
the sparsity preserving column transformation such that that defines the pivoting threshold can be given in which case it defines the then the first element defines the pivoting tolerance for the unsymmetric the values defined such that for full matrix
Definition: lu.cc:138
ComplexNDArray concat(NDArray &ra, ComplexNDArray &rb, const Array< octave_idx_type > &ra_idx)
Definition: CNDArray.cc:655
boolMatrix mx_el_le(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1519
FloatMatrix float_matrix_value(bool=false) const
FloatMatrix float_matrix_value(bool=false) const
Definition: ov-float.h:159
float float_value(bool=false) const
Definition: ov-float.h:149
#define DEFNDBINOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:294
#define INSTALL_ASSIGNOP(op, t1, t2, f)
Definition: ops.h:57
MatrixType matrix_type(void) const
Definition: ov-base-mat.h:129
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
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
boolMatrix mx_el_ge(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
Matrix x_el_div(double a, const SparseMatrix &b)
Definition: sparse-xdiv.cc:371
boolMatrix mx_el_gt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
octave_value elem_xpow(double a, const SparseMatrix &b)
Definition: sparse-xpow.cc:247
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)
Matrix xleftdiv(const SparseMatrix &a, const Matrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:464
boolMatrix mx_el_ne(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
FloatNDArray float_array_value(bool=false) const
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
const octave_char_matrix & v2
boolMatrix mx_el_or(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
#define DEFNDASSIGNOP_FN(name, t1, t2, e, f)
Definition: ops.h:105
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1527
void install_fm_fs_ops(void)
Definition: op-fm-fs.cc:134
boolMatrix mx_el_and(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
#define DEFNDASSIGNOP_OP(name, t1, t2, f, op)
Definition: ops.h:119
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
boolMatrix mx_el_lt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
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
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1525
#define DEFNDBINOP_OP(name, t1, t2, e1, e2, op)
Definition: ops.h:269
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:341
boolMatrix mx_el_eq(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
octave_value op_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1516
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())? '\'': '"'))