1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
// Template array classes with like-type math ops
/*
Copyright (C) 1996, 1997 John W. Eaton
This file is part of Octave.
Octave is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 2, or (at your option) any
later version.
Octave is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with Octave; see the file COPYING. If not, write to the Free
Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
02110-1301, USA.
*/
#if defined (__GNUG__) && defined (USE_PRAGMA_INTERFACE_IMPLEMENTATION)
#pragma interface
#endif
#if !defined (octave_MDiagArray2_h)
#define octave_MDiagArray2_h 1
#include "DiagArray2.h"
#include "MArray2.h"
// Two dimensional diagonal array with math ops.
// But first, some preprocessor abuse...
#include "MArray-defs.h"
MDIAGARRAY2_OPS_FORWARD_DECLS (MDiagArray2)
template <class T>
class
MDiagArray2 : public DiagArray2<T>
{
protected:
MDiagArray2 (T *d, int r, int c) : DiagArray2<T> (d, r, c) { }
public:
MDiagArray2 (void) : DiagArray2<T> () { }
MDiagArray2 (int r, int c) : DiagArray2<T> (r, c) { }
MDiagArray2 (int r, int c, const T& val) : DiagArray2<T> (r, c, val) { }
MDiagArray2 (const MDiagArray2<T>& a) : DiagArray2<T> (a) { }
MDiagArray2 (const DiagArray2<T>& a) : DiagArray2<T> (a) { }
explicit MDiagArray2 (const Array<T>& a) : DiagArray2<T> (a) { }
~MDiagArray2 (void) { }
MDiagArray2<T>& operator = (const MDiagArray2<T>& a)
{
DiagArray2<T>::operator = (a);
return *this;
}
operator MArray2<T> () const
{
int nr = DiagArray2<T>::dim1 ();
int nc = DiagArray2<T>::dim2 ();
MArray2<T> retval (nr, nc, T (0));
int len = nr < nc ? nr : nc;
for (int i = 0; i < len; i++)
retval.xelem (i, i) = this->xelem (i, i);
return retval;
}
static MDiagArray2<T> nil_array;
// Currently, the OPS functions don't need to be friends, but that
// may change.
// MDIAGARRAY2_OPS_FRIEND_DECLS (MDiagArray2)
};
#endif
/*
;;; Local Variables: ***
;;; mode: C++ ***
;;; End: ***
*/
|
