Re: dof with singular mesh

[Edouard Oudet]

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,
best,
Edouard.

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');
> gf_plot_mesh(tmesh);
> 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,
>
> Yves
>
>
>
>
>
>
>
> Le 11/05/2021 à 17:57, Edouard Oudet a écrit :
> > oups, sorry
> >
> > https://www-ljk.imag.fr/membres/Edouard.Oudet/download/dof_singular.png
> >
> > Best
> >
> > 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,
> > >
> > > Yves
> > >
> > >
> > > 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
> > > >
> > > > https://www-ljk.imag.fr/membres/Edouard.Oudet/download/singmesh.png
> > > >
> > > > This mesh can be downloaded at
> > > >
> > > > https://www-ljk.imag.fr/membres/Edouard.Oudet/download/savemesh.mesh
> > > >
> > > > 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
> > > >
> > > > https://www-ljk.imag.fr/membres/Edouard.Oudet/download/savemesh.mesh
> > > >
> > > > 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,
> > > >
> > > > best,
> > > >
> > > > Edouard.

--
Edouard Oudet : http://www-ljk.imag.fr/membres/Edouard.Oudet/
IMAG - Bureau 164
700 avenue Centrale
38400 Saint Martin d'Hères
+33 (0)4 57 42 17 71  (office LJK)