dune-localfunctions  2.2.0
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,         // position
00059                                           std::vector<typename Traits::JacobianType>& out) const      // return value
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