Issue to define FE space of MeshS type

General Discussion
1

Hello All,

I build a 2D mesh on the SALOME Plateform (a grid in the plane OXY). Then, I import this mesh in FreeFEM like this

meshS TH=readmeshS("coarse_grid.mesh");

However, it seems that I cannot do

fespace P0P0(TH,[P0, P0])

on this mesh (no matter the FE space in the bracket). I have a Segmentation Error.

Before, I was doing this with a 2D mesh of type mesh (no meshS) and it was working well. I don’t understand why.

Could someone help me with this problem ?

Otherwise, how could I change meshS type to mesh Type ?

Thank you in advance,

Best regards,

Loïc,

2

For more precisions,

This FE space VH=fespace(P0) works well, but this FE space VH=fespace([P0]) does not work.

It makes me think at the problem:

but I am on version 4.9.

Loïc,

3

I have try this

load "msh3"
meshS TH=square3(3,3);
//However, it seems that I cannot do
fespace P0P0(TH,[P0, P0]);

It seems Ok

4

Hello @frederichecht ,

Thank you for your reply,

This does not work with my version of FreeFEM (I am on version 4.9). I don’t understand why. It works also with @prj.

Do you have an idea about why this does not work on my machine ?

Thank you in advance,

Best regards,

Loïc,

5
  1. Can you send the “coarse_grid.mesh” , I want to try on my laptop.
  2. You run FreeFEM++ on with type of computer?
6

Hello @frederichecht ,

Thank you for your message, here my mesh :
coarse_grid.mesh (266 Bytes)

I have even tried this:

load "msh3"
meshS TH=square3(3,3);
//However, it seems that I cannot do
fespace P0P0(TH,[P0, P0]);

But it does not work too.

I run FreeFEM++ on UBUNTU 20.

Best regards,

Loïc,

7

I will try on this computer.

8

Thank you @frederichecht,

Let me know about the result,

Best regards,

Loïc,

9

I think the bug come form and correct in
commit baf1996585e5dbc09b591e2dcf62b9e5e8e88a31
Author: Frederic Hecht frederic.hecht@upmc.fr
Date: Wed Mar 24 14:52:06 2021 +0100

correct problem of initialisatio of static after usage on  ubuntu os,
remove the static of G (barycenter) use in PkLagrange.hpp had to find ...

So no way with version 4.9

10

Thank you for your reply @frederichecht .

So I have to compile the FreeFEM version v.10 ?

@prj already advise me to compile on dev branch.

Best regards,

Loïc,

11

Works well with FreeFEM version v10.

Loïc,