Thanks for the clarification, Dr. Bangerth.
I think I will use the zero boundary condition to test whether finite
element spaces are appropriate.
Best,
Zhuoran
On Tue, Jun 8, 2021 at 12:19 AM Wolfgang Bangerth
wrote:
> On 6/6/21 10:35 PM, Zhuoran Wang wrote:
> > I'm now trying to use WG(P1,P1)
On 6/6/21 10:35 PM, Zhuoran Wang wrote:
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
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
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
Hi,
I have a question about the Dirichlet boundary conditions.
I'm using fe(FE_DGP(1),1, FE_FaceP(1),1) as the finite element
space, FE_DGP is defined in interiors, FE_FaceP is defined on faces.
When I try to set the Dirichlet boundary conditions to the system by
*VectorTools::interpolate_boundar