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-m.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-re-mat.h"
33 #include "ov-cx-mat.h"
34 #include "ops.h"
35 #include "xdiv.h"
36 
37 #include "sparse-xpow.h"
38 #include "sparse-xdiv.h"
39 #include "smx-scm-m.h"
40 #include "smx-m-scm.h"
41 #include "ov-cx-sparse.h"
42 
43 // sparse complex matrix by matrix ops.
44 
45 DEFBINOP_OP (add, sparse_complex_matrix, matrix, +)
46 DEFBINOP_OP (sub, sparse_complex_matrix, matrix, -)
47 
48 DEFBINOP_OP (mul, sparse_complex_matrix, matrix, *)
49 
50 DEFBINOP (div, sparse_complex_matrix, matrix)
51 {
53  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
54  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
55  MatrixType typ = v2.matrix_type ();
56 
58  v2.matrix_value (), typ);
59 
60  v2.matrix_type (typ);
61  return ret;
62 }
63 
64 DEFBINOPX (pow, sparse_complex_matrix, matrix)
65 {
66  error ("can't do A ^ B for A and B both matrices");
67 }
68 
69 DEFBINOP (ldiv, sparse_complex_matrix, matrix)
70 {
72  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
73  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
74 
75  if (v1.rows () == 1 && v1.columns () == 1)
76  {
77  Complex d = v1.complex_value ();
78 
79  if (d == 0.0)
81 
82  return octave_value (v2.array_value () / d);
83  }
84  else
85  {
86  MatrixType typ = v1.matrix_type ();
87 
89  v2.matrix_value (), typ);
90 
91  v1.matrix_type (typ);
92  return ret;
93  }
94 }
95 
96 DEFBINOP_FN (lt, sparse_complex_matrix, matrix, mx_el_lt)
97 DEFBINOP_FN (le, sparse_complex_matrix, matrix, mx_el_le)
98 DEFBINOP_FN (eq, sparse_complex_matrix, matrix, mx_el_eq)
99 DEFBINOP_FN (ge, sparse_complex_matrix, matrix, mx_el_ge)
100 DEFBINOP_FN (gt, sparse_complex_matrix, matrix, mx_el_gt)
101 DEFBINOP_FN (ne, sparse_complex_matrix, matrix, mx_el_ne)
102 
103 DEFBINOP_FN (el_mul, sparse_complex_matrix, matrix, product)
104 DEFBINOP_FN (el_div, sparse_complex_matrix, matrix, quotient)
105 
106 DEFBINOP (el_pow, sparse_complex_matrix, matrix)
107 {
109  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
110  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
111 
112  return octave_value
114  (v2.matrix_value ())));
115 }
116 
117 DEFBINOP (el_ldiv, sparse_complex_matrix, matrix)
118 {
120  = dynamic_cast<const octave_sparse_complex_matrix&> (a1);
121  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
122 
123  return octave_value
125 }
126 
127 DEFBINOP_FN (el_and, sparse_complex_matrix, matrix, mx_el_and)
128 DEFBINOP_FN (el_or, sparse_complex_matrix, matrix, mx_el_or)
129 
130 DEFCATOP (scm_m, sparse_complex_matrix, matrix)
131 {
133  = dynamic_cast<octave_sparse_complex_matrix&> (a1);
134  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
135  SparseMatrix tmp (v2.matrix_value ());
136  return octave_value
138 }
139 
140 DEFASSIGNOP (assign, sparse_complex_matrix, matrix)
141 {
143  = dynamic_cast<octave_sparse_complex_matrix&> (a1);
144  const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2);
145 
147  v1.assign (idx, tmp);
148  return octave_value ();
149 }
150 
151 void
153 {
167  el_mul);
169  el_div);
171  el_pow);
173  el_ldiv);
175  el_and);
177  el_or);
178 
180 
182  assign);
183 }
ComplexColumnVector product(const ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:151
#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
NDArray array_value(bool=false) const
Definition: ov-re-mat.h:172
#define DEFBINOPX(name, t1, t2)
Definition: ops.h:218
const octave_base_value const Array< octave_idx_type > & ra_idx
octave_value op_eq(const octave_value &a1, const octave_value &a2)
Definition: ov.h:1524
octave_idx_type columns(void) const
Definition: ov-base.h:319
Complex complex_value(bool=false) const
void install_scm_m_ops(void)
Definition: op-scm-m.cc:152
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
boolMatrix mx_el_le(const boolMatrix &m1, const boolMatrix &m2)
Definition: boolMatrix.cc:90
void error(const char *fmt,...)
Definition: error.cc:570
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
Matrix matrix_value(bool=false) const
Definition: ov-re-mat.cc:141
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
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
#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
ComplexColumnVector quotient(const ComplexColumnVector &x, const ComplexColumnVector &y)
Definition: CColVector.h:151
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
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
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
ComplexMatrix complex_matrix_value(bool=false) const
Definition: ov-re-mat.cc:185
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
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
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
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())? '\'': '"'))