]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add two more pictures since they are so nice.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 24 Jan 2008 16:24:58 +0000 (16:24 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 24 Jan 2008 16:24:58 +0000 (16:24 +0000)
git-svn-id: https://svn.dealii.org/trunk@15667 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/include/numerics/vectors.h
deal.II/deal.II/include/numerics/vectors.templates.h

index 4d229dbf4add3e06ad23c1b8b48e4a5e31aa1335..1b446557680d8688bd35786691a0d0c298747082 100644 (file)
@@ -964,6 +964,32 @@ class VectorTools
                                      * face of the right cell belongs
                                      * to it. This case is not
                                      * currently implemented.
+                                     *
+                                     *
+                                     * <h4>Results</h4>
+                                     *
+                                     * Because it makes for good
+                                     * pictures, here are two images
+                                     * of vector fields on a circle
+                                     * and on a sphere to which the
+                                     * constraints computed by this
+                                     * function have been applied:
+                                     *
+                                     * <p ALIGN="center">
+                                     * @image html no_normal_flux_5.png
+                                     * @image html no_normal_flux_6.png
+                                     * </p>
+                                     *
+                                     * The vectors fields are not
+                                     * physically reasonable but the
+                                     * tangentiality constraint is
+                                     * clearly enforced. The fact
+                                     * that the vector fields are
+                                     * zero at some points on the
+                                     * boundary is an artifact of the
+                                     * way it is created, it is not
+                                     * constrained to be zero at
+                                     * these points.
                                      */
     template <int dim, template <int> class DH>
     static
index 3c071a17e61beaaa3957feca6507884a1adcf20a..e944621fa770b5a4430c98b89db66f4e6645155e 100644 (file)
@@ -2333,7 +2333,7 @@ VectorTools::compute_no_normal_flux_constraints (const DH<dim>         &dof_hand
                                             // then construct constraints
                                             // from this:
            const internal::VectorTools::VectorDoFTuple<dim> &
-               dof_indices = same_dof_range[0]->first;     
+             dof_indices = same_dof_range[0]->first;       
            internal::VectorTools::add_constraint (dof_indices, normal,
                                                   constraints);
 

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.