00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #ifndef SUNDANCE_LAGRANGE_H
00032 #define SUNDANCE_LAGRANGE_H
00033
00034 #include "SundanceDefs.hpp"
00035 #include "Teuchos_RefCountPtr.hpp"
00036 #include "SundanceBasisFamilyBase.hpp"
00037
00038 namespace Sundance
00039 {
00040
00041
00042
00043 class Lagrange : public ScalarBasis
00044 {
00045 public:
00046
00047 Lagrange(int order);
00048
00049
00050
00051
00052 bool supportsCellTypePair(
00053 const CellType& maximalCellType,
00054 const CellType& cellType
00055 ) const ;
00056
00057
00058 void print(std::ostream& os) const ;
00059
00060
00061 int order() const {return order_;}
00062
00063
00064 int nReferenceDOFsWithoutFacets(
00065 const CellType& maximalCellType,
00066 const CellType& cellType
00067 ) const ;
00068
00069
00070 void getReferenceDOFs(
00071 const CellType& maximalCellType,
00072 const CellType& cellType,
00073 Array<Array<Array<int> > >& dofs) const ;
00074
00075
00076 void refEval(
00077 const CellType& cellType,
00078 const Array<Point>& pts,
00079 const SpatialDerivSpecifier& deriv,
00080 Array<Array<Array<double> > >& result,
00081 int verbosity=0) const ;
00082
00083
00084 void getConstrainsForHNDoF(
00085 const int indexInParent,
00086 const int maxCellDim,
00087 const int maxNrChild,
00088 const int facetDim,
00089 const int facetIndex,
00090 const int nodeIndex,
00091 Array<int>& localDoFs,
00092 Array<int>& parentFacetDim,
00093 Array<int>& parentFacetIndex,
00094 Array<int>& parentFacetNode,
00095 Array<double>& coefs
00096 );
00097
00098
00099 GET_RCP(BasisFamilyBase);
00100
00101
00102 std::string description() const
00103 {return "Lagrange(" + Teuchos::toString(order()) + ")";}
00104
00105 private:
00106 int order_;
00107
00108 static Array<int> makeRange(int low, int high);
00109
00110
00111 void evalOnLine(const Point& pt,
00112 const MultiIndex& deriv,
00113 Array<double>& result) const ;
00114
00115
00116 void evalOnTriangle(const Point& pt,
00117 const MultiIndex& deriv,
00118 Array<double>& result) const ;
00119
00120
00121 void evalOnquad(const Point& pt,
00122 const MultiIndex& deriv,
00123 Array<double>& result) const ;
00124
00125
00126 void evalOnTet(const Point& pt,
00127 const MultiIndex& deriv,
00128 Array<double>& result) const ;
00129
00130
00131 void evalOnBrick(const Point& pt,
00132 const MultiIndex& deriv,
00133 Array<double>& result) const ;
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143 void returnDoFPositionOnRef(
00144 const int maxCellDim,
00145 const int facetDim,
00146 const int facetIndex,
00147 const int nodeIndex,
00148 Point& pos) const;
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158 void getDoFsOrdered(
00159 const CellType maxCellDim,
00160 int nrDoF,
00161 Array<int>& facetD,
00162 Array<int>& facetI,
00163 Array<int>& facetN);
00164
00165
00166 bool doFInfromationCalculated_;
00167
00168
00169 Array<int> facetD_;
00170
00171
00172 Array<int> facetI_;
00173
00174
00175 Array<int> facetN_;
00176 };
00177 }
00178
00179 #endif