Dear dealii community,
I am trying to use "DoFTools::map_dofs_to_support_points" function for a two-field problem. The components of the first (vectorial) field are approximated using "FE_Q" and the second (scalar) field with "FE_DGPMonomial". The "component_mask" parameter is used to select only the dofs associated with the first field. However I get the following error message (debug mode): * void dealii::DoFTools::internal::(anonymous namespace)::map_dofs_to_support_points(const hp::MappingCollection<dim, spacedim> &, const DoFHandler<dim, spacedim> &, std::map<types::global_dof_index, Point<spacedim>> &, const dealii::ComponentMask &) [dim = 2, spacedim = 2]* *The violated condition was: * * (fe_collection[fe_index].n_dofs_per_cell() == 0) || (fe_collection[fe_index].has_support_points())* *Additional information: * * You are trying to access the support points of a finite element that* * either has no support points at all, or for which the corresponding* * tables have not been implemented.* >From inspecting the method responsible for throwing the error is appears that all fields in the collection are checked via the Assert. Is the method applicable to problems where some fields have support points and others don't but are masked? The mask appears to only be appropriate to mask fields with support points. Please find attached a minimal code for a two-field formulation that demonstrates the problem. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/ef66c211-1290-43b4-9126-a9eabef00c55n%40googlegroups.com.
<<attachment: test_mdtsp.zip>>
