sparse_shell_matrix.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2014 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 #ifndef LIBMESH_SPARSE_SHELL_MATRIX_H
21 #define LIBMESH_SPARSE_SHELL_MATRIX_H
22 
23 
24 // Local includes
25 #include "libmesh/libmesh_common.h"
27 #include "libmesh/libmesh.h"
28 #include "libmesh/shell_matrix.h"
29 #include "libmesh/sparse_matrix.h"
30 
31 namespace libMesh
32 {
33 
34 
41 template <typename T>
42 class SparseShellMatrix : public ShellMatrix<T>
43 {
44 public:
49  explicit
50  SparseShellMatrix (const SparseMatrix<T>& new_m);
51 
55  virtual ~SparseShellMatrix ();
56 
61  virtual numeric_index_type m () const;
62 
67  virtual numeric_index_type n () const;
68 
73  virtual void vector_mult (NumericVector<T>& dest,
74  const NumericVector<T>& arg) const;
75 
79  virtual void vector_mult_add (NumericVector<T>& dest,
80  const NumericVector<T>& arg) const;
81 
85  virtual void get_diagonal (NumericVector<T>& dest) const;
86 
87 protected:
92 
93 };
94 
95 
96 
97 //-----------------------------------------------------------------------
98 // SparseShellMatrix inline members
99 template <typename T>
100 inline
102  ShellMatrix<T>(new_m.comm()),
103  _m(new_m)
104 {}
105 
106 
107 
108 template <typename T>
109 inline
111 {}
112 
113 
114 
115 template <typename T>
116 inline
118 {
119  return _m.m();
120 }
121 
122 
123 
124 template <typename T>
125 inline
127 {
128  return _m.n();
129 }
130 
131 
132 
133 template <typename T>
134 inline
136 {
137  _m.get_diagonal(dest);
138 }
139 
140 
141 } // namespace libMesh
142 
143 
144 #endif // LIBMESH_SPARSE_SHELL_MATRIX_H

Site Created By: libMesh Developers
Last modified: February 07 2014 16:57:07 UTC

Hosted By:
SourceForge.net Logo