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.

dofs of FacetFESpace

More
6 years 6 months ago #206 by rhebergens
dofs of FacetFESpace was created by rhebergens
Hello,

I'm running the following code:
Code:
from ngsolve import * from netgen.geom2d import SplineGeometry geo = SplineGeometry() geo.AddRectangle( (0.0, 0.0), (1.0, 1.0), bc=1 ) mesh = Mesh( geo.GenerateMesh(maxh=1.6)) mesh.Refine() mesh.Refine() mesh.Refine() Q = FacetFESpace(mesh, order=1) print("ndof qfe:", Q.ndof)

For the FacetFESpace, for order 1, I expect 2 DOFs per edge, but it seems I'm getting slightly more. There are a total of 208 edges, and I get a total of 493 DOFs for the whole mesh (I expect 416). For order 2, I expect to 624 DOFs, but get 701. As I increase the order, I realised that the total number of DOFs I get for the FacetFESpace, for the given mesh, is: DOFs = (p+1)*nrEdges + 77. Could you explain to me where the 77 extra DOFs comes from?

Thanks,
Sander
More
6 years 6 months ago #208 by schruste
Replied by schruste on topic dofs of FacetFESpace
Dear Sander,

This is related to facets from the former meshes that you refined from. To these older facets the low order dofs are still created. However, you may check the CouplingTypes to see that they are marked as unused as they actually never appear in the usual assembly (no non-zero couplings will be reserved for these dofs). If you don't like to have these 'virtual' facet dofs, you can save and reload the mesh which removes the refinement history and in turn also the virtual facets and hence the virtual dofs. You will get as many dofs as you expect.

Best,
Christoph
More
6 years 6 months ago #209 by rhebergens
Replied by rhebergens on topic dofs of FacetFESpace
Thanks for the explanation!
Sander
Time to create page: 0.104 seconds

Login