GNU Octave  3.8.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-scm-m.cc
Go to the documentation of this file.
1 /*
2 
3 Copyright (C) 2004-2013 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 #ifdef HAVE_CONFIG_H
25 #include <config.h>
26 #endif
27 
28 #include "gripes.h"
29 #include "oct-obj.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  const octave_matrix&);
54  MatrixType typ = v2.matrix_type ();
55 
57  v2.matrix_value (), typ);
58 
59  v2.matrix_type (typ);
60  return ret;
61 }
62 
63 DEFBINOPX (pow, sparse_complex_matrix, matrix)
64 {
65  error ("can't do A ^ B for A and B both matrices");
66  return octave_value ();
67 }
68 
69 DEFBINOP (ldiv, sparse_complex_matrix, matrix)
70 {
72 
73  if (v1.rows () == 1 && v1.columns () == 1)
74  {
76 
77  if (d == 0.0)
79 
80  return octave_value (v2.array_value () / d);
81  }
82  else
83  {
84  MatrixType typ = v1.matrix_type ();
85 
87  v2.matrix_value (), typ);
88 
89  v1.matrix_type (typ);
90  return ret;
91  }
92 }
93 
94 DEFBINOP_FN (lt, sparse_complex_matrix, matrix, mx_el_lt)
95 DEFBINOP_FN (le, sparse_complex_matrix, matrix, mx_el_le)
96 DEFBINOP_FN (eq, sparse_complex_matrix, matrix, mx_el_eq)
97 DEFBINOP_FN (ge, sparse_complex_matrix, matrix, mx_el_ge)
98 DEFBINOP_FN (gt, sparse_complex_matrix, matrix, mx_el_gt)
99 DEFBINOP_FN (ne, sparse_complex_matrix, matrix, mx_el_ne)
100 
101 DEFBINOP_FN (el_mul, sparse_complex_matrix, matrix, product)
102 DEFBINOP_FN (el_div, sparse_complex_matrix, matrix, quotient)
103 
104 DEFBINOP (el_pow, sparse_complex_matrix, matrix)
105 {
107  const octave_matrix&);
108 
109  return octave_value
111  (v2.matrix_value ())));
112 }
113 
114 DEFBINOP (el_ldiv, sparse_complex_matrix, matrix)
115 {
117  const octave_matrix&);
118 
119  return octave_value
121 }
122 
123 DEFBINOP_FN (el_and, sparse_complex_matrix, matrix, mx_el_and)
124 DEFBINOP_FN (el_or, sparse_complex_matrix, matrix, mx_el_or)
125 
126 DEFCATOP (scm_m, sparse_complex_matrix, matrix)
127 {
129  SparseMatrix tmp (v2.matrix_value ());
130  return octave_value
131  (v1.sparse_complex_matrix_value (). concat (tmp, ra_idx));
132 }
133 
134 DEFASSIGNOP (assign, sparse_complex_matrix, matrix)
135 {
137 
139  v1.assign (idx, tmp);
140  return octave_value ();
141 }
142 
143 void
145 {
159  el_mul);
161  el_div);
163  el_pow);
165  el_ldiv);
167  el_and);
169  el_or);
170 
172 
174  assign);
175 }