Hi Dr.Bangerth, Thank you for your quick reply! Yes, I'm still using deal.II to implement the WG. I'm now trying to use WG(P1,P1) to solve Darcy. For WG(P0,P0), I used FE_DGQ and FE_FaceQ which are the same as FE_DGP and FE_FaceP. But for higher elements, I need to use FE_FaceP. If the face dof is correctly figured out, I suppose I can assign the D.B.C to constraints by myself. Is that correct?
Thank you. Best, Zhuoran On Sunday, June 6, 2021 at 9:21:44 PM UTC+8 Wolfgang Bangerth wrote: > On 6/6/21 3:55 AM, Zhuoran Wang wrote: > > > > I thought that FE_FaceP has support points on boundaries. Does it mean > that I > > can't set the Dirichlet boundary conditions in this way? Or did I > > misunderstand anything? > > Hi Zhuoran, > nice to see that you're still using deal.II :-) > > I believe that the FaceP spaces don't have support points in 3d whereas > the > FaceQ spaces do. That's for the same reason as why a FE_DGP<2> can have no > support points: because they're not a complete tensor product polynomial. > > Do you need to use FE_FaceP? Could you use FE_FaceQ? > > Best > W. > > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@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/d6925ad8-a2de-4952-8902-af6174c16147n%40googlegroups.com.
