dune-common  2.3.1
dynmatrix.hh
Go to the documentation of this file.
1 // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=4 sw=2 sts=2:
3 // $Id: fmatrix.hh 6181 2010-10-13 18:53:40Z christi $
4 #ifndef DUNE_DYNMATRIX_HH
5 #define DUNE_DYNMATRIX_HH
6 
7 #include <cmath>
8 #include <cstddef>
9 #include <iostream>
10 
12 #include <dune/common/dynvector.hh>
15 
16 namespace Dune
17 {
18 
28  template< class K > class DynamicMatrix;
29 
30  template< class K >
32  {
34 
36 
38  typedef const row_type &const_row_reference;
39 
40  typedef std::vector<K> container_type;
41  typedef K value_type;
42  typedef typename container_type::size_type size_type;
43  };
44 
45  template< class K >
47  {
50  };
51 
56  template<class K>
57  class DynamicMatrix : public DenseMatrix< DynamicMatrix<K> >
58  {
59  std::vector< DynamicVector<K> > _data;
61  public:
62  typedef typename Base::size_type size_type;
63  typedef typename Base::value_type value_type;
64  typedef typename Base::row_type row_type;
65 
66  //===== constructors
69 
72  _data(r, row_type(c, v) )
73  {}
74 
75  //==== resize related methods
77  {
78  _data.resize(0);
79  _data.resize(r, row_type(c, v) );
80  }
81 
82  //===== assignment
83  using Base::operator=;
84 
85  // make this thing a matrix
86  size_type mat_rows() const { return _data.size(); }
87  size_type mat_cols() const {
88  assert(this->rows());
89  return _data.front().size();
90  }
91  row_type & mat_access(size_type i) { return _data[i]; }
92  const row_type & mat_access(size_type i) const { return _data[i]; }
93  };
94 
97 } // end namespace
98 
99 #endif