## esys-escript

Dear all,
How could I do, if I want to know how many gauss points in the circle region with radius r0.
What I have done is'
gx = Function(mydomain).getX()
region = r0-length(gx-ic)
nt = whereNegative(region)'
But I don't know how to print how many points where nt is equal to 0
Regards,
Lou Jichao

## Question information

Language:
English Edit question
Status:
Solved
For:
esys-escript Edit question
Assignee:
No assignee Edit question
Solved by:
Solved:
Last query:
 Revision history for this message Lou Jichao (ljc-1) said on 2019-09-29: #1

And what should I do if I want to know the nt value of with the number 5 gauss point.

 Revision history for this message  Adam Ellery (aellery) said on 2019-09-29: #2

If the circle is centered on (x0,y0) and has radius r0 then it is described by the equation

(x-x0)^2 + (y-y0)^2 = r0^2,

where x and y are the spatial coordinates, and all interior points satisfy the inequality

(x-x0)^2 + (y-y0)^2 - r0^2 < 0.

You can find these points using something like the following:

gx = Function(mydomain).getX()
tmp = (gx-ic)*(gx-ic)
tmpB = tmp+tmp-r0*r0
tmpC = whereNegative(tmpB)
 Revision history for this message Lou Jichao (ljc-1) said on 2019-09-30: #3