Added function to compute wheter the initial points are in or out the primaries' realms.

This improves the choice of the sign of the velocity, to better resemble manifold intersection computation
Correct double conversion from pointer to value
Updated version
included inorout.cpp