From: bangerth Date: Tue, 3 Dec 2013 13:59:47 +0000 (+0000) Subject: Include a suggestion by Denis Davydov. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=bc7ab2f88727abaf01cf9f9398de0904573164ed;p=dealii-svn.git Include a suggestion by Denis Davydov. git-svn-id: https://svn.dealii.org/trunk@31853 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-18/step-18.cc b/deal.II/examples/step-18/step-18.cc index f587d740be..1d6347f8ab 100644 --- a/deal.II/examples/step-18/step-18.cc +++ b/deal.II/examples/step-18/step-18.cc @@ -1718,6 +1718,17 @@ namespace Step18 // to the vertex we consider at present, as they may or may not be ordered // in the same order. // + // This inconvenience could be avoided if finite elements have support + // points on vertices (which the one here has; for the concept of support + // points, see @ref GlossSupport "support points"). For such a case, one + // could construct a custom quadrature rule using + // FiniteElement::get_unit_support_points(). The first + // GeometryInfo@::vertices_per_cell*fe.dofs_per_vertex + // quadrature points will then correspond to the vertices of the cell and + // are ordered consistent with cell-@>vertex(i), taking into + // account that support points for vector elements will be duplicated + // fe.dofs_per_vertex times. + // // Another point worth explaining about this short function is the way in // which the triangulation class exports information about its vertices: // through the Triangulation::n_vertices function, it