libMesh
tensor_shell_matrix.h
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2017 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_TENSOR_SHELL_MATRIX_H
21 #define LIBMESH_TENSOR_SHELL_MATRIX_H
22 
23 
24 // Local includes
25 #include "libmesh/libmesh_common.h"
26 #include "libmesh/reference_counted_object.h"
27 #include "libmesh/libmesh.h"
28 #include "libmesh/shell_matrix.h"
29 #include "libmesh/numeric_vector.h"
30 
31 namespace libMesh
32 {
33 
42 template <typename T>
43 class TensorShellMatrix : public ShellMatrix<T>
44 {
45 public:
51  const NumericVector<T> & w);
52 
56  virtual ~TensorShellMatrix ();
57 
58  virtual numeric_index_type m () const libmesh_override;
59 
60  virtual numeric_index_type n () const libmesh_override;
61 
62  virtual void vector_mult (NumericVector<T> & dest,
63  const NumericVector<T> & arg) const libmesh_override;
64 
65  virtual void vector_mult_add (NumericVector<T> & dest,
66  const NumericVector<T> & arg) const libmesh_override;
67 
68  virtual void get_diagonal (NumericVector<T> & dest) const libmesh_override;
69 
70 protected:
75 
80 };
81 
82 
83 
84 //-----------------------------------------------------------------------
85 // TensorShellMatrix inline members
86 template <typename T>
87 inline
89  const NumericVector<T> & w):
90  ShellMatrix<T>(v.comm()),
91  _v(v),
92  _w(w)
93 {}
94 
95 
96 
97 template <typename T>
98 inline
100 {}
101 
102 
103 
104 template <typename T>
105 inline
107 {
108  return _v.size();
109 }
110 
111 
112 
113 template <typename T>
114 inline
116 {
117  return _w.size();
118 }
119 
120 
121 } // namespace libMesh
122 
123 
124 #endif // LIBMESH_TENSOR_SHELL_MATRIX_H
virtual numeric_index_type m() const libmesh_override
const NumericVector< T > & _w
The row vector.
Numeric vector.
Definition: dof_map.h:66
virtual void get_diagonal(NumericVector< T > &dest) const libmesh_override
Copies the diagonal part of the matrix into dest.
The libMesh namespace provides an interface to certain functionality in the library.
virtual ~TensorShellMatrix()
Destructor.
Shell matrix that is given by a tensor product of two vectors, i.e.
virtual numeric_index_type n() const libmesh_override
dof_id_type numeric_index_type
Definition: id_types.h:92
virtual void vector_mult_add(NumericVector< T > &dest, const NumericVector< T > &arg) const libmesh_override
Multiplies the matrix with arg and adds the result to dest.
TensorShellMatrix(const NumericVector< T > &v, const NumericVector< T > &w)
Constructor; takes references to the two vectors as arguments.
virtual void vector_mult(NumericVector< T > &dest, const NumericVector< T > &arg) const libmesh_override
Multiplies the matrix with arg and stores the result in dest.
const Parallel::Communicator & comm() const
const NumericVector< T > & _v
The column vector.
Generic shell matrix, i.e.