fix(geospatial): increase epsilon tolerance to prevent particle loss

- Change default epsilon from 1e-6 to 1e-3 to ensure boundary points are included
- Update docstrings to avoid specifying parameter defaults
- Remove redundant test_geospatial_benchmark.py (covered by test_geospatial_performance.py)

Fixes particle tracking (PRT) test failures where 10% of particles were lost
at cell boundaries due to insufficient geometric tolerance.

🤖 Generated with [Claude Code](https://claude.com/claude-code)

Co-Authored-By: Claude <[email protected]>

Author: adacovsk