Dear all,
When I try to create a locally refined mesh by the following codes:
* GridGenerator::hyper_cube(triangulation, 0, 1); Point<dim>
center (0.5, 0.5, 0.5); triangulation.refine_global(5);
for (auto cell : triangulation.active_cell_iterators()) {
if (cell->is_locally_owned()) { if
(cell->point_inside(center)){
cell->set_refine_flag(); } }
}triangulation.execute_coarsening_and_refinement();*
The library return this error
* const dealii::FullMatrix<double>& dealii::FiniteElement<<anonymous>,
<anonymous> >::constraints(const dealii::internal::SubfaceCase<dim>&) const
[with int dim = 3; int spacedim = 3]The violated condition was:
(this->dofs_per_face == 0) || (interface_constraints.m() != 0)Additional
information: The finite element for which you try to obtain hanging
node constraints does not appear to implement them.*
Does this error comes from the AffineConstraints when I try to implement
the DoFTools::make_hanging_node_constraints?
Could you please help me to overcome this error?
Thank you very much.
Best regards,
Kien
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/fdcc4fd8-761d-4cc4-a954-8631e0cf117d%40googlegroups.com.
