Excluding grid in calculations

[Antoine Thomas]

Hello,

I'm trying to calculate HRTFs for various models with bempp. I'm using the sphere scattering tutorial as a reference. In the tutorial, we don't compute the results inside the sphere thanks to its equation on the plane. I've been looking all day for a way to do the same with models that are more complex, but haven't found any.

Is there a known way to do this ?

Thanks a lot,

Antoine THOMAS

[Matthew Scroggs]

Hello,

For shapes more complicated than a sphere or cube, you'll need a way to decide whether or not a point is inside (eg count the number of times a line from (-10,-10,-10) to the point crosses the surface. Once you've got this, something along the lines of:

```

def inside(x, y, z):

    # A function that returns True if (x, y, z) is inside shape

    

idx = []

for i, (x, y, z) in enumerate(points)

    if inside(x, y, z):

        idx.append(i)

```

Then using the potential operators as in the example.

Hope this helps,

Matthew