Mercurial > octave
annotate libinterp/operators/op-sm-m.cc @ 23083:e9a0469dedd9 stable
maint: strip extra trailing newlines from files.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 20 Jan 2017 12:19:08 -0500 |
parents | 34ce5be04942 |
children | ef4d915df748 3ac9f9ecfae5 |
rev | line source |
---|---|
5164 | 1 /* |
2 | |
22323
bac0d6f07a3e
maint: Update copyright notices for 2016.
John W. Eaton <jwe@octave.org>
parents:
21724
diff
changeset
|
3 Copyright (C) 2004-2016 David Bateman |
11523 | 4 Copyright (C) 1998-2004 Andy Adler |
7016 | 5 |
6 This file is part of Octave. | |
5164 | 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 | |
7016 | 10 Free Software Foundation; either version 3 of the License, or (at your |
11 option) any later version. | |
5164 | 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 | |
7016 | 19 along with Octave; see the file COPYING. If not, see |
20 <http://www.gnu.org/licenses/>. | |
5164 | 21 |
22 */ | |
23 | |
21724
aba2e6293dd8
use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents:
21647
diff
changeset
|
24 #if defined (HAVE_CONFIG_H) |
21301
40de9f8f23a6
Use '#include "config.h"' rather than <config.h>.
Rik <rik@octave.org>
parents:
21200
diff
changeset
|
25 # include "config.h" |
5164 | 26 #endif |
27 | |
21100
e39e05d90788
Switch gripe_XXX to either err_XXX or warn_XXX naming scheme.
Rik <rik@octave.org>
parents:
20940
diff
changeset
|
28 #include "errwarn.h" |
20940
48b2ad5ee801
maint: Rename oct-obj.[cc|h] to ovl.[cc|h] for clarity.
Rik <rik@octave.org>
parents:
20893
diff
changeset
|
29 #include "ovl.h" |
5164 | 30 #include "ov.h" |
31 #include "ov-typeinfo.h" | |
32 #include "ov-re-mat.h" | |
33 #include "ops.h" | |
34 #include "xdiv.h" | |
35 | |
36 #include "sparse-xpow.h" | |
37 #include "sparse-xdiv.h" | |
38 #include "smx-sm-m.h" | |
39 #include "smx-m-sm.h" | |
40 #include "ov-re-sparse.h" | |
41 | |
42 // sparse matrix by matrix ops. | |
43 | |
44 DEFBINOP_OP (add, sparse_matrix, matrix, +) | |
45 DEFBINOP_OP (sub, sparse_matrix, matrix, -) | |
46 | |
5429 | 47 DEFBINOP_OP (mul, sparse_matrix, matrix, *) |
5164 | 48 |
49 DEFBINOP (div, sparse_matrix, matrix) | |
50 { | |
21647
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
51 const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1); |
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
52 const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2); |
5785 | 53 MatrixType typ = v2.matrix_type (); |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
54 |
5785 | 55 Matrix ret = xdiv (v1.matrix_value (), v2.matrix_value (), typ); |
56 | |
57 v2.matrix_type (typ); | |
58 return ret; | |
5164 | 59 } |
60 | |
61 DEFBINOPX (pow, sparse_matrix, matrix) | |
62 { | |
63 error ("can't do A ^ B for A and B both matrices"); | |
64 } | |
65 | |
66 DEFBINOP (ldiv, sparse_matrix, matrix) | |
67 { | |
21647
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
68 const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1); |
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
69 const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2); |
6221 | 70 |
14846
460a3c6d8bf1
maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists.
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
71 if (v1.rows () == 1 && v1.columns () == 1) |
6221 | 72 { |
73 double d = v1.scalar_value (); | |
74 | |
75 if (d == 0.0) | |
21100
e39e05d90788
Switch gripe_XXX to either err_XXX or warn_XXX naming scheme.
Rik <rik@octave.org>
parents:
20940
diff
changeset
|
76 warn_divide_by_zero (); |
5322 | 77 |
6221 | 78 return octave_value (v2.array_value () / d); |
79 } | |
80 else | |
81 { | |
82 MatrixType typ = v1.matrix_type (); | |
5322 | 83 |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
84 Matrix ret = xleftdiv (v1.sparse_matrix_value (), |
9931
fb6b6fcafa62
untabify files in src/OPERATORS directory
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
85 v2.matrix_value (), typ); |
6221 | 86 |
87 v1.matrix_type (typ); | |
88 return ret; | |
89 } | |
5164 | 90 } |
91 | |
7802
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
92 DEFBINOP_FN (trans_mul, sparse_matrix, matrix, trans_mul); |
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
93 |
5164 | 94 DEFBINOP_FN (lt, sparse_matrix, matrix, mx_el_lt) |
95 DEFBINOP_FN (le, sparse_matrix, matrix, mx_el_le) | |
96 DEFBINOP_FN (eq, sparse_matrix, matrix, mx_el_eq) | |
97 DEFBINOP_FN (ge, sparse_matrix, matrix, mx_el_ge) | |
98 DEFBINOP_FN (gt, sparse_matrix, matrix, mx_el_gt) | |
99 DEFBINOP_FN (ne, sparse_matrix, matrix, mx_el_ne) | |
100 | |
101 DEFBINOP_FN (el_mul, sparse_matrix, matrix, product) | |
102 DEFBINOP_FN (el_div, sparse_matrix, matrix, quotient) | |
103 | |
104 DEFBINOP (el_pow, sparse_matrix, matrix) | |
105 { | |
21647
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
106 const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1); |
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
107 const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2); |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
108 |
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
109 return octave_value (elem_xpow (v1.sparse_matrix_value (), |
9931
fb6b6fcafa62
untabify files in src/OPERATORS directory
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
110 SparseMatrix (v2.matrix_value ()))); |
5164 | 111 } |
112 | |
113 DEFBINOP (el_ldiv, sparse_matrix, matrix) | |
114 { | |
21647
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
115 const octave_sparse_matrix& v1 = dynamic_cast<const octave_sparse_matrix&> (a1); |
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
116 const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2); |
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
117 |
5164 | 118 return octave_value |
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
119 (quotient (v2.matrix_value (), v1.sparse_matrix_value ())); |
5164 | 120 } |
121 | |
122 DEFBINOP_FN (el_and, sparse_matrix, matrix, mx_el_and) | |
123 DEFBINOP_FN (el_or, sparse_matrix, matrix, mx_el_or) | |
124 | |
125 DEFCATOP (sm_m, sparse_matrix, matrix) | |
126 { | |
21647
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
127 octave_sparse_matrix& v1 = dynamic_cast<octave_sparse_matrix&> (a1); |
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
128 const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2); |
5164 | 129 SparseMatrix tmp (v2.matrix_value ()); |
130 return octave_value (v1.sparse_matrix_value (). concat (tmp, ra_idx)); | |
131 } | |
132 | |
133 DEFASSIGNOP (assign, sparse_matrix, matrix) | |
134 { | |
21647
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
135 octave_sparse_matrix& v1 = dynamic_cast<octave_sparse_matrix&> (a1); |
66cae7a6dc47
eliminate some macros for operator definitions
John W. Eaton <jwe@octave.org>
parents:
21301
diff
changeset
|
136 const octave_matrix& v2 = dynamic_cast<const octave_matrix&> (a2); |
5164 | 137 |
138 SparseMatrix tmp (v2.matrix_value ()); | |
139 v1.assign (idx, tmp); | |
140 return octave_value (); | |
141 } | |
142 | |
143 void | |
144 install_sm_m_ops (void) | |
145 { | |
146 INSTALL_BINOP (op_add, octave_sparse_matrix, octave_matrix, add); | |
147 INSTALL_BINOP (op_sub, octave_sparse_matrix, octave_matrix, sub); | |
148 INSTALL_BINOP (op_mul, octave_sparse_matrix, octave_matrix, mul); | |
149 INSTALL_BINOP (op_div, octave_sparse_matrix, octave_matrix, div); | |
150 INSTALL_BINOP (op_pow, octave_sparse_matrix, octave_matrix, pow); | |
151 INSTALL_BINOP (op_ldiv, octave_sparse_matrix, octave_matrix, ldiv); | |
7802
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
152 INSTALL_BINOP (op_trans_mul, octave_sparse_matrix, octave_matrix, trans_mul); |
1a446f28ce68
implement optimized sparse-dense transposed multiplication
Jaroslav Hajek <highegg@gmail.com>
parents:
7017
diff
changeset
|
153 INSTALL_BINOP (op_herm_mul, octave_sparse_matrix, octave_matrix, trans_mul); |
5164 | 154 INSTALL_BINOP (op_lt, octave_sparse_matrix, octave_matrix, lt); |
155 INSTALL_BINOP (op_le, octave_sparse_matrix, octave_matrix, le); | |
156 INSTALL_BINOP (op_eq, octave_sparse_matrix, octave_matrix, eq); | |
157 INSTALL_BINOP (op_ge, octave_sparse_matrix, octave_matrix, ge); | |
158 INSTALL_BINOP (op_gt, octave_sparse_matrix, octave_matrix, gt); | |
159 INSTALL_BINOP (op_ne, octave_sparse_matrix, octave_matrix, ne); | |
160 INSTALL_BINOP (op_el_mul, octave_sparse_matrix, octave_matrix, el_mul); | |
161 INSTALL_BINOP (op_el_div, octave_sparse_matrix, octave_matrix, el_div); | |
162 INSTALL_BINOP (op_el_pow, octave_sparse_matrix, octave_matrix, el_pow); | |
163 INSTALL_BINOP (op_el_ldiv, octave_sparse_matrix, octave_matrix, el_ldiv); | |
164 INSTALL_BINOP (op_el_and, octave_sparse_matrix, octave_matrix, el_and); | |
165 INSTALL_BINOP (op_el_or, octave_sparse_matrix, octave_matrix, el_or); | |
166 | |
167 INSTALL_CATOP (octave_sparse_matrix, octave_matrix, sm_m); | |
168 | |
169 INSTALL_ASSIGNOP (op_asn_eq, octave_sparse_matrix, octave_matrix, assign); | |
170 } |