- Home
- About DUNE
- Download
- Documentation
- Community
- Development
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
Generated on Fri Apr 29 2011 with Doxygen (ver 1.7.1) [doxygen-log,error-log].