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_GAUSSIANQUADRATURE_H
00032 #define SUNDANCE_GAUSSIANQUADRATURE_H
00033
00034 #include "SundanceDefs.hpp"
00035 #include "Teuchos_RefCountPtr.hpp"
00036 #include "SundanceQuadratureFamilyBase.hpp"
00037
00038 namespace Sundance
00039 {
00040 using namespace Teuchos;
00041
00042
00043
00044
00045
00046
00047 class GaussianQuadrature : public QuadratureFamilyBase
00048 {
00049 public:
00050
00051 GaussianQuadrature(int order);
00052
00053
00054 virtual ~GaussianQuadrature(){;}
00055
00056
00057
00058 virtual XMLObject toXML() const ;
00059
00060
00061 virtual std::string description() const
00062 {return "GaussianQuadrature[order=" + Teuchos::toString(order())
00063 + "]";}
00064
00065
00066 GET_RCP(QuadratureFamilyStub);
00067
00068
00069 virtual void getAdaptedWeights(const CellType& cellType ,
00070 int cellDim,
00071 int celLID ,
00072 int facetIndex ,
00073 const Mesh& mesh ,
00074 const ParametrizedCurve& globalCurve ,
00075 Array<Point>& quadPoints ,
00076 Array<double>& quadWeights ,
00077 bool &isCut) const ;
00078
00079 protected:
00080
00081 virtual void getLineRule(Array<Point>& quadPoints,
00082 Array<double>& quadWeights) const ;
00083
00084
00085 virtual void getTriangleRule(Array<Point>& quadPoints,
00086 Array<double>& quadWeights) const ;
00087
00088
00089 virtual void getQuadRule(Array<Point>& quadPoints,
00090 Array<double>& quadWeights) const ;
00091
00092
00093 virtual void getTetRule(Array<Point>& quadPoints,
00094 Array<double>& quadWeights) const ;
00095
00096
00097 virtual void getBrickRule(Array<Point>& quadPoints,
00098 Array<double>& quadWeights) const ;
00099
00100 };
00101 }
00102
00103
00104 #endif