From: Wolfgang Bangerth Date: Wed, 18 Jan 2017 14:38:00 +0000 (-0700) Subject: Minor cleanups to the implementation of FETools::compute_node_matrix(). X-Git-Tag: v8.5.0-rc1~225^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9810c3b5b360814f409812c1e540b7a5df676415;p=dealii.git Minor cleanups to the implementation of FETools::compute_node_matrix(). --- diff --git a/source/fe/fe_tools.cc b/source/fe/fe_tools.cc index 4a054a7704..666ef8857f 100644 --- a/source/fe/fe_tools.cc +++ b/source/fe/fe_tools.cc @@ -1624,9 +1624,9 @@ namespace FETools const std::vector > &points = fe.get_generalized_support_points(); - // We need the values of the - // polynomials in all generalized - // support points. + // We need the values of the polynomials in all generalized support points. + // This function specifically works for the case where shape functions + // have 'dim' vector components, so allocate that much space std::vector > values (dim, std::vector(points.size())); @@ -1634,12 +1634,10 @@ namespace FETools // result of the interpolation std::vector local_dofs(n_dofs); - // One row per shape - // function. Remember that these - // are the 'raw' shape functions - // where the inverse node matrix is - // empty. Otherwise, this would - // yield identity. + // Get the values of each shape function in turn. Remember that these + // are the 'raw' shape functions (i.e., where the element has not yet + // computed the expansion coefficients with regard to the basis + // provided by the polynomial space). for (unsigned int i=0; i