From: bangerth Date: Thu, 24 Jan 2008 16:24:58 +0000 (+0000) Subject: Add two more pictures since they are so nice. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8fd3dfc2685875aba5c0d7be8b820f7cc32dc8a2;p=dealii-svn.git Add two more pictures since they are so nice. git-svn-id: https://svn.dealii.org/trunk@15667 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/include/numerics/vectors.h b/deal.II/deal.II/include/numerics/vectors.h index 4d229dbf4a..1b44655768 100644 --- a/deal.II/deal.II/include/numerics/vectors.h +++ b/deal.II/deal.II/include/numerics/vectors.h @@ -964,6 +964,32 @@ class VectorTools * face of the right cell belongs * to it. This case is not * currently implemented. + * + * + *

Results

+ * + * 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: + * + *

+ * @image html no_normal_flux_5.png + * @image html no_normal_flux_6.png + *

+ * + * 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 class DH> static diff --git a/deal.II/deal.II/include/numerics/vectors.templates.h b/deal.II/deal.II/include/numerics/vectors.templates.h index 3c071a17e6..e944621fa7 100644 --- a/deal.II/deal.II/include/numerics/vectors.templates.h +++ b/deal.II/deal.II/include/numerics/vectors.templates.h @@ -2333,7 +2333,7 @@ VectorTools::compute_no_normal_flux_constraints (const DH &dof_hand // then construct constraints // from this: const internal::VectorTools::VectorDoFTuple & - dof_indices = same_dof_range[0]->first; + dof_indices = same_dof_range[0]->first; internal::VectorTools::add_constraint (dof_indices, normal, constraints);