rename this to something more logical & give it a better return type
make a function that returns a vector1<pair<Size,Size> > containing all interacting pairs. This would be easy to implement.