Mercurial > octave-nkf

annotate liboctave/MArrayN.cc @ 9607:1be3c73ed7b5

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
reuse temporary arrays in nested expressions
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 03 Sep 2009 08:48:51 +0200
parents 3a1dd361f978
children 66970dd627f6
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
1 /*
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
2
9245
16f53d29049f update copyright notices
John W. Eaton <jwe@octave.org>
parents: 9200
diff changeset
3 Copyright (C) 1996, 1997, 2003, 2004, 2005, 2007, 2009 John W. Eaton
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
4
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
5 This file is part of Octave.
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
6
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
7 Octave is free software; you can redistribute it and/or modify it
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
8 under the terms of the GNU General Public License as published by the
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
10 option) any later version.
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
11
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
12 Octave is distributed in the hope that it will be useful, but WITHOUT
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
15 for more details.
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
16
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
17 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
19 <http://www.gnu.org/licenses/>.
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
20
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
21 */
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
22
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
23 #ifdef HAVE_CONFIG_H
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
24 #include <config.h>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
25 #endif
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
26
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
27 #include "MArrayN.h"
4669
334a27c8f453 [project @ 2003-11-26 07:02:42 by jwe]
jwe
parents: 4646
diff changeset
28 #include "Array-util.h"
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
29 #include "lo-error.h"
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
30
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
31 #include "MArray-defs.h"
9607
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
32 #include "mx-inlines.cc"
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
33
4513
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe]
jwe
parents: 4508
diff changeset
34 // N-dimensional array with math ops.
9607
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
35 template <class T>
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
36 void
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
37 MArrayN<T>::changesign (void)
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
38 {
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
39 if (Array<T>::is_shared ())
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
40 *this = - *this;
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
41 else
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
42 do_mx_inplace_op<MArrayN<T> > (*this, mx_inline_uminus2);
1be3c73ed7b5 reuse temporary arrays in nested expressions
Jaroslav Hajek <highegg@gmail.com>
parents: 9557
diff changeset
43 }
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
44
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
45 // Element by element MArrayN by scalar ops.
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
46
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
47 template <class T>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
48 MArrayN<T>&
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
49 operator += (MArrayN<T>& a, const T& s)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
50 {
9546
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
51 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
52 return a = a + s;
4646
eff8f977508c [project @ 2003-11-23 18:13:19 by jwe]
jwe
parents: 4645
diff changeset
53 DO_VS_OP2 (T, a, +=, s)
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
54 return a;
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
55 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
56
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
57 template <class T>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
58 MArrayN<T>&
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
59 operator -= (MArrayN<T>& a, const T& s)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
60 {
9546
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
61 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
62 return a = a - s;
4646
eff8f977508c [project @ 2003-11-23 18:13:19 by jwe]
jwe
parents: 4645
diff changeset
63 DO_VS_OP2 (T, a, -=, s)
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
64 return a;
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
65 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
66
9546
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
67 template <class T>
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
68 MArrayN<T>&
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
69 operator *= (MArrayN<T>& a, const T& s)
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
70 {
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
71 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
72 return a = a * s;
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
73 DO_VS_OP2 (T, a, *=, s)
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
74 return a;
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
75 }
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
76
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
77 template <class T>
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
78 MArrayN<T>&
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
79 operator /= (MArrayN<T>& a, const T& s)
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
80 {
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
81 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
82 return a = a / s;
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
83 DO_VS_OP2 (T, a, /=, s)
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
84 return a;
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
85 }
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
86
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
87 // Element by element MArrayN by MArrayN ops.
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
88
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
89 template <class T>
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
90 MArrayN<T>&
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
91 operator += (MArrayN<T>& a, const MArrayN<T>& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
92 {
9546
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
93 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
94 return a = a + b;
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
95
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
96 dim_vector a_dims = a.dims ();
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
97 dim_vector b_dims = b.dims ();
4645
bd2067547b40 [project @ 2003-11-23 08:07:52 by jwe]
jwe
parents: 4543
diff changeset
98
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
99 if (a_dims != b_dims)
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
100 gripe_nonconformant ("operator +=", a_dims, b_dims);
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
101 else
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
102 {
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
103 octave_idx_type l = a.length ();
4645
bd2067547b40 [project @ 2003-11-23 08:07:52 by jwe]
jwe
parents: 4543
diff changeset
104
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
105 if (l > 0)
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
106 DO_VV_OP2 (T, a, +=, b);
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
107 }
4645
bd2067547b40 [project @ 2003-11-23 08:07:52 by jwe]
jwe
parents: 4543
diff changeset
108
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
109 return a;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
110 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
111
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
112 template <class T>
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
113 MArrayN<T>&
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
114 operator -= (MArrayN<T>& a, const MArrayN<T>& b)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
115 {
9546
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
116 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases
Jaroslav Hajek <highegg@gmail.com>
parents: 9245
diff changeset
117 return a = a - b;
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
118
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
119 dim_vector a_dims = a.dims ();
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
120 dim_vector b_dims = b.dims ();
4645
bd2067547b40 [project @ 2003-11-23 08:07:52 by jwe]
jwe
parents: 4543
diff changeset
121
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
122 if (a_dims != b_dims)
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
123 gripe_nonconformant ("operator -=", a_dims, b_dims);
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
124 else
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
125 {
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
126 octave_idx_type l = a.length ();
4645
bd2067547b40 [project @ 2003-11-23 08:07:52 by jwe]
jwe
parents: 4543
diff changeset
127
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
128 if (l > 0)
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
129 DO_VV_OP2 (T, a, -=, b);
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
130 }
9548
e5f7aee2ab8c optimize &=, |= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9546
diff changeset
131
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
132 return a;
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
133 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
134
9557
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
135 template <class T>
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
136 MArrayN<T>&
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
137 product_eq (MArrayN<T>& a, const MArrayN<T>& b)
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
138 {
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
139 if (a.is_shared ())
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
140 return a = product (a, b);
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
141
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
142 dim_vector a_dims = a.dims ();
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
143 dim_vector b_dims = b.dims ();
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
144
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
145 if (a_dims != b_dims)
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
146 gripe_nonconformant ("operator .*=", a_dims, b_dims);
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
147 else
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
148 {
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
149 octave_idx_type l = a.length ();
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
150
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
151 if (l > 0)
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
152 DO_VV_OP2 (T, a, *=, b);
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
153 }
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
154
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
155 return a;
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
156 }
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
157
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
158 template <class T>
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
159 MArrayN<T>&
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
160 quotient_eq (MArrayN<T>& a, const MArrayN<T>& b)
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
161 {
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
162 if (a.is_shared ())
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
163 return a = quotient (a, b);
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
164
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
165 dim_vector a_dims = a.dims ();
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
166 dim_vector b_dims = b.dims ();
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
167
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
168 if (a_dims != b_dims)
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
169 gripe_nonconformant ("operator ./=", a_dims, b_dims);
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
170 else
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
171 {
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
172 octave_idx_type l = a.length ();
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
173
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
174 if (l > 0)
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
175 DO_VV_OP2 (T, a, /=, b);
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
176 }
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
177
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
178 return a;
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
179 }
3a1dd361f978 optimize .*=, ./= operators
Jaroslav Hajek <highegg@gmail.com>
parents: 9548
diff changeset
180
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
181 // Element by element MArrayN by scalar ops.
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
182
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
183 #define MARRAYN_NDS_OP(OP) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
184 template <class T> \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
185 MArrayN<T> \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
186 operator OP (const MArrayN<T>& a, const T& s) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
187 { \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
188 MArrayN<T> result (a.dims ()); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
189 T *r = result.fortran_vec (); \
5275
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe]
jwe
parents: 5105
diff changeset
190 octave_idx_type l = a.length (); \
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
191 const T *v = a.data (); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
192 DO_VS_OP (r, l, v, OP, s); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
193 return result; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
194 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
195
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
196 MARRAYN_NDS_OP (+)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
197 MARRAYN_NDS_OP (-)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
198 MARRAYN_NDS_OP (*)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
199 MARRAYN_NDS_OP (/)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
200
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
201 // Element by element MArrayN by scalar ops.
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
202
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
203 #define MARRAYN_SND_OP(OP) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
204 template <class T> \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
205 MArrayN<T> \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
206 operator OP (const T& s, const MArrayN<T>& a) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
207 { \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
208 MArrayN<T> result (a.dims ()); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
209 T *r = result.fortran_vec (); \
5275
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe]
jwe
parents: 5105
diff changeset
210 octave_idx_type l = a.length (); \
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
211 const T *v = a.data (); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
212 DO_SV_OP (r, l, s, OP, v); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
213 return result; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
214 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
215
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
216 MARRAYN_SND_OP (+)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
217 MARRAYN_SND_OP (-)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
218 MARRAYN_SND_OP (*)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
219 MARRAYN_SND_OP (/)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
220
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
221 #define MARRAY_NDND_OP(FCN, OP) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
222 template <class T> \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
223 MArrayN<T> \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
224 FCN (const MArrayN<T>& a, const MArrayN<T>& b) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
225 { \
4513
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe]
jwe
parents: 4508
diff changeset
226 dim_vector a_dims = a.dims (); \
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe]
jwe
parents: 4508
diff changeset
227 dim_vector b_dims = b.dims (); \
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
228 int dims_ok = 1; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
229 int any_dims_zero = 0; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
230 if (a_dims.length () != b_dims.length ()) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
231 dims_ok = 0; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
232 else \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
233 { \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
234 for (int i = 0; i < a_dims.length (); i++) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
235 { \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
236 if (a_dims (i) != b_dims (i)) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
237 { dims_ok = 0; break; } \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
238 if (a_dims (i) == 0) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
239 any_dims_zero = 1; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
240 } \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
241 } \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
242 if (!dims_ok) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
243 { \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
244 gripe_nonconformant (#FCN, a_dims, b_dims); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
245 return MArrayN<T> (); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
246 } \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
247 if (any_dims_zero) \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
248 return MArrayN<T> (a_dims); \
5275
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe]
jwe
parents: 5105
diff changeset
249 octave_idx_type l = a.length (); \
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
250 MArrayN<T> result (a_dims); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
251 T* r = result.fortran_vec (); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
252 const T *x = a.data (); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
253 const T *y = b.data (); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
254 DO_VV_OP (r, l, x, OP, y); \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
255 return result; \
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
256 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
257
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
258 MARRAY_NDND_OP (operator +, +)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
259 MARRAY_NDND_OP (operator -, -)
4543
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
260 MARRAY_NDND_OP (product, *)
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe]
jwe
parents: 4513
diff changeset
261 MARRAY_NDND_OP (quotient, /)
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
262
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
263 template <class T>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
264 MArrayN<T>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
265 operator + (const MArrayN<T>& a)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
266 {
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
267 return a;
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
268 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
269
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
270 template <class T>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
271 MArrayN<T>
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
272 operator - (const MArrayN<T>& a)
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
273 {
5275
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe]
jwe
parents: 5105
diff changeset
274 octave_idx_type l = a.length ();
4508
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
275 MArrayN<T> result (a.dims ());
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
276 T *r = result.fortran_vec ();
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
277 const T *x = a.data ();
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
278 NEG_V (r, l, x);
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
279 return result;
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
280 }
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
281
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
282 /*
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
283 ;;; Local Variables: ***
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
284 ;;; mode: C++ ***
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
285 ;;; End: ***
2bc437206787 [project @ 2003-09-11 01:06:43 by jwe]
jwe
parents:
diff changeset
286 */