I have modified PointInPolygon#point_location to use random_ray directly instead
Then it becomes very fast for me
And the result is acceptable too (observed on a map with markers and polygons)
So I wonder if it's possible to add an option to pass when calculating
I might make a PR later