Mercurial > hg > octave-lyh
annotate liboctave/MSparse.h @ 10485:b4e14e628fc9
Truncate argv() for scripts used without command line parameters. Bug #29423
author | Judd Storrs <jstorrs@gmail.com> |
---|---|
date | Fri, 02 Apr 2010 13:28:02 -0400 |
parents | 19e1e4470e01 |
children | f0266ee4aabe |
rev | line source |
---|---|
5164 | 1 /* |
2 | |
8920 | 3 Copyright (C) 2004, 2005, 2007, 2008 David Bateman |
7016 | 4 Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004 Andy Adler |
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 | |
24 #if !defined (octave_MSparse_h) | |
25 #define octave_MSparse_h 1 | |
26 | |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10312
diff
changeset
|
27 #include "MArray.h" |
5164 | 28 |
29 #include "Sparse.h" | |
30 | |
31 // Two dimensional sparse array with math ops. | |
32 | |
33 // But first, some preprocessor abuse... | |
34 | |
35 #include "MSparse-defs.h" | |
36 | |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10312
diff
changeset
|
37 SPARSE_OPS_FORWARD_DECLS (MSparse, MArray, ) |
5164 | 38 |
39 template <class T> | |
40 class | |
41 MSparse : public Sparse<T> | |
42 { | |
43 public: | |
44 | |
45 MSparse (void) : Sparse<T> () { } | |
46 | |
5275 | 47 MSparse (octave_idx_type n, octave_idx_type m) : Sparse<T> (n, m) { } |
5164 | 48 |
6823 | 49 MSparse (const dim_vector& dv, octave_idx_type nz = 0) : |
50 Sparse<T> (dv, nz) { } | |
51 | |
5164 | 52 MSparse (const MSparse<T>& a) : Sparse<T> (a) { } |
53 | |
54 MSparse (const MSparse<T>& a, const dim_vector& dv) : Sparse<T> (a, dv) { } | |
55 | |
56 MSparse (const Sparse<T>& a) : Sparse<T> (a) { } | |
57 | |
10479
ded9beac7582
optimize sparse matrix assembly
Jaroslav Hajek <highegg@gmail.com>
parents:
10421
diff
changeset
|
58 MSparse (const Array<T>& a, const idx_vector& r, const idx_vector& c, |
ded9beac7582
optimize sparse matrix assembly
Jaroslav Hajek <highegg@gmail.com>
parents:
10421
diff
changeset
|
59 octave_idx_type nr = -1, octave_idx_type nc = -1, bool sum_terms = true) |
ded9beac7582
optimize sparse matrix assembly
Jaroslav Hajek <highegg@gmail.com>
parents:
10421
diff
changeset
|
60 : Sparse<T> (a, r, c, nr, nc, sum_terms) { } |
ded9beac7582
optimize sparse matrix assembly
Jaroslav Hajek <highegg@gmail.com>
parents:
10421
diff
changeset
|
61 |
5275 | 62 explicit MSparse (octave_idx_type r, octave_idx_type c, T val) : Sparse<T> (r, c, val) { } |
5164 | 63 |
5275 | 64 MSparse (octave_idx_type r, octave_idx_type c, octave_idx_type num_nz) : Sparse<T> (r, c, num_nz) { } |
5164 | 65 |
66 ~MSparse (void) { } | |
67 | |
68 MSparse<T>& operator = (const MSparse<T>& a) | |
69 { | |
70 Sparse<T>::operator = (a); | |
71 return *this; | |
72 } | |
73 | |
5275 | 74 MSparse<T>& insert (const Sparse<T>& a, octave_idx_type r, octave_idx_type c) |
5164 | 75 { |
76 Sparse<T>::insert (a, r, c); | |
77 return *this; | |
78 } | |
79 | |
6823 | 80 MSparse<T>& insert (const Sparse<T>& a, const Array<octave_idx_type>& indx) |
81 { | |
82 Sparse<T>::insert (a, indx); | |
83 return *this; | |
84 } | |
85 | |
5164 | 86 MSparse<T> transpose (void) const { return Sparse<T>::transpose (); } |
87 | |
88 MSparse<T> squeeze (void) const { return Sparse<T>::squeeze (); } | |
89 | |
10421
99e9bae2d81e
improve sparse indexing interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10350
diff
changeset
|
90 MSparse<T> index (const idx_vector& i, bool resize_ok) const |
5164 | 91 { return Sparse<T>::index (i, resize_ok); } |
92 | |
10421
99e9bae2d81e
improve sparse indexing interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10350
diff
changeset
|
93 MSparse<T> index (const idx_vector& i, const idx_vector& j, bool resize_ok) const |
5164 | 94 { return Sparse<T>::index (i, j, resize_ok); } |
95 | |
96 MSparse<T> reshape (const dim_vector& new_dims) const | |
97 { return Sparse<T>::reshape (new_dims); } | |
98 | |
5275 | 99 MSparse<T> permute (const Array<octave_idx_type>& vec, bool inv = false) const |
5164 | 100 { return Sparse<T>::permute (vec, inv); } |
101 | |
5275 | 102 MSparse<T> ipermute (const Array<octave_idx_type>& vec) const |
5164 | 103 { return Sparse<T>::ipermute (vec); } |
104 | |
105 | |
7620
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7602
diff
changeset
|
106 MSparse<T> diag (octave_idx_type k = 0) const |
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7602
diff
changeset
|
107 { |
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7602
diff
changeset
|
108 return Sparse<T>::diag (k); |
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7602
diff
changeset
|
109 } |
36594d5bbe13
Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents:
7602
diff
changeset
|
110 |
9812
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
111 // FIXME: should go away. |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
112 template <class U> |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
113 MSparse<U> |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
114 map (U (&fcn) (T)) const |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
115 { return Sparse<T>::template map<U> (fcn); } |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
116 |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
117 template <class U> |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
118 MSparse<U> |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
119 map (U (&fcn) (const T&)) const |
f80c566bc751
improve unary mapper system
Jaroslav Hajek <highegg@gmail.com>
parents:
8920
diff
changeset
|
120 { return Sparse<T>::template map<U> (fcn); } |
7602
7bfaa9611558
Rewrite sparse mappers in terms of a functor template function
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
121 |
5164 | 122 // Currently, the OPS functions don't need to be friends, but that |
123 // may change. | |
124 | |
10350
12884915a8e4
merge MArray classes & improve Array interface
Jaroslav Hajek <highegg@gmail.com>
parents:
10312
diff
changeset
|
125 // SPARSE_OPS_FRIEND_DECLS (MSparse, MArray) |
5164 | 126 }; |
127 | |
128 #endif |