From: bangerth Date: Sat, 11 Aug 2012 14:05:27 +0000 (+0000) Subject: Add some comments. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c7d313635d3355b60af7603e3e05054aafb909f0;p=dealii-svn.git Add some comments. git-svn-id: https://svn.dealii.org/trunk@25901 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/dofs/dof_tools.cc b/deal.II/source/dofs/dof_tools.cc index c1204f1f43..b385dd908b 100644 --- a/deal.II/source/dofs/dof_tools.cc +++ b/deal.II/source/dofs/dof_tools.cc @@ -3453,12 +3453,20 @@ namespace DoFTools // assignment is done by blocks, not by // components, as specified by // component_select. The additional argument - // component_select only is used for + // component_select is only used for // non-primitive FEs, where we need it since // more components couple, and no unique // component can be assigned. Then, we sort // them to the first selected component of the // vector system. + // + // the output array dofs_by_component + // lists for each dof the corresponding + // vector component. if the DoFHandler is + // based on a parallel distributed + // triangulation then the output array is + // index by + // dof.locally_owned_dofs().index_within_set(indices[i]) template inline void @@ -3474,10 +3482,20 @@ namespace DoFTools ExcDimensionMismatch(dofs_by_component.size(), dof.n_locally_owned_dofs())); - // next set up a table for the - // degrees of freedom on each of - // the cells whether it is - // something interesting or not + // next set up a table for the degrees + // of freedom on each of the cells + // (regardless of the fact whether it + // is listed in the component_select + // argument or not) + // + // for each element 'f' of the + // FECollection, + // local_component_association[f][d] + // then returns the vector component + // that degree of freedom 'd' belongs + // to (or, in the case of + // sort_by_blocks, the block that it + // corresponds to) std::vector > local_component_association (fe_collection.size()); for (unsigned int f=0; f