I have a rectangular mesh of 500000x500000. I extracted the coordinates of the mesh elements. Now I would like to remove elements in the semicircular form from this rectangular mesh. I have the function to remove elements. But I am struggling to set the condition for the semicircular region to determine all the elements inside the semicircular region. For example if ex and ey are the coordinates of elements in mesh, then I am using
find(ex<=R*cos(theta) & ey <=R*sin(theta)).
How to set this condition according to the grid points in the rectangular mesh. I have also attached the image for more clarity.