How to use notaregion (identify if a point is within the domain)

[julienG (Julien)]

You can do this

mesh Th;
Th = square(10,10,[0.+1.*x,0.+1.*y]);
fespace V0h(Th,P0);
V0h Inside=1;
//
real xx=0.5; real yy=0.5;
if(!Inside(xx,yy)) cout << “outside” << endl;
else cout << “inside” << endl;