00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024 #if defined (__GNUG__) && defined (USE_PRAGMA_INTERFACE_IMPLEMENTATION)
00025 #pragma interface
00026 #endif
00027
00028 #if !defined (octave_MDiagArray2_h)
00029 #define octave_MDiagArray2_h 1
00030
00031 #include "DiagArray2.h"
00032 #include "MArray2.h"
00033
00034
00035
00036
00037
00038 #include "MArray-defs.h"
00039
00040 MDIAGARRAY2_OPS_FORWARD_DECLS (MDiagArray2)
00041
00042 template <class T>
00043 class
00044 MDiagArray2 : public DiagArray2<T>
00045 {
00046 protected:
00047
00048 MDiagArray2 (T *d, int r, int c) : DiagArray2<T> (d, r, c) { }
00049
00050 public:
00051
00052 MDiagArray2 (void) : DiagArray2<T> () { }
00053
00054 MDiagArray2 (int r, int c) : DiagArray2<T> (r, c) { }
00055
00056 MDiagArray2 (int r, int c, const T& val) : DiagArray2<T> (r, c, val) { }
00057
00058 MDiagArray2 (const MDiagArray2<T>& a) : DiagArray2<T> (a) { }
00059
00060 MDiagArray2 (const DiagArray2<T>& a) : DiagArray2<T> (a) { }
00061
00062 explicit MDiagArray2 (const Array<T>& a) : DiagArray2<T> (a) { }
00063
00064 ~MDiagArray2 (void) { }
00065
00066 MDiagArray2<T>& operator = (const MDiagArray2<T>& a)
00067 {
00068 DiagArray2<T>::operator = (a);
00069 return *this;
00070 }
00071
00072 operator MArray2<T> () const
00073 {
00074 int nr = DiagArray2<T>::dim1 ();
00075 int nc = DiagArray2<T>::dim2 ();
00076
00077 MArray2<T> retval (nr, nc, T (0));
00078
00079 int len = nr < nc ? nr : nc;
00080
00081 for (int i = 0; i < len; i++)
00082 retval.xelem (i, i) = this->xelem (i, i);
00083
00084 return retval;
00085 }
00086
00087 static MDiagArray2<T> nil_array;
00088
00089
00090
00091
00092
00093
00094 };
00095
00096 #endif
00097
00098
00099
00100
00101
00102