GitLab

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:

1. compute the bounding box
2. pick one corner and move from there away from the opposite corner by an arbitrary length
3. pick an arbitrary triangle from the surface mesh
4. compute its balance point and move a small distance from it in normal direction
5. 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