Re: [deal.II] make hanging node constraint for locally refined mesh
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 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 > nodes. > > Best > W. > > > -- > > Wolfgang Bangerth email: bange...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- > 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. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/44ffc21e-c3e2-4f48-a4ec-1537f37e8b42%40colostate.edu > . > -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAAo%2BSZcyBae4j_-0ahdsQSrVCbTCrrVOHpixTGiraB8_zONUjg%40mail.gmail.com.
Re: [deal.II] make hanging node constraint for locally refined mesh
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 nodes. Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/44ffc21e-c3e2-4f48-a4ec-1537f37e8b42%40colostate.edu.
Re: [deal.II] make hanging node constraint for locally refined mesh
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 5, 2020 vào lúc 05:42 Wolfgang Bangerth < bange...@colostate.edu> đã viết: > 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) > > 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? > > Well, the error message says pretty clearly what is happening: The element > you > use does not provide hanging node information. What element are you using? > > Best > W. > > -- > > Wolfgang Bangerth email: bange...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- > 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. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/fc63875b-df35-eb80-02cd-c503ca81d82e%40colostate.edu > . > -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAAo%2BSZfXpVTE1NW1gMpbsSsng08KXFaJfHqO%2B%2BUVDiMoqFjGXA%40mail.gmail.com.
Re: [deal.II] make hanging node constraint for locally refined mesh
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) 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? Well, the error message says pretty clearly what is happening: The element you use does not provide hanging node information. What element are you using? Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/fc63875b-df35-eb80-02cd-c503ca81d82e%40colostate.edu.
[deal.II] make hanging node constraint for locally refined mesh
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()) { 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& 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)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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/fdcc4fd8-761d-4cc4-a954-8631e0cf117d%40googlegroups.com.