Forum de mathématiques - Bibm@th.net

ccapucine

Membre

Inscription : 19-05-2018

Messages : 180

FreeFem++

Bonjour
est-ce qu'il y a parmi vous quelqu'un qui connaît bien Freefem++ pour une petite question.

Bien cordialement

Roro

Membre expert

Inscription : 07-10-2007

Messages : 1 565

Re : FreeFem++

Bonsoir,

On ne peut pas dire que je connaisse bien, mais pose quand même la question !

Roro.

ccapucine

Membre

Inscription : 19-05-2018

Messages : 180

Re : FreeFem++

Bonjour
ma question est comment programmer en freefem++ la fnonction suivante

$$
\begin{cases}
1& \mbox{si} (\dfrac{x}{\epsilon},\dfrac{y}{\epsilon}) \in B + \mathbb{Z}^2\\
0& \mbox{sinon}
\end{cases}
$$

où $B$ est une boule de mesure $1/w$.

Cordialement

Roro

Membre expert

Inscription : 07-10-2007

Messages : 1 565

Re : FreeFem++

Bonsoir,

Tu peux t'inspirer de la fonction suivante :

func f=1*(x^2+y^2<1);

qui doit (sauf erreur) retourner $1$ si $(x,y)$ est dans la boule de centre $0$ et de rayon $1$...

Concernant la périodicité, en général on travaille sur une "cellule" et on impose des conditions au bord de type "périodique". Ce qui en Freefem ressemble à :

mesh Th=square(10,10,[2*x*pi,2*y*pi]);
fespace Vh(Th,P2,periodic=[[2,y],[4,y],[1,x],[3,x]]);

Roro.

ccapucine

Membre

Inscription : 19-05-2018

Messages : 180

Re : FreeFem++

Merci Roro

rab1920

Invité

Re : FreeFem++

Est ce quelqu un peut m'aider sur la fonction "convect" de freefem car j'arrive pas à comprendre son fonctionnement?