From: wolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 21 Apr 2003 16:15:40 +0000 (+0000)
Subject: Clarify one comment.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4a1befef88ae5870cb6d5e2df76d2c8c531a1b85;p=dealii-svn.git

Clarify one comment.


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

diff --git a/deal.II/deal.II/source/fe/fe_q_hierarchical.cc b/deal.II/deal.II/source/fe/fe_q_hierarchical.cc
index 8783bf4c9c..e161c60731 100644
--- a/deal.II/deal.II/source/fe/fe_q_hierarchical.cc
+++ b/deal.II/deal.II/source/fe/fe_q_hierarchical.cc
@@ -421,6 +421,25 @@ void FE_Q_Hierarchical<dim>::initialize_unit_support_points ()
                                    // support point. There is
                                    // only one support point per
                                    // vertex, line, quad, hex, etc.
+                                   //
+                                   // note, however, that the support
+                                   // points thus associated with
+                                   // shape functions are not unique:
+                                   // the linear shape functions are
+                                   // associated with the vertices,
+                                   // but all others are associated
+                                   // with either line, quad, or hex
+                                   // midpoints, and there may be
+                                   // multiple shape functions
+                                   // associated with them. there
+                                   // really is no other useful
+                                   // numbering, since the
+                                   // hierarchical shape functions do
+                                   // not vanish at all-but-one
+                                   // interpolation points (like the
+                                   // Lagrange functions used in
+                                   // FE_Q), so there's not much we
+                                   // can do here.
   unsigned int k=0;
   for (unsigned int iz=0; iz <= ((dim>2) ? degree : 0) ; ++iz)
     for (unsigned int iy=0; iy <= ((dim>1) ? degree : 0) ; ++iy)