5164
|
1 /* |
|
2 |
|
3 Copyright (C) 2004 David Bateman |
|
4 Copyright (C) 1998-2004 Andy Adler |
|
5 |
|
6 Octave is free software; you can redistribute it and/or modify it |
|
7 under the terms of the GNU General Public License as published by the |
|
8 Free Software Foundation; either version 2, or (at your option) any |
|
9 later version. |
|
10 |
|
11 Octave is distributed in the hope that it will be useful, but WITHOUT |
|
12 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
|
13 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
|
14 for more details. |
|
15 |
|
16 You should have received a copy of the GNU General Public License |
|
17 along with this program; see the file COPYING. If not, write to the Free |
|
18 Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
|
19 |
|
20 */ |
|
21 |
|
22 |
|
23 #if !defined (octave_sparse_base_lu_h) |
|
24 #define octave_sparse_base_lu_h 1 |
|
25 |
|
26 #include "MArray.h" |
|
27 |
|
28 template <class lu_type, class lu_elt_type, class p_type, class p_elt_type> |
|
29 class |
|
30 sparse_base_lu |
|
31 { |
|
32 public: |
|
33 |
|
34 sparse_base_lu (void) { } |
|
35 |
|
36 sparse_base_lu (const sparse_base_lu& a) |
|
37 : Lfact (a.Lfact), Ufact (a.Ufact), cond (a.cond), P (a.P), Q (a.Q) { } |
|
38 |
|
39 sparse_base_lu& operator = (const sparse_base_lu& a) |
|
40 { |
|
41 if (this != &a) |
|
42 { |
|
43 Lfact = a.Lfact; |
|
44 Ufact = a.Ufact; |
|
45 cond = a.cond; |
|
46 P = a.P; |
|
47 Q = a.Q; |
|
48 } |
|
49 return *this; |
|
50 } |
|
51 |
|
52 ~sparse_base_lu (void) { } |
|
53 |
|
54 lu_type L (void) const { return Lfact; } |
|
55 |
|
56 lu_type U (void) const { return Ufact; } |
|
57 |
|
58 p_type Pc (void) const; |
|
59 |
|
60 p_type Pr (void) const; |
|
61 |
|
62 MArray<int> row_perm (void) const { return P; } |
|
63 |
|
64 MArray<int> col_perm (void) const { return Q; } |
|
65 |
|
66 double rcond (void) const { return cond; } |
|
67 |
|
68 protected: |
|
69 |
|
70 lu_type Lfact; |
|
71 lu_type Ufact; |
|
72 |
|
73 double cond; |
|
74 |
|
75 MArray<int> P; |
|
76 MArray<int> Q; |
|
77 }; |
|
78 |
|
79 #endif |
|
80 |
|
81 /* |
|
82 ;;; Local Variables: *** |
|
83 ;;; mode: C++ *** |
|
84 ;;; End: *** |
|
85 */ |