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-fm-fm.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-re-mat.h"
31 #include "ov-flt-re-mat.h"
32 #include "ov-typeinfo.h"
33 #include "ov-null-mat.h"
34 #include "ops.h"
35 #include "xdiv.h"
36 #include "xpow.h"
37 
38 // matrix unary ops.
39 
40 DEFNDUNOP_OP (not, float_matrix, float_array, !)
41 DEFNDUNOP_OP (uplus, float_matrix, float_array, /* no-op */)
42 DEFNDUNOP_OP (uminus, float_matrix, float_array, -)
43 
44 DEFUNOP (transpose, float_matrix)
45 {
47 
48  if (v.ndims () > 2)
49  {
50  error ("transpose not defined for N-d objects");
51  return octave_value ();
52  }
53  else
54  return octave_value (v.float_matrix_value ().transpose ());
55 }
56 
57 DEFNCUNOP_METHOD (incr, float_matrix, increment)
58 DEFNCUNOP_METHOD (decr, float_matrix, decrement)
59 DEFNCUNOP_METHOD (changesign, float_matrix, changesign)
60 
61 // matrix by matrix ops.
62 
63 DEFNDBINOP_OP (add, float_matrix, float_matrix, float_array, float_array, +)
64 DEFNDBINOP_OP (sub, float_matrix, float_matrix, float_array, float_array, -)
65 
66 DEFBINOP_OP (mul, float_matrix, float_matrix, *)
67 
68 DEFBINOP (div, float_matrix, float_matrix)
69 {
71  MatrixType typ = v2.matrix_type ();
72 
74  v2.float_matrix_value (), typ);
75 
76  v2.matrix_type (typ);
77  return ret;
78 }
79 
80 DEFBINOPX (pow, float_matrix, float_matrix)
81 {
82  error ("can't do A ^ B for A and B both matrices");
83  return octave_value ();
84 }
85 
86 DEFBINOP (ldiv, float_matrix, float_matrix)
87 {
89  MatrixType typ = v1.matrix_type ();
90 
92  v2.float_matrix_value (), typ);
93 
94  v1.matrix_type (typ);
95  return ret;
96 }
97 
98 DEFBINOP (trans_mul, float_matrix, float_matrix)
99 {
104 }
105 
106 DEFBINOP (mul_trans, float_matrix, float_matrix)
107 {
112 }
113 
114 DEFBINOP (trans_ldiv, float_matrix, float_matrix)
115 {
117  MatrixType typ = v1.matrix_type ();
118 
120  v2.float_matrix_value (), typ, blas_trans);
121 
122  v1.matrix_type (typ);
123  return ret;
124 }
125 
126 DEFNDBINOP_FN (lt, float_matrix, float_matrix, float_array,
127  float_array, mx_el_lt)
128 DEFNDBINOP_FN (le, float_matrix, float_matrix, float_array,
129  float_array, mx_el_le)
130 DEFNDBINOP_FN (eq, float_matrix, float_matrix, float_array,
131  float_array, mx_el_eq)
132 DEFNDBINOP_FN (ge, float_matrix, float_matrix, float_array,
133  float_array, mx_el_ge)
134 DEFNDBINOP_FN (gt, float_matrix, float_matrix, float_array,
135  float_array, mx_el_gt)
136 DEFNDBINOP_FN (ne, float_matrix, float_matrix, float_array,
137  float_array, mx_el_ne)
138 
139 DEFNDBINOP_FN (el_mul, float_matrix, float_matrix, float_array,
140  float_array, product)
141 DEFNDBINOP_FN (el_div, float_matrix, float_matrix, float_array,
142  float_array, quotient)
143 DEFNDBINOP_FN (el_pow, float_matrix, float_matrix, float_array,
144  float_array, elem_xpow)
145 
146 DEFBINOP (el_ldiv, float_matrix, float_matrix)
147 {
149 
151  v1.float_array_value ()));
152 }
153 
154 DEFNDBINOP_FN (el_and, float_matrix, float_matrix, float_array,
155  float_array, mx_el_and)
156 DEFNDBINOP_FN (el_or, float_matrix, float_matrix, float_array,
157  float_array, mx_el_or)
158 DEFNDBINOP_FN (el_not_and, float_matrix, float_matrix, float_array,
159  float_array, mx_el_not_and)
160 DEFNDBINOP_FN (el_not_or, float_matrix, float_matrix, float_array,
161  float_array, mx_el_not_or)
162 DEFNDBINOP_FN (el_and_not, float_matrix, float_matrix, float_array,
163  float_array, mx_el_and_not)
164 DEFNDBINOP_FN (el_or_not, float_matrix, float_matrix, float_array,
165  float_array, mx_el_or_not)
166 
167 
168 
169 DEFNDCATOP_FN (fm_fm, float_matrix, float_matrix, float_array,
170  float_array, concat)
171 
172 DEFNDCATOP_FN (m_fm, matrix, float_matrix, float_array, float_array, concat)
173 
174 DEFNDCATOP_FN (fm_m, float_matrix, matrix, float_array, float_array, concat)
175 
176 DEFNDASSIGNOP_FN (assign, float_matrix, float_matrix, float_array, assign)
177 
178 DEFNDASSIGNOP_FN (dbl_assign, matrix, float_matrix, array, assign)
179 
180 DEFNULLASSIGNOP_FN (null_assign, float_matrix, delete_elements)
181 
182 DEFNDASSIGNOP_OP (assign_add, float_matrix, float_matrix, float_array, +=)
183 DEFNDASSIGNOP_OP (assign_sub, float_matrix, float_matrix, float_array, -=)
184 DEFNDASSIGNOP_FNOP (assign_el_mul, float_matrix, float_matrix, float_array,
186 DEFNDASSIGNOP_FNOP (assign_el_div, float_matrix, float_matrix, float_array,
188 
189 CONVDECL (float_matrix_to_matrix)
190 {
192 
193  return new octave_matrix (v.array_value ());
194 }
195 
196 void
198 {
204 
205  INSTALL_NCUNOP (op_incr, octave_float_matrix, incr);
206  INSTALL_NCUNOP (op_decr, octave_float_matrix, decr);
208 
228  el_and_not);
230  el_or_not);
232  el_not_and);
234  el_not_or);
236  trans_mul);
238  mul_trans);
240  trans_mul);
242  mul_trans);
244  trans_ldiv);
246  trans_ldiv);
247 
251 
254  INSTALL_ASSIGNOP (op_asn_eq, octave_matrix,
255  octave_float_matrix, dbl_assign);
256 
258  null_assign);
260  null_assign);
262  null_assign);
263 
265  assign_add);
267  assign_sub);
269  assign_el_mul);
271  assign_el_div);
272 
273  INSTALL_CONVOP (octave_float_matrix, octave_matrix, float_matrix_to_matrix);
274 }
ComplexColumnVector quotient_eq(ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:149
ComplexColumnVector product(const ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:149
FloatMatrix float_matrix_value(bool=false) const
Definition: ov-ch-mat.h:111
octave_value op_mul_trans(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1357
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
virtual MatrixType matrix_type(void) const
Definition: ov-base.cc:347
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
#define DEFBINOPX(name, t1, t2)
Definition: ops.h:276
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1335
boolNDArray mx_el_and_not(const boolNDArray &m1, const boolNDArray &m2)
Definition: boolNDArray.cc:137
#define DEFNDASSIGNOP_FNOP(name, t1, t2, f, fnop)
Definition: ops.h:146
boolNDArray mx_el_not_or(const boolNDArray &m1, const boolNDArray &m2)
Definition: boolNDArray.cc:137
#define DEFNULLASSIGNOP_FN(name, t, f)
Definition: ops.h:116
boolNDArray mx_el_not_and(const boolNDArray &m1, const boolNDArray &m2)
Definition: boolNDArray.cc:137
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
boolMatrix mx_el_le(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
void error(const char *fmt,...)
Definition: error.cc:476
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 DEFNDBINOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:330
#define INSTALL_ASSIGNOP(op, t1, t2, f)
Definition: ops.h:55
MatrixType matrix_type(void) const
Definition: ov-base-mat.h:120
#define CAST_BINOP_ARGS(t1, t2)
Definition: ops.h:79
ComplexMatrix mul_trans(const ComplexMatrix &m, const SparseComplexMatrix &a)
Definition: CSparse.cc:7415
#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 quotient(const ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:149
#define INSTALL_BINOP(op, t1, t2, f)
Definition: ops.h:46
boolMatrix mx_el_ge(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
#define CAST_CONV_ARG(t)
Definition: ops.h:83
octave_value op_not(const octave_value &a)
Definition: ov.h:1293
boolMatrix mx_el_gt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
octave_value op_trans_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1356
octave_value op_transpose(const octave_value &a)
Definition: ov.h:1297
octave_value elem_xpow(double a, const SparseMatrix &b)
Definition: sparse-xpow.cc:258
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_mul_herm(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1359
Matrix xleftdiv(const SparseMatrix &a, const Matrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:466
boolMatrix mx_el_ne(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
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
virtual FloatMatrix float_matrix_value(bool=false) const
Definition: ov-base.cc:549
#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
ComplexMatrix xgemm(const ComplexMatrix &a, const ComplexMatrix &b, blas_trans_type transa, blas_trans_type transb)
Definition: CMatrix.cc:3686
const octave_char_matrix & v2
boolMatrix mx_el_or(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
void assign(octave_value::assign_op, const octave_value &)
Definition: oct-lvalue.cc:33
#define DEFNDASSIGNOP_FN(name, t1, t2, e, f)
Definition: ops.h:125
virtual FloatNDArray float_array_value(bool=false) const
Definition: ov-base.cc:566
Matrix xdiv(const Matrix &a, const SparseMatrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:135
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1338
ComplexMatrix trans_mul(const SparseComplexMatrix &m, const ComplexMatrix &a)
Definition: CSparse.cc:7445
boolMatrix mx_el_and(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
boolNDArray mx_el_or_not(const boolNDArray &m1, const boolNDArray &m2)
Definition: boolNDArray.cc:137
#define DEFNDASSIGNOP_OP(name, t1, t2, f, op)
Definition: ops.h:135
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
boolMatrix mx_el_lt(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
#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
octave_value op_herm_mul(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1358
#define DEFNDUNOP_OP(name, t, e, op)
Definition: ops.h:241
octave_value op_hermitian(const octave_value &a)
Definition: ov.h:1298
void install_fm_fm_ops(void)
Definition: op-fm-fm.cc:197
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1336
#define DEFNDBINOP_OP(name, t1, t2, e1, e2, op)
Definition: ops.h:313
#define DEFNDCATOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:364
boolMatrix mx_el_eq(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
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
ComplexColumnVector product_eq(ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:149
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())? '\'': '"'))
FloatNDArray float_array_value(bool=false) const
Definition: ov-ch-mat.h:117
#define DEFNCUNOP_METHOD(name, t, method)
Definition: ops.h:264