Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

Issue working in FES for DG

More
4 years 6 months ago #1920 by Jon
I am trying to implement discontinuous Galerkin method for steady state Stokes equations using the following finite element space:

fes = L2(mesh,order=order_p,dim=3,dgjumps=True)

However I am running into issues where my matrix is singular. I have the code working correctly for this finite element space:

X = L2(mesh,order=order_p)
fes = FESpace([X,X,X], dgjumps=True)

When I check the number of degrees of freedom is 3 times as many, am I wrong to think they should be the same as the first finite element space? While I do have it working for this space, I would really like to have code using the first listed finite element space to more easily allow me to incorporate into other code.

I've attached both my working and non-working code.
Time to create page: 0.141 seconds