From 224f4ba4835487b1dc08c114078aa57cd3597be5 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Fri, 31 Jul 2009 21:24:45 +0000
Subject: [PATCH] Provide a shortcut to get at the boundary object of a
 face/edge.

git-svn-id: https://svn.dealii.org/trunk@19152 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/deal.II/include/grid/tria_accessor.h       | 14 ++++++++++++++
 .../deal.II/include/grid/tria_accessor.templates.h | 13 +++++++++++++
 deal.II/doc/news/changes.h                         | 11 +++++++++++
 3 files changed, 38 insertions(+)

diff --git a/deal.II/deal.II/include/grid/tria_accessor.h b/deal.II/deal.II/include/grid/tria_accessor.h
index 0c09f3af59..47c97d1127 100644
--- a/deal.II/deal.II/include/grid/tria_accessor.h
+++ b/deal.II/deal.II/include/grid/tria_accessor.h
@@ -36,6 +36,8 @@ template <typename Accessor> class TriaRawIterator;
 template <typename Accessor> class TriaIterator;
 template <typename Accessor> class TriaActiveIterator;
 
+template <int dim, int spacedim> class Boundary;
+
 
 namespace internal
 {
@@ -1217,6 +1219,18 @@ class TriaAccessor : public TriaAccessorBase<structdim, dim, spacedim>
 				      */
     bool at_boundary () const;
 
+				     /**
+				      * Return a constant reference to a
+				      * boundary object used for this
+				      * object. This function is a shortcut to
+				      * retrieving the boundary indicator
+				      * using boundary_indicator() and then
+				      * asking the
+				      * Triangulation::get_boundary() function
+				      * for the boundary object.
+				      */
+    const Boundary<dim,spacedim> & get_boundary () const;
+    
 				     /**
 				      * @}
 				      */
diff --git a/deal.II/deal.II/include/grid/tria_accessor.templates.h b/deal.II/deal.II/include/grid/tria_accessor.templates.h
index 0fd98bff7c..f83877e04d 100644
--- a/deal.II/deal.II/include/grid/tria_accessor.templates.h
+++ b/deal.II/deal.II/include/grid/tria_accessor.templates.h
@@ -1862,6 +1862,19 @@ TriaAccessor<structdim, dim, spacedim>::at_boundary () const
 
 
 
+template <int structdim, int dim, int spacedim>
+const Boundary<dim,spacedim> &
+TriaAccessor<structdim, dim, spacedim>::get_boundary () const
+{
+  Assert (structdim<dim, ExcImpossibleInDim(dim));
+  Assert (this->used(), TriaAccessorExceptions::ExcCellNotUsed());
+
+  return this->tria->get_boundary(this->objects()
+				  .material_id[this->present_index]);
+}
+
+
+
 template <int structdim, int dim, int spacedim>
 double
 TriaAccessor<structdim, dim, spacedim>::diameter () const
diff --git a/deal.II/doc/news/changes.h b/deal.II/doc/news/changes.h
index 2e32cc1c11..6fc71ba726 100644
--- a/deal.II/doc/news/changes.h
+++ b/deal.II/doc/news/changes.h
@@ -276,6 +276,17 @@ inconvenience this causes.
 <h3>deal.II</h3>
 
 <ol>
+  <li>
+  <p>
+  New: Instead of asking face or edge iterators for their boundary indicator
+  using TriaAccessor::boundary_indicator() and then the triangulation for
+  the boundary object using Triangulation::get_boundary(), you can now directly
+  ask the iterator for the boundary object using TriaAccessor::get_boundary().
+  <br>
+  (WB 2009/07/31)
+  </p>
+  </li>
+
   <li>
   <p>
   Fixed: The CellAccessor::recursively_set_material_id function did not
-- 
2.39.5