Dear Praveen, > > When I use FE_TraceQ and am setting initial condition with > VectorTools::interpolate, I get this error > > -------------------------------------------------------- > > An error occurred in line <127> of file > </Users/praveen/Applications/deal.II/git/include/deal.II/numerics/vector_tools.templates.h> > in function > > void dealii::VectorTools::interpolate(const Mapping<dim, spacedim> > &, const DoFHandlerType<dim, spacedim> &, const Function<spacedim, > typename VectorType::value_type> &, VectorType &, const > dealii::ComponentMask &) [VectorType = > dealii::TrilinosWrappers::MPI::Vector, dim = 2, spacedim = 2, > DoFHandlerType = DoFHandler] > > The violated condition was: > > > (fe[fe_index].base_element(fe[fe_index].component_to_base_index(component_index).first).has_support_points()) > || > (fe[fe_index].base_element(fe[fe_index].component_to_base_index(component_index).first).dofs_per_cell > == 0) > > The name and call sequence of the exception was: > > ExcNonInterpolatingFE() > > Additional Information: > > You are attempting an operation that requires the finite element > involved to be 'interpolating', i.e., it needs to have support points. > The finite element you are using here does not appear to have those > for the required components. > > > which shouldn't occur since this is interpolating space. You are right, we should allow for interpolation of initial condition. We also allow this for FE_FaceQ. I will soon open a PR. On the other hand, we will probably not be able to project onto FE_TraceQ with our tools because this element does not define basis functions inside the elements, only on the faces. (One can write manual project() functions for FE_FaceQ, but I wouldn't know how to do this for FE_TraceQ.)
> Also in FE_TraceQ, there is a private member FE_Q<dim,spacedim> fe_q. > What is the use of this ? It is strictly for internal use in order to simplify the implementation. We use it for querying all kind of information, given that the degrees of freedom of FE_TraceQ are a subset of the ones in FE_Q. The last degrees of freedom of FE_Q, belonging to the interior of the element, are left out, and polynomial spaces are interpreted as the restriction of the FE_Q spaces to the surfaces. Best, Martin -- 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. For more options, visit https://groups.google.com/d/optout.
