From 9810c3b5b360814f409812c1e540b7a5df676415 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 18 Jan 2017 07:38:00 -0700 Subject: [PATCH] Minor cleanups to the implementation of FETools::compute_node_matrix(). --- source/fe/fe_tools.cc | 21 +++++++++------------ 1 file changed, 9 insertions(+), 12 deletions(-) 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