raviartthomas02dlocalbasis.hh
Go to the documentation of this file.00001 #ifndef DUNE_RT0TRIANGLELOCALBASIS_HH
00002 #define DUNE_RT0TRIANGLELOCALBASIS_HH
00003
00004 #include <dune/common/fmatrix.hh>
00005
00006 #include <dune/localfunctions/common/localbasis.hh>
00007
00008 namespace Dune
00009 {
00018 template<class D, class R>
00019 class RT02DLocalBasis
00020 {
00021 public:
00022 typedef LocalBasisTraits<D,2,Dune::FieldVector<D,2>,R,2,Dune::FieldVector<R,2>,
00023 Dune::FieldMatrix<R,2,2> > Traits;
00024
00026 RT02DLocalBasis ()
00027 {
00028 sign0 = sign1 = sign2 = 1.0;
00029 }
00030
00032 RT02DLocalBasis (unsigned int s)
00033 {
00034 sign0 = sign1 = sign2 = 1.0;
00035 if (s&1) sign0 = -1.0;
00036 if (s&2) sign1 = -1.0;
00037 if (s&4) sign2 = -1.0;
00038 }
00039
00041 unsigned int size () const
00042 {
00043 return 3;
00044 }
00045
00047 inline void evaluateFunction (const typename Traits::DomainType& in,
00048 std::vector<typename Traits::RangeType>& out) const
00049 {
00050 out.resize(3);
00051 out[0][0] = sign0*in[0]; out[0][1]=sign0*(in[1]-1.0);
00052 out[1][0] = sign1*(in[0]-1.0); out[1][1]=sign1*in[1];
00053 out[2][0] = sign2*in[0]; out[2][1]=sign2*in[1];
00054 }
00055
00057 inline void
00058 evaluateJacobian (const typename Traits::DomainType& in,
00059 std::vector<typename Traits::JacobianType>& out) const
00060 {
00061 out.resize(3);
00062 out[0][0][0] = sign0; out[0][0][1] = 0;
00063 out[0][1][0] = 0; out[0][1][1] = sign0;
00064 out[1][0][0] = sign1; out[1][0][1] = 0;
00065 out[1][1][0] = 0; out[1][1][1] = sign1;
00066 out[2][0][0] = sign2; out[2][0][1] = 0;
00067 out[2][1][0] = 0; out[2][1][1] = sign2;
00068 }
00069
00071 unsigned int order () const
00072 {
00073 return 1;
00074 }
00075
00076 private:
00077 R sign0, sign1, sign2;
00078 };
00079 }
00080 #endif