Update: It looks as if this issue: https://github.com/dealii/dealii/issues/2237 might be related to the currently occuring problem?
Am Donnerstag, 17. Januar 2019 08:42:01 UTC+1 schrieb Maxi Miller: > > For convergence tests I wanted to replace the commonly used elements FE_Q > with FE_Bernstein. Unfortunately apparently it is not possible to simply > replace all occurences. When running my program, I get the following error: > -------------------------------------------------------- > An error occurred in line <1126> of file <~/Downloads/git-files/dealii/ > source/fe/fe.cc> in function > dealii::Point<(dim - 1)> dealii::FiniteElement<<anonymous>, > <anonymous> >::unit_face_support_point(unsigned int) const [with int dim = > 2; int spacedim = 2] > The violated condition was: > unit_face_support_points.size() == this->dofs_per_face > 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. > > Stacktrace: > ----------- > #0 /opt/dealii/lib/libdeal_II.g.so.9.1.0-pre: dealii::FiniteElement<2, > 2>::unit_face_support_point(unsigned int) const > #1 /opt/dealii/lib/libdeal_II.g.so.9.1.0-pre: dealii::FESystem<2, > 2>::unit_face_support_point(unsigned int) const > #2 ./main: ) [0x8953ea] > #3 ./main: void dealii::VectorTools::interpolate_boundary_values<2, 2, > dealii::DoFHandler, double>(dealii::Mapping<2, 2> const&, > dealii::DoFHandler<2, 2> const&, std::map<unsigned int, dealii::Function<2, > double> const*, std::less<unsigned int>, std::allocator<std::pair<unsigned > int const, dealii::Function<2, double> const*> > > const&, > std::map<unsigned int, double, std::less<unsigned int>, > std::allocator<std::pair<unsigned int const, double> > >&, > dealii::ComponentMask const&) > #4 ./main: void dealii::VectorTools::interpolate_boundary_values<2, 2, > dealii::DoFHandler, double>(dealii::Mapping<2, 2> const&, > dealii::DoFHandler<2, 2> const&, std::map<unsigned int, dealii::Function<2, > double> const*, std::less<unsigned int>, std::allocator<std::pair<unsigned > int const, dealii::Function<2, double> const*> > > const&, > dealii::AffineConstraints<double>&, dealii::ComponentMask const&) > #5 ./main: void dealii::VectorTools::interpolate_boundary_values<2, 2, > dealii::DoFHandler, double>(dealii::Mapping<2, 2> const&, > dealii::DoFHandler<2, 2> const&, unsigned int, dealii::Function<2, double> > const&, dealii::AffineConstraints<double>&, dealii::ComponentMask const&) > #6 ./main: void dealii::VectorTools::interpolate_boundary_values<2, 2, > dealii::DoFHandler, double>(dealii::DoFHandler<2, 2> const&, unsigned int, > dealii::Function<2, double> const&, dealii::AffineConstraints<double>&, > dealii::ComponentMask const&) > #7 ./main: ) [0x7f32a3] > #8 ./main: ) [0x7f0d58] > #9 ./main: main > > > > with interpolate_boundary_values() called as > > VectorTools::interpolate_boundary_values(dof_handler, 0, > BoundaryValuesCarriers<dim, EQU>(n_components), boundary_constraints, > A_mask); > VectorTools::interpolate_boundary_values(dof_handler, 0, > BoundaryValuesB<dim, EQU>(n_components), boundary_constraints, B_mask); > VectorTools::interpolate_boundary_values(dof_handler, 0, > BoundaryValuesC<dim, EQU>(n_components), boundary_constraints, C_mask); > > boundary_constraints.close(); > > > Is there another way how I could replace FE_Q-elements with > FE_Bernstein-elements? > Thanks! > -- 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.
