Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Phạm Ngọc Kiên
Dear Prof. Wolfgang Bangerth, Thank you very much. I will change the codes as your guidance. Yours sincerely, Kien Vào Th 4, 27 thg 5, 2020 vào lúc 12:58 Wolfgang Bangerth < bange...@colostate.edu> đã viết: > On 5/26/20 7:30 PM, Phạm Ngọc Kiên wrote: > > > > Thus, the conforming FE_Nedelec

Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Wolfgang Bangerth
On 5/26/20 7:30 PM, Phạm Ngọc Kiên wrote: Thus, the conforming FE_Nedelec and FE_NedelecSZ can only work with the conforming mesh. The non-conforming mesh, which has hanging nodes, does not work with these elements. The FE_NedelecSZ element doesnt, but FE_Nedelec does work with hanging

Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Phạm Ngọc Kiên
Dear Prof. Wolfgang Bangerth, I used FE_NedelecSZ element. Thus, the conforming FE_Nedelec and FE_NedelecSZ can only work with the conforming mesh. The non-conforming mesh, which has hanging nodes, does not work with these elements. Thank you very much. Yours sincerely, Kien Vào Th 4, 27 thg

Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Wolfgang Bangerth
On 5/26/20 1:34 AM, Phạm Ngọc Kiên wrote: /*  const dealii::FullMatrix& dealii::FiniteElement<, >::constraints(const dealii::internal::SubfaceCase&) const [with int dim = 3; int spacedim = 3] The violated condition was:     (this->dofs_per_face == 0) || (interface_constraints.m() != 0)

[deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Phạm Ngọc Kiên
Dear all, When I try to create a locally refined mesh by the following codes: * GridGenerator::hyper_cube(triangulation, 0, 1);Point center (0.5, 0.5, 0.5);triangulation.refine_global(5); for (auto cell : triangulation.active_cell_iterators()) {