Hello Martin
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. Also in FE_TraceQ,
there is a private member FE_Q<dim,spacedim> fe_q. What is the use of this ?
Thanks
praveen
--
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.
