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-scm-s.cc
Go to the documentation of this file.
1 /*
2 
3 Copyright (C) 2004-2017 David Bateman
4 Copyright (C) 1998-2004 Andy Adler
5 
6 This file is part of Octave.
7 
8 Octave is free software; you can redistribute it and/or modify it
9 under the terms of the GNU General Public License as published by the
10 Free Software Foundation; either version 3 of the License, or (at your
11 option) any later version.
12 
13 Octave is distributed in the hope that it will be useful, but WITHOUT
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
16 for more details.
17 
18 You should have received a copy of the GNU General Public License
19 along with Octave; see the file COPYING. If not, see
20 <http://www.gnu.org/licenses/>.
21 
22 */
23 
24 #if defined (HAVE_CONFIG_H)
25 # include "config.h"
26 #endif
27 
28 #include "errwarn.h"
29 #include "ovl.h"
30 #include "ov.h"
31 #include "ov-typeinfo.h"
32 #include "ov-cx-mat.h"
33 #include "ov-scalar.h"
34 #include "ops.h"
35 #include "xpow.h"
36 
37 #include "sparse-xpow.h"
38 #include "sparse-xdiv.h"
39 #include "smx-scm-s.h"
40 #include "smx-s-scm.h"
41 #include "ov-re-sparse.h"
42 #include "ov-cx-sparse.h"
43 
44 // sparse complex matrix by scalar ops.
45 
46 DEFBINOP_OP (add, sparse_complex_matrix, scalar, +)
47 DEFBINOP_OP (sub, sparse_complex_matrix, scalar, -)
48 DEFBINOP_OP (mul, sparse_complex_matrix, scalar, *)
49 
50 DEFBINOP (div, sparse_complex_matrix, scalar)
51 {
53  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
54  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
55 
56  double d = v2.double_value ();
58 
59  if (d == 0.0)
61 
62  retval = octave_value (v1.sparse_complex_matrix_value () / d);
63 
64  return retval;
65 }
66 
67 DEFBINOP (pow, sparse_complex_matrix, scalar)
68 {
70  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
71  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
72 
73  double tmp = v2.scalar_value ();
74  if (static_cast<int> (tmp) == tmp)
75  return xpow (v1.sparse_complex_matrix_value (), tmp);
76  else
77  return xpow (v1.complex_matrix_value (), tmp);
78 }
79 
80 DEFBINOP (ldiv, sparse_complex_matrix, scalar)
81 {
83  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
84  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
85 
86  if (v1.rows () == 1 && v1.columns () == 1)
87  {
88  Complex d = v1.complex_value ();
89 
90  if (d == 0.0)
92 
93  return octave_value (SparseComplexMatrix (1, 1, v2.scalar_value () / d));
94  }
95  else
96  {
97  MatrixType typ = v1.matrix_type ();
99  Matrix m2 = Matrix (1, 1, v2.scalar_value ());
100  ComplexMatrix ret = xleftdiv (m1, m2, typ);
101  v1.matrix_type (typ);
102  return ret;
103  }
104 }
105 
106 DEFBINOP_FN (lt, sparse_complex_matrix, scalar, mx_el_lt)
107 DEFBINOP_FN (le, sparse_complex_matrix, scalar, mx_el_le)
108 DEFBINOP_FN (eq, sparse_complex_matrix, scalar, mx_el_eq)
109 DEFBINOP_FN (ge, sparse_complex_matrix, scalar, mx_el_ge)
110 DEFBINOP_FN (gt, sparse_complex_matrix, scalar, mx_el_gt)
111 DEFBINOP_FN (ne, sparse_complex_matrix, scalar, mx_el_ne)
112 
113 DEFBINOP_OP (el_mul, sparse_complex_matrix, scalar, *)
114 
115 DEFBINOP (el_div, sparse_complex_matrix, scalar)
116 {
118  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
119  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
120 
121  double d = v2.double_value ();
123 
124  if (d == 0.0)
126 
127  retval = octave_value (v1.sparse_complex_matrix_value () / d);
128 
129  return retval;
130 }
131 
132 DEFBINOP_FN (el_pow, sparse_complex_matrix, scalar, elem_xpow)
133 
134 DEFBINOP (el_ldiv, sparse_complex_matrix, scalar)
135 {
137  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
138  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
139 
140  return octave_value
142 }
143 
144 DEFBINOP_FN (el_and, sparse_complex_matrix, scalar, mx_el_and)
145 DEFBINOP_FN (el_or, sparse_complex_matrix, scalar, mx_el_or)
146 
147 DEFCATOP (scm_s, sparse_complex_matrix, scalar)
148 {
150  = dynamic_cast<octave_sparse_complex_matrix&> (a1);
151  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
152  SparseComplexMatrix tmp (1, 1, v2.complex_value ());
153  return octave_value
155 }
156 
157 DEFASSIGNOP (assign, sparse_complex_matrix, scalar)
158 {
160  = dynamic_cast<octave_sparse_complex_matrix&> (a1);
161  const octave_scalar& v2 = dynamic_cast<const octave_scalar&> (a2);
162 
163  SparseComplexMatrix tmp (1, 1, v2.complex_value ());
164  v1.assign (idx, tmp);
165  return octave_value ();
166 }
167 
168 void
170 {
184  el_mul);
186  el_div);
188  el_pow);
190  el_ldiv);
192  el_and);
194  el_or);
195 
197 
199  assign);
200 }
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
const octave_base_value const Array< octave_idx_type > & ra_idx
Complex complex_value(bool=false) const
Definition: ov-scalar.h:176
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_idx_type columns(void) const
Definition: ov-base.h:319
Complex complex_value(bool=false) const
double scalar_value(bool=false) const
Definition: ov-scalar.h:150
octave_value op_el_ldiv(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1532
boolMatrix mx_el_le(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
SparseComplexMatrix concat(const SparseComplexMatrix &rb, const Array< octave_idx_type > &ra_idx)
Definition: CSparse.cc:581
octave_value op_pow(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1519
SparseComplexMatrix sparse_complex_matrix_value(bool=false) const
Definition: ov-cx-sparse.h:125
#define INSTALL_ASSIGNOP(op, t1, t2, f)
Definition: ops.h:57
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
#define DEFBINOP_OP(name, t1, t2, op)
Definition: ops.h:228
const octave_base_value & a2
void assign(const octave_value_list &idx, const T &rhs)
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
ComplexMatrix complex_matrix_value(bool=false) const
#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
#define DEFASSIGNOP(name, t1, t2)
Definition: ops.h:74
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)
double tmp
Definition: data.cc:6300
octave_value retval
Definition: data.cc:6294
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
void warn_divide_by_zero(void)
Definition: errwarn.cc:320
Definition: dMatrix.h:37
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
#define DEFCATOP(name, t1, t2)
Definition: ops.h:321
double double_value(bool=false) const
Definition: ov-scalar.h:145
const octave_char_matrix & v2
boolMatrix mx_el_or(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:87
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
MatrixType matrix_type(void) const
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
void install_scm_s_ops(void)
Definition: op-scm-s.cc:169
std::complex< double > Complex
Definition: oct-cmplx.h:31
octave_value op_ge(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1525
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_idx_type rows(void) const
Definition: ov-base.h:312
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())? '\'': '"'))