From: wolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 16 Dec 1999 14:27:00 +0000 (+0000)
Subject: Add HalfHyperBallBoundary and move implementation of constructors of HyperBallBoundar... 
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2b302ed4a6265371ef95490741d892312cb8c874;p=dealii-svn.git

Add HalfHyperBallBoundary and move implementation of constructors of HyperBallBoundary to .cc file.


git-svn-id: https://svn.dealii.org/trunk@2073 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/deal.II/include/grid/tria_boundary_lib.h b/deal.II/deal.II/include/grid/tria_boundary_lib.h
index 5f4839d907..1145323449 100644
--- a/deal.II/deal.II/include/grid/tria_boundary_lib.h
+++ b/deal.II/deal.II/include/grid/tria_boundary_lib.h
@@ -31,8 +31,8 @@ class HyperBallBoundary : public StraightBoundary<dim> {
 				     /**
 				      * Constructor
 				      */
-    HyperBallBoundary (const Point<dim> p=Point<dim>(), const double radius=1.0) :
-		    center(p), radius(radius) {};
+    HyperBallBoundary (const Point<dim> p      = Point<dim>(),
+		       const double     radius = 1.0);
 
 				     /**
 				      * Refer to the general documentation of
@@ -51,7 +51,7 @@ class HyperBallBoundary : public StraightBoundary<dim> {
     get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const;
 
 
-  private:
+  protected:
 				     /**
 				      * Center point of the hyperball.
 				      */
@@ -63,6 +63,42 @@ class HyperBallBoundary : public StraightBoundary<dim> {
     const double radius;
 };
 
+
+
+/**
+ * Variant of #HyperBallBoundary# which denotes a half hyper ball
+ * where the first coordinate is restricted to the range $x>=0$ (or
+ * $x>=center(0)$). In two dimensions, this equals the right half
+ * circle, in three space dimensions it is a half ball.
+ *
+ * @author Wolfgang Bangerth, 1999
+ */
+template <int dim>
+class HalfHyperBallBoundary : public HyperBallBoundary<dim> {
+  public:
+				     /**
+				      * Constructor
+				      */
+    HalfHyperBallBoundary (const Point<dim> p      = Point<dim>(),
+			   const double     radius = 1.0);
+
+				     /**
+				      * Check if on the line #x==0#,
+				      * otherwise pass to the base
+				      * class.
+				      */
+    virtual Point<dim>
+    get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const;
+
+				     /**
+				      * Check if on the line #x==0#,
+				      * otherwise pass to the base
+				      * class.
+				      */
+    virtual Point<dim>
+    get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const;
+};
+
     
 
 
diff --git a/deal.II/deal.II/source/grid/tria_boundary_lib.cc b/deal.II/deal.II/source/grid/tria_boundary_lib.cc
index d2d012b320..2aed965181 100644
--- a/deal.II/deal.II/source/grid/tria_boundary_lib.cc
+++ b/deal.II/deal.II/source/grid/tria_boundary_lib.cc
@@ -10,6 +10,14 @@
 
 
 
+template <int dim>
+HyperBallBoundary<dim>::HyperBallBoundary (const Point<dim> p,
+					   const double     radius) :
+		center(p), radius(radius)
+{};
+
+
+
 template <int dim>
 Point<dim>
 HyperBallBoundary<dim>::get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const
@@ -59,6 +67,42 @@ get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) c
 
 
 
+template <int dim>
+HalfHyperBallBoundary<dim>::HalfHyperBallBoundary (const Point<dim> center,
+						   const double     radius) :
+		HyperBallBoundary<dim> (center, radius)
+{};
+
+
+
+template <int dim>
+Point<dim>
+HalfHyperBallBoundary<dim>::
+get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const 
+{
+  const Point<dim> line_center = line->center();
+  if (line_center(0) == center(0))
+    return line_center;
+  else
+    return HyperBallBoundary<dim>::get_new_point_on_line (line);
+};
+
+
+
+
+template <int dim>
+Point<dim>
+HalfHyperBallBoundary<dim>::
+get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const 
+{
+  const Point<dim> quad_center = quad->center();
+  if (quad_center(0) == center(0))
+    return quad_center;
+  else
+    return HyperBallBoundary<dim>::get_new_point_on_quad (quad);
+};
+
+
 
 template <int dim>
 HyperShellBoundary<dim>::HyperShellBoundary (const Point<dim> &center) :
@@ -129,4 +173,5 @@ get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) c
 
 // explicit instantiations
 template class HyperBallBoundary<deal_II_dimension>;
+template class HalfHyperBallBoundary<deal_II_dimension>;
 template class HyperShellBoundary<deal_II_dimension>;