Why is my code note working ?( simple code)

General Discussion
1

Good eveneing everyone , I am a begginer on freemfem++. My programm doesnt compile and i just cant understand why. The problems seems to be due to the int2d function. Thank you for your help! Error message : "int2d(Vh) error operator ( <PP5v_fes>
List of (hoices <20CDomainOfIntegration> : … )
( <20CDomainOfIntegration> : … )
( <20CDomainOfIntegration> : … )

Error line number 21, "
: the programm is as follow ( the comments are in french ).
real a = 1; // Demi axe vertical

real r0 = 1.5; // Demi axe horizontal

real R = 5; // Taille de la boîte englobante

real f=20; //precision du maillage

border omega(t=0, pi){x =r0sin(t); y = acos(t);label=0;};

border axis2(t=-a,-R){x=0;y=t;label=1;};

border side3(t=0,R){x=t;y=-R;label=1;};

border side4(t=-R,R){x=R;y=t;label=1;};

border side5(t=R,0){x=t;y=R;label=1;};

border axis1(t=R,a){x=0;y=t;label=1;};

plot(omega(f/2)+axis2(f)+side3(f)+side4(2*f)+side5(f)+axis1(f));

mesh maillage =buildmesh (omega(f/2)+axis2(f)+side3(f)+side4(2*f)+side5(f)+axis1(f));

//plot(maillage);

//mesh maillageup=adaptmesh(maillage,50);

//plot(maillageup);

fespace Vh(maillage ,P1);

Vh uh,vh;

problem MonProbleme(uh,vh)=

int2d(Vh)(dx(uh)*dx(vh) + dy(uh)*dy(vh))

+on(1,uh=0)

+on(0,uh=1);

MonProbleme;

2

Integration must be on the mesh, not the FE-space, so: int2d(maillage)(...);

1 Like