i’m assuming you don’t wanna use actual sdf volumes but you have closed geometry? you could call intersect_all and if the number of intersections is odd you’re inside
this is the crossing number
geomalgorithms.com/a03-_inclusion.html
if you can precalculate, use winding number sop