[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: dof with singular mesh

From: Edouard Oudet
Subject: Re: dof with singular mesh
Date: Tue, 11 May 2021 23:36:05 +0200
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:78.0) Gecko/20100101 Thunderbird/78.7.1

Dear Yves,
Thank you so much for your fast answer: it saved my time!
I realized that I was using an old 5.2 version of getfem. I upgraded my library to 5.4.1, recompiled and now everything is ok: the dof points list seems to be complete.
Thanks a lot,

On 5/11/21 7:38 PM, Yves Renard wrote:

Dear Edouard,

I tested with the Matlab interface a similar experience :

tmesh = gf_mesh('load', 'savemesh.mesh')
tmeshfem = gf_mesh_fem(tmesh, 1)
gf_mesh_fem_set(tmeshfem, 'fem', gf_fem('FEM_PK(2, 4)'))
pdof = gf_mesh_fem_get(tmeshfem, 'basic dof nodes');
hold on; scatter(pdof(1,:), pdof(2,:), '+'); hold off
And I do not have the problem (see attached image of the small circle)

I do not see where your result comes from.

Best regards,


Le 11/05/2021 à 17:57, Edouard Oudet a écrit :
oups, sorry


On 5/11/21 5:46 PM, Yves Renard wrote:

Dear Edouard,

Your second picture is missing. Can you send the link again ?

Best regards,


Le 11/05/2021 à 16:11, Edouard Oudet a écrit :
Dear getfem users,

Using the python interface of Getfem++, I obtained surprising results related to the representation of degrees of freedom for a singular mesh associated to a square with two circular holes. One of the hole has a very small radius see picture

This mesh can be downloaded at

Using the following instructions

filemesh = '/tmp/savemesh.mesh'
tmesh = getfem.Mesh("load", filemesh)
tmeshfem = getfem.MeshFem(tmesh, 1)
tmeshfem.set_fem(getfem.Fem('FEM_PK(2, 4)'))
pdof = getfem.MeshFem.get(tmeshfem, 'basic dof nodes')

and plotting these pdof points I obtained the following picture around the small circle

As it can be observed some dof are missing.. I do not get this problem for a degree smaller than 4.
Any suggestions are welcome,



Edouard Oudet :
IMAG - Bureau 164
700 avenue Centrale
38400 Saint Martin d'Hères
+33 (0)4 57 42 17 71  (office LJK)

reply via email to

[Prev in Thread] Current Thread [Next in Thread]