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 #include "SundanceReducedQuadrature.hpp"
00032 #include "SundanceGaussianQuadrature.hpp"
00033
00034 namespace Sundance
00035 {
00036 using namespace Teuchos;
00037
00038
00039 ReducedQuadrature::ReducedQuadrature()
00040 : QuadratureFamilyBase(1)
00041 {}
00042
00043 XMLObject ReducedQuadrature::toXML() const
00044 {
00045 XMLObject rtn("ReducedQuadrature");
00046 return rtn;
00047 }
00048
00049
00050
00051 void ReducedQuadrature::getPoints(const CellType& cellType,
00052 Array<Point>& quadPoints,
00053 Array<double>& quadWeights) const
00054 {
00055 GaussianQuadrature basic(1);
00056 basic.getPoints(cellType, quadPoints, quadWeights);
00057 }
00058
00059 int ReducedQuadrature::getNumPoints( const CellType &cellType ) const
00060 {
00061 return 1;
00062 }
00063
00064 void ReducedQuadrature::getAdaptedWeights(const CellType& cellType ,
00065 int cellDim,
00066 int celLID ,
00067 int facetIndex ,
00068 const Mesh& mesh ,
00069 const ParametrizedCurve& globalCurve ,
00070 Array<Point>& quadPoints ,
00071 Array<double>& quadWeights ,
00072 bool &isCut) const
00073 {
00074 TEUCHOS_TEST_FOR_EXCEPT(true);
00075 }
00076
00077 }