00001 #ifndef SERENDIPITY_H
00002 #define SERENDIPITY_H
00003
00004 #include "SundanceDefs.h"
00005
00006 #include "ScalarBasis.h"
00007 #include "XMLObject.h"
00008
00009
00010 namespace Sundance
00011 {
00012
00013 using namespace TSF;
00014 using std::string;
00015
00016 using std::ostream;
00017
00018
00019
00020
00021
00022 class Serendipity : public ScalarBasis
00023 {
00024 public:
00025
00026
00027
00028 Serendipity(int order);
00029
00030
00031
00032
00033
00034 Serendipity(const XMLObject& object);
00035
00036
00037
00038
00039 virtual void setNodeNumbers(const Cell& cell,
00040 const TSFArray<int>& funcID,
00041 TemporaryMapper& tmpMap) const ;
00042
00043
00044
00045
00046 virtual void getNodeNumbers(const Cell& cell,
00047 const TSFArray<int>& funcID,
00048 TemporaryMapper& tmpMap,
00049 TSFArray<TSFArray<int> >& indices) const ;
00050
00051
00052
00053
00054 virtual void getFacetDOFs(const Cell& c, int facetDim, int facetIndex,
00055 const TSFArray<int>& maximalDOFs,
00056 TSFArray<int>& facetDOFs) const ;
00057
00058
00059
00060
00061
00062 virtual void getNodalPoints(const Cell& cell,
00063 const TSFArray<bool>& activeNodes,
00064 TSFArray<Point>& evalPoints) const ;
00065
00066 protected:
00067
00068 virtual int nNodesOnLine() const ;
00069
00070
00071 virtual int nNodesOnTriangle() const ;
00072
00073
00074 virtual int nNodesOnQuad() const ;
00075
00076
00077 virtual int nNodesOnTet() const ;
00078
00079
00080 virtual int nNodesOnBrick() const ;
00081
00082
00083 virtual void evalOnLine(const Point& pt,
00084 const MultiIndex& deriv,
00085 DenseSerialVector& result) const ;
00086
00087
00088 virtual void evalOnTriangle(const Point& pt,
00089 const MultiIndex& deriv,
00090 DenseSerialVector& result) const ;
00091
00092
00093 virtual void evalOnQuad(const Point& pt,
00094 const MultiIndex& deriv,
00095 DenseSerialVector& result) const ;
00096
00097
00098 virtual void evalOnTet(const Point& pt,
00099 const MultiIndex& deriv,
00100 DenseSerialVector& result) const ;
00101
00102
00103 virtual void evalOnBrick(const Point& pt,
00104 const MultiIndex& deriv,
00105 DenseSerialVector& result) const ;
00106 };
00107
00108 extern "C"
00109 {
00110
00111 BasisFamilyBase* createXMLSerendipity(const XMLObject& xml);
00112 }
00113
00114 }
00115 #endif