* weak formulation of your discontinuous Galerkin scheme.
* @author Ralf Hartmann, 1998
*/
-
-
template <int dim>
class FEDGLinear : public FELinear<dim>{
public:
*/
FEDGLinear();
+ /**
+ * Refer to the base class for detailed
+ * information on this function.
+ */
+ virtual void get_face_support_points (const DoFHandler<dim>::face_iterator &face,
+ const Boundary<dim> &boundary,
+ vector<Point<dim> > &support_points) const;
+
+
/**
* This function returns an error since the
* correct use of the restriction
* weak formulation of your discontinuous Galerkin scheme.
* @author Ralf Hartmann, 1998
*/
-
-
template <int dim>
class FEDGQuadraticSub : public FEQuadraticSub<dim>{
public:
* Constructor
*/
FEDGQuadraticSub();
+
+ /**
+ * Refer to the base class for detailed
+ * information on this function.
+ */
+ virtual void get_face_support_points (const DoFHandler<dim>::face_iterator &face,
+ const Boundary<dim> &boundary,
+ vector<Point<dim> > &support_points) const;
/**
* This function returns an error since the
* weak formulation of your discontinuous Galerkin scheme.
* @author Ralf Hartmann, 1998
*/
-
-
template <int dim>
class FEDGCubicSub : public FECubicSub<dim>{
public:
*/
FEDGCubicSub();
+ /**
+ * Refer to the base class for detailed
+ * information on this function.
+ */
+ virtual void get_face_support_points (const DoFHandler<dim>::face_iterator &face,
+ const Boundary<dim> &boundary,
+ vector<Point<dim> > &support_points) const;
+
/**
* This function returns an error since the
* correct use of the restriction
* weak formulation of your discontinuous Galerkin scheme.
* @author Ralf Hartmann, 1998
*/
-
-
template <int dim>
class FEDGQuarticSub : public FEQuarticSub<dim>{
public:
*/
FEDGQuarticSub();
+ /**
+ * Refer to the base class for detailed
+ * information on this function.
+ */
+ virtual void get_face_support_points (const DoFHandler<dim>::face_iterator &face,
+ const Boundary<dim> &boundary,
+ vector<Point<dim> > &support_points) const;
+
/**
* This function returns an error since the
* correct use of the restriction
+template <int dim>
+void
+FEDGLinear<dim>::get_face_support_points (const typename DoFHandler<dim>::face_iterator &,
+ const Boundary<dim> &,
+ vector<Point<dim> > &support_points) const {
+ Assert ((support_points.size() == 0),
+ ExcWrongFieldDimension (support_points.size(),0));
+};
+
+
template <int dim>
const FullMatrix<double> &
FEDGLinear<dim>::restrict (const unsigned int child) const {
};
+
+template <int dim>
+void
+FEDGQuadraticSub<dim>::get_face_support_points (const typename DoFHandler<dim>::face_iterator &,
+ const Boundary<dim> &,
+ vector<Point<dim> > &support_points) const {
+ Assert ((support_points.size() == 0),
+ ExcWrongFieldDimension (support_points.size(),0));
+};
+
+
template <int dim>
const FullMatrix<double> &
FEDGQuadraticSub<dim>::restrict (const unsigned int child) const {
+template <int dim>
+void
+FEDGCubicSub<dim>::get_face_support_points (const typename DoFHandler<dim>::face_iterator &,
+ const Boundary<dim> &,
+ vector<Point<dim> > &support_points) const {
+ Assert ((support_points.size() == 0),
+ ExcWrongFieldDimension (support_points.size(),0));
+};
+
+
template <int dim>
const FullMatrix<double> &
FEDGCubicSub<dim>::restrict (const unsigned int child) const {
+
+template <int dim>
+void
+FEDGQuarticSub<dim>::get_face_support_points (const typename DoFHandler<dim>::face_iterator &,
+ const Boundary<dim> &,
+ vector<Point<dim> > &support_points) const {
+ Assert ((support_points.size() == 0),
+ ExcWrongFieldDimension (support_points.size(),0));
+};
+
+
template <int dim>
const FullMatrix<double> &
FEDGQuarticSub<dim>::restrict (const unsigned int child) const {
};
-
template <int dim>
void
FEDGConstant<dim>::get_face_support_points (const typename DoFHandler<dim>::face_iterator &,
- const Boundary<dim> &,
- vector<Point<dim> > &support_points) const {
- Assert (false, ExcNotImplemented());
-
- Assert ((support_points.size() == dofs_per_face) &&
- (support_points.size() == GeometryInfo<dim>::vertices_per_face),
- ExcWrongFieldDimension (support_points.size(),
- GeometryInfo<dim>::vertices_per_face));
+ const Boundary<dim> &,
+ vector<Point<dim> > &support_points) const {
+ Assert ((support_points.size() == 0),
+ ExcWrongFieldDimension (support_points.size(),0));
};
+
template <int dim>
const FullMatrix<double> &
FEDGConstant<dim>::restrict (const unsigned int child) const {