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-fcs-fm.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 "mx-fcs-fnda.h"
28 #include "mx-fnda-fcs.h"
29 
30 #include "errwarn.h"
31 #include "ovl.h"
32 #include "ov.h"
33 #include "ov-complex.h"
34 #include "ov-flt-complex.h"
35 #include "ov-cx-mat.h"
36 #include "ov-flt-cx-mat.h"
37 #include "ov-flt-re-mat.h"
38 #include "ov-re-mat.h"
39 #include "ov-typeinfo.h"
40 #include "ops.h"
41 #include "xdiv.h"
42 #include "xpow.h"
43 
44 // complex scalar by matrix ops.
45 
46 DEFNDBINOP_OP (add, float_complex, float_matrix, float_complex, float_array, +)
47 DEFNDBINOP_OP (sub, float_complex, float_matrix, float_complex, float_array, -)
48 DEFNDBINOP_OP (mul, float_complex, float_matrix, float_complex, float_array, *)
49 
50 DEFBINOP (div, float_complex, float_matrix)
51 {
52  const octave_float_complex& v1 = dynamic_cast<const octave_float_complex&> (a1);
53  const octave_float_matrix& v2 = dynamic_cast<const octave_float_matrix&> (a2);
54 
57  MatrixType typ = v2.matrix_type ();
58 
59  FloatComplexMatrix ret = xdiv (m1, m2, typ);
60 
61  v2.matrix_type (typ);
62  return ret;
63 }
64 
65 DEFBINOP_FN (pow, float_complex, float_matrix, xpow)
66 
67 DEFBINOP (ldiv, float_complex, float_matrix)
68 {
69  const octave_float_complex& v1 = dynamic_cast<const octave_float_complex&> (a1);
70  const octave_float_matrix& v2 = dynamic_cast<const octave_float_matrix&> (a2);
71 
73 
74  if (d == 0.0f)
76 
77  return octave_value (v2.float_array_value () / d);
78 }
79 
80 DEFNDCMPLXCMPOP_FN (lt, float_complex, float_matrix, float_complex,
81  float_array, mx_el_lt)
82 DEFNDCMPLXCMPOP_FN (le, float_complex, float_matrix, float_complex,
83  float_array, mx_el_le)
84 DEFNDCMPLXCMPOP_FN (eq, float_complex, float_matrix, float_complex,
85  float_array, mx_el_eq)
86 DEFNDCMPLXCMPOP_FN (ge, float_complex, float_matrix, float_complex,
87  float_array, mx_el_ge)
88 DEFNDCMPLXCMPOP_FN (gt, float_complex, float_matrix, float_complex,
89  float_array, mx_el_gt)
90 DEFNDCMPLXCMPOP_FN (ne, float_complex, float_matrix, float_complex,
91  float_array, mx_el_ne)
92 
93 DEFNDBINOP_OP (el_mul, float_complex, float_matrix, float_complex,
94  float_array, *)
95 DEFNDBINOP_FN (el_div, float_complex, float_matrix, float_complex,
96  float_array, x_el_div)
97 DEFNDBINOP_FN (el_pow, float_complex, float_matrix, float_complex,
98  float_array, elem_xpow)
99 
100 DEFBINOP (el_ldiv, float_complex, float_matrix)
101 {
102  const octave_float_complex& v1 = dynamic_cast<const octave_float_complex&> (a1);
103  const octave_float_matrix& v2 = dynamic_cast<const octave_float_matrix&> (a2);
104 
106 
107  if (d == 0.0f)
109 
110  return octave_value (v2.float_array_value () / d);
111 }
112 
113 DEFNDBINOP_FN (el_and, float_complex, float_matrix, float_complex,
114  float_array, mx_el_and)
115 DEFNDBINOP_FN (el_or, float_complex, float_matrix, float_complex,
116  float_array, mx_el_or)
117 
118 DEFNDCATOP_FN (fcs_fm, float_complex, float_matrix, float_complex_array,
119  float_array, concat)
120 
121 DEFNDCATOP_FN (cs_fm, complex, float_matrix, float_complex_array,
122  float_array, concat)
123 
124 DEFNDCATOP_FN (fcs_m, float_complex, matrix, float_complex_array,
125  float_array, concat)
126 
127 void
129 {
146  el_ldiv);
149 
153 
158 }
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
void install_fcs_fm_ops(void)
Definition: op-fcs-fm.cc:128
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 * f
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
#define DEFNDBINOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:294
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
FloatComplexMatrix float_complex_matrix_value(bool=false) const
octave_int< T > pow(const octave_int< T > &a, const octave_int< T > &b)
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
FloatComplex float_complex_value(bool=false) const
const octave_char_matrix & v2
boolMatrix mx_el_or(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
Matrix xdiv(const Matrix &a, const SparseMatrix &b, MatrixType &typ)
Definition: sparse-xdiv.cc:133
octave_value op_ne(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1527
boolMatrix mx_el_and(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
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
#define DEFNDCMPLXCMPOP_FN(name, t1, t2, e1, e2, f)
Definition: ops.h:305
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
std::complex< float > FloatComplex
Definition: oct-cmplx.h:32
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
#define INSTALL_ASSIGNCONV(t1, t2, tr)
Definition: ops.h:66
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())? '\'': '"'))