3 #ifndef DUNE_RAVIARTTHOMASPREBASIS_HH
4 #define DUNE_RAVIARTTHOMASPREBASIS_HH
8 #include <dune/geometry/genericgeometry/topologytypes.hh>
14 template <
unsigned int dim,
class Field>
15 struct RTPreBasisFactory;
16 template <
unsigned int dim,
class Field>
22 typedef typename MBasisFactory::Object
MBasis;
27 typedef unsigned int Key;
31 template <
class Topology,
class Field >
34 template <
unsigned int dim,
class Field>
42 template <
unsigned int dd,
class FF>
47 template<
class Topology >
51 typename Traits::MBasis *mbasis = Traits::MBasisFactory::template create<Topology>(order+1);
52 typename remove_const<Object>::type *tmBasis =
53 new typename remove_const<Object>::type(*mbasis);
54 tmBasis->fill(vecMatrix);
58 template <
class Topology,
class Field>
61 static const unsigned int dim = Topology::dimension;
67 FieldVector< MI, dim > x;
68 for(
unsigned int i = 0; i <
dim; ++i )
70 std::vector< MI > val( basis.
size() );
74 unsigned int notHomogen = 0;
76 notHomogen = basis.
sizes()[order-1];
77 unsigned int homogen = basis.
sizes()[order]-notHomogen;
78 row_ = (notHomogen*dim+homogen*(dim+1))*
dim;
82 for (
unsigned int i=0; i<notHomogen+homogen; ++i)
84 for (
unsigned int r=0; r<
dim; ++r)
86 for (
unsigned int rr=0; rr<
dim; ++rr)
89 for (
unsigned int j=0; j<
col_; ++j)
99 for (
unsigned int i=0; i<homogen; ++i)
101 for (
unsigned int r=0; r<
dim; ++r)
104 for (
unsigned int j=0; j<
col_; ++j)
109 MI xval = val[notHomogen+i];
111 for (w=homogen+notHomogen; w<val.size(); ++w)
119 assert(w<val.size());
126 for (
unsigned int i=0; i<
rows(); ++i) {
137 template <
class Vector>
138 void row(
const unsigned int row, Vector &vec )
const
140 const unsigned int N =
cols();
141 assert( vec.size() == N );
142 for (
unsigned int i=0; i<N; ++i)
151 #endif // DUNE_RAVIARTTHOMASPREBASIS_HH