WIP: Optimise inside bounds checks

Closes #666 (closed)

This is rationalising:

• the calls to transform(...) should be kept low (reuse if possible)
• the information that a position is actually on the surface is used in the intersect(...) method