dune-common  2.3.0
identitymatrix.hh
Go to the documentation of this file.
1 #ifndef DUNE_COMMON_IDENTITYMATRIX_HH
2 #define DUNE_COMMON_IDENTITYMATRIX_HH
3 
4 #include <dune/common/fmatrix.hh>
5 #include <dune/common/ftraits.hh>
6 #include <dune/common/math.hh>
8 
17 namespace Dune
18 {
19 
20  // IdentityMatrix
21  // --------------
22 
34  template< class K, int N >
36  {
38  typedef K field_type;
40  typedef std::size_t size_type;
41 
43  DUNE_CONSTEXPR size_type rows () const { return N; }
45  DUNE_CONSTEXPR size_type cols () const { return N; }
46 
48  template< class X, class Y >
49  void mv ( const X &x, Y &y ) const
50  {
51  y = x;
52  }
53 
55  template< class X, class Y >
56  void mtv ( const X &x, Y &y ) const
57  {
58  y = x;
59  }
60 
62  template< class X, class Y >
63  void umv ( const X &x, Y &y ) const
64  {
65  y += x;
66  }
67 
69  template< class X, class Y >
70  void umtv ( const X &x, Y &y ) const
71  {
72  y += x;
73  }
74 
76  template< class X, class Y >
77  void umhv ( const X &x, Y &y ) const
78  {
79  y += x;
80  }
81 
83  template< class X, class Y >
84  void mmv ( const X &x, Y &y ) const
85  {
86  y -= x;
87  }
88 
90  template< class X, class Y >
91  void mmtv ( const X &x, Y &y ) const
92  {
93  y -= x;
94  }
95 
97  template< class X, class Y >
98  void mmhv ( const X &x, Y &y ) const
99  {
100  y -= x;
101  }
102 
104  template< class X, class Y >
105  void usmv ( const field_type &alpha, const X &x, Y &y ) const
106  {
107  y.axpy( alpha, x );
108  }
109 
111  template< class X, class Y >
112  void usmtv ( const field_type &alpha, const X &x, Y &y ) const
113  {
114  y.axpy( alpha, x );
115  }
116 
118  template< class X, class Y >
119  void usmhv ( const field_type &alpha, const X &x, Y &y ) const
120  {
121  y.axpy( alpha, x );
122  }
123 
126  {
127  return std::sqrt( frobenius_norm2() );
128  }
129 
132  {
134  }
135 
138  {
140  }
141 
144  {
146  }
147 
150  {
151  FieldMatrix< field_type, N, N > fieldMatrix( 0 );
152  for( int i = 0; i < N; ++i )
153  fieldMatrix[ i ][ i ] = field_type( 1 );
154  return fieldMatrix;
155  }
156  };
157 
158 } // namespace Dune
159 
160 #endif // #ifndef DUNE_COMMON_IDENTITYMATRIX_HH