00001
00002
00003
00004
00005
00006
00007
00008 #ifndef SUNDANCETRAPESOIDQUADRATURE_HPP_
00009 #define SUNDANCETRAPESOIDQUADRATURE_HPP_
00010
00011 #include "SundanceDefs.hpp"
00012 #include "SundanceQuadratureFamilyBase.hpp"
00013
00014 namespace Sundance {
00015
00016 class TrapesoidQuadrature : public QuadratureFamilyBase {
00017
00018 public:
00019
00020
00021 TrapesoidQuadrature(int resolution);
00022
00023
00024 virtual ~TrapesoidQuadrature()
00025 {;}
00026
00027
00028 virtual XMLObject toXML() const;
00029
00030
00031 virtual std::string description() const
00032 {
00033 return "TrapesoidQuadrature[order=" + Teuchos::toString(order()) + "]";
00034 }
00035
00036
00037 GET_RCP(QuadratureFamilyStub);
00038
00039 protected:
00040
00041 virtual void getLineRule(Array<Point>& quadPoints,
00042 Array<double>& quadWeights) const;
00043
00044
00045 virtual void getTriangleRule(Array<Point>& quadPoints,
00046 Array<double>& quadWeights) const;
00047
00048
00049 virtual void getQuadRule(Array<Point>& quadPoints,
00050 Array<double>& quadWeights) const;
00051
00052
00053 virtual void getTetRule(Array<Point>& quadPoints,
00054 Array<double>& quadWeights) const ;
00055
00056
00057 virtual void getBrickRule(Array<Point>& quadPoints,
00058 Array<double>& quadWeights) const;
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073 virtual void getAdaptedWeights(const CellType& cellType, int cellDim,
00074 int cellLID, int facetIndex, const Mesh& mesh,
00075 const ParametrizedCurve& globalCurve,
00076 Array<Point>& quadPoints, Array<double>& quadWeights,
00077 bool &weightsChanged) const;
00078
00079 private:
00080
00081
00082 int resolution_;
00083
00084 };
00085
00086 }
00087
00088 #endif