Fixed reference point for inside check
In the class ArbitraryDomain globalIsInside_m is hard coded. This should be computed automatically for any domain and shouldn't be hard coded or needed from the input file. Possibility:
- compute the bounding box
- pick one corner and move from there away from the opposite corner by an arbitrary length
- pick an arbitrary triangle from the surface mesh
- compute its balance point and move a small distance from it in normal direction
- test whether the connecting line crosses an odd number of times
- if the test is positive then you're done
- otherwise move the same distance from the triangle to the other side
- compute number crosses with surface
- if the test is still negative then diminish distance from the triangle by factor 2 and try on both sides if necessary