]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix the behaviour of the new boundary a bit and update the documentation.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 18 Feb 1999 16:06:53 +0000 (16:06 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 18 Feb 1999 16:06:53 +0000 (16:06 +0000)
git-svn-id: https://svn.dealii.org/trunk@838 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/grid/tria_boundary.h
deal.II/deal.II/source/grid/tria_boundary.cc

index ec61dc0a1c03f6095745718b41a02d6df87b8cc9..d6153543fea4681389b55a078a9aadae56879771 100644 (file)
@@ -132,17 +132,26 @@ template <int dim>
 class StraightBoundary : public Boundary<dim> {
   public:
                                     /**
+                                     * Let the new point be the arithmetic
+                                     * mean of the two vertices of the line.
+                                     *
                                      * Refer to the general documentation of
                                      * this class and the documentation of the
-                                     * base class.
+                                     * base class for more information.
                                      */
     virtual Point<dim>
     get_new_point_on_line (const typename Triangulation<dim>::line_iterator &line) const;
 
                                     /**
+                                     * Let the new point be the arithmetic mean
+                                     * of the four vertices of this quad and
+                                     * the four midpoints of the lines, which
+                                     * are already created at the time of calling
+                                     * this function.
+                                     *
                                      * Refer to the general documentation of
                                      * this class and the documentation of the
-                                     * base class.
+                                     * base class for more information.
                                      */
     virtual Point<dim>
     get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const;
index 6fa6609c6639c0390e372335f831bd88fb8c712d..7a076d862ecba3502bb9104b0bcbcecb26c6ef7f 100644 (file)
@@ -55,7 +55,11 @@ Point<dim>
 StraightBoundary<dim>::get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const 
 {
   return (quad->vertex(0) + quad->vertex(1) +
-         quad->vertex(2) + quad->vertex(3)) / 4;
+         quad->vertex(2) + quad->vertex(3) +
+         quad->line(0)->child(0)->vertex(1) +
+         quad->line(1)->child(0)->vertex(1) +
+         quad->line(2)->child(0)->vertex(1) +
+         quad->line(3)->child(0)->vertex(1)) / 8;
 };
 
 #endif

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.